Question

matlab fsolve函数应用

我想求出x^2+y^2+30*y=175,x^2+y^2-90*x=-800中x,y,的解请问matlab的程序应该怎么呀？

Answer 1

fsolve是MATLAB软件的一条函数。它采用最小二乘法来求解非线性方程。它的一般求解方式为:X=FSOLVE(FUN,X0,OPTIONS)。
matlab fsolve函数的应用举例如下，

求下列非线性方程组在(0.5,0.5) 附近的数值解。
(1) 建立函数文件myfun.m。
function q=myfun(p)
x=p(1);
y=p(2);
q(1)=x-0.6*sin(x)-0.3*cos(y);
q(2)=y-0.6*cos(x)+0.3*sin(y);
(2) 在给定的初值x0=0.5,y0=0.5下，调用fsolve函数求方程的根。
x=fsolve('myfun',[0.5,0.5],optimset('Display','off'))
x =
0.6354
0.3734
将求得的解代回原方程，可以检验结果是否正确，命令如下:
q=myfun(x)
q =
1.0e-009 *
0.2375 0.2957
可见得到了较高精度的结果。

Answer 2

function h=afang(xy) % 保存为 afang.m 文件
x=xy(1);y=xy(2);
h=[x^2+y^2+30*y-175;
x^2+y^2-90*x+800];
---------------------------------------------------------
[x,fval,exitflag]=fsolve('afang',[1;1])
% x 为所求的解
%fval为解代入方程的验证结果
% exitflag=1表示函数收敛于解；反之，=0

追问

是哪里出错了呢？

追答

function h=afang(xy) % 保存为 afang.m 文件
x=xy(1);y=xy(2);
h=[x^2+y^2+30*y-175;
x^2+y^2-90*x+800];