数字信号处理实验教案

（1112341用） 数字信号解决实验规 预备实验用于信号解决的MATLAB基本操 实验 系统响应及系统稳定 实验 时域采样与频域采 实验 用FFT对信号作频谱分 实验 IIR数字滤波器设计及软件实 实验 FIR数字滤波器设计与软件实 附录1：MATLAB编程及绘图办 附录2：MATLAB矩阵及矩阵的运 实验基本规定（以软件设计为重要体现形式预备实验MATLAB（3实验 系统响应及系统稳定性（3学时实验 时域采样与频域采样（3学时实验 用FFT对信号作频谱分析（3学时实验 IIR数字滤波器设计及软件实现（3学时实验 FIR数字滤波器设计与软件实现（3学时姓名 学号 实验日期 [1]高西全，丁玉美.数字信号解决（第三版）西电出版社张志勇等.MATLABRa.北京航空航天大学出版社.MATLAB[3]程佩清.数字信号解决教程[M].清华大学出版社.预备实验用于信号解决的MATLABMATLAB②MATLAB①熟悉MATLAB运用冒号（：）X1=[1234X2=[1.0001.5002.000X3=[5432C）分别输入矩阵123456789D）分别计算A+B，A-B，A+3，A-4，A*B，A.*B，C=inv（A）,A/B,A./BE）分别计算sin(x1),cos(x1),exp(x1),lg(x2),sqrt(x2)运用冒号（：）X=j:kX=[j,j+1,j+2,…,k-I>0j<kI<0j>k,则生成向量I<0j>k,则生成向量例X1=[1234X2=[1.0001.5002.0002.500X3=5：-432 0B= n*n0B= m*n0B=Zeros(Size(A)) A0注：m,n二矩阵的算术运算A-对于矩阵和标量（一种数）A-

Sin(A),cos(A),exp(A),log(A)( 式中A能够是标量也能够是矩阵练习xlabel('0-四．绘图（二维绘图1plot(y);%绘持续曲线图。stem(y);%绘离散序列波形yyy中元2plot(x,y);绘持续曲线图。stem(x,y);%绘离散序列波形x,yxy元素的值为纵坐标，一一MATLAB会自动设定所画曲线的颜色和线型。可自动设立，在输入变量组后，加一叉标，-.点划线，下面给出窗口分割语句x，yx*ysubplot(3,2,2)3*21MATLAB2MATLAB1：MATLABMATLABMATLABC语言：以字母开头，由字母、下31个。MATLAB惯用运算符（1＋—矩阵右除(矩阵元素右除矩阵左除（矩阵元素左除矩阵乘（矩阵元素乘∧（.矩阵求幂（矩阵元素求幂[%‘><&|~MATLABlog(x)、log2(x0、e（2、10）[1,2,…]=函数名（1MATLAB3n1*n2title(‘字符串xlable(‘字符串xxlable(‘字符串y1yyy中元2x,yxy元素的值为纵坐标，一一MATLAB>>Y=[0.7,0.5,0.2,->>Y=[0.7,0.5,0.2,- %画连线图, %画散点图,图1.1持续波形 图1.2散点波形;..]holdon2holdholdon3plot 4plotyy命令，设有两个纵坐标1：figure（n）：MATLABplotplotfigure命令打开新的图形窗口。格式2：clf：去除现在图形窗的内容3：close:4：closeall:subplot(2,2,4);fill(t,y,’r’);subplot(2,2,4);fill(t,y,’r’);MATLAB与大多数计算机语言同样，MATLABMATLAB的次序构造就是由体现式构成的语句，语句之间用逗号或分号隔开，这些MATLABforwhilefor-endforfor循环变量＝1212313之间，如果条件满（DFT

jNX(k)x(n)Wnk NN

k

,N fork=0:N-1forn=0:N-while-endwhilewhile3if-else-endififswitchswitchcase体现式的值＝1为真时的语句组case值2体现式的值＝23forM文献的运行时间,>>elapsed_time>>elapsed_time=00.38forfork=1:100在命令窗口运行,>>elapsed_time >>elapsed_time

2：MATLAB矩阵须用中括号[]>>x=[12345678];%一维1×8>>x=0:0.02:1%以起始值=0、增量值=0.02、终止值=1的矩阵(阵列)%linspace函数，起始值=0,终止值=1,以等间隔（1-0）/100100A（i,j）Aij例1．A>>>>>>A162738495MATLABA（u,v）A的子矩阵，u,vA([1,3],:)%A1ansA(:,[2,4,5]A2,4,5ans A（L1，：）表达提取A的L1指定的行、全部列；A（：，L2A([1,3],:)%A1ansA(:,[2,4,5]A2,4,5ans 123457930.5>>rand(>>rand(‘seed’,0);>>x=rand(1,10);%1X10>>xColumns1through Columns7throughx4A2,4,6>>>>>>A(:,[2,4,6])=[]A= 5>>A 35>>A 35079000000使用，这时使用*.mat>>n=1:;x=cos(0.2*pi*n);>>savec:\temp.matx>>load save途径\文献名 >>n=1:;x=cos(0.2*pi*n);>>savec:\temp.matx>>loadMATLABMATLAB命令窗口定义的变量、运算成果和有关数据，此内存空间就是workspaceworkspace*.matFilesaveworkspaceas….项即可矩阵 矩阵 a-ABAB矩阵的除法从逆矩阵的概念引申，