Question

matlab 对矩阵各元素进行二元二次线性方程组的求解，用到for循环，但一直报错出不来结果，why？

已知ar br Sr Hr为3*3矩阵，对矩阵中的每一个元素有一组二元二次线性方程组的对应关系，要求出aAr和bAr.程序代码如下
ar=[2,5,8;
-5,-6,-9;
-9,-7,4];
br=[8 9 2;
-2 2 -6;
9 -44 -15];
Hr=[30 31 32;
50 14 28;
71 48 39];
Sr=[10 25 36;
45 88 15;
26 88 95];
Si=size(ar);
aAr=zeros(Si);
bAr=zeros(Si);
for i=1:Si(2)
for j=1:Si(1)
if ar(i,j)<0 & abs(Hr(i,j)-41.6)<=15
[aAr(i,j),bAr(i,j)]=solve('(aAr(i,j)-125)^2-(bAr(i,j)-172)^2=125^2+172^2','aAr(i,j)^2+bAr(i,j)^2=Sr(i,j)^2*(Hr(i,j)-41.6)^2/（15^2)');
aAr(i,j)=-abs(aAr(i,j));
bAr(i,j)=abs(bAr(i,j));
end
end
end
最后需要的是矩阵aAr和bAr中的每个元素的值。麻烦高手指正一下问题在哪儿，一直调试报错aAr和bAr都还是零矩阵，无法对元素赋值。拜托了比较着急！

Answer 1

原来的方程有问题，这个OK，不过还是要把变量名进行修改！
[aAr(i,j),bAr(i,j)]=solve('(aAr(i,j)-125)^2-(bAr(i,j)-172)^2=125^2+172^2','aAr(i,j)^2+bAr(i,j)^2-Sr(i,j)^2*(Hr(i,j)-41.6)^2/(15^2)=0');

追问

你好，我照你上面的进行了修改还是报错 警告为2 equations in 1 variables 和Subscripted assignment dimension mismatch.，且solve这一行还是有错，问题在哪儿呢？你们说的这个字符串处理是什么意思咧？aAr和bAr是先申请为3*3的零矩阵，通过求解方差组再对矩阵中的每个元素进行赋值，但一直都赋值不上，麻烦你可以再帮我看看么?非常谢谢

追答

Subscripted assignment dimension 这个是指维数不匹配，警告的话说明你某个方程没有被列入，但比较奇怪的是为何变为1个变量了。这个觉得你还是需要把你的问题描述记录清楚，不然别人没法解决的。

追问

追问限制字数 我一直提交不了 所以发了私信，先谢谢你啦

Answer 2

[aAr(i,j),bAr(i,j)]=solve('(aAr(i,j)-125)^2(bAr(i,j)-172)^2=125^2+172^2','aAr(i,j)^2+bAr(i,j)^2=Sr(i,j)^2*(Hr(i,j)-41.6)^2/（15^2)'); 中 '(aAr(i,j)-125)^2(bAr(i,j)-172)^2=125^2+172^2' aAr(i,j)只是当做一个字符串处理了，没有代入值

Answer 3

能，两两组合，取交集