关于matlab潮流程序数据导入

xiaotiejiang523

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

milesnail

等待高人解答

xiaotiejiang523

data=zeros(1,4)
/ n& y1 ]: T5 z+ G4 s3 \$ q6 S* Kdata=load('d:\MATLAB7\work\data.txt')
4 V! Z/ a/ g' `& Q7 [disp('节点数：')
n=data(1)
) P8 S6 V3 z2 d8 ?9 r9 ydisp('支路数：')
nl=data(2)
disp('平衡节点编号：')
isb=data(3)
7 V0 H1 R9 R4 }$ D( {& l" Fdisp('误差精度：')
/ Q; ?; }; C4 j; Q& `* W3 d, Upr=data(4)
disp('由支路参数形成的矩阵：')
; M3 m( @7 H4 \$ D6 f/ {2 OB1=load('d:\MATLAB7\work\B1data.txt')
disp('各节点参数形成的矩阵：')
# w& ]' }  s, e2 g+ `B2=load('d:\MATLAB7\work\B2data.txt')
! B) H3 y# G: U0 cY=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);
) P: {$ x4 E; c' a' eO=zeros(1,n);S1=zeros(nl);
for i=1:nl
    if B1(i,6)==0;
6 G. F% F8 }/ c; N, H3 E$ ]         p=B1(i,1);q=B1(i,2);
5 {4 t1 y% ?1 I* O+ H- b3 i! u     else p=B1(i,2);q=B1(i,1);
     end
& _2 n/ X4 Q! D5 x. G3 n8 }+ n$ Q, s     Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
" g' h" B& u. ?9 F+ u     Y(q,p)=Y(p,q);
* Q8 `- E$ S/ q0 p" m) p     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
# u2 E( D, C  C9 r: p$ b%求导纳矩阵
disp('导纳矩阵Y=');
disp(Y);
这是我的程序，txt文件里都是矩阵形式的数据

xiaotiejiang523

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

sparkyan

牛人啊

just1do2it

多谢楼主

mutouzhang

学习了

小虎

真厉害。

csusttang2006

太厉害了！1

杰克上尉

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