关于matlab潮流程序数据导入

xiaotiejiang523

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

milesnail

等待高人解答

xiaotiejiang523

data=zeros(1,4)
data=load('d:\MATLAB7\work\data.txt')
" d: I& K5 z5 r5 Bdisp('节点数：')
. I; A) T0 `6 f& an=data(1)
, s) H$ N2 j3 @) G+ ]' _disp('支路数：')
nl=data(2)
% J/ Z* P2 _4 Vdisp('平衡节点编号：')
) P4 _) }! T4 H( |isb=data(3)
disp('误差精度：')
pr=data(4)
disp('由支路参数形成的矩阵：')
B1=load('d:\MATLAB7\work\B1data.txt')
disp('各节点参数形成的矩阵：')
7 B3 y4 F* @+ ^) eB2=load('d:\MATLAB7\work\B2data.txt')
Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);
5 L. Q& J- l; n, b) R  QO=zeros(1,n);S1=zeros(nl);
1 d9 _$ l( s' N, }; F9 U# xfor i=1:nl
. q6 T( x& `' C. l    if B1(i,6)==0;
! D; S* _! W* u0 D% c8 ^         p=B1(i,1);q=B1(i,2);
     else p=B1(i,2);q=B1(i,1);
/ U) b2 ~' v+ a/ ~8 F: _" W0 F     end
     Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
7 E* C* N' r0 V0 u/ U     Y(q,p)=Y(p,q);
  Y4 E7 c3 _8 e8 l! _) ]     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
7 \+ ^6 l) A  `+ N4 N%求导纳矩阵
disp('导纳矩阵Y=');
disp(Y);
这是我的程序，txt文件里都是矩阵形式的数据

xiaotiejiang523

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

sparkyan

牛人啊

just1do2it

多谢楼主

mutouzhang

学习了

小虎

真厉害。

csusttang2006

太厉害了！1

杰克上尉

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