Hi,
I checked Matlab code of 'case 5' with the data of 'Small test systems for
power system economic studies' but in the part '*generator data*' I could
not get how to write *'Pg'* , *'Qg'*,* 'Qmax',* *'Qmin'*, *'Pmax,'* and *'*
*Qmin'*.
For example, if I use Table 1 (Generation parameters) the
gt;
> Hope this helps
>
> --
> *Od:* bounce-121948758-78499...@list.cornell.edu [
> bounce-121948758-78499...@list.cornell.edu] v imenu naime ahmadi [
> naimeahm...@gmail.com]
> *Poslano:* 16. oktober 2017 11:02
> *Za:* matpowe...@list.cornell.e
I found the data in "small test system for power system economic studies".
The problem is for getting this data do I need to run load flow?
If yes how I can run it for 5 case buses.
On Mon, Oct 16, 2017 at 8:59 PM, naime ahmadi <naimeahm...@gmail.com> wrote:
> Sorry, you mea
Hi,
I could not understand what PF, QF, PT, QT means?
For example if I have a p1(power which is injected to bus one) or P12(power
flow between 1 and 2) what should I have to define it with PF.
Also in the part test_se in state estimation, I could not understand what
idx.idx_zPF = [1;2] means. is
al or
> README <https://github.com/MATPOWER/matpower/#running-matpower>.
>
> runpf('case5')
>
> Does that answer your question?
>
> Ray
>
>
>> On Oct 17, 2017, at 4:57 AM, naime ahmadi <naimeahm...@gmail.com> wrote:
>>
>> I found the da
her than bus injections then you want to use Bf rather than Bbus from
> the return values of makeBdc().
>
> Ray
>
>
> On Apr 19, 2018, at 10:16 PM, naime ahmadi <naimeahm...@gmail.com> wrote:
>
> Dear Ray,
> Thank you. I want to calculate this Jacobin.
> P(i,j)=
ng of theory...
>
> Regards
> Saeed Ahmed
>
>
> On 9 April 2018 at 15:46, naime ahmadi <naimeahm...@gmail.com> wrote:
>
>> Hi,
>> How can I use matpower to calculate DC state estimation?
>> Kind Regards,
>> Naime
>>
>
>
Hi,
How can I use matpower to calculate DC state estimation?
Kind Regards,
Naime
edu/matpower/docs/ref/matpower6.0/makeBdc.html>.
>
> Ray
>
>
> On Apr 10, 2018, at 5:04 AM, naime ahmadi <naimeahm...@gmail.com> wrote:
>
> Thank you.
> Actually I read that chapter. My problem is getting DC Jacobian matrix
> from matpower.
> Does j11 be the DC J
Dear Ray,
Thank you. I want to calculate this Jacobin.
P(i,j)=B(i,j)*(theta(i)-theta(j))
PP=[P(1,2);P(1,4);P(1,5);P(2,3);P(2,4);P(4,5)](measurement vector)
J=double(jacobian(PP,theta)) getting derivatives of power flow between
line i and j in respect to all other theta
J = J(:,2:end);
Hi,
The result of voltage and angle in AC state estimation should be the same
as AC power flow.
Is it similar for DC case?
Kind Regards,
Naime
What should I do?
I have a question.
On Mon, Mar 26, 2018 at 5:28 PM, naime ahmadi <naimeahm...@gmail.com> wrote:
> Dear all,
> I have a two questions.
> I want to change Jacobin Matrix for IEEE 14 bus system and after that get
> the state estimation.
> I do not know whe
Hi,
I have 2 questions.
1-I run AC state estimation in matpower (use runse and state_est).
When I see these plots
subplot(3,2,1), plot(180/pi*(angle(Vlf)-angle(V)),'.'), title('Voltage
Angle (deg)');
subplot(3,2,2), plot(abs(Vlf)-abs(V),'.'), title('Voltage Magnitude
(p.u.)');
Error in angle is
ation of buses.this picture of document of
> Math power .
> beast regards.
>
> On Thu, Dec 20, 2018 at 7:17 AM naime ahmadi
> wrote:
>
>> I want to plot IEEE 14 and 30 bus system.
>> how can I do this with this bus_position?
>>
>> On Thu, Dec 20, 2018 at 12
Hi,
I have question about branch data in case 57 bus system in matpower.
When I look at diagram of IEEE 57 bus system, I can see it have some
crossing in the network between lines.
Are there any crossing in power system?
Thanks,
Naime
Hi,
How can I get position of buses in IEEE 14 or 30 bus system in matpower?
Thanks,
Naime
I want to plot IEEE 14 and 30 bus system.
how can I do this with this bus_position?
On Thu, Dec 20, 2018 at 12:11 PM Meraj Fouad wrote:
> Hi, this will give the bus position.
>
> file='case30';
>
> mpc=loadcase(file);
>
> bus_position=mpc.bus(:,1)
>
> On Wed, De
Hi,
I have a question about 300 bus system. When I opened the case300, in bus
data I saw 1190 as bus number. or in branch data, bus 9001 is connected to
9005. As case300 has only 300 bus so what is 1190 and 9001?
How can I fix this case?
Thanks,
Naime
Hi
I have a question about state estimation in MATPOWER.
In 118 bus system the slack bus is 69. The voltage and angle on bus 69
should be 1 and 0.
But when I check Vlf (power flow voltage of bus 69) it is 0.9083 + 0.3784i.
So, when I check the error in function runse
Hi
I have a question about state estimation part ''index vector for state
variables to be updated''
ww = [ nref; nb+nref ];
In runse, angle and voltage of slack bus is known. If I only know angle of
slack bus, how I can change the code? I want to calculate voltage magnitude
of slack bus.
Thanks,
Hi friends,
I have a question about the covariance matrix in 'State_est' function.
Below is the sigma which is used in Matpower.
sigma = [
0.02 * abs(Sf) + 0.0052 * fullscale * ones(nbr,1);
0.02 * abs(St) + 0.0052 * fullscale * ones(nbr,1);
0.02 * abs(Sbus)+ 0.0052 *
Hello Matpower users,
I want to remove duplicate branches in IEEE 118 bus system and put an
equivalent one. For example:
fbus tbus rxb rateA rateB rateC ratio angle status angmin angmax
duplicate branches:
42 49 0.0715 0.323 0.086 0 0 0 0 0 1 -360 360;
42 49 0.0715 0.323
t; However, a quick correction … while the total impedance (*r* and *x*) are
> halved, the total line charging (*b*) is doubled.
>
> Ray
>
>
> On Apr 19, 2021, at 9:59 AM, naime ahmadi wrote:
>
> Hello Matpower users,
>
> I want to remove duplicate branches in
>
> You can also easily check by running a power flow before and after
> consolidating the lines. You should get *exactly* the same solution if
> you did it correctly.
>
> Hope this clarifies,
>
>Ray
>
> On Apr 19, 2021, at 11:37 AM, naime ahmadi wrote:
>
>
Dear Professor Ray,
I have a 3 bus system and want to calculate the DC Jacobian. I attached the
result that I want. My problem is that I do not know which function in
MATPOWER does exactly this calculation.
I also attached the 3 bus system in MATPOWER format.
function mpc =
Dear Dr. Sanchez,
Thank you for your time and help.
Kind Regards,
Naime
On Thu, Nov 25, 2021 at 7:18 AM Carlos E Murillo-Sanchez <
ce.murillosanc...@gmail.com> wrote:
> naime ahmadi wrote:
> > Dear Professor Ray,
> > I have a 3 bus system and want to calculate the DC Jac
Hello,
How can I change the power injection to flow measurement?
For example, in the IEEE 14-bus system, bus 2 is connected to buses 1, 3,
4, and 5. We have injection measurement at bus 2. From here:
Sbus_lf = V .* conj(Ybus * V);
P_inj= real(Sbus_lf);
P_inj(2)
After this, I do not know how to
Hello,
I want to get PMU data for IEEE 30 bus system. According to my
understanding, this data is generated using the Simscape Electrical toolbox.
Consider that two PMUs are located at buses j and k where both buses are
adjacent to each other.
For every time stamp t, I want to obtain arriving PMU
Hello,
How can I use *fmincopf_solver in *MATPOWER* for *IEEE 30 bus system?
case_test = '30';
mpc = loadcase(['case',case_test]);
I do not understand how to get 'om'
* [results, success, raw] = fmincopf_solver(om, mpopt) *
Thank you in advance,
Naime
Hello,
I run the below code
[MVAbase, bus, gen, branch, success]=runpf('case14');
and I got this error:
MATPOWER Version 6.0, 16-Dec-2016 -- AC Power Flow (Newton)
Operator '-' is not supported for operands of type 'function_handle'.
Error in newtonpf (line 52)
mis = V .* conj(Ybus * V) - Sbus;
else
Q_inj_max(i,:) = Q_bus(i,:)-0.05*Q_bus(i,:);
Q_inj_min(i,:) = Q_bus(i,:)+0.05*Q_bus(i,:);
end
end
On Thu, Oct 13, 2022 at 1:33 AM Ray Daniel Zimmerman
wrote:
> I’m afraid I don’t understand what you are trying to calculate.
>
> Ray
>
>
> On Oct 11, 2022, a
(Ybus * V);
> Sf = V .* conj(Yf * V);
> St = V .* conj(Yt * V);
>
> Here Sbus are the bus injections into the branches (equal to the
> generation – load), and Sf and St are the branch power flows at the *from*
> and * to* ends, respectively.
>
> Ray
>
>
>
>
=branch(:,QT)/baseMVA;
P_inj_max = Pfe+0.05*Pfe;
P_inj_min = Pfe-0.05*Pfe;
Q_inj_max = Qfe+0.05*Qfe;
Q_inj_min = Qfe-0.05*Qfe;
Is this way correct?
Kind regards,
Naime
On Tue, Oct 11, 2022 at 1:45 PM naime ahmadi wrote:
> Thank you for your time
>
> On Tue, Oct 11, 2022 at 7:14 AM Carlos
Hello,
I have a question regarding equations (2) and (3) in the attached picture.
I am not sure equation (3) is correct. But *sj* is the power injection and
*f{i,j}* is the power flow.
I am a littel bit confused, how should I calculate y*{j,k} and *Yjk* in
Matpower?
According to the attached
I have forgotten to attach the picture.
On Mon, Oct 10, 2022 at 3:40 PM naime ahmadi wrote:
> Hello,
> I have a question regarding equations (2) and (3) in the attached picture.
> I am not sure equation (3) is correct. But *sj* is the power injection
> and *f{i,j}* is the power
Thank you
On Wed, Nov 16, 2022 at 4:32 AM Carlos E Murillo-Sanchez <
ce.murillosanc...@gmail.com> wrote:
> naime ahmadi wrote:
>
> Hello,
> I want to calculate shunt capacitance (*c*) in the attached figure.
> In the branch data of MATPOWER we have *b*.
> Is it correct
Hello,
I want to calculate shunt capacitance (*c*) in the attached figure.
In the branch data of MATPOWER we have *b*.
Is it correct to use this formulation to calculate *c*
*c=b/(2*pi*f)*
Kind regards,
Naime
Hello,
How can I calculate the nominal value of active and reactive power
injection in MATPOWER?
Kind regards,
Naime
Hello,
I want to add PV to bus 5 in the 14 bus system.
I used the code below
pv_bus = 5;
load_5 = sum(mpc.bus(mpc.bus(:,1)==pv_bus, 3));
pv_capacity = load_5 * 1.2; mpc.bus(mpc.bus(:,1)==pv_bus, 3) = 0;
in the generator data, I added this row
5 0 0 0 0 1 100 1 pv_capacity 0 0 0 0 0 0 0 0 0 0 0
you
> remove ";" at the end of each line, the command window will change the bass
> from type 1 to type 2. But in the same way, the voltage and power values
> should be changed according to the problem.
>
> On Sun, Feb 26, 2023 at 7:18 PM naime ahmadi
> wrote:
>
>>
40 matches
Mail list logo
