Question

matlab 求代数方程组 a*x^2+b*x+c=0 x+y=0 关于x,y的解，并分别绘制x和y关于b和c的图像（a视作常数）

Answer 1

1. 求代数方程组的解：

>> [x,y]=solve('a*x^2+b*x+c=0','x+y=0','x,y')

 

x =

  1/2/a*(-b+(b^2-4*a*c)^(1/2))

 1/2/a*(-b-(b^2-4*a*c)^(1/2))

 y =

  -1/2/a*(-b+(b^2-4*a*c)^(1/2))

 -1/2/a*(-b-(b^2-4*a*c)^(1/2))

 >> 

2. 从上面的解可以看出，x，y都有两组解且x，y互为相反数。

假设a=1，这里有两种方法绘制x，y关于b，c的图像：

（1）隐函数绘图

x1=subs(x(1),'a',1);

x2=subs(x(2),'a',1);

y1=subs(y(1),'a',1);

y2=subs(y(2),'a',1);

figure

po=get(gcf,'position');

set(gcf,'position',[po(1)-0.5*po(3) po(2) 2*po(3) po(4)]);

subplot(121)

ezsurf(x1,[-10 10])

hold on

ezsurf(x2,[-10 10])

subplot(122)

ezsurf(y1,[-10 10])

hold on

ezsurf(y2,[-10 10])

（2）根据方程式直接绘图

>> a=1;

>> [b,c]=meshgrid(-10:0.5:10);

>> delta=b.^2-4*a*c;

>> delta(delta<0)=NaN;

>> x1=0.5/a*(-b+sqrt(delta));

>> x2=0.5/a*(-b-sqrt(delta));

>> y1=-x1;

>> y2=-x2;

>> figure

>> po=get(gcf,'position');

>> set(gcf,'position',[po(1)-0.5*po(3) po(2) 2*po(3) po(4)]);

>> subplot(121)

>> surf(b,c,x1)

>> hold on

>> surf(b,c,x2)

>> xlabel('b');ylabel('c');zlabel('x')

>> subplot(122)

>> surf(b,c,y1)

>> hold on

>> surf(b,c,y2)

>> xlabel('b');ylabel('c');zlabel('y')