现代信号处理及其应用仿真.docVIP

仿真作业 周子超 200820003043 3.17 (1) clear all close all f1=0.15; f2=0.17; f3=0.26; v=randn(1,256); A1=10^(30/20); A2=10^(30/20); A3=10^(30/27); t=0:255; x1=A1*exp(j*2*pi*f1*t); x2=A2*exp(j*2*pi*f2*t); x3=A3*exp(j*2*pi*f3*t); u=v+x1+x2+x3; u=u(:,1:32); u1=[u,zeros(1,32)]; uk=fft(u1); for i=1:64 uuk(i)=(1/32)*uk(i)*conj(uk(i)); end r0=ifft(uuk); rr1=[r0(1,34:64),r0(1,1:32)]; r=0; r10=0; r1=zeros(1,63); u=[u,zeros(1,32)]; for n=1:32 r=(1/32)*u(n)*conj(u(n)); r10=r+r10; end for m=1:31 for n=m+1:32 r=(1/32)*u(n)*conj(u(n-m)); r1(1,m)=r1(1,m)+r; end end for m=32:62 for n=1:31 r=(1/32)*u(n)*conj(u(n+m-31)); r1(1,m)=r1(1,m)+r; end end rr2=[r1(1,62:-1:32),r10,r1(1,1:31)]; 通过上面的计算可以发现基于FFT的自相关函数快速算法估计的自相关函数和式(3.1.2)估计出的自相关函数相等。 (2) clear all close all f1=0.15; f2=0.17; f3=0.26; v=randn(1,256); A1=10^(30/20); A2=10^(30/20); A3=10^(27/20); t=0:255; x1=A1*exp(j*2*pi*f1*t); x2=A2*exp(j*2*pi*f2*t); x3=A3*exp(j*2*pi*f3*t); u=v+x1+x2+x3; U=fft(u); for i=1:256 SPER(i)=abs(U(i)*conj(U(i)))/256; end SPER=10*log(SPER); n=0:length(u)-1; k=n/(length(u)-1); figure plot(k,SPER) xlabel(归一化频率) ylabel(功率谱(dB)) title(周期图法) u=u(:,1:256); r=0; r10=0; r1=zeros(1,511); u=[u,zeros(1,256)]; for n=1:256 r=(1/256)*u(n)*conj(u(n)); r10=r+r10; end for m=1:255 for n=m+1:256 r=(1/256)*u(n)*conj(u(n-m)); r1(1,m)=r1(1,m)+r; end end for m=256:510 for n=1:255 r=(1/256)*u(n)*conj(u(n+m-255)); r1(1,m)=r1(1,m)+r; end end rr=[r1(1,510:-1:256),r10,r1(1,1:255)]; rr=rr(1,193:319); k=-63:63; w=pi*k/63+pi; Sbtw=rr*(exp(-j)).^(k*w); SBTW=abs(Sbtw); P=10*log(SBTW); figure plot(w/(2*pi),P) xlabel(归一化频率) ylabel(功率谱(dB)) title(BT法) (3) clear all close all f1=0.15; f2=0.17; f3=0.26; v=randn(1,256); A1=10^(30/20); A2=10^(30/20); A3=10^(27/20); t=0:255; x1=A1*exp(j*2*pi*f1*t); x2=A2*exp(j*2*pi*f2*t); x3=A3*exp(j*2*pi*f3*t); u=v+x1+x2+x3; r=0; r10=0; r1=zeros(1,511); u=[u,zeros(1,256)]; for n=1:256 r=(1/256)*u(n)*conj(u(n)); r10=r+r10; end for m=