%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%　　　對VSP地震資料進行模擬，并利用譜比法進行Q值估計　%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear?all
close?all
%%??參數設置??%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Fs=1000;????%%%采樣頻率
Ts=1/Fs;
t=（0:999）*Ts;
Lt=length（t）;
vp（1）=1500;??%%%%%%%%%?速度
vp（2）=2000;
vp（3）=2500;?
vp（4）=3000;?
h（1）=300;???%%%%%%%%%?深度
h（2）=700;?
h（3）=1200;?
Q（1）=40;????%%%%%%%%%?品質因子
Q（2）=70;
Q（3）=100;
Q（4）=50;
Dx=20;???????%%%%%%%%%?道間距
w0=100*pi;??
size=Lt;
a（1）=h（1）/Dx;
a（2）=h（2）/Dx;
a（3）=h（3）/Dx;
Ds=80;???????%%%%%%%總道數
%%?求各道的到時??%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for?j=1:1:Ds
????if?j<=a（1）
???????T（j）=j*Dx/vp（1）;
????elseif?j>a

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2832??2012-03-07?21:51??譜比法\pubifa.m
?????文件????????2846??2011-07-01?09:35??譜比法\wigb.m
?????目錄???????????0??2014-05-25?08:09??譜比法\