数字信号处理试验源程序详解.doc

试验程序 试验一 n=[-3:3];xn=[(n-0)==0]; %序列的起点为-3，终点为3，在n=0处生成一个单位脉冲； subplot(1,1,1);stem(n,xn,.); %stem：离散序列绘图指令； Title(单位脉冲序列); line([-3,3],[0,0]) %创建直线； axis([-3,3,0,1]); %控制坐标轴比例和外观； xlabel(n);ylabel(x(n)) %X轴和Y轴标注； n=[-10:10];xn=[(n-0)=0]; %序列的起点为-10，终点为10，在n=0处生成一个单位阶跃； subplot(1,1,1);stem(n,xn,.); Title(单位阶跃序列); line([-10,10],[0,0]) axis([-10,10,0,1]); xlabel(n);ylabel(x(n)) clear,close all %清除所有程序； n=[-10:10];xn1=[(n-0)=0]; xn2=[(n-5)=0]; %定义两个阶跃序列； xn=xn1-xn2; %两个阶跃序列之差得到矩形序列； subplot(1,1,1);stem(n,xn,.); Title(矩形序列); line([-10,10],[0,0]) axis([-10,10,0,1]); xlabel(n);ylabel(x(n)) T=5,t=-2*T:0.01:2*T;duty=50; %duty—信号为正的区域在一个周期内所占的百分比； x=square(t,duty); subplot(1,1,1) plot(t,x) title(连续时间函数-周期方波脉冲) %title—标出图名； axis([-10,10,-1.2,1.2]) line([-10,10],[0,0]) T=4;t=-2.5*T:0.01:2.5*T; x=sawtooth(t-2*T/3,0.5); %sawtooth（x,width）：三角波函数，width为该位置的横坐标和周期的比值（0-1的标量）。width=1时产生锯齿波，width=0.5时产生标准对称三角波； subplot(1,1,1); plot(t,x); title(周期三角波脉冲，T=4) axis([-10,10,-1.2,1.2]); line([-10,10],[0,0]) 实验二 clear;close all; N=64;n=[0:N-1]; A=[1,-0.9];B=[0.05,0.05]; %系统差分方程系数向量； x=[n==0]; %产生deta(n)； x1n=[1 1 1 1 1 1 1 1 zeros(1,56)]; %产生 x2n=ones(1,64); %产生 hn=filter(B,A,x); %filter求系统对序列x的响应； subplot(2,2,1); stem(n,hn,.); line([0,64],[0,0]); axis([0,64,0,0.1]); title((a)系统单位脉冲响应h(n)) y1n=filter(B,A,x1n) %求系统对序列x1n的响应； subplot(2,2,2); stem(n,y1n,.);line([0,64],[0,0]); axis([0,64,0,0.6]); title((b)系统对R8（n）的响应y1(n)) y2n=filter(B,A,x2n); %求系统对序列x2n的响应； subplot(2,2,4); stem(n,y2n,.);line([0,64],[0,0]); axis([0,64,0,1]); title((c)系统对u(n)的响应y2(n)) close all;clear all x1n=[1 1 1 1 1 1 1 1];N=8;M1=10;M2=4; h1n=ones(1,10) ; h2n=[1 2.5 2.5 1]; y1n=conv(x1n,h1n); %conv求卷积； y2n=conv(x1n,h2n);