|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
请教各位高手: 我在用matlab求线性方程组的时候,方程组的解输出到文件中,为什么输出不完全(应该有225个数据,文件却只有216个)?
以下是程序代码
clear
fid=fopen('a.data','r');
a=fscanf(fid,'%f',[225,225]);
r=rank(a)
s=svd(a)
[l,a1,m]=svd(a);
co=cond(a)
%a=a-0.01*eye(size(a));
fid3=fopen('a0.data','w');
fprintf(fid3,'%f\n',a);
n=det(a)
c=pinv(a);
plot(c)
d=norm(a*c-eye(size(a)))
%d=norm(a*c*a-a)
e=eig(a);
fid4=fopen('b.data','r');
b=fscanf(fid4,'%f',[225]);
x=c*b;
fid2=fopen('eig.data','w');
fprintf(fid2,'%f\n',e);
fidd=fopen('ia.data','w');
fprintf(fidd,'%18.15f\n',c);
fid5=fopen('x.data','w');
fprintf(fid5,'%18.15f\n',x
fclose(fid);
fclose(fid4); |
|