View License

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

Highlights from
Newton-Raphson Loadflow

Join the 15-year community celebration.

Play games and win prizes!

» Learn more

4.26
4.3 | 53 ratings Rate this file 179 Downloads (last 30 days) File Size: 7.54 KB File ID: #21059 Version: 1.5

Newton-Raphson Loadflow

by

Praviraj PG (view profile)

 

13 Aug 2008 (Updated )

These matlab m files are used to calculate bus voltages and angles using Newton Raphson iterative me

| Watch this File

File Information
Description

These matlab m files are used to calculate bus voltages and angles, power flows using Newton Raphson iterative method

Acknowledgements

This file inspired Ieee 30 Bus Newton Raphson Load Flow Analysis.

MATLAB release MATLAB 7 (R14)
Other requirements Matlab 6 or higher
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (107)
22 Aug 2016 ogochukwu Chuks-Ezike

can you please send a matlab code for newton raphson method with UPFC incorporated to my email- ogochukwuchuks@gmail.com. Thanks for your anticipated cooperation

Comment only
19 Aug 2016 Katarina Gajic

Hi, I am calculeting PTDFs for 30 bus system, both with DC and AC method and now I am stack in formation of transaction matrix, apparently for both of these methods they are the same size 41*870 but they are not suppposed (in theory) equally filled and I can't understand why, and I am wandering if PQ bus can also inject some power in the network?

Comment only
19 Aug 2016 madiha elnagar

hi....i need the IEEE 33-bus system line data and bus data... Ive been searching for them for a while and i can't find them....could anyone help...thanx

Comment only
25 Jun 2016 Saurav Raj

Thankyou sir.....

27 May 2016 BAVA Robert

Good job

Comment only
08 May 2016 ahmet tan  
04 May 2016 Berny Rox  
26 Apr 2016 amanda tadros

please can you send me matlab newton raphson power flow with UPFC
tadros.amanda@hotmail.com

Comment only
25 Apr 2016 pranay prakash

sir , i am getting function undefined error in this code .where should i define this

Comment only
22 Apr 2016 Florian Rousseau  
22 Apr 2016 PANKAJ UPADHYAY

Thank you it worked.
Thanks Erik for correction.

25 Mar 2016 Kamsundaram Sree Latha

wError in ==> nrlfppg at 160

loadflow(nbus,V,del,BMva); %

Calling Loadflow.m..
after running all i am getting this error
can anyone help please?

Comment only
25 Mar 2016 Kamsundaram Sree Latha

which program we need to run at first i am getting error if i am running busdat

Comment only
17 Mar 2016 Mohammad Anas Anees  
10 Mar 2016 tabia

tabia (view profile)

Error in loadflow.m
Attempted to access Vm(2); index out of bounds because numel(Vm)=1.

Error in loadflow (line 28)
Iij(p,q) = -(Vm(p) - Vm(q))*Y(p,q); % Y(m,n) = -y(m,n)..
How to correct this?

07 Mar 2016 t s

t s (view profile)

 
03 Mar 2016 mohammed Alsaraj

It work successfully.
Thank you very much Praviraj.
Thank you Erik for correction.

Comment only
16 Feb 2016 pratik chavadekar

I want load flow analysis program for 15 bus.

Comment only
08 Dec 2015 mitali chakravorty

i want the newton raphson code with L-index formulation for voltage stability analysis

Comment only
08 Dec 2015 amanpreet kaur

plz give 33 IEEE bus data and line data plzzzzzzz

Comment only
25 Nov 2015 Parminder Gill

function Y = ybusppg(num) % Returns Y
|
Error: Function definitions are not permitted in this context.

sir plz tell how to overcome this error

Comment only
20 Nov 2015 amir hossein aminian

thank you Thomas, thanks a looooooot.

Comment only
24 Oct 2015 Zahid Abbas  
16 Sep 2015 Thomas

Thomas (view profile)

For anyone getting 'Function is not defined for sparse' type errors:

You need to converter spare matrix element to full matrix so 'fprintf' can display it. eg:

