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

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图像来

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) 结果：
再问： 通过二维画图 可以近似得到下图中的曲线我感觉你的程序得修改一下吧 我没怎么看懂