实验一 Matlab软件的使用.docVIP

实验一 Matlab软件的使用 1、实验目的和要求 ㈠ 掌握Matlab软件的使用方法； ㈡ 常用离散时间信号的产生、显示和运算。 2、实验内容 ㈠ Matlab软件常用命令和工具箱的使用，学会简单的编程； ㈡ 编程实现常用离散时间信号； ㈢ 离散时间信号的叠加、移位、线性卷积等基本运算。 3、实验原理 利用软件生成数字信号处理系统中所涉及的信号及信号的加、移位、卷积等运算。 4、主要仪器设备 计算机，Matlab软件。 实验步骤： 程序1： n=-5:30; x=[zeros(1,5),1,zeros(1,30)]; stem(n,x,fill); grid on 程序2： n=-5:30; x=[zeros(1,5),ones(1,8),zeros(1,23)]; stem(n,x,fill); grid on； 程序3： n=-15:20; x=[zeros(1,15),1,2,3,4,5,zeros(1,16)]; stem(n,x,fill); grid on; 程序4： N=10; M=10; L=N+M-1; x=[1,2,3,4,5,zeros(1,5)]; h=[1,2,1,2,zeros(1,6)]; y=conv(x,h); n=0:L-1; stem(n,y,*k); grid on ; 实验结果： 1、 2、 3、 4、 实验二 DTFT和Z变换的Matlab实现 1、实验目的和要求 ㈠ 掌握DTFT和Z变换的Matlab实现的方法； ㈡ 进一步了解离散时域系统的时频特性。 2、实验内容 ㈠ DTFT和z变换的计算和基本性质； ㈡ 离散时域系统的仿真； ㈢ 利用函数实现离散时域系统的频率响应，分析零、极点的分布，进一步了解离散时域系统的时域特性。 3、实验原理 利用Matlab软件计算出系统函数的零极点分布、单位脉冲响应和系统频率响应等的图像并于笔算结果进行比较，找出异同。编译合适程序能计算取值范围不同的离散卷积。 4、主要仪器设备 计算机，Matlab软件。 实验步骤： 程序1： n=0:30; x=zeros(1,length(n)); h=zeros(1,length(n)); x([find((n=0)(n=4))])=1; h([find((n=0)(n=8))])=0.5; figure(1) subplot(3,1,1); stem(n,x); axis([0,30,0,2]); title(输入序列); xlabel(n); ylabel(x(n)); subplot(3,1,2); stem(n,h); axis([0,30,0,2]); title(冲激响应序列); xlabel(n); ylabel(h(n)); y=conv(x,h); subplot(3,1,3); n=0:length(y)-1; stem(n,y); title(输出响应); xlabel(n); ylabel(y(n)); 程序2： b=[0.0181,0.0543, 0.0543,0.0181]; a=[1.000,-1.76, 1.1829,-0.2781]; w=pi*freqspace(500); H=freqz(b,a,w); MH=abs(H); AH=angle(H); subplot(2,1,1); plot(w/pi,MH); grid; axis([0,1,0,1]); xlabel(w(pi)); ylabel(|H|); title(幅度、相位响应); subplot(2,1,2); plot(w/pi,AH); grid; xlabel(w(pi)); ylabel(angle(H));  实验结果： 1、 2、 实验三 DFS和DFT的Matlab实现 1、实验目的和要求 ㈠ 加深DFS和DFT的算法原理和基本性质的理解； ㈡ 熟悉利用FFT对典型数字信号进行频谱分析。 2、实验内容 ㈠ 使用Matlab实现程序周期序列的DFS； ㈡ 利用FFT和IFFT编程实现两个有限长序列的循环卷积； 3、实验原理 对有限长序列使用离散Fouier变换(DFT)可以很好的反映序列的频谱特性，而且易于用快速算法在计算机上实现。 4、主要仪器设备 计算机，Matlab软件。  实验步骤： 程序1： function shiyan22() F=50; N=64; T=0.000625; n=1:N; x=cos(2*pi*F*n*T); subplot(2,1,1); plot(n,x); hold on title(x(n)); xlabel(n); X=fft(x); subplot(2,1,2); plot(n,X); grid title(FFT|X