fprintf('%8.3f', full(Pij(p,q)))

A similar update to the code is need for all sparse matrix elements being displayed by 'fprintf'.

Alternatively, you can go further back in the code and change the sparse matrix to full versions, thought this is at expense of increased memory requirement.

Comment only
23 Aug 2015 Aakanksha kumawat

hey..i m working on 33 bus test system and i want the ideal bus data...plz help me

Comment only
27 Jul 2015 Naem

Naem (view profile)

thus is the error that I've faced with:

??? Error using ==> fprintf

Function is not defined for sparse

inputs.

Error in ==> loadflow at 97

fprintf('%8.3f', Pij(p,q));

Error in ==> nrlfppg at 160

loadflow(nbus,V,del,BMva); %

Calling Loadflow.m..

can anyone help please?

Comment only
12 Jun 2015 shijith pt

sir after running the program it showing the following errors

Error in ==> loadflow at 2
Y = ybusppg(nb); % Calling Ybus progra
please help me correct it...

shijithpala@gmail.com

11 Jun 2015 wajahat hussain

hi everybody
can anybody send me code for generating
gui for 30 bus bar ieee through newton raphson method at
s.ahil@kimo.com

Comment only
29 May 2015 juliana Klas

Very nice, congrats and many thanks!

12 May 2015 B V NAGESWARUDU

i am doing project on the congestion management power system using differential evolution algorithm.the equality constraints are handled by NR method,but the program showing the error as follows
Line FLow and Losses
-------------------------------------------------------------------------------------
|From|To | P | Q | From| To | P | Q | Line Loss |
|Bus |Bus| MW | MVar | Bus | Bus| MW | MVar | MW | MVar |
so please send me the matlab coding,which is helpful for my project.send coding to this mail nageshbv13@gmail.com

08 May 2015 Zaid Yousaf

Hello everyone;
Anyone who have experience to work on "Power System Improvement using UPFC" (Newton Raphson algorithm used in it and MATLAB used as a Tool).
If we take 3 bus system and find the power flow using Newton Raphson Method, and again take this system by improve power system stability by using UPFC with same algorithm (Newton Raphson Method) used.
Kindly guide me about following points;
a) Parameters used for UPFC in this?
b) What important things should be kept in mind while load flow calculations using UPFC??
Please guide me guys.....urgent please

03 May 2015 Ahmad Farhan

hello sir... after executing the program the command shows the following error:

-------------------------------------------------------------------------------------
Line FLow and Losses
-------------------------------------------------------------------------------------
|From|To | P | Q | From| To | P | Q | Line Loss |
|Bus |Bus| MW | MVar | Bus | Bus| MW | MVar | MW | MVar |
-------------------------------------------------------------------------------------
1 2Error using fprintf
Function is not defined for sparse inputs.

Error in loadflow (line 95)
fprintf('%4g', p); fprintf('%4g', q); fprintf(' %8.3f', Pij(p,q)); fprintf(' %8.3f', Qij(p,q));

can you help me to solve this problem?

Comment only
17 Apr 2015 killu bhashna

good work mam...
plz.. can any 1 help me how to install dg (distributed generation) in load flow solution by NR method ..sonia.bhashu@gmail.com

Comment only
15 Apr 2015 bhawana

??? Input argument "nb" is undefined.

Error in ==> loadflow at 2
Y = ybusppg(nb); % Calling Ybus program..
please help me..

Comment only
12 Apr 2015 CHAN ROY

Can you send me matlab code?
email: benmilous@hotmail.com

Comment only
06 Apr 2015 Bhaskar Pulimamidi

can you please send me matlab code for newton raphson method with UPFC incorporated.please kindly help me at bhasker.31121993@gmail.com

Comment only
18 Mar 2015 ajay konathala

can i have the newton raphson program that has been successfully compiled,i am using STATCOM in the load flow analysis for a IEEE 30 bus system to reduce reactive power,i need the code of newton raphson method.
if anyone has it kindly reply to konathalaajay57@gmail.com

