声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1366|回复: 5

[FFT] 关于白噪声的两个问题:频谱估计 和 与其他信号的卷积的频谱估计

[复制链接]
发表于 2008-6-12 17:30 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
问题描述:
某信号x具有频谱X,x和白噪声wn作卷积后生成x_conv,再作频谱估计得到X_CONV,理论上应该和X一样,但是我做出来的却不一样:1、最大值不同;2、有泄漏现象。
关于问题1,我单独写了一段程序,发现当对白噪声作分段频谱估计的时候,却与生成白噪声的频谱WN的不一样。
关于问题2,这种泄漏现象是怎样造成的,怎样才能生成下图

卷积后功率谱

卷积后功率谱

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

 楼主| 发表于 2008-6-12 17:31 | 显示全部楼层
白噪声卷积MATLAB脚本
clc
clear all
close all

lg = inline('10*log10(abs(squeeze(x)))','x');
%% 生成伪随机信号
N = 1024;      %数据点数
X = ones(N,1)*10^(-10);
X(51:200) = ones(150,1);
RandPhas = exp(j*(rand(N,1)*2*pi));    %随机相位
X = X.*RandPhas;
X(N/2+2:end) = conj(X(N/2:-1:2));      %频谱作“共轭对称”处理

x = ifft(X,'symmetric')*N;             %伪随机信号生成

figure(1)
subplot(2,1,1)
plot(lg(X))
title('伪随机信号')
ylabel('Amp dB');
axis tight
subplot(2,1,2)
plot(x)
axis tight
%% 生成白噪声
N1 = 100*N;
WN = ones(N1,1);
RandPhas = exp(j*(rand(N1,1)*2*pi));
WN = WN.*RandPhas;
WN(N1/2+2:end) = conj(WN(N1/2:-1:2));
wn = ifft(WN,'symmetric')*N1;
%wn = wn-mean(wn);

figure(2)
subplot(2,1,1)
plot(lg(WN))
title('白噪声')
ylabel('Amp dB');
axis tight
subplot(2,1,2)
plot(wn)
axis tight

x_conv = conv(x,wn);
len = length(x_conv);
ave = floor(len/N);
X_CONV = zeros(N,1);

for i = 1:ave
    tmp = x_conv((i-1)*N+(1:N));
    TMP = fft(tmp)/N;
    X_CONV = X_CONV+abs(TMP);
end
X_CONV = X_CONV/ave;

figure(3)
plot(lg(X_CONV));
title('卷积信号的频谱')
ylabel('Amp dB');
axis tight
 楼主| 发表于 2008-6-12 17:35 | 显示全部楼层
上面脚本的结果:

fig 1

fig 1

fig 2

fig 2

fig 3

fig 3
 楼主| 发表于 2008-6-12 17:37 | 显示全部楼层
关于白噪声分段频谱估计的程序及结果:
clc
clear all
close all
lg = inline('10*log10(abs(squeeze(x)))','x');
N = 1024;
N1 = 100*N;
WN = ones(N1,1);
RandPhas = exp(j*(rand(N1,1)*2*pi));
WN = WN.*RandPhas;
WN(N1/2+2:end) = conj(WN(N1/2:-1:2));
wn = ifft(WN,'symmetric')*N1;
ave = N1/N;
WN2 = zeros(N,1);
for i = 1:ave
    tmp = wn((i-1)*N+(1:N));
    TMP = fft(tmp)/N;
    WN2 = WN2+abs(TMP);
end
WN2 = WN2/ave;
figure(1)
plot((1:N1)/100,lg(WN),'g',1:N,lg(WN2),'b');
ylabel('Amp dB');
axis tight;
4.jpg
 楼主| 发表于 2008-6-12 17:41 | 显示全部楼层
小弟的信号处理是自学的,学的不到家,在上面的问题上停步了很久:@Q ,请高手指点迷津,小弟在此先谢过了:@)
 楼主| 发表于 2008-6-16 09:40 | 显示全部楼层
:@( 请大家给点提示吧,我自己也在想,但终不得其解:@Q :@Q :@Q :@Q
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-13 22:42 , Processed in 0.067735 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表