最近遇到了求二元二次方程组的问题,故写此文以记录解决方法。
首先创建如下两个文件,并输入对应的代码
1.myfun1
function q=myfun1(p)
x=p(1);
y=p(2);
q(1)=(x-a)/(y-b)-(x-c)/(y-d)-e;
q(2)=(x-A)/(y-B)-(x-C)/(y-D)-E;
end
q(1)和q(2)后面就是要求的方程组,我这里假设是这两个方程。
2.untitle4(这个文件名字随意,上面那个不能改)
x=fsolve('myfun1',[0.5,0.5],optimset('Display','off'));
x(1)
x(2)
完成这一步后效果如下(未展示untit4文件)
然后使用时
将a,b,c等等参数,输入到myfun中,就是说直接就a等字母改成值
类似这样
然后点击untitled4文件,注意是untitled4文件,点击运行,即可得到两个值,第一个是x,第二个是y
类似这样