18 Mar 2015 ajay konathala

can i have the newton raphson program that has been successfully compiled,i am using STATCOM in the load flow analysis for a IEEE 30 bus system to reduce reactive power,i need the code of newton raphson method.

11 Mar 2015 nur ashiqin sabarudin

Hello,can you give me matlab code for newton raphson load flow with UPFC at pieceyo_xoxo@yahoo.com

Tq 4 the reply.

Comment only
22 Feb 2015 Ahmed Hussain

can you please give me code for nr method using UPFC

at
AHMEDSSUET@hotmail.com

05 Feb 2015 boat bodin

thankyou

12 Jan 2015 Vanka DEGA RAJAJI

your work is very nice , but am in need of calculation of switching angles alpha1,alpha2 etc in 5 level multilevel inverters using newton- raphson method. please send me this file

Comment only
11 Jan 2015 mohamed

very nice work that help me
thank you

06 Jan 2015 Kumar Mahesh

hello dear i need Matlab IEEE 33 code, if anyone can help

Comment only
15 Dec 2014 Aakash

Aakash (view profile)

Can i get the matlab code for incorporating upfc in newton raphson algorithm for ieee 30 bus system...
plz mail me it at aakash274@gmail.com

Comment only
23 Nov 2014 Ishan Purohit

Hey can anybody tell em how to copy the complete Newton Raphson Loadflow Analysis and line losses data in excel sheet

Comment only
21 Nov 2014 khalaf

khalaf (view profile)

I have changed the 33-bus code (mesh system) to 34-bus(radial system) and i have added the 34th bus data and changed the line data.

when i run, it works but never close to converge. it means the tolerance is very high in kilos and the results are unstable( goes from zeros to kilos and never stop).

anybody can help. is it ok to run the code for radial system is staid of mesh system? please help me if you can

15 Nov 2014 sandeep

hello mr. hossein. I read your paper on distributed generation. I am also working on same topic is it possible for you to help me?

Comment only
14 Nov 2014 Erik

Erik (view profile)

Hanumant, change lines 44 and 45 of loadflow.m to read:

Pij = full(real(Sij));
Qij = full(imag(Sij));

Comment only
09 Nov 2014 mohammad hossein

hi my freinds...
i need to your help for using matlab to run newton raphson load flow by using genetic algorithm...
my problem is in 30 & 33 bus IEEE ...!!
thanks for your participation for solving the problem, with me...!!
god bless yo...

07 Nov 2014 J Gunda

Sorry I just had a rough look and did not checked the program completely. I understood you have rearranged the equations while calculating Jacobian elements. Please ignore my last comment.

Comment only
06 Nov 2014 J Gunda

I think there was a mistake in calculation of Jacobian elements in your program. When m=n, if you do the differentiation you will get 2*Vm but you have used only Vm. Anyway good job....

Comment only
03 Sep 2014 Zaid Yousaf

I want to perform
"LOAD FLOW STUDIES WITH UPFC POWER
INJECTION MODEL FOR 30 BUSES",
Kindly help me coding(MATLAB simulation).
Can i get complete code for 30 buses.

PLEASE HELP ME .
SINCERE REGARDS;
Zaid

03 Sep 2014 Zaid Yousaf

how can i watch PRAVI posts......i have already made account....kindly guide me ...i am new user

Comment only
02 Sep 2014 Mohsen

Mohsen (view profile)

Thank you very much

01 Sep 2014 mourad

mourad (view profile)

How to insert a photovoltaic source and estimat its impact in the distribution network ?

04 Aug 2014 Merah ahmede

Nice Joooob.

09 Apr 2014 Hanumant Sarde

hello sir after executing the program the command shows the following error

1 2Error using fprintf
Function is not defined for sparse inputs.

Error in loadflow (line 95)
fprintf('%4g', p); fprintf('%4g', q); fprintf(' %8.3f', Pij(p,q)); fprintf(' %8.3f',
Qij(p,q));

