1、数字信号处理 实验实验1xn=1 1 1 1 ; %输入时域序列向量xn=R4(n)Xk8=fft(xn,8); %计算xn的8点DFTXk16=fft(xn,16); %计算xn的16点DFTXk32=fft(xn,32); %计算xn的32点DFTk=0:7;wk=2*k/8; %产生8点DFT对应的采样点频率(关于归一化值)subplot(3,2,1);stem(wk,abs(Xk8),.); %绘制8点DFT的幅频特性图title(a) 8点DFT的幅频特性图);xlabel(/);ylabel(幅度);grid subplot(3,2,2);stem(wk,angle(Xk8),.)
2、; %绘制8点DFT的相频特性图title(b) 8点DFT的相频特性图);grid on;xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5)k=0:15;wk=2*k/16; %产生16点DFT对应的采样点频率(关于归一化值)subplot(3,2,3);stem(wk,abs(Xk16),.); %绘制16点DFT的幅频特性图title(c) 16点DFT的幅频特性图);xlabel(/);ylabel(幅度);grid subplot(3,2,4);stem(wk,angle(Xk16),.); %绘制16点DFT的相频特性图title(d) 16点DFT的相
3、频特性图) ;grid on;xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5)k=0:31;wk=2*k/32; %产生32点DFT对应的采样点频率(关于归一化值)subplot(3,2,5);stem(wk,abs(Xk32),.); %绘制32点DFT的幅频特性图title(e) 32点DFT的幅频特性图);xlabel(/);ylabel(幅度);gridsubplot(3,2,6);stem(wk,angle(Xk32),.); %绘制32点DFT的相频特性图title(f) 32点DFT的相频特性图); grid on;xlabel(/);ylabel(
4、相位);axis(0,2,-3.5,3.5);实验2xn=1 2 3 4 4 3 2 1; %输入时域序列向量x2Xk8=fft(xn,8); %计算xn的8点DFTXk16=fft(xn,16); %计算xn的16点DFT%以下为绘图部分k=0:7;wk=2*k/8; %产生8点DFT对应的采样点频率(关于归一化值)subplot(2,2,1);stem(wk,abs(Xk8),.); %绘制8点DFT的幅频特性图title(a) 8点DFT的幅频特性图);xlabel(/);ylabel(幅度);grid subplot(2,2,2);stem(wk,angle(Xk8), .); %绘制
5、8点DFT的相频特性图title(b) 8点DFT的相频特性图);grid xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5)k=0:15;wk=2*k/16; %产生16点DFT对应的采样点频率(关于归一化值)subplot(2,2,3);stem(wk,abs(Xk16), .); %绘制16点DFT的幅频特性图title(c) 16点DFT的幅频特性图);xlabel(/);ylabel(幅度);grid subplot(2,2,4);stem(wk,angle(Xk16), .); %绘制16点DFT的相频特性图title(d) 16点DFT的相频特性图);
6、grid xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5)实验3xn=4 3 2 1 1 2 3 4; %输入时域序列向量x2Xk8=fft(xn,8); %计算xn的8点DFTXk16=fft(xn,16); %计算xn的16点DFT%以下为绘图部分k=0:7;wk=2*k/8; %产生8点DFT对应的采样点频率(关于归一化值)subplot(2,2,1);stem(wk,abs(Xk8),.); %绘制8点DFT的幅频特性图title(a) 8点DFT的幅频特性图);xlabel(/);ylabel(幅度);grid subplot(2,2,2);stem(w
7、k,angle(Xk8), .); %绘制8点DFT的相频特性图title(b) 8点DFT的相频特性图);grid xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5)k=0:15;wk=2*k/16; %产生16点DFT对应的采样点频率(关于归一化值)subplot(2,2,3);stem(wk,abs(Xk16), .); %绘制16点DFT的幅频特性图title(c) 16点DFT的幅频特性图);xlabel(/);ylabel(幅度);grid subplot(2,2,4);stem(wk,angle(Xk16), .); %绘制16点DFT的相频特性图ti
8、tle(d) 16点DFT的相频特性图);grid xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5)实验3xn=4 3 2 1 1 2 3 4; %输入时域序列向量x2Xk8=fft(xn,8); %计算xn的8点DFTXk16=fft(xn,16); %计算xn的16点DFT%以下为绘图部分k=0:7;wk=2*k/8; %产生8点DFT对应的采样点频率(关于归一化值)subplot(2,2,1);stem(wk,abs(Xk8),.); %绘制8点DFT的幅频特性图title(a) 8点DFT的幅频特性图);xlabel(/);ylabel(幅度);grid
9、subplot(2,2,2);stem(wk,angle(Xk8), .); %绘制8点DFT的相频特性图title(b) 8点DFT的相频特性图);grid xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5)k=0:15;wk=2*k/16; %产生16点DFT对应的采样点频率(关于归一化值)subplot(2,2,3);stem(wk,abs(Xk16), .); %绘制16点DFT的幅频特性图title(c) 16点DFT的幅频特性图);xlabel(/);ylabel(幅度);grid subplot(2,2,4);stem(wk,angle(Xk16),
10、.); %绘制16点DFT的相频特性图title(d) 16点DFT的相频特性图);grid xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5)实验4n=0:50;xn=cos(pi/4*n); %输入时域序列向量xn=cos(pi/4.*n)Xk4=fft(xn,4); %计算xn的4点DFTXk8=fft(xn,8); %计算xn的8点DFTXk16=fft(xn,16); %计算xn的16点DFT k=0:3;wk=2*k/4; %产生4点DFT对应的采样点频率(关于归一化值)subplot(3,2,1);stem(wk,abs(Xk4),.); %绘制4点D
11、FT的幅频特性图title(a) 4点DFT的幅频特性图);xlabel(/);ylabel(幅度);grid subplot(3,2,2);stem(wk,angle(Xk4), .); %绘制4点DFT的相频特性图title(b) 4点DFT的相频特性图);grid on;xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5) k=0:7;wk=2*k/8; %产生8点DFT对应的采样点频率(关于归一化值)subplot(3,2,3);stem(wk,abs(Xk8),.); %绘制8点DFT的幅频特性图title(c) 8点DFT的幅频特性图);xlabel(/)
12、;ylabel(幅度);grid subplot(3,2,4);stem(wk,angle(Xk8), .); %绘制8点DFT的相频特性图title(d) 8点DFT的相频特性图) ;grid on;xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5)k=0:15;wk=2*k/16; %产生16点DFT对应的采样点频率(关于归一化值)subplot(3,2,5);stem(wk,abs(Xk16),.); %绘制16点DFT的幅频特性图title(e) 16点DFT的幅频特性图);xlabel(/);ylabel(幅度);gridsubplot(3,2,6);st
13、em(wk,angle(Xk16), .); %绘制16点DFT的相频特性图title(f) 16点DFT的相频特性图); grid on;xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5);实验6n=0:50;xn=10*(0.8.n); %输入时域序列向量xn1=10*(0.8.(-n);xep=(xn+xn1)/2;xop=(xn-xn1)/2;Xk32=fft(xn,32); Xkep=fft(xep,32); Xkop=fft(xop,32); rexk=real(Xk32);imxk=j.*imag(Xk32);k=0:31;wk=2*k/32; sub
14、plot(6,2,1);stem(n,xep,.);title(xep(n)时域);subplot(6,2,2);stem(n,xop,.);title(xop(n)时域);subplot(6,2,3);stem(wk,abs(Xkep),.); title(xep(n) 32点DFT的幅频特性图);xlabel(/);ylabel(幅度);gridsubplot(6,2,4);stem(wk,angle(Xkep), .); %绘制16点DFT的相频特性图title(xep(n) 32点DFT的相频特性图); grid on;xlabel(/);ylabel(相位);axis(0,2,-3.
15、5,3.5);k=0:31;wk=2*k/32; %产生16点DFT对应的采样点频率(关于归一化值)subplot(6,2,5);stem(wk,abs(Xkop),.); %绘制16点DFT的幅频特性图title(xkop 32点DFT的幅频特性图);xlabel(/);ylabel(幅度);gridsubplot(6,2,6);stem(wk,angle(Xkop), .); %绘制16点DFT的相频特性图title(xkop 32点DFT的相频特性图); grid on;xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5);k=0:31;wk=2*k/32; %
16、产生32点DFT对应的采样点频率(关于归一化值)subplot(6,2,7);stem(wk,abs(Xk32),.); %绘制32点DFT的幅频特性图title(Xk32 32点DFT的幅频特性图);xlabel(/);ylabel(幅度);gridsubplot(6,2,8);stem(wk,angle(Xk32),.); %绘制32点DFT的相频特性图title(Xk32 32点DFT的相频特性图); grid on;xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5); subplot(6,2,9);stem(wk,abs(rexk),.);title(shi
17、pu 32点DFT的幅频特性图);xlabel(/);ylabel(幅度);gridsubplot(6,2,10);stem(wk,angle(rexk),.); %绘制32点DFT的相频特性图title(shipu 32点DFT的相频特性图); grid on;xlabel(/);ylabel(相位);axis(0,2,-3.5,3.5); subplot(6,2,11);stem(wk,abs(imxk),.);title(xupu32点DFT的幅频特性图);xlabel(/);ylabel(幅度);gridsubplot(6,2,12);stem(wk,angle(imxk),.); %绘制32点DFT的相频特性图title(xupu 32点DFT的相频特性图); grid on;