电力系统潮流程序

电力系统的潮流程序，经典方法，高斯法，牛拉法，PQ法

format long;
u1=1.0;u2=0.9;e=0.0001;
Fun=[11];n=1;
while （max（Fun）>=e | min（Fun）<=-e）
U0=[u1;u2];
J=[2+u2 u1;-u22-u1];%J怎么计算?表达式手算;
Jni=pinv（J）;%求J的逆;
Fun=[2*u1+u1*u2-3.3;2*u2-u1*u2-0.9];%函数
Ux=U0-Jni*Fun;%第x次迭代后的值，公式10.26，PAGE340
u1=Ux（11）;u2=Ux（21）;
D=‘迭代次数:‘;
disp（D）;disp（n）;
E=‘值：‘;
disp（E）;
disp（Ux）;
n=n+1;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1485  2008-12-03 15:06  no3.m

     文件        370  2008-12-03 14:55  no1.m

     文件       1539  2008-12-03 14:42  no2.m

----------- ---------  ---------- -----  ----

                 3394                    3