关于matlab潮流程序数据导入

xiaotiejiang523

各位大侠，本人在编写牛拉潮流程序时，想导入数据，这样就不用一步步输入了，单看某个矩阵，赋值是对的，但是最后结果不对啊。请高人指点。谢谢。

milesnail

等待高人解答

xiaotiejiang523

data=zeros(1,4)
# l% \; e/ [' x5 O+ V9 e; pdata=load('d:\MATLAB7\work\data.txt')
+ q6 o1 W' u2 Ldisp('节点数：')
n=data(1)
9 L+ O* V& ?9 P7 Tdisp('支路数：')
4 B( \; c/ e  z4 r: Pnl=data(2)
disp('平衡节点编号：')
1 N9 O$ \- r' d: o. ~+ \isb=data(3)
disp('误差精度：')
' O7 K% a8 b7 K, Z# npr=data(4)
% d$ S2 ]5 }4 f& W  Tdisp('由支路参数形成的矩阵：')
B1=load('d:\MATLAB7\work\B1data.txt')
. L9 V" K$ e% b) q! wdisp('各节点参数形成的矩阵：')
B2=load('d:\MATLAB7\work\B2data.txt')
) k- @9 D$ o- \! ?) \5 J6 q; HY=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);
O=zeros(1,n);S1=zeros(nl);
for i=1:nl
+ `  _$ z4 Y4 `7 n. `: M    if B1(i,6)==0;
" V% |7 d( J5 ]# j7 v, x         p=B1(i,1);q=B1(i,2);
% c" ?: w) l* t0 J" E( e3 D     else p=B1(i,2);q=B1(i,1);
     end
* U/ E2 [9 R) T, L* A7 _0 z& S; X     Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
     Y(q,p)=Y(p,q);
5 ?/ y( [5 w6 }4 w     Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;
     Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
end
%求导纳矩阵
disp('导纳矩阵Y=');
2 O8 {! `, b; c" ldisp(Y);
/ |5 x% p9 u, I2 N这是我的程序，txt文件里都是矩阵形式的数据

xiaotiejiang523

我知道了，我的数据里不只是数据，还有字母，也就是复数，用了dlmread导入就好了。

sparkyan

牛人啊

just1do2it

多谢楼主

mutouzhang

学习了

小虎

真厉害。

csusttang2006

太厉害了！1

杰克上尉

楼主很强大啊，有人做过关于某种算法的无功优化没有，需要指点。