clc;
clear;
Ap=0.3;fp=300;Np=300;h=-150;dt=0.0005;
As=1;fs=200;Ns=450;
for?i=300:1000
?A1（i）=Ap*exp（h*（i*dt-Np*dt））?*sin（2*pi*fp*（i*dt-Np*dt））;
end

for?i=450:1000
?????A2（i）=-As*exp（h*（i*dt-Ns*dt））?*sin（2*pi*fs*（i*dt-Ns*dt））;
end

for?i=1:1000
s（i）=A1（i）+A2（i）;
end



?

NOISE=randn（size（s））;

Y=s+0.03*NOISE;


A=abs（Y）;



for?i=180:1000
????LTA1（i）=0;
????for?n=i-60+1:i
????LTA1（i）=LTA1（i）+A（n）*A（n）;
????end?
end


for?i=180:1000
????LTA3（i）=0;
????for?n=i-180+1:i
????LTA3（i）=LTA3（i）+A（n）*A（n）;
????end?
end





for?i=180:1000
????STA（i）=0;
????for?n=i-20+1:i
???STA（i）=STA（i）+A（n）*A（n）;
????end?
end
sta=STA./20;
r1=3.*STA./LTA1;
lta1=LTA1./60;
r3=9.*STA./LTA3;
lta3=LTA3./180;
THR（1:1000）=3;
figure（1）
subplot?511?;plot（Y）;xlim（

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1436??2016-04-25?21:08??duanshic.m
?????文件?????????933??2016-04-25?21:21??jianquanmer.m
?????文件?????????588??2016-04-25?21:20??jiaquanxishu（1）.m
?????文件????????1228??2016-04-25?17:28??lta.m
?????文件????????1224??2016-04-25?21:34??CHANG.m