答案:clear; close; t=0:pi/100:2*pi; x1=2 sqrt(5)*cos(t); y1=3-2*x1 sqrt(5)*sin(t); x2=3 sqrt(2)*cos(t); y2=6*sin(t); plot(x1,y1,x2,y2); grid on; %作图发现4个解的大致位置,然后分别求解 y1=fsolve('[(x(1)-2)^2 (x(2)-3 2*x(1))^2-5,2*(x(1)-3)^2 (x(2)/3)^2-4]',[1.5,2]) y2=fsolve('[(x(1)-2)^2 (x(2)-3 2*x(1))^2-5,2*(x(1)-3)^2 (x(2)/3)^2-4]',[1.8,-2]) y3=fsolve('[(x(1)-2)^2 (x(2)-3 2*x(1))^2-5,2*(x(1)-3)^2 (x(2)/3)^2-4]',[3.5,-5]) y4=fsolve('[(x(1)-2)^2 (x(2)-3 2*x(1))^2-5,2*(x(1)-3)^2 (x(2)/3)^2-4]',[4,-4])