作业帮 > 综合 > 作业

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

来源:学生作业帮 编辑:拍题作业网作业帮 分类:综合作业 时间:2024/04/30 01:45:05
matlab 求代数方程组 a*x^2+b*x+c=0 x+y=0 关于x,y的解,并分别绘制x和y关于b和c的图像(a视作常数)
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')