仿真实验3完整版本.docxVIP

3.17：

在计算机上用如下方法产生随机信号u(n)观测样本：首先产生一段零均值、方差为σ2的复高斯白噪声序列v(n)；然后在v(n)上叠加三个复正弦信号，它们的归一化频率分别是f1=0.15、f2=0.17和f3=0.26。调整σ2和正弦信号的幅度，使在f1、f2和f3处的信噪比分别为30dB、30dB和27dB。

（1）令信号观测样本长度N=32，试用3.1.1节讨论的基于FFT的自相关函数快速计算方法估计出自相关函数r?0(m)，并与教材式（3.1.2）估计出的自相关函数r?(m)作比较。

（2）令信号观测样本长度N=256，试分别用BT法和周期图法估计u(n)的功率谱，这里设BT法中所用自相关函数的单边长度M=64。

（3）令信号观测样本长度N=256，试用Levinson-Durbin迭代算法求解AR模型的系数并估计u(n)的功率谱，模型的阶数取为p=16。

(1)程序如下：

N=32;

noise=(randn(1,N)+1i*randn(1,N))/sqrt(2);

f1=0.15;f2=0.17;f3=0.26;SNR1=30;SNR2=30;SNR3=27;

A1=10^(SNR1/20);A2=10^(SNR2/20);A3=10^(SNR3/20);

signal1=A1*exp(1i*2*pi*f1*(0:N-1));

signal2=A2*exp(1i*2*pi*f2*(0:N-1));

signal3=A3*exp(1i*2*pi*f3*(0:N-1));

un=signal1+signal2+signal3+noise;

Uk=fft(un,2*N);Sk=(1/N)*abs(Uk).^2;r0=ifft(Sk);

r1=[r0(N+2:2*N),r0(1:N)];r=xcorr(un,N-1,biased);

r11=real(r1);r12=imag(r1);r1=real(r);r2=imag(r);

m=1-N:N-1;subplot(2,2,1);stem(m,r11,o);

xlabel(m);ylabel(实部);title(基于FFT的自相关函数);

subplot(2,2,2);stem(m,r12,o);xlabel(m);ylabel(虚部);

subplot(2,2,3);stem(m,r1);xlabel(m);ylabel(实部);

title(基于直接计算的自相关函数);

subplot(2,2,4);stem(m,r2);xlabel(m);ylabel(虚部);

(2)

N=256;noise=(randn(1,N)+1i*randn(1,N))/sqrt(2);

f1=0.15;f2=0.17;f3=0.26;

SNR1=30;SNR2=30;SNR3=27;

A1=10^(SNR1/20);A2=10^(SNR2/20);A3=10^(SNR3/20);

signal1=A1*exp(1i*2*pi*f1*(0:N-1));

signal2=A2*exp(1i*2*pi*f2*(0:N-1));

signal3=A3*exp(1i*2*pi*f3*(0:N-1));

un=signal1+signal2+signal3+noise;

NF=1024;Spr=fftshift((1/NF)*abs(fft(un,NF)).^2);

m=-511:512;

subplot(2,1,1);plot(m/1024,10*log10(Spr)/max(10*log10(Spr)));

axis([-0.5,0.5,0,1]);

xlabel(\omega/2*pi);ylabel(功率谱);title(周期图法);

M=64;r=xcorr(un,M,biased);NF=1024;

BT=fftshift(fft(r,NF));

subplot(2,1,2);plot(m/1024,10*log10(BT)/max(10*log10(BT)));

axis([-0.5,0.5,0,1]);

xlabel(\omega/2*pi);ylabel(功率谱);title(BT法);

(3)

N=256;

noise=(randn(1,N)+1i*randn(1,N))/sqrt(2);

f1=0.15;f2=0.17;f3=0.26;

SNR1=30;SNR2=30;SNR3=27;

A1=10^(SNR1/20);A2=10^(SNR2/20);A3=10^(SNR3/20);

signal1=A1*exp(1i*2*pi*f1*(0:N-1));

signal2=A2*exp(1i*2*pi*f2*(0:N-1));

signal3=A3*exp(1i*2*pi