上学时搞得,粘到m文件中执行就可以了，必须要明白基本原理哈，否则看不懂的。
clear
syms x t dr A pd w1 w f1 acf
A=0.01;dr=0.8;w1=1000;acf=(200*pi)/30;f1=0.6;
x=A.*exp(-dr.*w1.*t).*(sin(sqrt(1-dr^2).*w1.*t).*(1./sqrt((w1^2-(acf.*t).^2).^2+(2.*dr.*w1.*(acf.*t)).^2)).*(w1^2-(acf.*t).^2)+cos(sqrt(1-dr^2).*w1.*t).*(1./sqrt((w1^2-(acf.*t)^2).^2+(2.*dr.*w1.*(acf.*t)).^2)).*2.*dr.*w1.*(acf.*t))+((w1^2-(acf.*t).^2).*f1.*sin((acf.*t).*t)-2.*dr.*w1.*(acf.*t).*f1.*cos((acf.*t).*t))./((w1^2-(acf.*t).^2).^2+(2.*dr.*w1.*(acf.*t)).^2);
t=linspace(0,12,12000);
subs(x);
plot(t,subs(x))