《西北工业大学Matlab实验报告》由会员分享，可在线阅读，更多相关《西北工业大学Matlab实验报告（34页珍藏版）》请在金锄头文库上搜索。
1、word西北工业大学基于MATLAB的数字信号处理实验报告学 院： 计算机学院 学号：姓名：专 业： 计算机科学与技术 西北工业大学2017年 07月实验一 MATLAB基本编程实验一、实验目的及要求1熟悉MATLAB运行环境；2. 掌握MATLAB的基本语法和函数；3. 掌握MATLAB的基本绘图功能二、实验设备（环境）及要求1 计算机2 Matlab软件编程实验平台三、实验容与步骤1 设求使2 . 分别用for或while循环结构编写程序，求出 。并考虑一种避免循环语句的程序设计算法实现同样的运算。3. 将图形窗口4等分，按顺序绘制出下列函数的曲线图及曲面图并标记。四、设计思想 1.，则X
2、=B*inv（A） 2.编写循环程序，循环105次，每次循环中计算的值对做数学运算，可以得到3.使用subplot将图形窗口分为四个子窗口，在子窗口中进行绘图五、程序代码及注释% question 1*% 计算矩阵B*inv(A)A= 0 2 1 2 -1 3 -3 3 -4;B= 1 2 3 2 -3 1;X=B*inv(A)%*% question2*% 分别使用for及while循环计算矩阵K的值K=0;for i=1:105 K=K+sqrt(3)/2(i);endK%-K=0;i=1;while i106 K=K+sqrt(3)/2(i); i=i+1;endK%*% questio
3、n 3*% 首先计算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=meshgrid(-2:0.1:2,-4:0.4:4);Z=X.2/22+Y.2/42;subplot(2,2,1)plot(x,y1)axis(0 6.5 min(y1)-0.1 max(y1)+0.1);xlabel(x)ylabel(y1)title(x-y1 graph)subplot(2,2,2)plot(x,y2)axis(0 6.5 min(y2)-0.1 max(y2)+0.1);xlabel(x)ylab
4、el(y2)title(x-y2 graph)subplot(2,2,3)plot(x,y3)axis(0 6.5 min(y3)-0.1 max(y3)+0.1);xlabel(x)ylabel(y3)title(x-y3 graph)subplot(2,2,4)surf(X,Y,Z);axis(-2 2 -4 4 0 2);xlabel(X)ylabel(Y)zlabel(Z)title(X-Y-Z graph)%*六、实验结果1.3.实验二 MATLAB基本编程实验一、实验目的及要求1. 回顾数字信号处理的主要容；2. 掌握利用MATLAB进行信号处理的方法；3. 了解信号处理工具箱中一
5、些函数的功能；二、实验设备（环境）及要求1 计算机2 Matlab软件编程实验平台三、实验容1 .任何实数序列x(n)都能分解成为它的偶部分量和奇部分量之和，编写一个MATLAB函数文件，该函数的功能是将一给定序列分解成为它的偶部分量和奇部分量。并通过M文件或命令行调用该函数文件将以下序列分解成为偶部分量和奇部分量。2. 一数字滤波器由 频率响应函数描述 1）写出其差分方程表示； 2）画出上面滤波器的幅频和相频图； 3）产生信号 的200个样本，通过该滤波器得到输出y(n), 试将输出y(n)的稳态部分与x(n)作比较，说明这两个正弦信号的幅度和相位是如何受该滤波器影响的。3. 设计带通滤波器
6、（Butterworth、椭圆、窗函数），采样率fs2000Hz, 通带围为300 -600Hz，阶数自选, 画出频率响应并分析比较。四、设计思想 1.首选对原始序列做折叠处理，得到折叠后的序列，再利用计算公式求的序列的奇偶分量 2.差分方程如下根据频率响应函数，使用freqz函数计算幅频和相频图取样出x(n)的200个样本，通过filter函数求出滤波器作用后的y(n),对比输入输出信号，进行分析 3.首先设定带通滤波器阶数，本实验设定为10阶，窗函数部分使用了矩形窗，然后根据通带围频率求得各个滤波器设计所需的参数，然后调用matlab滤波器生成函数，生成相应的滤波器。五、程序代码及注释代码
7、1function xe,xo,m=Calculatexexo(x,n)%输入一个序列，输出分解后的奇偶序列%2017-8-19 by Binwangm=-fliplr(n);n2=min(n,m);n3=max(n,m);m=n2:n3; %计算折叠处理之后的横轴坐标I=find(m=n(1);x1=zeros(1,length(m);x1(I:I+length(n)-1) = x; %在新的横轴坐标域表示原始序列x=x1;xe=0.5*(x+fliplr(x);xo=0.5*(x-fliplr(x); %计算奇偶分量end%计算序列的奇偶分量 n=0:100; x=exp(-0.05*n)
8、.*sin(0.1*pi*n+pi/3); %生成序列 xe,xo,m=Calculatexexo(x,n); %计算序列奇偶分量 subplot(2,1,1) stem(n,x,ro); %原始序列 title(original graph) subplot(2,2,3) stem(m,xe,*); %偶分量 title(xe graph) subplot(2,2,4) stem(m,xo,*); %奇分量 title(xo graph)fs=1000;b=1 0 0 0 1; %表达式分子系数a=1 0 0 0 -0.8145 %表达式分母系数H,w=freqz(b,a,256,fs);
9、mag=abs(H); %计算幅值ph=angle(H); %计算相位ph=ph*180/pi;subplot(2,1,1)plot(w,mag);grid; xlabel(frequency Hz);ylabel(magnitude);subplot(2,1,2)plot(w,ph);grid;xlabel(frequency Hz);ylabel(phase);n=1:0.025:5; x=sin(pi*n/2)+5*cos(pi*n); stem(n,x) b=1 0 0 0 1; %表达式分子系数 a=1 0 0 0 -0.8145; %表达式分母系数 y=filter(b,a,x);
10、 hold on stem(n,real(y),r*) legend(x(n),y(n); title(输入输出序列对比)代码3fs=2000; fc1=300; fc2=600; N=10; %阶数设置为10 wlp=fc1/(fs/2); whp=fc2/(fs/2); wn=wlp,whp; %计算带通围%矩形窗带通滤波器 w1=boxcar(N); %矩形窗时域响应 b1=fir1(N-1,wn,w1); h,w=freqz(b1,1,256,fs); subplot(2,3,1) plot(w,20*log10(abs(h); xlabel(frequency Hz); ylabel
11、(magnitude /dB); title(矩形窗函数的幅频特性); subplot(2,3,4) plot(w,180/pi*unwrap(angle(h); xlabel(frequency Hz); ylabel(phase); title(矩形窗函数的相频特性);%巴特沃斯带通滤波器 b,a=butter(N,wn); h1,w1=freqz(b,a,256,fs); subplot(2,3,2) plot(w1,20*log10(abs(h1)/max(abs(h1); xlabel(frequency Hz); ylabel(magnitude /dB); title(巴特沃斯的
12、幅频特性); subplot(2,3,5) plot(w1,180/pi*(angle(h1); xlabel(frequency Hz); ylabel(phase); title(巴特沃斯的相频特性);%椭圆带通滤波器 rs=60;rp=0.1; %椭圆带通滤波器的指标 B,A=ellip(N,rp,rs,wn) h2,w2=freqz(B,A,256,fs); subplot(2,3,3) plot(w2,20*log10(abs(h2)/max(abs(h2); xlabel(frequency Hz); ylabel(magnitude /dB); title(椭圆滤波器的幅频特性); subplot(2,3,6) plot(w2,180/pi*(angle(h2); xlabel(frequency Hz); ylabel(phase); title(椭圆滤波器的相频特性);六、实验结果1原始序列及得到的奇偶序列如下图所示计算出的幅频和相频图如下2.3 输入序列与滤波器的输出序列如下图所示分析结果可知，输出相对于输入在相位上有一定的延迟，幅度上有了提升。两个正弦信号在对应的频率处均有一定给上的幅度增加和相位延迟3带通滤波器设计根据上图结构可以看出：巴特沃斯带通滤波器的通带无波动，且设计出来
