《基于MATLAB的数字信号处理》实验报告课件.docVIP

西北工业大学 《基于MATLAB的数字信号处理》 学 院： 学　　号： 姓　　名： 专 业： 实验时间： 2014年7月 指导教师： 西北工业大学 2014年 07 月  实验一 MATLAB基本1．熟悉MATLAB；1．2． 1． 2. 分别用for或while循环结构编写程序，求出 。并考虑一种避免循环语句的程序设计算法实现同样的运算。 3. 在同一坐标系下绘制以下3条曲线,并作标记。 四、设计思想 设计思想1.1： 对于AX=B类型的方程组求解，需要先对非奇异矩阵A求逆，若存在逆矩阵，然后与B相乘即为所求 设计思想1.2： 用for或while循环结构编写程序时，首先对K初始化为0,确定循环的起始条件、步长和终值条件；用非循环语句时，先确定首项、公比，然后带入到等比数列求和公式中 设计思想1.3： 首先确定绘制曲线的每两点之间的间隔，然后在一幅图上绘制出所有曲线，并进行标注 五、程序代码及注释 程序代码1： A=[6 3 4; -2 5 7; 8 -4 -3]; %矩阵A X=inv(A) X = 0.1300 -0.0700 0.0100 0.5000 -0.5000 -0.5000 -0.3200 0.4800 0.3600 B=[3 -4 -7]; X=X*B 程序代码2： 方法1：for循环结构实现 %结果初始化0 K= 0.0; for i = 1: 106 %默认步长为1 K=K+sqrt(3)/(2^i); end; fprintf(K=%f\n,K); 方法2：while循环结构实现 K = 0.0; i = 1; while i 107 K = K +sqrt(3)/(2^i); i = i + 1; end; fprintf(K=%f\n,K); 方法3：不使用循环，利用等比数列的前106项和公式实现 K = sqrt(3)/2.0; %数列首项 K = K *(1-0.5^106)/(1-0.5);%公比为0.5 fprintf(K=%f\n,K); 程序代码3： x=0:0.01:pi; %步长0.01 y1=sin(x); y2=sin(x).*sin(10*x); y3=-cos(x); grid on%显示网格 hold on %在一幅图中绘制多条曲线 plot(x, y1, r,LineWidth,2); %y1红线表示，线宽为2 plot(x, y2, g,LineWidth,2); %y2绿线表示，线宽为2 plot(x, y3, b,LineWidth,2); %y3蓝线表示，线宽为2 %x轴范围0~3.5，y轴范围-1.2~1.2 axis([0,3.5,-1.2,1.2]); %图例注释自动放在最佳位置 legend(sin(x),sin(x)sin(10x),-cos(x)); xlabel(x坐标轴); %标注横坐标 ylabel(y坐标轴); %标注纵坐标 title(三曲线对比图); %标注标题 六、实验结果 实验结果1： X = 0.6000 7.0000 -5.4000 所以，X=[]=[ 0.6000, 7.0000, -5.4000]’ 实验结果2： K=1.732051 实验结果3： 三曲线的对比图如下所示： 图1.1 三曲线的对比 实验二 基于MATLAB信号处理实验 一、实验目的及要求 1. 回顾数字信号处理的主要内容； 2. 掌握利用MATLAB进行信号处理的方法； 3. 了解信号处理工具箱中一些函数的功能； 二、实验设备（环境）及要求 1．2． 1． 设序列 其中， 2 . 一数字滤波器由 频率响应函数描述  1）写出其差分方程表示； 2）画出上面滤波器的幅频和相频图； 3）产生信号 的200个样本，通过该滤波器得到输出y(n), 试将输出y(n)的稳态部分与x(n)作比较，说明这两个正弦信号的幅度和相位是如何受该滤波器影响的。 3. 设计F