matlab差分阶跃响应,四、离散LTI冲激和阶跃响应实验报告.docx

实验题目离散LTI系统的冲激响应与阶跃响应　实验时间　实验台号　实验目的差分方程；单位冲激响应；单位阶跃响应；　实验仪器PC台式机，软件MATLAB 实验预习实验预习：1、LTI差分方程描述及系统响应；2、差分方程三种解法；　3、系统函数的描述方法；实验原理：1、系统响应： 令： 则：2、LTI系统响应：可以改写成：：系统幅度响应；：系统相位响应；3、差分方程： 4、差分方程解法： 卷积和计算法； 递推法； 经典解法； 5、MATLAB对LTI系统描述：1、传递函数模型(tf) 2、零--极点增益模型(zp)： 3、二阶分时模型(sos)： 4、状态空间模型(ss)： 其中A,B,C,D表示系统状态空间矢量实验步骤应用函数ss2tf();zp2tf();sos2tf()；sos2zp();zp2sos;题目4-1：求离散时间系统的零、极点向量和增益系数？b=[2,3];a=[1,0.4,1];[b,a]=eqtflength(b,a);[z,p,k]=tf2zp(b,a) %结果4-1题目4-2：二、书上例题：求二阶系统对100点噪声的响应曲线？b=[2,-3.4,5.5];a=[1,-1.2,0.8];x=randn(1,100);y=dlsim(b,a,x);plot(y); title('随机噪声响应曲线');xlabel('t/s');ylabel('幅度');b=[2,-3.4,5.5];a=[1,-1.2,0.8];h=impz(b,a);subplot(2,1,1);stem(h);title('系统单位冲激响应');xlabel('n');ylabel('h(n)');x=randn(1,100);y=conv(x,h);subplot(2,1,2);plot(y(1,1:100));title('随机噪声响应曲线');xlabel('t/s');ylabel('幅度'); %图像4-2题目4-3：三、已知离散线性时不变系统的差分方程，请分别用impz和dstep、filtic和filter子函数两种方法求解系统的冲激响应和阶跃响应。①x(n)＋x(n－6)＝y(n) ②2y(n)－3y(n－1)＋y(n－2)＝x(n－1)a=[1 0 0 0 0 0 0];b=[1 0 0 0 0 0 1];N=32;%作图点数n=0:N-1;hn=impz(b,a,N);%时域单位冲激函数gn=dstep(b,a,N);%时域单位阶跃函数subplot(1,2,1),stem(n,hn);title('系统单位冲激函数')%axis([0 N 1.1*min(hn) 1.1*max(hn)])subplot(1,2,2),stem(n,gn);title('系统单位阶跃函数')%axis([0 N 1.1*min(gn) 1.1*max(gn)]) %图像4-3a=[1 0 0 0 0 0 0];b=[1 0 0 0 0 0 1];tmp=a(1);a=a./tmp;%归一化b=b./tmp;N=32;%作图点数n=0:N-1;x1= n==0;%建立输入单位冲激信号x1(n)%xi=filtic(b,a,0);%求等效初始条件的输入序列hn=filter(b,a,x1);%对输入单位冲激信号进行滤波，求单位冲激响应x2= n>=0;%建立输入单位阶跃信号x2(n)gn=filter(b,a,x2);%对输入单位冲激信号进行滤波，求单位冲激响应subplot(1,2,1),stem(n,hn);title('系统单位冲激函数')%axis([0 N 1.1*min(hn) 1.1*max(hn)])subplot(1,2,2),stem(n,gn);title('系统单位阶跃函数')%axis([0 N 1.1*min(gn) 1.1*max(gn)]) %图像4-4a=[2 -3 1];b=[0 1 0];tmp=a(1);a=a./tmp;%归一化b=b./tmp;N=32;%作图点数n=0:N-1;hn=impz(b,a,N);%时域单位冲激函数gn=dstep(b,a,N);%时域单位阶跃函数subplot(1,2,1),stem(n,hn);title('系统单位冲激函数')%axis([0 N 1.1*min(hn) 1.1*max(hn)])subplot(1,2,2),stem(n,gn);title('系统单位阶跃函数')%axis([0 N 1.1*min(gn) 1.1*max(gn)]) %图像4-5a=[2 -3 1];b=[0 1 0];tmp=a(1);a=a./tmp;%归一化b=b./tmp