M=40;
T=0.025;%%취
A=1;
h=0.15;
wr=1000*pi;
Fs=2000;
N=5000;
t=(1:N+1)/(2*N)
st=A.*exp(-h.*wr.*t).*sin(wr.*t).*Step(t);
for ii=1:N+1
xt=0;
tt=t(ii);
for m=-M:M
% xt=xt+A.*exp(-h*wr*(t-m*T-Sum_tao(T,m,M)))*sin(wr.*(t-m*T-Sum_tao(T,m,M)))*Step(t-m*T-Sum_tao(T,m,M));
% xt=xt+A.*exp(-h*wr.*(tt-m*T-Sum_tao(T,m,M))).*sin(wr.*(tt-m*T-Sum_tao(T,m,M))).*Step(tt-m*T-Sum_tao(T,m,M));
xt=Sum_tao(T,m,M);
end;
xx(ii)=xt;
end;
%
figure(1);
plot(t,xx);
%%
%%-------------------------------------------%%
%%
yyi=abs(fft(yi))*2/N;
yys=abs(fft(ys))*2/N;
yyi=yyi(1:N/2);
yys=yys(1:N/2);
ff=0:fs/N:(N/2-1)*fs/N;
figure(2);
plot(ff,yyi);
hold on;
plot(ff,yys,'r');
%%
%------------------------------------------%%
yyei=envelope_h(yi);
yyei=yyei(1:N/2);
yyes=envelope_h(ys);
yyes=yyes(1:N/2);
figure(3);
plot(ff,yyei);
hold on;
plot(ff,yyes,'r');
figure(4);
plot(t,yyf);
figure(5);
yyef=envelope_h(yyf);
yyef=yyef(1:N/2);
plot(ff,yyef);