Error in nrlfppg (line 160)
loadflow(nbus,V,del,BMva); % Calling Loadflow.m..

can u please make me understand why the error occured

11 Mar 2014 pruthvi

hello sir after executing the program the command shows the following error

1 2Error using fprintf
Function is not defined for sparse inputs.

Error in loadflow (line 95)
fprintf('%4g', p); fprintf('%4g', q); fprintf(' %8.3f', Pij(p,q)); fprintf(' %8.3f',
Qij(p,q));

Error in nrlfppg (line 160)
loadflow(nbus,V,del,BMva); % Calling Loadflow.m..

can u please make me understand why the error occured

24 Jan 2014 Rohit kumar Bharadia  
23 Jan 2014 SRINIVAS GAMPA  
23 Dec 2013 Priyadharshini M

Hai im in need of Newton-Raphson load flow analysis Matlab code for contingency analysis in power system..I want to apply the line outages,load outages, generation outages randomly..i need your help..please help me..

21 Sep 2013 mangeshkumar rupnar  
12 Sep 2013 azri zawawi

ello..i'm azri...I want to do power system load flow analysis using matlab by fast decoupled method....can u help me....???please help me...

Comment only
17 Jun 2013 Mahdi

Mahdi (view profile)

Hi every body
I want to do load flow analysis for IEEE 13 bus (unbalance distribution network).
How can I perform it?
Please help me.
Thanks a lot.

28 Apr 2013 vit university  
27 Apr 2013 Waleed

Waleed (view profile)

Hi all,

I am running this code for fixed real powers at generator buses but it changes the values after execution, from those values that I have input. PV (generator) buses are supposed to have a fixed P and V, and Q and angle get calculated from them by Newton Raphson. Now I don't know what to do. Kindly look into my issue urgently, I have a project that is due in a day.

Comment only
06 Apr 2013 Meshal AM  
06 Mar 2013 JMI ENGG COLLEGE

In NRLF simulation, I got most of the results as NaN. How can I solve the problem?
Does anyone has solution to this problem?

Comment only
01 Jan 2013 arman ghaderi

do you consider the limit of reactive power of generators?

Comment only
01 Jan 2013 umakanth ramineedi

Error in ==> at 81
G = real(Y); please help me i you help me it really use full to my project

Comment only
13 Nov 2012 JAIME MORONI

