[数字信号处理实验四.docVIP

实验一：DFS、DFT与FFT 一、实验内容 2、已知某周期序列的主值序列为x(n)=[0,1,2,3,2,1,0]，编程显示2个周期的序列波形。要求： ① 用傅里叶级数求信号的幅度谱和相位谱，并画出图形 ② 求傅里叶级数逆变换的图形，并与原序列进行比较。 N=7; xn=[0,1,2,3,2,1,0]; xn=[xn,xn]; n=0:2*N-1; k=0:2*N-1; Xk=xn*exp(-j*2*pi/N).^(n*k); x=(Xk*exp(j*2*pi/N).^(n*k))/N; subplot(2,2,1);stem(n,xn); title(x(n));axis([-1,2*N,1.1*min(xn),1.1*max(xn)]); subplot(2,2,2);stem(n,abs(x)); title(IDFS|X(k)|);axis([-1,2*N,1.1*min(x),1.1*max(x)]); subplot(2,2,3),stem(k,abs(Xk)); title(|X(k)|);axis([-1,2*N,1.1*min(abs(Xk)),1.1*max(abs(Xk))]); subplot(2,2,4),stem(k,angle(Xk)); title(arg|X(k)|);axis([-1,2*N,1.1*min(angle(Xk)),1.1*max(angle(Xk))]); 3、已知有限长序列x(n)=[1,0.5,0,0.5,1,1,0.5,0]，要求： ① 求该序列的DFT、IDFT的图形； xn=[1,0.5,0,0.5,1,1,0.5,0]; N=length(xn); n=0:N-1; k=0:N-1; Xk=xn*exp(-1i*2*pi/N).^(n*k); x=(Xk*exp(1i*2*pi/N).^(n*k))/N; subplot(2,2,1);stem(n,xn); title(x(n));axis([-1,N,1.1*min(xn),1.1*max(xn)]); subplot(2,2,2);stem(n,abs(x)); title(IDFT|X(k)|);axis([-1,N,1.1*min(x),1.1*max(x)]); subplot(2,2,3),stem(k,abs(Xk)); title(|X(k)|);axis([-1,N,1.1*min(abs(Xk)),1.1*max(abs(Xk))]); subplot(2,2,4),stem(k,angle(Xk)); title(arg|X(k)|);axis([-1,N,1.1*min(angle(Xk)),1.1*max(angle(Xk))]); ② 用FFT算法求该序列的DFT、IDFT的图形； xn=[1,0.5,0,0.5,1,1,0.5,0]; N=length(xn); subplot(2,2,1);stem(n,xn); title(x(n)); k=0:N-1; Xk=fft(xn,N); subplot(2,1,2);stem(k,abs(Xk)); title(Xk=DFT(xn)); xn1=ifft(Xk,N); subplot(2,2,2);stem(n,xn1); title(x(n)=IDFT(Xk)); ③ 假定采用频率Fs=20Hz，序列长度N分别取8、32和64，用FFT计算其幅度谱和相位谱。 Ts=0.05;C=[32,64,8]; for r=0:2; N=C(r+1); n=0:N-1; xn=[1,0.5,0,0.5,1,1,0.5,0,zeros(1,N-8)]; D=2*pi/(N*Ts); k=floor(-(N-1)/2:(N-1)/2); X=fftshift(fft(xn,N)); subplot(3,2,2*r+1);plot(k*D,abs(X)); subplot(3,2,2*r+2);stairs(k*D,angle(X)); end 4、已知一个无限长序列x(n)=0.5n（n≥0），采