数字信号处理实验指书钟黎萍.docVIP

实验1 Matlab上机操作实践 一、实验目的 1．Matlab环境，掌握基本指令 2．学习数组运算、赋值、绘图及程序的编写 二、实验涉及的Matlab函数 1．abs ：求模值（绝对值） 2．plot：绘图 3．stem：绘脉冲杆图 4．subplot：建立子图 5．title：标注图名 6．xlabel：横坐标标注 7．ylabel：纵坐标说明 8．sawtooth： 锯齿波产生函数 9．square： 方波产生函数 三、实验内容 1．数组赋值 2．数组运算 3．常用函数及波形显示 显示曲线：f（t）=2sin（2*pi*t） 例程： t=0:0.05:3; f=2*sin(2*pi*t); plot(t,f); title(‘f(t)-t曲线’); xlabel(‘t’); ylabel(‘f(t)’); subplot(2,1,1) stem(t,f) subplot（2,1,2） stairs(t,f) 四、实验报告 产生并显示下列波形：三角波，锯齿波，方波，sin（t^2）,exp(j*t^2) 实验二 序列的时域特性 一、实验目的 通过对序列进行简单的运算处理，生成具有所需特性的输出信号，进而观察序列的线性、时延、卷积等特性。 二、实验原理 对离散系统y(n)=T[x(n)]，若满足： 当x(n)=ax1(n)+bx2(n),时，输出为y(n)= ay1(n)+by2(n), 其中：y1(n)= T[x1(n)], y2(n)= T[x2(n)]。 则此系统称为线性系统。 对离散系统y(n)=T[x(n)]，若满足： 当x(n)=x1(n-n0)时，输出为y(n)= y1(n- n0)， 其中：y1(n)= T[x1(n)]。 则此系统称为时不变系统 既满足线性又满足时不变特性的系统称为线性时不变系统（LTI）。 任何输入序列的输出可以用输入序列与系统单位冲击响应的卷积来表示： y(n)= ∑h(n-k)x(k)。 有界输入产生有界输出的系统称为稳定系统。 MATLAB相关函数 filter: 功能：产生离散系统的输出。 调用格式：y=filter(num,den,x) 其中：num: 系统函数的分子向量，num=[a0, a1, ... aN] den: 系统函数的分母向量，den=[b0, b1, ... bM] x: 输入序列 impz: 功能：求系统的冲激响应 调用格式：y=impz(num,den,N),N为样本数。 conv: 功能：求卷积 调用格式：y=conv(x1,x2) 三、实验内容 1.滑动平均系统 y(n)=[x(n)+x(n-1)+...+x(n-M+1)]/M: 表示求输入序列的平均值，实际上起到低通滤波的作用。 例程： n=0:100； x1=cos(2*pi*0.05*n); x2=cos(2*pi*0.5*n) x=x1+x2; M=length(x1) num=ones(1,M); y=filter(num,1,x)/M; subplot(2,2,1) plot(n,x1) axis([0,100,-2.2]) xlabel(‘时间序号n’)；ylabel（‘振幅’）； title（‘信号x1’） subplot(2,2,2) plot(n,x2) axis([0,100,-2.2]) xlabel(‘时间序号n’)；ylabel（‘振幅’）； title（‘信号x2’） subplot(2,2,3) plot(n,x) axis([0,100,-2.2]) xlabel(‘时间序号n’)；ylabel（‘振幅’）； title（‘信号x’） subplot(2,2,4) plot(n,y) axis([0,100,-2.2]) xlabel(‘时间序号n’)；ylabel（‘振幅’）； title（‘信号y’） 比较输出序列y与输入序列x，看看发生了什么变化。 2.线性特性 系统差分方程为：y(n)-0.4y(n-1)+0.75y(n-2)=2.24x(n)+2.49x(n-1)+2.24x(n-2),验证此系统 的线性特性。 例程： n=0:40; a=3;b=3; x1=cos(2*pi*0.1*n); x2=cos(2*pi*0.2*n); x=a*x1+b*x2; num=[2.24 2.49 2.24]; den=[1 -0.4 0.75]; ic=[0 0];%初始值 y1=filter(num,den,