《基于matlab的数字信号处理》实验报告

《《基于matlab的数字信号处理》实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、西北工业大学《基于MATLAB的数字信号处理》实验报告学院：学　　号：姓　　名：专业：实验时间：2014年7月指导教师：西北工业大学2014年07月57实验一MATLAB基本编程实验一、实验目的及要求1．熟悉MATLAB运行环境；2.掌握MATLAB的基本语法和函数；3.掌握MATLAB的基本绘图功能二、实验设备（环境）及要求1．计算机2．Matlab软件编程实验平台三、实验内容与步骤1．求下列线性方程组的解2.分别用for或while循环结构编写程序，求出。并考虑一种避免循环语句的程序设计算法实现同样的运算。3.在同一坐标系下绘制以下3条曲线

2、,并作标记。四、设计思想设计思想1.1：对于AX=B类型的方程组求解，需要先对非奇异矩阵A求逆，若存在逆矩阵，然后与B相乘即为所求设计思想1.2：57用for或while循环结构编写程序时，首先对K初始化为0,确定循环的起始条件、步长和终值条件；用非循环语句时，先确定首项、公比，然后带入到等比数列求和公式中设计思想1.3：首先确定绘制曲线的每两点之间的间隔，然后在一幅图上绘制出所有曲线，并进行标注五、程序代码及注释程序代码1：>>A=[634;-257;8-4-3];%矩阵A>>X=inv(A)X=0.1300-0.07000.01000.50

3、00-0.5000-0.5000-0.32000.48000.3600>>B=[3-4-7]';X=X*B程序代码2：方法1：for循环结构实现%结果初始化0K=0.0;fori=1:106%默认步长为1K=K+sqrt(3)/(2^i);end;fprintf('K=%f',K);方法2：while循环结构实现K=0.0;57i=1;whilei<107K=K+sqrt(3)/(2^i);i=i+1;end;fprintf('K=%f',K);方法3：不使用循环，利用等比数列的前106项和公式实现K=sqrt(3)/2.0;%数列首项

4、K=K*(1-0.5^106)/(1-0.5);%公比为0.5fprintf('K=%f',K);程序代码3：x=0:0.01:pi;%步长0.01y1=sin(x);y2=sin(x).*sin(10*x);y3=-cos(x);gridon%显示网格holdon%在一幅图中绘制多条曲线plot(x,y1,'r','LineWidth',2);%y1红线表示，线宽为2plot(x,y2,'g','LineWidth',2);%y2绿线表示，线宽为2plot(x,y3,'b','LineWidth',2);%y3蓝线表示，线宽为2%x轴范围

5、0~3.5，y轴范围-1.2~1.2axis([0,3.5,-1.2,1.2]);%图例注释自动放在最佳位置legend('sin(x)','sin(x)sin(10x)','-cos(x)');xlabel('x坐标轴');%标注横坐标ylabel('y坐标轴');%标注纵坐标title('三曲线对比图');%标注标题六、实验结果57实验结果1：X=0.60007.0000-5.4000所以，X=[]=[0.6000,7.0000,-5.4000]’实验结果2：K=1.732051实验结果3：三曲线的对比图如下所示：图1.1三曲线的对比57实

6、验二基于MATLAB信号处理实验一、实验目的及要求1.回顾数字信号处理的主要内容；2.掌握利用MATLAB进行信号处理的方法；3.了解信号处理工具箱中一些函数的功能；二、实验设备（环境）及要求1．计算机2．Matlab软件编程实验平台三、实验内容1．设序列其中，2.一数字滤波器由频率响应函数描述1）写出其差分方程表示；2）画出上面滤波器的幅频和相频图；3）产生信号的200个样本，通过该滤波器得到输出y(n),试将输出y(n)的稳态部分与x(n)作比较，说明这两个正弦信号的幅度和相位是如何受该滤波器影响的。3.设计FIR滤波器（带通滤波器，四种窗

7、函数，fs＝2000Hz,带通频率300Hz—600Hz，阶数自选,画频率特性并分析比较）。四、设计思想设计思想2.1：将求互相关问题转化为求卷积的问题，互相关的区间大小为两信号的区间端点分别相加，区间长度为两信号区间长度之和设计思想2.2:由系统的差分方程57，得到方程的各个项的系数，进而可求出系统的幅频和相频响应设计思想2.3：首先选定带通滤波器的阶数，根据采样频率、上、下限截止频率求得滤波器的各个参数，然后得出窗函数的频率特性五、程序代码及注释程序代码2.1：clc;clearall;x=[3,1,7,0,-1,4,2];%冲击信号x的幅

8、度nx=[-3:3];%冲击信号x的取值范围%对x进行右移2个单位，得到x(n-2)ny=nx+2;y=x;w=normrnd(0,1,[1,7]);