数字信号处理Matlab课后实验(吴镇扬)课案.doc

数字信号处理实验报告实验一熟悉 MATLAB 环境实验二信号的采样与重建实验三快速变换及其应用实验四 IIR 数字滤波器的设计实验五 FIR 数字滤波器的设计实验一熟悉 MATLAB 环境一、实验目的(1) 熟悉 MATLAB 的主要操作命令。(2) 学会简单的矩阵输入和数据读写。(3) 掌握简单的绘图命令。(4) 用 MATLAB 编程并学会创建函数。(5) 观察离散系统的频率响应。二、实验内容认真阅读本章附录,在 MATLAB 环境下重新做一遍附录中的例子,体会各条命令的含义。在熟悉了 MATLAB 基本命令的基础上,完成以下实验。上机实验内容: (1) 数组的加、减、乘、除和乘方运算。输入 A=[1 23 4], B=[3 45 6],求 C=A+B , D=A-B , E=A.*B , F=A./B , G=A.^B 并用 stem 语句画出 A、B、C、D、E、F、G。实验程序: A=[1 23 4]; B=[3 45 6]; n=1:4; C=A+B;D=A-B;E=A.*B;F=A./B;G=A.^B; subplot(4,2,1);stem(n,A, 'fill' );xlabel (' 时间序列 n' );ylabel( 'A' ); subplot(4,2,2);stem(n,B, 'fill' );xlabel (' 时间序列 n' );ylabel( 'B' ); subplot(4,2,3);stem(n,C, 'fill' );xlabel (' 时间序列 n' );ylabel( 'A+B' ); subplot(4,2,4);stem(n,D, 'fill' );xlabel (' 时间序列 n' );ylabel( 'A-B' ); subplot(4,2,5);stem(n,E, 'fill' );xlabel (' 时间序列 n' );ylabel( 'A.*B' ); subplot(4,2,6);stem(n,F, 'fill' );xlabel (' 时间序列 n' );ylabel( 'A./B' ); subplot(4,2,7);stem(n,G, 'fill' );xlabel (' 时间序列 n' );ylabel( 'A.^B' ); 运行结果: (2) 用 MATLAB 实现以下序列。 a) x(n)= n0≤n≤ 15 实验程序: n=0:15;x=.^n; stem(n,x, 'fill' ); xlabel (' 时间序列 n' );ylabel( 'x(n)=^n' ); b) x(n)=e (+3j)n 0≤n≤ 15 实验程序: n=0:15;x=exp((+3*j)*n); stem(n,x, 'fill' ); xlabel (' 时间序列 n' );ylabel( 'x(n)=exp((+3*j)*n)' ); 运行结果: a )的时间序列 b )的时间序列 c) x(n)=3cos( π n+ π)+2sin( π n+ π)0≤n≤ 15 实验程序: n=0:1:15; x=3*cos(*pi*n+*pi)+2*sin(*pi*n+*pi); stem(n,x, 'fill' ); xlabel( ' 时间序列 n' ); ylabel( 'x(n)=3*cos(*pi*n+*pi)+2*sin(*pi*n+*pi)' ); 运行结果: d)将c)中的 x(n) 扩展为以 16为周期的函数 x 16 (n)=x(n+16), 绘出四个周期实验程序: n=0:1:63; x=3*cos(*pi*rem(n,16)+*pi)+2*sin(*pi*rem(n,16)+*pi); stem(n,x, 'fill' ); xlabel (' 时间序列 n' );ylabel( 'x16(n)' ); e)将c)中的 x(n) 扩展为以 10为周期的函数 x 10 (n)=x(n+10) ,绘出四个周期实验程序: n=0:1:39; x=3*cos(*pi*rem(n,10)+*pi)+2*sin(*pi*rem(n,10)+*pi); stem(n,x, 'fill' ); xlabel (' 时间序列 n' );ylabel( 'x10(n)' ); 运行结果: d )的时间序列 e )的时间序列(3) x(n)=[1,-1,3,5] ,产生并绘出下列序列的样本。 a)x 1 (n)=2x(n+2)-x(n-1)-2x(n) 实验程序: n=0:3; x=[1 -1 3 5