数字信号处理实验.zip

用matlab实现fft,fir,iir的数学建模，设计滤波器的参数指标。

x1n=[ones（14）];     %产生R4（n）序列向量
X1k8=fft（x1n8）;     %计算x1n的8点DFT
X1k16=fft（x1n16）;   %计算x1n的16点DFT
%以下绘制幅频特性曲线
N=8;
f=2/N*（0:N-1）;
figure（1）;
subplot（121）;stem（fabs（X1k8）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（1a） 8点DFT[x_1（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
N=16;
f=2/N*（0:N-1）;
subplot（122）;stem（fabs（X1k16）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（1a） 16点DFT[x_1（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
%x2n 和 x3n
M=8;xa=1:（M/2）;  xb=（M/2）:-1:1; 
x2n=[xaxb];    %产生长度为8的三角波序列x2（n）
x3n=[xbxa];
X2k8=fft（x2n8）;
X2k16=fft（x2n16）;
X3k8=fft（x3n8）;
X3k16=fft（x3n16）;
figure（2）;
N=8;
f=2/N*（0:N-1）;
subplot（221）;stem（fabs（X2k8）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（2a） 8点DFT[x_2（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
subplot（223）;stem（fabs（X3k8）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（3a） 8点DFT[x_3（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
N=16;
f=2/N*（0:N-1）;
subplot（222）;stem（fabs（X2k16）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（2a） 16点DFT[x_2（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
subplot（224）;stem（fabs（X3k16）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（3a） 16点DFT[x_3（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
%x4n 和 x5n
N=8;n=0:N-1;
x4n=cos（pi*n/4）;
x5n=cos（pi*n/4）+cos（pi*n/8）;
X4k8=fft（x4n8）;
X4k16=fft（x4n16）;
X5k8=fft（x5n8）;
X5k16=fft（x5n16）;
figure（3）;
N=8;
f=2/N*（0:N-1）;
subplot（221）;stem（fabs（X4k8）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（4a） 8点DFT[x_4（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
subplot（223）;stem（fabs（X5k8）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（5a） 8点DFT[x_5（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
N=16;
f=2/N*（0:N-1）;
subplot（222）;stem（fabs（X4k16）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（4a） 16点DFT[x_4（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
subplot（224）;stem（fabs（X5k16）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（5a） 16点DFT[x_5（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;

%x8n
Fs=64; T=1/Fs; 
N=16;n=0:N-1; %对于N=16的情况
nT = n*T;
x8n=cos（8*pi*nT）+cos（16*pi*nT）+cos（20*pi*nT）
X8k16=fft（x8n16）;
N=16;
f=2/N*（0:N-1）;
figure（4）;
subplot（221）;stem（fabs（X8k16）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（8a） 16点DFT[x_8（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
N=32;n=0:N-1; %对于N=16的情况
nT = n*T;
x8n=cos（8*pi*nT）+cos（16*pi*nT）+cos（20*pi*nT）
X8k32=fft（x8n32）;
N=32;
f=2/N*（0:N-1）;
subplot（222）;stem（fabs（X8k32）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（8a） 32点DFT[x_8（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;
N=64;n=0:N-1; %对于N=16的情况
nT = n*T;
x8n=cos（8*pi*nT）+cos（16*pi*nT）+cos（20*pi*nT）
X8k64=fft（x8n64）;
N=64;
f=2/N*（0:N-1）;
subplot（223）;stem（fabs（X8k64）‘.‘）; %绘制8点DFT的幅频特性图
title（‘（8a） 64点DFT[x_8（n）]‘）;xlabel（‘ω/π‘）;ylabel（‘幅度‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-06-16 13:52  数字信号处理实验\
     目录           0  2018-12-31 10:14  数字信号处理实验\实验一\
     文件        1390  2018-12-16 19:15  数字信号处理实验\实验一\shiyan1_1.asv
     文件        2718  2018-12-16 19:22  数字信号处理实验\实验一\shiyan1_1.m
     文件      184320  2018-12-16 19:40  数字信号处理实验\实验一\实验一截图.doc
     文件      193024  2018-12-16 19:21  数字信号处理实验\实验一\实验一：用FFT对信号作频谱分析_实验报告.doc
     目录           0  2018-12-31 10:14  数字信号处理实验\实验三\
     文件         269  2018-12-16 19:56  数字信号处理实验\实验三\myplot.m
     文件         132  2018-12-27 11:14  数字信号处理实验\实验三\shiyan3_1.m
     文件         149  2018-12-27 13:57  数字信号处理实验\实验三\shiyan3_2.m
     文件         368  2018-12-27 14:35  数字信号处理实验\实验三\shiyan3_3.m
     文件         224  2018-12-16 19:54  数字信号处理实验\实验三\tplot.m
     文件        1317  2018-12-17 13:22  数字信号处理实验\实验三\xtg.m
     文件       73728  2018-12-17 13:17  数字信号处理实验\实验三\实验三--FIR数字滤波器设计与软件实现.doc
     文件      201216  2018-12-27 14:37  数字信号处理实验\实验三\实验三截图.doc
     目录           0  2018-12-31 10:14  数字信号处理实验\实验二\
     文件         559  2018-12-16 19:54  数字信号处理实验\实验二\mstg.m
     文件         269  2018-12-16 19:56  数字信号处理实验\实验二\myplot.m
     文件         873  2018-12-16 20:25  数字信号处理实验\实验二\shiyan2.m
     文件         224  2018-12-16 19:54  数字信号处理实验\实验二\tplot.m
     文件      203776  2018-12-16 20:26  数字信号处理实验\实验二\实验二截图.doc
     文件      228238  2018-12-16 20:26  数字信号处理实验\实验二\实验二：IIR数字滤波器设计及软件实现.doc
     目录           0  2018-12-31 10:14  数字信号处理实验\实验四\
     文件         287  2018-12-27 15:57  数字信号处理实验\实验四\shiyan4_1.m
     文件         149  2018-12-27 16:22  数字信号处理实验\实验四\shiyan4_2.m
     文件       79872  2018-12-27 16:23  数字信号处理实验\实验四\实验四截图.doc