Question

matlab求解下列方程的解，画出（x,y）的曲线图

(2*x-1200*sin(t)-400)^2+(2*y+1200*cos(t)-1200-sqrt(1280000))^2=1440000;
(2*x-2400*sin(t))^2+(2*y+2400*cos(t)-2400-2*sqrt(1280000))^2=160000+(l-sqrt(1280000))^2;
其中给t一个范围t=-20/180*pi:pi/1800:20/180*pi; 能否画出xy图像来

Answer 1

t在这个范围内只有一个解。程序如下
t=-20/180*pi:pi/1800:20/180*pi;
f=@(x)([(2*x(1)-1200*sin(t)-400).^2+(2*x(2)+1200*cos(t)-1200-sqrt(1280000)).^2-1440000;(2*x(1)-2400*sin(t)).^2+(2*x(2)+2400*cos(t)-2400-2*sqrt(1280000)).^2-160000-(1-sqrt(1280000))^2]);
x1=fsolve(f,[0 0])；
x=x1(:,1)
y=x1(:,2)
结果：
x =
-294.7252

y =
727.1720

追问

给t是一个范围（一个矩阵） 那么得到的xy也应有对应的矩阵啊 怎么会只有一个解呢？xy不应可以会出一条曲线来吗？求解

追答

程序改下：

x=[];
y=[];
i=1;
for t=-20/180*pi:pi/1800:20/180*pi
   
f=@(x)([(2*x(1)-1200*sin(t)-400).^2+(2*x(2)+1200*cos(t)-1200-sqrt(1280000)).^2-1440000;(2*x(1)-2400*sin(t)).^2+(2*x(2)+2400*cos(t)-2400-2*sqrt(1280000)).^2-160000-(1-sqrt(1280000))^2]);
x1=fsolve(f,[0 0]);
x(:,i)=x1(:,1);
y(:,i)=x1(:,2);
i=i+1;
end
plot(x,y)

 

结果：

追问

通过二维画图  可以近似得到下图中的曲线

我感觉你的程序得修改一下吧  我没怎么看懂