系统响应及系统稳定性

精品.实验一：系统响应及系统稳定性1、在时域求系统响应的方法有两种：第一种是通过解差分方程求得系统输2、检验系统的稳定性，其方法是在输入端加入单位阶跃序列，观察输出波形，如果波形稳定在一个常数值(包括零)上，系统稳定，否则不稳定。—幅频、相频)、系统函数/特性(系统单位脉冲响应的Z变换)、零极点分布结n精品.附录1.1、实验结果见(1)用实验方法检查系统是否稳定，输入信号为u(n)时，画出系统输出波(2)给定输入信号为精品.1015202530351520253035404520304050500001005响应输出信号向量Y，B、A是差分方程的系数向量。即B=[a1,a2……am]精品.((d)系统单位脉冲响应h1(n)2.50.80.60.40.202468101214(e)h1(n)与R8(n)的卷积y21(n)86654432210510152025(g)h2(n)与R8(n)的卷积y22(n)(f)系统单位脉冲响应h2(n)510152025510150.50000002711(d)(f)单位脉冲响应点数与程序要求一致；(e)(g)卷积点数满足M+N-1的要求，图形也满足要求。精品.hun)的响应y31(n)0.050-0.05050100150200250300i信号的响应y32(n)10.50-0.5-1050100150200250300在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数(包括零)，就可以断定系统是稳定的，(h)中的输出显然趋近于零，所以系统是稳在(i)中，谐振器具有对某个频率进行谐振的性质，本实验中的谐振器的的单位脉冲响应是有限长序列，可否精品.冲响应h(n)与各段的卷积；精品.号的形，如果波形稳定在一个常数值(包括零)上，系统稳定，否则不稳定。在此次试验中，通过课堂所留的三个例子，温习了关于MATLAB软件的操MATLAB软件的了解，体会到了MATLAB具有完备的图形处理功能，实现计算精品.精品.本里面知识点的应用，更深一层的加深了对MATLAB软件的使用。这对自己以六、附件：MATLAB原程序清单精品stept=2*pi/n;w=stept:stept:2*pi;y=sin(2.5*w)./sin(0.5*w);axis([stept2*pi-26]);ylabel('y=sin(2.5*pi)/sin(0.5*pi)');xlabel('w=0~2*pi');gridon;课堂作业stept=2*pi/n;w=stept:stept:2*pi;y=sin(2.5*w)./sin(0.5*w);izewaxis([stept2*pi06]);ylabel('y=sin(2.5*pi)/sin(0.5*pi)');xlabel('w=0~2*pi');gridon;.6665544pi3ns*5ns=y01-134w=0~2*pi34w=0~2*pi5656s=yni(sni()(i/302222(i))L=N+M-1;x=[1,2,3,4,5];h=[6,2,3,6,4,2];y=conv(x,h);nx0:N-1;nh=0:M-1;ny:L-1;subplot(231);stem(nx,x,'.k');xlabel('n');ylabel('x(n)');gridon;subplot(232);stem(nh,h,'.k');xlabel('n');ylabel('h(n)');gridon;subplotstemnyy,'.k');xlabel('n');ylabel('y(n)');gridon;精品5568046043222010n3、分析频谱(a，b，c保证幅频特性的最大值为1)。5n040n000几4pp=0.8;r=0.85;alpha=pi/4;N=25;b1=[1,1];a1=[1-p];a=(1-p)/2;b1=b1*a;b2=[1,-1];a2=[1-p];b=(1+p)/2;b2=b2*b;b3=[10-1];a3=[1-2*r*cos(alpha)r*r];c1=exp(j*2*alpha);c=abs(1-c1)/((1-r)*abs(1-r*c1));b3=b3/c;h1=impz(b1,a1,N);subplot(331);stem(h1,'.k');subplot(332);zplane(b1,a1)[H1,P]=freqz(b1,a1,256,'whole',1);subplot(333);plot(P,abs(H1));gridon;h2=impz(b2,a2,N);subplot(334);stem(h2,'.k')holdon;plot(zeros(size(h2)));subplot(335);zplane(b2,a2)[H2,P]=freqz(b2,a2,256,'whole',1);subplot(336);plot(P,abs(H2));gridon;h3=impz(b3,a3,N);subplot(337);stem(h3,'.k')holdon;plot(zeros(size(h3)));subplot(338);zplane(b3,a3)[H3,P]=freqz(b3,a3,256,'whole',1);subplot(339);plot(P,abs(H3精))品;gridon;精品yyy.yyy.215.105001.50.2.10-0.1555101510151015PnamaIPn