数字信号处理 实验.docx

数字信号处理 实验.docx

《数字信号处理 实验.docx》由会员分享，可在线阅读，更多相关《数字信号处理 实验.docx（9页珍藏版）》请在冰豆网上搜索。

数字信号处理实验

实验1

xn=[1111];%输入时域序列向量xn=R4（n）

Xk8=fft（xn,8）;%计算xn的8点DFT

Xk16=fft（xn,16）;%计算xn的16点DFT

Xk32=fft（xn,32）;%计算xn的32点DFT

k=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）,'.'）;%绘制8点DFT的相频特性图

title（'（b）8点DFT的相频特性图'）;gridon;

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的相频特性图'）;gridon;

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（'幅度'）;grid

subplot（3,2,6）;stem（wk,angle（Xk32）,'.'）;%绘制32点DFT的相频特性图

title（'（f）32点DFT的相频特性图'）;gridon;

xlabel（'ω/π'）;ylabel（'相位'）;axis（[0,2,-3.5,3.5]）;

实验2

xn=[12344321];%输入时域序列向量x2

Xk8=fft（xn,8）;%计算xn的8点DFT

Xk16=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）,'.'）;%绘制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的相频特性图'）;grid

xlabel（'ω/π'）;ylabel（'相位'）;axis（[0,2,-3.5,3.5]）

实验3

xn=[43211234];%输入时域序列向量x2

Xk8=fft（xn,8）;%计算xn的8点DFT

Xk16=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）,'.'）;%绘制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的相频特性图'）;grid

xlabel（'ω/π'）;ylabel（'相位'）;axis（[0,2,-3.5,3.5]）

实验3

xn=[43211234];%输入时域序列向量x2

Xk8=fft（xn,8）;%计算xn的8点DFT

Xk16=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）,'.'）;%绘制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的相频特性图'）;grid

xlabel（'ω/π'）;ylabel（'相位'）;axis（[0,2,-3.5,3.5]）

实验4

n=0:

50;

xn=cos（pi/4*n）;%输入时域序列向量xn=cos（pi/4.*n）

Xk4=fft（xn,4）;%计算xn的4点DFT

Xk8=fft（xn,8）;%计算xn的8点DFT

Xk16=fft（xn,16）;%计算xn的16点DFT

k=0:

3;wk=2*k/4;%产生4点DFT对应的采样点频率（关于π归一化值）

subplot（3,2,1）;stem（wk,abs（Xk4）,'.'）;%绘制4点DFT的幅频特性图

title（'（a）4点DFT的幅频特性图'）;xlabel（'ω/π'）;ylabel（'幅度'）;grid

subplot（3,2,2）;stem（wk,angle（Xk4）,'.'）;%绘制4点DFT的相频特性图

title（'（b）4点DFT的相频特性图'）;gridon;

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（'ω/π'）;ylabel（'幅度'）;grid

subplot（3,2,4）;stem（wk,angle（Xk8）,'.'）;%绘制8点DFT的相频特性图

title（'（d）8点DFT的相频特性图'）;gridon;

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（'幅度'）;grid

subplot（3,2,6）;stem（wk,angle（Xk16）,'.'）;%绘制16点DFT的相频特性图

title（'（f）16点DFT的相频特性图'）;gridon;

xlabel（'ω/π'）;ylabel（'相位'）;axis（[0,2,-3.5,3.5]）;

实验6

n=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;

subplot（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（'幅度'）;grid

subplot（6,2,4）;stem（wk,angle（Xkep）,'.'）;%绘制16点DFT的相频特性图

title（'xep（n）32点DFT的相频特性图'）;gridon;

xlabel（'ω/π'）;ylabel（'相位'）;axis（[0,2,-3.5,3.5]）;

k=0:

31;wk=2*k/32;%产生16点DFT对应的采样点频率（关于π归一化值）

subplot（6,2,5）;stem（wk,abs（Xkop）,'.'）;%绘制16点DFT的幅频特性图

title（'xkop32点DFT的幅频特性图'）;xlabel（'ω/π'）;ylabel（'幅度'）;grid

subplot（6,2,6）;stem（wk,angle（Xkop）,'.'）;%绘制16点DFT的相频特性图

title（'xkop32点DFT的相频特性图'）;gridon;

xlabel（'ω/π'）;ylabel（'相位'）;axis（[0,2,-3.5,3.5]）;

k=0:

31;wk=2*k/32;%产生32点DFT对应的采样点频率（关于π归一化值）

subplot（6,2,7）;stem（wk,abs（Xk32）,'.'）;%绘制32点DFT的幅频特性图

title（'Xk3232点DFT的幅频特性图'）;xlabel（'ω/π'）;ylabel（'幅度'）;grid

subplot（6,2,8）;stem（wk,angle（Xk32）,'.'）;%绘制32点DFT的相频特性图

title（'Xk3232点DFT的相频特性图'）;gridon;

xlabel（'ω/π'）;ylabel（'相位'）;axis（[0,2,-3.5,3.5]）;

subplot（6,2,9）;stem（wk,abs（rexk）,'.'）;

title（'shipu32点DFT的幅频特性图'）;xlabel（'ω/π'）;ylabel（'幅度'）;grid

subplot（6,2,10）;stem（wk,angle（rexk）,'.'）;%绘制32点DFT的相频特性图

title（'shipu32点DFT的相频特性图'）;gridon;

xlabel（'ω/π'）;ylabel（'相位'）;axis（[0,2,-3.5,3.5]）;

subplot（6,2,11）;stem（wk,abs（imxk）,'.'）;

title（'xupu32点DFT的幅频特性图'）;xlabel（'ω/π'）;ylabel（'幅度'）;grid

subplot（6,2,12）;stem（wk,angle（imxk）,'.'）;%绘制32点DFT的相频特性图

title（'xupu32点DFT的相频特性图'）;gridon;