the solution for the fprintf...
fprintf(' %8.3f', Pij(p,q)); fprintf('
%8.3f', Qij(p,q));....

you have to change fprintf('
%8.3f', Qij(p,q)) by fprintf('
%8.3f', full(Qij(p,q))) ....

you must do the same for Pij Qij any vector...

i have a cuestion the results are a little diferents. that a flow

http://catedrasempresa.esi.us.es/endesared/documentos/Tesis_Marano.pdf

30 Oct 2012 Mamta

Mamta (view profile)

 
30 Oct 2012 Mamta

Mamta (view profile)

 
01 Aug 2012 Shahariar Kabir

Hi, I am having the the same problem as
Munibullah. Please help us urgently.

Comment only
12 Jul 2012 Munibullah

hey im having the same problem.

1 2??? Error using ==> fprintf
Function is not defined for sparse inputs.

Error in ==> loadflow at 96
fprintf('%4g', p); fprintf('%4g', q);
fprintf(' %8.3f', Pij(p,q)); fprintf('
%8.3f', Qij(p,q));

Error in ==> nrlfppg at 160
loadflow(nbus,V,del,BMva); %
Calling Loadflow.m..

HOw to fix this? i need to get the results its very imp.

Comment only
08 Jul 2012 Umar

Umar (view profile)

How to determine line current flow, the commands are given in loadflow.m but this m file does not execute. Please Help!!!!!!!!

08 Jun 2012 mary

mary (view profile)

hi, thanks for it,
i need to calculate jacobian matrix 1 time in distribution system but i dont know how do it, can u help me?

18 May 2012 vikram singh

thanks

16 Apr 2012 Houssemus

@leonardo: yes it does!!

19 Feb 2012 sumit verma

v.helpful program ......thanx pravi...

10 Feb 2012 jayendra javi  
19 Jan 2012 leonardo

I wonder if this script calculates the number of iterations performed, and if they do not calculate whether there is any way to calculate. Thanks in advance.

Comment only
19 Jan 2012 leonardo

first time thank for your attention.

Comment only
19 Jan 2012 leonardo

hi, i'm have some problem with this script, you could help me?
Error in ==> Untitled at 4
Y = ybusppg(nbus); % Calling
ybusppg.m to get Y-Bus Matrix..

Comment only
01 Nov 2011 meysam shamshiri

Thanks

Comment only
23 Oct 2011 reygen oraiz

hi..i have the same mistake with omar and as you said, make this "Sij = sparse(Sij);" be commented at line 60..(its line 43) but an error occur again, it says:

??? Input argument "nb" is undefined.

Error in ==> loadflow at 5
Y = ybusppg(nb);

---------please help------------
i really need to have this right! thanks

Comment only
13 Aug 2011 Fabián

Why do you add .01 to voltage magnitude when a pv bus is out of mvar limits??

Comment only
09 Jul 2011 balaji

balaji (view profile)

When i run for 69 bus i get the error ""??? Attempted to access V(66); index out of bounds because numel(V)=65."
can anyone help

06 Jul 2011 Mehdi Davoudi

it works properly but why the results are far from the real ones even with too many iterations?

15 Jun 2011 rasoul bayani

thanks .very much

Comment only
11 May 2011 Ali

Ali (view profile)

 
23 Feb 2011 gmrit

gmrit (view profile)

PLEASE SEND NR METHOD LOADFLOW STUDY WITH UPFC. I LIKE THIS PROGRMME SOMUCH.IT IS FANTASTICK.PLEASE SEND UPFC PROGRAM TO ramesh.eee44@gmail.com

08 Oct 2010 Praviraj PG

Praviraj PG (view profile)

@Omar: I guess, this is error is bcs of the different versions of matlab. Keep the line 60 (loadflow.m) commented.

%Sij = sparse(Sij);

And try again.

Comment only
17 May 2010 Omar

Omar (view profile)

i have one mistake in the display Line Flow and losses:
error in loadflow at 95(Function is not defined for sparse inputs.) and nrlfppg at 160 , i need help

Comment only
23 Mar 2010 senthil

please provide me the FAST DECOUPLED load flow method

Comment only
23 Mar 2010 senthil

thank you to provide this

Comment only
25 May 2009 Erdal Bizkevelci  
11 Dec 2008 Raju C V

Simple steps,Good program for novice programmers. Thanks

22 Nov 2008 Husam Aldahiyat  
19 Oct 2008 prasert limsakun

very good for program loadflow

27 Sep 2008 S F

best

24 Sep 2008 Randhir Soni

if u ve made this prog with including Tap changing transformer then it will be better for someone to have complete practical idea about load flow problem

Comment only
12 Sep 2008 sougata mandal

we need the programme such a simple steps,which is better 2 understand for the new learners. thank u........

04 Sep 2008 Jayanti Bhattacharjee

needs excellent speed

Comment only
18 Aug 2008 saloman danaraj

This program does not include the effect of tap changing transformers on powerflow which very much essential for modelling the problem.
The IEEE 30 bus system has 3 transformers .The author simply ignores that.

Comment only
14 Aug 2008 B. Roossien

- C-type of programming instead of Matlab-style (vectorised) programming.

- Script files instead of functions

- Very little documentation is provided making this contribution only useful for the author.

Author please remove this submission as it has no value to anyone but you. Why submit a file that no one else can understand?

Updates
21 Nov 2008 1.1

Updated for Tap changing transformers and Q-limit violations

23 Nov 2008 1.2

Updated busdata

23 Nov 2008 1.3

Updated busdata

03 Feb 2009 1.4

corrected the error with loadflow.m

08 May 2009 1.5

Included files for calculating power flows

Contact us