现代数字信号处理仿真作业.docxVIP

现代数字信号处理仿真作业 1.仿真题3.17 仿真结果及图形： 图 SEQ 图 \* ARABIC 1 基于FFT的自相关函数计算 图 图 SEQ 图 \* ARABIC 2 基于式3.1.2的自相关函数的计算 图 SEQ 图 \* ARABIC 3 周期图法和BT法估计信号的功率谱 图 SEQ 图 \* ARABIC 4 利用LD迭代对16阶AR模型的功率谱估计 16阶AR模型的系数为： a1=-0.402637623107952-0.919787323662670i; a2=-0.013530139693503+0.024214641171318i; a3=-0.074241889634714-0.088834852915013i; a4=0.027881022353997-0.040734794506749i; a5=0.042128517350786+0.068932699075038i; a6=-0.0042799971761507 + 0.028686095385146i; a7=-0.048427890183189 - 0.019713457742372i; a8=0.0028768633718672 - 0.047990801912420i a9=0.023971346213842+ 0.046436389191530i; a10=0.026025963987732 + 0.046882756497113i; a11= -0.033929397784767 - 0.0053437929619510i; a12=0.0082735406293574 - 0.016133618316269i; a13=0.031893903622978 - 0.013709547028453i ; a14=0.0099274520678052 + 0.022233240051564i; a15=-0.0064643069578642 + 0.014130696335881i; a16=-0.061704614407581- 0.077423818476583i. 仿真程序（3_17）: clear all clc %% 产生噪声序列 N=32; %基于FFT的样本长度 %N=256; %周期图法，BT法，AR模型功率谱估计的长度 vn=(randn(1,N)+1i*randn(1,N))/sqrt(2); %%产生复正弦信号 f=[0.15 0.17 0.26]; %归一化频率 SNR=[30 30 27]; %信噪比 A=10.^(SNR./20); %幅度 signal=[A(1)*exp(1i*2*pi*f(1)*(0:N-1)); %复正弦信号 A(2)*exp(1i*2*pi*f(2)*(0:N-1)); A(3)*exp(1i*2*pi*f(3)*(0:N-1))]; %% 产生观察样本 un=sum(signal)+vn; %% 利用3.1.1的FFT估计 Uk=fft(un,2*N); Sk=(1/N)*abs(Uk).^2; r0=ifft(Sk); r1=[r0(N+2:2*N),r0(1:N)]; %% 利用3.1.2估计R r2=xcorr(un,N-1,biased); % 画图 k=-N+1:N-1; figure(1) subplot(1,2,1) stem(k,real(r1)) xlabel(m);ylabel(实部); subplot(1,2,2) stem(k,imag(r1)) xlabel(m);ylabel(虚部); figure(2) subplot(1,2,1) stem(k,real(r2)) xlabel(m);ylabel(实部); subplot(1,2,2) stem(k,imag(r2)) xlabel(m);ylabel(虚部); %% 周期图法 NF=1024; Spr=fftshift((1/NF)*abs(fft(un,NF)).^2); kk=-0.5+(0:NF-1)*(1/(NF-1)); Spr_norm=10*log10(abs(Spr)/max(abs(Spr))); %% BT法 M=64; r3=xcorr(un,M,biased); BT=fftshift(fft(r3,NF)); BT_norm=10*log10(abs(BT)/max(abs(BT))); figure(3) subplot(1,2,1) plot(kk,Spr_norm) xlabel(w/2pi);ylabel(归一化功率谱/DB); title(周期图法) subplot(1,2