matlab 解方程组 算多次
a是系数矩阵,i是之前得到的坐标41对 我加了for循环 可最终只返还了我一组解
呃 还不是i=1的解.该怎么编.求c的41个值
a=[A1(i,1)*(-0.0913)-0.3280 A1(i,1)*0.0228+6.7123 A1(i,1)*(-0.0060)+1.2329;
A1(i,2)*(-0.0913)+3.3811 A1(i,2)*(-0.0228)+4.7034 A1(i,2)*(-0.0060)+5.2896;
B1(i,1)*0.0039-12.8350 B1(i,1)*0.0087+3.5069 B1(i,1)*(-0.0021)+0.2487];
b=[197.0020-A1(i,1);197.0020-A1(i,2);139.4985-B1(i,1)];
c=inv(a)*b
A1=load('sift1.txt')
B1=load('sift2.txt')
fid=fopen('reslut.txt','w')
[m n]=size(A1)
c=zeros(m,3)
for i=1:m;
a=[A1(i,1)*(-0.0913)-0.3280 A1(i,1)*0.0228+6.7123 A1(i,1)*(-0.0060)+1.2329;
A1(i,2)*(-0.0913)+3.3811 A1(i,2)*(-0.0228)+4.7034 A1(i,2)*(-0.0060)+5.2896;
B1(i,1)*0.0039-12.8350 B1(i,1)*0.0087+3.5069 B1(i,1)*(-0.0021)+0.2487];
for i=1:m;
b=[197.0020-A1(i,1);197.0020-A1(i,2);139.4985-B1(i,1)];
c1=inv(a)*b;
c(i,:)=c1;
end
fprintf(fid,'%f %frn',c);
fclose(fid)
end
全貌 呃