数字信号处理MATLAB仿真.doc

数字信号处理MATLAB仿真————————————————————————————————作者：————————————————————————————————日期：个人收集整理勿做商业用途个人收集整理勿做商业用途个人收集整理勿做商业用途实验一数字信号处理的Matlab仿真一、实验目的1、掌握连续信号及其MATLAB实现方法；2、掌握离散信号及其MATLAB实现方法3、掌握离散信号的基本运算方法，以及MATLAB实现4、了解离散傅里叶变换的MATLAB实现5、了解IIR数字滤波器设计6、了解FIR数字滤波器设计1二、实验设备计算机，Matlab软件三、实验内容(一）、连续信号及其MATLAB实现单位冲击信号例1.1:t=1/A=50时，单位脉冲序列的MATLAB实现程序如下：clearall;t1=—0.5:0。001：0；A=50；A1=1/A；n1=length(t1）；u1=zeros(1，n1);t2=0：0。001：A1；t0=0；u2=A*stepfun（t2，t0);t3=A1：0。001：1；n3=length（t3)；u3=zeros(1,n3)；t=[t1t2t3］;u=[u1u2u3］；plot(t，u)axis(［—0。510A+2］）任意函数例1.2：用MATLAB画出如下表达式的脉冲序列clearall;t=—2：1:3；N=length（t)；x=zeros（1，N);x（1）=0。4；x(2)=0。8x（3）=1。2;x(4）=1.5；x(5）=1.0；x(6)=0.7；stem（t,x);axis(［-2。23。201.7]）单位阶跃函数例1。3：用MATLAB实现单位阶跃函数clearall;t=—0。5：0.001：1;t0=0；u=stepfun（t,t0）；plot(t，u）axis(［-0.51—0.21。2］）斜坡函数例1.4:用MATLAB实现g(t）=3（t—1)clearall;t=0:0。01:3;B=3；t0=1；u=stepfun(t，t0）;n=length（t）;fori=1:nu（i）=B*u（i)＊(t(i）-t0）;endplot（t,u)axis（[—0。23。1—0。26.2］）实指数函数例1.5:用MATLAB实现clearall;t=0:0.001:3；A=3;a=0.5；u=A*exp（a＊t);plot(t，u)axis(［-0。23.1-0。214］)正弦函数例1。6:用MATLAB实现正弦函数f（t)=3cos（10πt+1）clearall；t=-0。5：0.001：1；A=3；f=5;fai=1；u=A＊sin（2*pi*f＊t+fai）;plot（t,u）axis(［—0.51-3。23.2]）（二）、离散信号及其MATLAB实现单位冲激序列例2.1：用MATLAB产生64点的单位冲激序列clearall;N=64;x=zeros（1,N)；x(1)=1;xn=0：N—1；stem(xn,x)axis([—16501。1])任意序列例2.2：用MATLAB画出如下表达式的脉冲序列clearall;N=8；x=zeros(1,N）；x（1）=8。0;x（2）=3.4x（3）=1。8;x(4)=5.6;x(5)=2。9;x（6）=0.7;xn=0：N—1;stem（xn，x)axis(［—1808.2])单位阶跃序列例2.3:用MATLAB实现单位阶跃函数clearall；N=32；x=ones(1，N);xn=0:N-1；stem（xn，x)axis（[—13201.1］)斜坡序列例2。4：用MATLAB实现g（n)=3（n-4)点数为32的斜坡序列clearall;N=32；k=4B=3；t0=1；x=［zeros(1,k)ones(1,N-k)]；fori=1：Nx（i)=B*x（i）＊（i-k）；endxn=0：N-1;stem（xn，x)axis([—132090]）正弦序列例2。5:用MATLAB实现幅度A=3，频率f=100，初始相位Φ=1.2,点数为32的正弦信号clearall；N=32；A=3;f=100；fai=1.2;xn=0：N-1;x=A＊sin（2＊pi＊f*(xn/N)+fai）；stem(xn,x）axis([-132—3.23.2]）实指数序列例2。6：用MATLAB实现,点数为32的实指数序列clearall；N=32;A=3；a=0。7；xn=0:N—1；x=A＊a.^xn；stem（xn,x）复指数序列例2.7:用MATLAB实现幅度A=3,a=0。7，角频率ω=314，点数为32的实指数序列clearall；N=32；A=3；a=0。7;w=314；xn=0：N—1