th.m

来自「利用FFT与IFFT制造满足平滑设计谱之人造地震波」· M 代码 · 共 17 行

clear all;

load d:\EQX.txt;                            
u=EQX(4000:14000,1);			

Ug=u/max(abs(u));                              
N=length(Ug);
dt=0.005;                                       
r=0.5;                                         

S_D1=0.52;                                      
S_DS=0.7;                                      
T_0D=S_D1/S_DS;                                 

S_aD=zeros(N,1);
T=zeros(1,N);
%劍﹚壉瞯㎝㏄戳