|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clc;
%预定参数;
node=5;
filepath='D:\history\';
msgbox('正在计算,请稍候...','运行程序','help','replace');
for i=1:node
A=load([filepath 'X-' num2str(i) '.txt']);
eval(['P',num2str(i),'=','A(:,2)',';']);%分别保存在P1,P2...
eval(['p_his','(:,i)','=','P',num2str(i),';']);%%
end %这里只是将信号组放在一个矩阵里面
msgbox('正在计算!','运行程序','help','replace');
for i=1:5
Fs=1/0.18; %采样频率
n=0:1/Fs:2952;
nfft=16400;
window2=hanning(16400); %汉宁窗
noverlap=20; %数据无重叠
range='half'; %频率间隔为[0 Fs/2],只计算一半的频率
[sp,f]=pwelch(p_his(:,i),window2,noverlap,nfft,Fs,range);
spp(:,i)=sp;ff(:,i)=f;
plot_spp=log10(abs(spp(:,i)));
figure(1);
grid on;
subplot(3,2,i);
plot(log10(abs(ff(:,i))),plot_spp);
end
根据上面的方法我求的的功率谱打出的图如下,怎么有负值啊?曲线好像是合理的,但是我要的是正值啊,幅值也不知道对不对,没有把握。
请大侠答疑啊。谢谢。
|
|