x.m

来自「应用经典prony法和AR参数建模进行模态参数识别」· M 代码 · 共 5 行

function y=x(t)
Wn=10;D=0.05;n=Wn*D;Wd=Wn*sqrt(1-D^2);X0=0.001;
A=sqrt(X0^2+(D*Wn*X0/Wd)^2);alp=atan(Wd*X0/(D*Wn*X0));
y=A*exp(-n*t)*sin(Wd*t+alp);