clc
close
clear
taxis=100;
yaxis=2;
syms t;
Am=1;
omega=2;
beta=0;
A=0;
B=1;
C=2;
i1=2;
i2=0;
xt1=Am*cos(omega*t+beta);

funca1=sym('k1*cos(omega*t)+k2*sin(omega*t)');
funca1=subs(funca1);
x1=A*diff(funca1,'t',2)+B*diff(funca1,'t',1)+C*funca1;
b='x1=xt1';
b=char(subs(b));
ff1=subs(b,t,0);
ff2=subs(b,t,pi/2/omega);
s=solve(ff1,ff2,'k1','k2');
k1=s.k1;
k2=s.k2;
yf=subs(funca1)

f='A*D2y+B*Dy+C*y=xt1';
f=char(subs(f));
ini1='y(0)=i1';
ini1=char(subs(ini1));
ini2='Dy(0)=i2';
ini2=char(subs(ini2));
if i2~=0&A~=0
y=dsolve(f,ini1,ini2);
else y=dsolve(f,ini1);
end
y=subs(y)
 
yn=y-yf
 
 f0='A*D2y+B*Dy+C*y=0';
f0=char(subs(f0));
ini1='y(0)=i1';
ini1=char(subs(ini1));
ini2='Dy(0)=i2';
ini2=char(subs(ini2));
if i2~=0&A~=0
yzi=dsolve(f0,ini1,ini2);
else
    yzi=dsolve(f0,ini1);
end

yzi=subs(yzi)
yzs=y-yzi


t=0:0.1:8;
y=subs(y);subplot(2,3,1)
plot(y);
grid 
title('全響應(yīng)')
subplot(2,3,2)
t=0:0.1:8;
yf=subs(yf);
plot(yf);
grid 
title('受迫響應(yīng)')
subplot(2,3,3)
t=0:0.1:8;
yn=subs(yn);
plot(yn);
grid 
title('自然響應(yīng)')
subplot(2,3,4)
t=0:0.1:80;
y=subs(yzi);
plot(yzi);
grid 
title('零輸入響應(yīng)')
subplot(2,3,5)
t=0:0.1:8;
y=subs(y);
plot(y);
grid 
title('零狀態(tài)響應(yīng)')