用FFT对信号作频谱分析(4)

X6k64=fftshift(X6k64); %将零频率移到频谱中心 Tp=N*T;F=1/Tp; %频率分辨率F

k=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心） subplot(3,1,3);stem(fk,abs(X6k64),'.'); box on%绘制8点DFT的幅频特性图 title('(6a) 64点|DFT[x_6(nT)]|');xlabel('f(Hz)');ylabel('幅度'); axis([-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k64))])

3.4对x9(n)作FFT变换区间N为8和16时的频谱分析

M=8;xa=1:(M/2); xb=(M/2):-1:1;

x2n=[xa,xb]; %产生长度为8的三角波序列x2(n) x3n=[xb,xa]; x9n=x2n+x3n*j X9k8=fft(x9n,8); X9k16=fft(x9n,16); figure(1); N=8;

f=2/N*(0:N-1);

subplot(2,2,1);stem(f,abs(X9k8),'.'); %绘制8点DFT的幅频特性图 title('(9a) 8点DFT[x_9(n)]');xlabel('ω/π');ylabel('幅度'); N=16;

f=2/N*(0:N-1);

subplot(2,2,2);stem(f,abs(X9k16),'.'); %绘制16点DFT的幅频特性图 title('(9b) 16点DFT[x_9(n)]');xlabel('ω/π');ylabel('幅度');

1,2,3,3,2,1?的有关问题 3.5 序列x(n)??w=2*pi*(0:255)/256;

xw=1+2*exp(-j*w)+3*exp(-j*2*w)+3*exp(-j*3*w)+2*exp(-j*4*w)+exp(-j*5*w);

figure(1); subplot(2,1,1);

plot(w/pi,abs(xw),'.'); title('x(n)的幅频曲线');

subplot(2,1,2); plot(w/pi,angle(xw)); line([0,2],[0,0]); title('x(n)的相频曲线'); figure(2); N1=6;K1=0:N1-1; N2=18;K2=0:N2-1; N3=36;K3=0:N3-1; xn=[1,2,3,3,2,1]; xk1=fft(xn,N1);

subplot(3,2,1);stem(K1,abs(xk1),'.'); title('N=6点时的DFT[x(n)]=x1(k)');hold on; plot(N1/2*w/pi,abs(xw),'r');

subplot(3,2,2);stem(K1,angle(xk1),'.'); title('x1(k)的相位'); xk2=fft(xn,N2);

subplot(3,2,3);stem(K2,abs(xk2),'.'); title('N=18点的DFT[x(n)]=X2(k)');hold on; plot(N2/2*w/pi,abs(xw),'.');

subplot(3,2,4);stem(K2,angle(xk2),'.'); title('X2(k)的相位');xk3=fft(xn,N3); subplot(3,2,5);stem(K3,abs(xk3),'.'); title('N=36点时的DFT[x(n)]=X3(k)');hold on; plot(N3/2*w/pi,abs(xw),'r');

subplot(3,2,6);stem(K3,angle(xk3),'.'); title('X3(k)的相位');hold on; figure(3);

xn1=ifft(xk1,N1);stem(K1,xn1);

title('x1(k)作2DFT,得到的xl(n)与原序列x(n)一致');

3.6 用FFT对信号x2(n)、x3(n)选择合适的变换区间N进行谱分析

clear all;close

all; %==================================================================

n2=0:50;n3=-10:10; ;N2=20;N3a=32;N3b=64;

x2n=2*sin(0.45*pi*n2).*sin(0.55*pi*n2); %计算序列x2(n) x3n=0.5.^abs(n3); %计算序列x3(n)

x3anp=zeros(1,N3a); %构造x3(n)的周期延拓序列,周期为N3a

for m=1:10,

x3anp(m)=x3n(m+10);x3anp(N3a+1-m)=x3n(11-m); end

x3bnp=zeros(1,N3b); %构造x3(n)的周期延拓序列,周期为N3b

for m=1:10,

x3bnp(m)=x3n(m+10);x3bnp(N3b+1-m)=x3n(11-m); end

X2k=fft(x2n,N2); %计算序列x2(n)的N2点DFT X3ak=fft(x3anp,N3a); %计算序列x3(n)的N3a点DFT X3bk=fft(x3bnp,N3b); %计算序列x3(n)的N3b点DFT %以下为绘图部分

%==============绘制x2(n)的频谱特性图============================

k=0:N2-1;wk=2*k/N2; %产生N2点DFT对应的采样点频率(关于π归一化值) subplot(3,2,2);stem(wk,abs(X2k),'.');grid on;box on %绘制x2(n)的N2点DFT的幅频特性图

title('(c) x2(n)的幅频特性图');xlabel('ω/π');ylabel('幅度')

subplot(3,2,4);stem(wk,angle(X2k),'.');grid on;box on %绘制x2(n)的N2点DFT的相频特性图

line([0,2],[0,0]) %画横坐标轴线 title('(d) x2(n)的相频特性图'); xlabel('ω/π');ylabel('相位'); figure(2)

%=============绘制32点周期延拓序列和32点DFT[x3(n)]的频谱特性图===============

n=0:N3a-1;subplot(3,2,1);stem(n,x3anp,'.');box on title('(e) x3(n)的32点周期延拓序列');

k=0:N3a-1;wk=2*k/N3a; %产生N3a点DFT对应的采样点频率(关于π归一化值) subplot(3,2,3);plot(wk,abs(X3ak)); %绘制x3(n)的N3a点DFT的幅频特性图 title('(f) DFT[x3(n)]_3_2的幅频特性图');xlabel('ω/π');ylabel('幅度') subplot(3,2,5);plot(wk,angle(X3ak)); %绘制x3(n)的N3点DFT的相频特性图 line([0,2],[0,0]) %画横坐标轴线 title(' (g) DFT[x3(n)]_3_2的相位'); xlabel('ω/π');ylabel('相位');

%=============绘制64点周期延拓序列和64点DFT[x3(n)]的频谱特性图===============

n=0:N3b-1;subplot(3,2,2);stem(n,x3bnp,'.');box on title('(h) x3(n)的64点周期延拓序列');

k=0:N3b-1;wk=2*k/N3b; %产生N3a点DFT对应的采样点频率(关于π归一化值) subplot(3,2,4);plot(wk,abs(X3bk)); %绘制x3(n)的N3a点DFT的幅频特性图 title('(i) DFT[x3(n)]_6_4的幅频特性图');xlabel('ω/π');ylabel('幅度') subplot(3,2,6);plot(wk,angle(X3bk)); %绘制x3(n)的N3点DFT的相频特性图 line([0,2],[0,0]) %画横坐标轴线

title(' (j) DFT[x3(n)]_3_2的相位'); xlabel('ω/π');ylabel('相位');