MATLAB求解方程组

来源:百度知道 编辑:UC知道 时间:2024/06/06 13:07:34
x1=6474.6;
y1=1816.08;
a=6378;
b=6370;

[x,y]=solve('((x1-x)*b^2*x+(y1-y)*a^2*y)/sqrt(((x1-x)^2+(y1-y)^2)*((b^2*x)^2+(a^2*y)^2))=cos(87*pi/180)','x^2/a^2+y^2/b^2=1');
x_y=vpa([x,y])
求解这方程组
这是我列举的程序。。我把原题发来帮我解决下..希望给我个程序..
求MATLAB高手解下面方程组
已知a= 6378 b=6370 x=6474.6 y= 1816
求解
(x1-x)*b^2*x+(y1-y)*a^2*y)/(((x1-x)^2+(y1-y)^2)((b^2*x)^2+(a^2*y)^2))^(1/2)=cos(87*pi/180)
x^2/a^2+y^2/b^2=1

你这都是啥玩意啊。。。语法有问题。

你定义好了几个参数,[x, y] = solve(....)就解决了啊。。。只要你把方程定义正确就行了。。。

我试了一下:

[x, y] = solve('((x1-x)*b^2*x+(y1-y)*a^2*y)/(((x1-x)^2+(y1-y)^2)*((b^2*x)^2+(a^2*y)^2))^(1/2)=cos(87*pi/180)','x^2/a^2+y^2/b^2=1');

这个语法应该是对的,但运行起来很慢。 你要仔细看一下方程的定义,能够简化的就简化,比如第一个分式,能不能先化简一下,这样算起来就比较快。

语法确实有点细节不够完善,但是也并不影响程序
这个方程组好像不大能解出来,我运行时,界面左下角一直显示“busy”中,没能解出来,这个方程太过于复杂了。
_________________________________-
还有一种很原始的方法,在一个图像上做出两条曲线,然后在图上点击交点,便会显示出X的数值