%Appl9_1
load wmq318.VHZ     %加載第一天的數(shù)據(jù)
load wmq319.VHZ    %加載第二天的數(shù)據(jù)
load wmq320.VHZ  %加載第三天的數(shù)據(jù)
load wmq321.VHZ  %加載第四天的數(shù)據(jù)
[m,n]=size(wmq318);   %將數(shù)據(jù)格式轉(zhuǎn)換成該程序使用的格式，以下同
nn=m*n;
b318=reshape(wmq318',nn,1);  
[m,n]=size(wmq319);
nn=m*n;
b319=reshape(wmq319',nn,1);
[m,n]=size(wmq320);
nn=m*n;
b320=reshape(wmq320',nn,1);
[m,n]=size(wmq321);
nn=m*n;
b321=reshape(wmq321',nn,1);
b=[b318;b319;b320;b321];   %將幾天的數(shù)據(jù)連接起來組成一個(gè)矩陣
[WMQZPxx,f] = psd(b,16384,1/10);   %功率譜分析
plot(f,WMQZPxx);   %繪圖
ylim([0 5.5e9])   %將圖放大到合適比例
xlim([1.2 6]*1.0E-3)   %為了清楚地觀看某一頻段的功率譜，可改變此間隔
xlabel('頻率/Hz');ylabel('功率譜')