一元多次方程在MATLAB中的解以及图形显示

• 一元二次方程

% 求解 ax2+bx+c=0   2次方程有两个个解

syms a b c x;
[x] = solve(a*x^2+b*x+c==0,x,'Real',false)
syms x
[x] = solve(x^2+x+1==0,x,'Real',false)  %false 求复数解 true 求实数解
%图形
x = -2:0.01:1;
y = x.^2+x+1;
y1 = zeros(size(x));
plot(x,y,'b');
hold on 
plot(x,y1,'r--')

        

• 一元三次方程

% 求解 ax^3+bx^2+cx+d=0
syms a b c d x
[x] = solve(a*x^3+b*x^2+c*x+d==0,x) 
%  root(az3+bz2+cz+d,z,1)
syms x;
[x] = solve(x^3+x^2+x+1==0,x)  % 解x = -1,-i,i
syms x;
[x] = solve(x^3+x^2+x+1==0,x,'real',true) % x = -1
syms x;
[x] = solve(2*x^3+5*x^2+x+1==3,x) % x = -1

%图形
x = -4:0.01:3;
y = 2*x.^3+5*x.^2+x+1;
y1 = zeros(size(x));
plot(x,y,'b');
hold on 
plot(x,y1,'r--')

x = 

 

• 一元四次方程

% 求解 ax^4+bx^3+cx^2+dx+e=0
syms a b c d e x;
[x] = solve(a*x^4+b*x^3+c*x^2+d*x+e==0,x)
syms  x;
[x] = solve(x^4+x^3+x^2+x+1==1,x,'Real',true) %实数解

%图形
x = -8:0.01:4;
y = 2*x.^4+17*x.^3+7*x.^2+8*x+1;
y1 = zeros(size(x));
plot(x,y,'b');
hold on 
plot(x,y1,'r--')

x =

x =

• 一元五次方程

% 求解 ax^5+bx^4+cx^3+dx^2+ex+f==0
syms a b c d e f x
[x] = solve(a*x^5+b*x^4+c*x^3+d*x^2+e*x+f==0,x)
syms  x
[x] = solve(x^5+x^4+x^3+x^2+x+1==0,x,'Real',true)

%图形
x = -6:0.01:4;
y = 2*x.^5+12*x.^4+8*x.^3+20*x.^2+42*x+1;
y1 = zeros(size(x));
plot(x,y,'b');
hold on 
plot(x,y1,'r--')

x =

x =

• 一元六次方程

% 求解 ax^6+bx^5+cx^4+dx^3+ex^2+fx+g=0
syms a b c d e f g x
[x] = solve(a*x^6+b*x^5+c*x^4+d*x^3+e*x^2+f*x+g==0,x)
syms x
[x] = solve(x^6+x^5+x^4+x^3+x^2+x+1==1,x,'Real',true)

%图形
x = -30:0.01:20;
y = x.^6+25*x.^5+20*x.^4+30*x.^3+40*x.^2+50*x+1;
y1 = zeros(size(x));
plot(x,y,'b');
hold on 
plot(x,y1,'r--')

x =

x =