关于matlab潮流程序数据导入

xiaotiejiang523 · 发表于 2010-12-8 12:04:43

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

您需要登录才可以下载或查看，没有账号？立即加入

×

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

milesnail · 发表于 2010-12-8 13:09:41

等待高人解答

xiaotiejiang523 · 发表于 2010-12-8 13:15:56

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文件里都是矩阵形式的数据

xiaotiejiang523 · 发表于 2010-12-8 13:25:43

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

sparkyan · 发表于 2011-5-5 17:56:50

牛人啊

just1do2it · 发表于 2011-5-8 10:46:41

多谢楼主

mutouzhang · 发表于 2011-5-9 20:24:57

学习了

小虎 · 发表于 2011-5-30 02:55:38

真厉害。

csusttang2006 · 发表于 2011-5-31 10:24:59

太厉害了！1

杰克上尉 · 发表于 2011-8-11 10:08:10

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

账号		自动登录	找回密码
密码			立即加入

关于matlab潮流程序数据导入

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

相关帖子