appl9_1.m

来自「MATLAB7.x数字信号处理 光盘内容」· M 代码 · 共 23 行

%Appl9_1
load wmq318.VHZ     %加载第一天的数据
load wmq319.VHZ    %加载第二天的数据
load wmq320.VHZ  %加载第三天的数据
load wmq321.VHZ  %加载第四天的数据
[m,n]=size(wmq318);   %将数据格式转换成该程序使用的格式，以下同
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];   %将几天的数据连接起来组成一个矩阵
[WMQZPxx,f] = psd(b,16384,1/10);   %功率谱分析
plot(f,WMQZPxx);   %绘图
ylim([0 5.5e9])   %将图放大到合适比例
xlim([1.2 6]*1.0E-3)   %为了清楚地观看某一频段的功率谱，可改变此间隔
xlabel('频率/Hz');ylabel('功率谱')