matlab 汉明窗求功率谱 错在哪里?求help
来源:学生作业帮 编辑:大师作文网作业帮 分类:综合作业 时间:2024/11/12 18:29:21
matlab 汉明窗求功率谱 错在哪里?求help
%汉明窗求功率谱
window=hamming(20); % 采用hanmming 窗,长度为20
noverlap=10; % 重叠的点数
Nfft=512; % 做FFT 的点数
Fs=1000; % 采样频率,为1000Hz
[Px,f]=pwelch(x,window,noverlap,Nfft,Fs,'onesided'); % 估计功率谱密度
f=[-fliplr(f) f(2:end)]; % 构造一个对称的频率,范围是[-Fs/2,Fs/2]
Py=[-fliplr(Py) Py(2:end)]; % 对称的功率谱
plot(f,10*log10(Py),'b');
错误信息:Error in ==> aaaa at 15
f=[-fliplr(f) f(2:end)]; % 构造一个对称的频率,范围是[-Fs/2,Fs/2]
%汉明窗求功率谱
window=hamming(20); % 采用hanmming 窗,长度为20
noverlap=10; % 重叠的点数
Nfft=512; % 做FFT 的点数
Fs=1000; % 采样频率,为1000Hz
[Px,f]=pwelch(x,window,noverlap,Nfft,Fs,'onesided'); % 估计功率谱密度
f=[-fliplr(f) f(2:end)]; % 构造一个对称的频率,范围是[-Fs/2,Fs/2]
Py=[-fliplr(Py) Py(2:end)]; % 对称的功率谱
plot(f,10*log10(Py),'b');
错误信息:Error in ==> aaaa at 15
f=[-fliplr(f) f(2:end)]; % 构造一个对称的频率,范围是[-Fs/2,Fs/2]
f=[-fliplr(f') (f(2:end))'];
Py=[-fliplr(Px') (Px(2:end))'];
%这两行改成上面的
Py=[-fliplr(Px') (Px(2:end))'];
%这两行改成上面的