华北电力大学数字信号处理实验报告.docVIP

华 北 电 力 大 学 实 验 报 告 实验环境 MATLAB 6.5 实验名称 实 验 目 的 MATLAB仿真软件的使用在掌握数字信号处理相关理论的基础上根据数字信号处理课程所学内容利用MATLAB对常系数差分方程求解分析不同初始条件对解的影响 2.由MATLAB软件实现常用序列的编程并会编写函数且调用MATLAB内置函数求解差分方程掌握差分方程的求解方法 实 验 原 理 一阶前项差分与后项差分方程并无本质差别。 通常所说的差分方程形式为： y(k)及其各移位序列均为常数，就称其为常系数差分方程；如果某些系数是变量k的函数，就称其为变系数差分方程。描述LTI离散系统的是常系数线性差分方程。 已知系统的输入序列，通过求解差分方程可以求出输出序列。求解差分方程的基本方法有以下三种： （1）经典解法。这种方法类似于模拟系统中求解微分方程的方法，它包括齐次解与特解，由边界条件求待定系数。 （2）递推解法。这种方法简单，且适合用计算机求解，但只能得到数值解，对于阶次较高的线性常系数差分方程不容易得到封闭式（公式）解答。 （3）变换域方法。这种方法是将差分方程变换到z域进行求解，方法简便有效。 当然还可以不直接求解差分方程，而是先由差分方程求出系统的单位脉冲响应，再与已知的输入序列进行卷积运算，得到系统输出。但是系统的单位脉冲响应如果不是预先知道，仍然需要求解差分方程，求其零状态响应解。 （4）卷h(n)，再与已知的x(n) 进行卷积，得到 y(n)。 实 验 内 容 一实验内容2、已知一个二阶线性常系数差分方程用下式表示： y(n)+a1y(n-1)+a2y(n-2)= b0x(n)+b1x(n-1)+b2x(n-2)， 要求： （1）参数a1、a2、b0、b1、b2由运行时输入； （2）出； （3）已知输入，求出并分别画出时域波形图； ①实现三种典型序列 打开MATLAB软件编写程序设计三种典型序列分别是单位取样序列单位冲击响应序列和矩形序列编写完毕后保存以待后续使用根据三种序列的特编写程序 ②解常系数差分方程 首先写出传参数的内置函数=fc(a,b),编写程序求差分方程的MATLAB内置滤波器函数 分别用线性卷积法和利用FFT的方法求输入为时系统的响应 在MATLAB信号处理工具中提供了利用n),对其求k),再求出n)的X（k），两者相乘，再求FFT反变换即可得系统的响应，再画出其时域波形图。 三、实验程序 ①、单位冲激序列 function[x,n]=impseq(n0,n1,n2) n=[n1:n2];x=[(n-n0)==0]; ⑵、单位阶跃序列 function[x,n]=stepseq(n0,n1,n2) n=[n1:n2];x=[(n-n0)=0]; ③、矩形系列 function[x,n]=rectseq(n0,ns,nf,N) n=[ns:nf]; x=[(n-n0)=0((n0+N-1)-n=0]; ④、线性卷积函数function[y,ny]=conv_m(x,nx,h,nh) nyb=nx(1)+nh(1);nye=nx(length(x))+nh(length(h)); ny=nyb:nye; y=conv(x,h); ⑤、求系统响应 function[]=fc(a,b) %求冲激响应 x=impseq(0,0,20);n=[0:20]; 实 验 h=filter(b,a,x); subplot(2,2,1);stem(n,h);title(冲激响应); xlabel(n);ylabel(h(n)); %线性卷积 n=0:20;x=(0.5).^n;nx=[0:20];nh=[0:20]; y=conv_m(x,nx,h,nh);ny=length(y); n=0:ny-1; subplot(2,2,2);stem(n,y);title(利用线性卷积求); xlabel(n);ylabel(y1(n)); %利用FFT求 L=length(x)+lemhth(h)-1; X=fft(x,L); H=fft(H,L); y2=ifft(X.*H) subplot(2,2,3);stem(n,y2);title(利用FFT求); xlabel(n);ylabel(y2(n)); hold on; 实验 结 果 及 分 析 ①、单位冲激序列 function[x,n]=impseq(n0,n1,n2) n=[n1:n2];x=[(n-n0)==0]; ②、单位阶跃序列 function[x,n]=stepseq(n0,n1,n2) n=[n1:n2];x=[(n-n0)=0]; ③、矩形系列 functio