带有脉冲响应的xcorr函数

Question

我想在Matlab中设计一个Wiener滤波器来解卷积问题，但是我遇到了很多问题。我有一个方差为1.2的高斯白噪声过程和一个长度为2的脉冲响应。它的值是g（0）= 5和g（1）= 4。稍后在这个过程中，我尝试确定Rxx（m）。为此，我需要计算g（m）* g（-m）（卷积），并建议在Matlab中使用xcorr函数，但我的结果没有意义。任何人都可以通过这个xcorr函数来帮助我，并提供关于如何使用这个脉冲响应的建议吗？我曾尝试过g的fourrier变换，但这一直没有帮助。

Answer 1

下面的代码只实现了我可以在描述中看到的一部分内容。它会产生噪音过程并执行第一部分所述的内容。自相关不是用滤波器系数计算的，而是用实际信号计算的。

% generate noise process y 
y = randn(1,N) * sqrt(1.2); 

% filter y with the fir-filter g 
g = [2, 0.6]; 
r = filter(g,1,y); 

% generate noise process d 
d = randn(1,N) * sqrt(0.2); 

% x is the sum of r and d 
x = r + d; 

% autocorrelation of x 
[Rxx,lagRxx] = xcorr(x); 

% plot autocorrelation 
figure; grid on; 
plot(lagRxx,Rxx); 
title('Biased Autocorrelation of Signal x'); 
xlabel('Lag'); 

% cross correlation between x and y 
[Rxy,lagRxy] = xcorr(x,y); 

% plot crosscorrelation 
figure; grid on; 
plot(lagRxy,Rxy); 
title('Biased Crosscorrelation of Signal x and y'); 
xlabel('Lag');