关于matlab潮流程序数据导入
各位大侠,本人在编写牛拉潮流程序时,想导入数据,这样就不用一步步输入了,单看某个矩阵,赋值是对的,但是最后结果不对啊。请高人指点。谢谢。 等待高人解答 data=zeros(1,4)data=load('d:\MATLAB7\work\data.txt')
disp('节点数:')
n=data(1)
disp('支路数:')
nl=data(2)
disp('平衡节点编号:')
isb=data(3)
disp('误差精度:')
pr=data(4)
disp('由支路参数形成的矩阵:')
B1=load('d:\MATLAB7\work\B1data.txt')
disp('各节点参数形成的矩阵:')
B2=load('d:\MATLAB7\work\B2data.txt')
Y=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
if B1(i,6)==0;
p=B1(i,1);q=B1(i,2);
else p=B1(i,2);q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
Y(q,p)=Y(p,q);
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=');
disp(Y);
这是我的程序,txt文件里都是矩阵形式的数据 我知道了,我的数据里不只是数据,还有字母,也就是复数,用了dlmread导入就好了。 牛人啊 多谢楼主 学习了 真厉害。 太厉害了!1 楼主很强大啊,有人做过关于某种算法的无功优化没有,需要指点。
页:
[1]
2
