资源描述：
《西北工业大学Matlab实验资料报告材料.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、西北工业大学《基于MATLAB的数字信号处理》实验报告学院：计算机学院学 号：姓 名：专业：计算机科学与技术西北工业大学2017年07月实验一MATLAB基本编程实验一、实验目的及要求1．熟悉MATLAB运行环境；2.掌握MATLAB的基本语法和函数；3.掌握MATLAB的基本绘图功能二、实验设备（环境）及要求1．计算机2．Matlab软件编程实验平台三、实验容与步骤1．设求使2.分别用for或while循环结构编写程序，求出。并考虑一种避免循环语句的程序设计算法实现同样的运算。3.将图形窗口4等分，
2、按顺序绘制出下列函数的曲线图及曲面图并标记。四、设计思想1.，则X=B*inv（A）2.编写循环程序，循环105次，每次循环中计算的值对做数学运算，可以得到3.使用subplot将图形窗口分为四个子窗口，在子窗口中进行绘图五、程序代码及注释%question1*************%计算矩阵B*inv(A)A=[0212-13-33-4];B=[1232-31];X=B*inv(A)%*************************%question2***************%分别使用for及
3、while循环计算矩阵K的值K=0;fori=1:105K=K+sqrt(3)/2^(i);endK%-------------------------K=0;i=1;whilei<106K=K+sqrt(3)/2^(i);i=i+1;endK%*************************%question3*************%首先计算x,y1,y2,y3,X,Y,Z的值，然后画图x=0:0.1:2*pi;y1=cos(x);y2=2*sin(2*x);y3=x.*cos(x);[X,Y]=
4、meshgrid(-2:0.1:2,-4:0.4:4);Z=X.^2/2^2+Y.^2/4^2;subplot(2,2,1)plot(x,y1)axis([06.5min(y1)-0.1max(y1)+0.1]);xlabel('x')ylabel('y1')title('x-y1graph')subplot(2,2,2)plot(x,y2)axis([06.5min(y2)-0.1max(y2)+0.1]);xlabel('x')ylabel('y2')title('x-y2graph')subplot
5、(2,2,3)plot(x,y3)axis([06.5min(y3)-0.1max(y3)+0.1]);xlabel('x')ylabel('y3')title('x-y3graph')subplot(2,2,4)surf(X,Y,Z);axis([-22-4402]);xlabel('X')ylabel('Y')zlabel('Z')title('X-Y-Zgraph')%*************************六、实验结果1.2.K=1.73213.实验二MATLAB基本编程实验一、实验目的
6、及要求1.回顾数字信号处理的主要容；2.掌握利用MATLAB进行信号处理的方法；3.了解信号处理工具箱中一些函数的功能；二、实验设备（环境）及要求1．计算机2．Matlab软件编程实验平台三、实验容1.任何实数序列x(n)都能分解成为它的偶部分量和奇部分量之和，编写一个MATLAB函数文件，该函数的功能是将一给定序列分解成为它的偶部分量和奇部分量。并通过M文件或命令行调用该函数文件将以下序列分解成为偶部分量和奇部分量。2.一数字滤波器由频率响应函数描述1）写出其差分方程表示；2）画出上面滤波器的幅频和相频
7、图；3）产生信号的200个样本，通过该滤波器得到输出y(n),试将输出y(n)的稳态部分与x(n)作比较，说明这两个正弦信号的幅度和相位是如何受该滤波器影响的。3.设计带通滤波器（Butterworth、椭圆、窗函数），采样率fs＝2000Hz,通带围为300-600Hz，阶数自选,画出频率响应并分析比较。四、设计思想1.首选对原始序列做折叠处理，得到折叠后的序列，再利用计算公式求的序列的奇偶分量2.差分方程如下根据频率响应函数，使用freqz函数计算幅频和相频图取样出x(n)的200个样本，通过filt
8、er函数求出滤波器作用后的y(n),对比输入输出信号，进行分析3.首先设定带通滤波器阶数，本实验设定为10阶，窗函数部分使用了矩形窗，然后根据通带围频率求得各个滤波器设计所需的参数，然后调用matlab滤波器生成函数，生成相应的滤波器。五、程序代码及注释代码1function[xe,xo,m]=Calculatexexo(x,n)%输入一个序列，输出分解后的奇偶序列%2017-8-19byBinwangm=-fliplr(n)
