信号实验报告

1、信号与系统实验报告课程名称：信号与系统实验实验时间：10 至 11 月任课教师：王群实验地点：4-427实验教师：何冰松, 范哲意实验类型： 原理验证 综合设计 创新学生尹迪学号/班级：1120131225/05111361组号：学院：信息与电子学院同组搭档：专业：信息工程成绩：信号与系统实验报告实验名称：信号的时域描述与运算课程名称：信号与系统实验实验时间：10 月 15 日任课教师：王群实验地点：4-427实验教师：何冰松, 范哲意实验类型： 原理验证 综合设计 创新学生尹迪学号/班级：1120131225/05111361组号：学院：信息与电子学院同组搭档：专业：信息工程成绩：一1.2.

2、3.二1.2.实验目的掌握信号的表达及其可视化方法掌握信号基本时域运算的实现方法。利用实验原理利用利用分析常用信号，加深对信号时域特性的理解表达连续时间信号进行连续时间信号的时域运算（相加，相乘，微分， 尺度变换）表达离散时间信号进行离散时间信号的时域运算,移位利用利用反转3.4.三1.相应程序的代码及实验结果1.1t=0:0.01:10;x=(1-exp(-0.5*t).*heaviside(t); plot(t,x)1.2t=0:0.01:10;x=cos(pi*t).*(heaviside(t)-heaviside(t-2); plot(t,x)1.3t=-2:0.01:10;x=(ab

3、s(t)/2).*cos(pi*t).*(heaviside(t+2)-heaviside(t-2); plot(t,x)1.4t=0:0.01:10;x=exp(-t).*sin(2*pi*t).*(heaviside(t)-heaviside(t-3);plot(t,x)2.2.1n=0:10;x=heaviside(n-3); stem(n,x,'filled');xlabel('n');title('x');2.2n=-10:10;x=(-0.5).(n).*heaviside(n); stem(n,x,'filled')

4、;xlabel('n');title('x');2.3n=-10:10;x=n.*(heaviside(n)-heaviside(n-5); stem(n,x,'filled');xlabel('n');title('x');2.4n=-10:10;x=sin(n.*pi/2).*heaviside(n); stem(n,x,'filled');xlabel('n');title('x');3.t=0:0.01:6;x=square(3*t); plot(t,x)4.

5、4.1t=0:0.01:4;x1=-t+4;x2=sin(2*pi*t); x3=x1+x2; plot(x3)4.2t=0:0.01:4;x1=-t+4;x2=sin(2*pi*t); x3=(x1).*(x2);plot(x3)4.3syms t;x1=(-t+4)*(heaviside(t)-heaviside(t-4); x2=subs(x1,t,-t);x3=(x1)+(x2);ezplot(x3,-4,4)4.4syms t;x1=(-t+4)*(heaviside(t)-heaviside(t-4); x2=sin(2*pi*t);x3=(x1)+(x2);x31=subs(x3

6、,t,t-1); x6=x2*x31; ezplot(x6,-4,10)5.5.1n=-3:4;x=0 1 2 3 3 3 3 0;stem(n,x,'filled');xlabel('n');title('x');5.2n=-3:4;x=0 1 2 3 3 3 3 0;stem(-n,x,'filled'); xlabel('n');title('x');5. 3n=-3:4;x=0 1 2 3 3 3 3 0;stem(n-2,x,'filled'); xlabel('n

7、');title('x');5.4n=-3:4;x=0 1 2 3 3 3 3 0;stem(n+2,x,'filled');xlabel('n');title('x');6.6.1t=0:0.1:20;x=1+cos(pi/4*t-pi/3)+2*cos(pi/2*t-pi/4)+cos(2*pi*t); plot(t,x)周 期 函 数6.2t=0:0.01:20;x=sin(t)+2*sin(pi*t); plot(t,x)非 周 期 函 数6.3t=0:10;x=2+3*sin(2*t*pi/3-pi/8); st

8、em(t,x,'filled');xlabel('t');title('x');周 期 函 数6.4t=0:10;x=cos(t*pi/6)+sin(t*pi/3)+cos(t*pi/2); stem(t,x,'filled');xlabel('t');title('x');周 期 函 数四试验中所遇到的问题及解决方法实验中首次运用进行实际编程，遇到了不熟练，操作不明确等问1.题，经过的努力摸索和查阅书籍，逐渐克服了此类问题。在计算过程中，发现 U（n）函数在零点的取值为 0.5，与实际情况不符，2

9、.通过与交流得知，在中，零点的值等于前值和后值之和的一半，虽然与书本不同，但是不影响实验结果。在进行反转，移位等操作时，出现定义域不明确的问题，通过仔细阅读书3.籍改正了心得体会的错误。五1.2.是一款非常强大的软件的时候，由于其所保有的函数库非常大，基本涵盖了我们在使用所有的应用，所以要勤于查阅，多使用 HELP 等帮助性功能，有助于我们学习。信号与系统实验报告实验名称：LTI 系统的时域分析课程名称：信号与系统实验实验时间：10 月 22 日任课教师：王群实验地点：4-427实验教师：何冰松, 范哲意实验类型： 原理验证 综合设计 创新学生尹迪学号/班级：1120131225/051113

10、61组号：学院：信息与电子学院同组搭档：专业：信息工程成绩：一1.2.3.4.实验目的掌握利用对系统进行时域分析的方法掌握连续时间系统零状态响应，冲击响应和阶跃响应的求解方法掌握求解离散时间系统响应，抽样响应的方法加深对卷积 和计算的方法实验原理和卷积和的理解。掌握利用计算机进行卷积和卷积二1.2.3.三利用利用实现连续时间系统的时域分析实现离散时间系统的时域分析卷积和与卷积相应程序的代码及实验结果第一题：1.1 b=1;a=1 sqrt(2) 1;sys=tf(b,a);t=0:0.01:10;impulse(sys)Impulse Response0.60.50.40.30.20.10-0

11、.1012345678910Time (seconds)1.2b=1 0 0;a=1 sqrt(2) 1;sys=tf(b,a);t=0:0.01:10;impulse(sys)AmplitudeImpulse Response0.20-0.2-0.4-0.6-0.8-1-1.2-1.4-1.60123456789Time (seconds)1.3b=1 0;a=1 1 1;sys=tf(b,a);t=0:0.01:10;impulse(sys)Impulse Response10.80.60.40.20-0.2-0.40246Time (seconds)810121.4a=1 1 1;b=2

12、 0 1;sys=tf(b,a);t=0:0.01:10;impulse(sys)AmplitudeAmplitudeImpulse Response0.50-0.5-1-1.5-20246Time (seconds)81012第二题:2.1（阶跃响应） a=1 1 6;b=1;sys=tf(b,a); step(sys)Step Response0.350.30.250.20.150.10.05002468101214Time (seconds)(冲激响应) a=1 1 6;b=1;sys=tf(b,a); impulse(sys)AmplitudeAmplitudeImpulse Resp

13、onse0.50.40.30.20.10-0.1-0.202468101214Time (seconds)2.2冲激响应的波形在在负无穷到正无穷的求和不是无穷， 所以系统是稳定的。2.3a=1 1 6;b=1;sys=tf(b,a);t=0:0.01:10;x=exp(-t);lsim(sys,x,t)Linear Simulation Results1.210.80.60.40.20-0.2012345678910Time (seconds)第三题3.1a=1 3 2;b=1;impz(b,a,0:10)AmplitudeAmplitudeImpulse Response2500200015

14、0010005000-500-1000-1500012345n (samples)678910不稳定系统3.2a=1 -0.5 0.8;b=1 -3;impz(b,a,0:10)Impulse Response2 521 510 50-0 5-1-1 5-2-2 5012345n (samples)678910稳定系统第四题抽样响应：a=1 1 0.25;b=1;impz(b,a,0:10)AmplitudeAmplitudeImpulse Response10 80 60.40 20-0 2-0.4-0 6-0 8-1012345n (samples)678910阶跃响应：n=0:10;a=

15、1 1 0.25;b=1;x=heaviside(n); y=filter(b,a,x);stem(n,y,'filled');xlabel('n');title('y');y0 50.450.40 350 30 250 20.150.10 050012345n678910第五题n=-3:4;x1=1,2,1,1; x2=1 1 1 1 1;x=conv(x1,x2); stem(n,x,'filled');xlabel('n');title('x')Amplitude0 50-3第六题n=0:40

16、;n1=0:80;h=sin(0.5*n);x=sin(0.2*n);y=conv(x,h); stem(n1,y,'filled');xlabel('n1');title('y');0-1-2-3-40第七题t1=-1:0.001:1;t2=-2:0.001:2;t3=-3:0.001:3; x1=2+t1-t1; x2=1+t2-t2;y=conv(x1,x2)./2000; plot(t3,y);2.521.510.50-3-2-10123四1.试验中所遇到的问题及解决方法实验中运用进行实际编程，遇到了不熟练，操作不明确等问题，经过的努力摸

17、索和查阅书籍，逐渐克服了此类问题。2.在进行卷积运算的时候，会出现卷积结果定义域不明确的问题，需要自己注意。心得体会是一款非常强大的软件五1.2.在使用的时候，由于其所保有的函数库非常大，基本涵盖了我们所有的应用，所以要勤于查阅，多使用 HELP 等帮助性功能，有助于我们学习。3.在卷积部分，对于卷积结果定义域的确定是非常重要的。信号与系统实验报告实验名称：信号的频域分析课程名称：信号与系统实验实验时间：10 月 23 日任课教师：王群实验地点：4-427实验教师：何冰松, 范哲意实验类型： 原理验证 综合设计 创新学生尹迪学号/班级：1120131225/05111361组号：学院：信息与电

18、子学院同组搭档：专业：信息工程成绩：一实验目的1. 深入理解信号频谱的概念，掌握信号的频域分析方法。2. 观察典型周期信号和非周期信号的频谱，掌握其频谱特性。实验原理与方法二1.2.3.4.利用利用利用利用分析连续周期信号的频谱 分析连续非周期信号的频谱分析离散周期时间信号的频域情况 分析离散非周期时间信号的频域情况三相应程序代码及实验结果1.1.2t=-10:0.01:10;N=input('N='); x=zeros(size(t); x=0.5;for n=1:1:N;x=x+(2/(pi*n).*sin(n*pi/2).*cos(0.5*pi*n*t);end plot

19、(t,x);当T=4;A=1;宽度为2时当T=6;A=1;宽度为4，N=10的时候1.210.80.60.40.20-0.2-10-8-6-4-20246810当 T=6;A=1;宽度为 4，N=100 的时候1.210.80.60.40.20-0.2-10-8-6-4-20246810当 T=6;A=1;宽度为 4，N=1000 的时候1.210.80.60.40.20-0.2-10-8-6-4-20246810当 T=6;A=1;宽度为 4，N=1000 的时候1.210.80.60.40.20-0.2-10-8-6-4-202468101.3N = 10;T = input('T

20、=');tt = input('tt=');A = 1;w = 2*pi/T; n1 = -N:-1;c1 = (2*A*l/T)*sin(n1*w*l/2)/w/l./n1*2; c0 = E*l/T;n2 = 1:N;c2 = (2*A*l/T)*sin(n1*w*l/2)/w/l./n1*2; cn = c1 c0 c2;n = -N:N;subplot(211); stem(n,abs(cn),'filled'); xlabel('omega/omega_0'); title('Magnitude of ck');

21、 subplot(212); stem(n,angle(cn),'filled'); xlabel('omega/omega_0'); title('Phase of ck');N = 10，A = 1，T=1tt = 0.5Magnitude of ck0.80.60.40.20-10-8-6-4-20w/w0Phase of ck24681043210-10-8-6-4-20w/w0246810tt = 0.3Magnitude of ck0.80.60.40.20-10-8-6-4-20w/w0Phase of ck24681043210-

22、10-8-6-4-20w/w0246810tt = 0.7Magnitude of ck0.80.60.40.20-10-8-6-4-20w/w0Phase of ck24681043210-10-8-6-4-20w/w0246810思考题：1.什么是吉伯斯现象？产生吉伯斯现象的是什么？答：将具有不连续点的周期函数（如矩形脉冲）进行傅立叶级数展开后，选取有限项进行。当选取的项数越多，在所的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时，该峰起值趋于一个常数，大约等于总跳变值的 9%。这种现象称为吉伯斯现象。是在不连续点附近所有的正弦信号均具有相同的变化趋势，该趋势在有限项内无法被消

23、除。2. 以周期矩形脉冲信号为例，说明周期信号的频谱有什么特点？ 答：周期信号的频谱是具有周期性的一系列的脉冲信号。3. 周期矩形脉冲信号的有效频带宽度与信号的时域宽度之间有什么关系？答：时域宽度越大，有效频带宽度越小。4. 随着矩形脉冲信号参数的变化，其频谱结构如何变化？答：频谱包络形状不变，过零点不变，普贤间隔随着 T 变大而缩小。2.2.2syms t A tt; tt=input('tt=');A=input('A=');x=A*( heaviside(t+ (tt/2) )-heaviside(t- (tt/2) ); X=fourier(x)ezpl

24、ot(abs(X),-20,20);A=1;tt=1;A=1;tt=2;A=1;tt=5;思考题：1).比较矩形脉冲信号和周期矩形脉冲信号的频谱，两者之间有何异同？答：矩形脉冲信号频谱为连续函数，周期矩形脉冲信号的频谱为一系列的脉冲。2).根据矩形脉冲宽度变化时频谱的变化规律，说明信号的有效频带宽度与其时域宽度之间有什么关系？答：信号有效频带宽度越大，时域宽度越小。3.3.3Double N; double N1; N=input('N=');N1=input('N1=');n=-N+1:N-1;x(n>=-N1 & n<= N1)=1; x

25、(n>N1 & n<(N-N1)=0; x(n>-N+N1 & n<-N1)=0; x(n>=N-N1 & n<=N+N1)=1;subplot(211); stem(n,x,'filled');xlabel('n');title('x(n)'); subplot(212); X=fft(x); stem(n,X,'filled');xlabel('k');title('X(k)');N1=5,N=20 的运行结果x(n)10.50-20-

26、15-10-50n X(k)510152020100-10-20-15-10-50k5101520x(n)10.50-20-15-10-50n X(k)5101520151050-5-20-15-10-50k5101520思考题：1).以周期方波序列为例，说明周期序列与连续周期信号的频谱有何异同。答：周期序列的频谱向外越来越大，连续周期信号频谱则是中间向两边越来越小。2).随着周期方波序列占空比的变化，其频谱如何随之变化？ 答：随着占空比越来越大，频谱密度也越来越大。4.4.1double N1; N1=input('N1='); w=-pi:0.01*pi:pi; n=-N1

27、:N1;x(n>=(-N1) & n<=N1)=1;X=x*exp(-j*(n'*w) ); subplot(211); stem(n,x,'filled');xlabel('n');title('x(n)'); subplot(212); plot(w/pi,abs(X);N=10x(n)10.50-10-8-6-4-20n2468103020100-1-0.8-0.6-0.4-0.200 20.40.60.81N=5x(n)10.50-5-4-3-2-10n12345151050-1-0.8-0.6-0.4-0.2

28、00 20.40.60.81思考题：1).随着矩形脉冲序列宽度的变化，其频谱如何随之变化？其宽度与频谱的有效宽度有何关系？答：宽度越大，频谱变化就越密集。矩形脉冲信号序列宽度越大，频谱有效宽度越大。试验中所遇到的问题及解决方法1. 此次实验难度与前两次相比较有明显的提升，花费的时间也明显增加。很多程序都是在与同学讨论之后才写出来的。本次实验中我遇到的问题就是容易误把频谱也做了累加，得到了错误的波形和结果， 在经过和上的范例对比之后，我找出了错误，得到了正确的波形和结果。心得体会四五1.主要收获就是从实验的角度了解了吉伯斯现象，同时也掌握了对连续时间信号以及离散时间信号的频域分析方法，很使用多时

29、域上难以分析的信号在频域上是非常直观和间接的。2. 信号的频域分析方法对于信号与系统这门课程的学习乃至是今后要从事的有关信号处理的工作而言都是关键、有用的。利用将十分抽象的表达式转化为了直观的图像，并且从侧面印证了很多课上学习到的东西，加深了对课程的理解。或许这个实验的意义就是印证课上知识加深理解，并掌握分析方法，为以后工作与学习打下坚实的基础。信号与系统实验报告实验名称：LTI 系统的频域分析课程名称：信号与系统实验实验时间：10 月 28 日任课教师：王群实验地点：4-427实验教师：何冰松, 范哲意实验类型： 原理验证 综合设计 创新学生尹迪学号/班级：1120131225/051113

30、61组号：学院：信息与电子学院同组搭档：专业：信息工程成绩：一、1.2.二、1.2.三、1.实验目的加深对 LTI 系统频率响应概念的掌握和理解学习和掌握 LTI 系统频率特性的分析方法实验原理与方法利用利用计算连续时间系统的频率响应计算离散时间系统的频率响应相应程序的代码及实验结果1.2b=0.04 0 0;a=0.04 0.4 2;freqs(b,a); H,w=freqs(b,a); subplot(211); plot(w,abs(H);axis(0,30,0,1);set(gca,'xtick',0:2:30);set(gca,'ytick',0 0.

31、4 0.707 1); xlabel('omega(rad/s)'); ylabel('Magnitude'); title('|H(jomega)|');grid on; subplot(212); plot(w,angle(H);xlabel('omega(rad/s)'); ylabel('Phase');title('phi(omega)'); grid on;|H(jw)|10.7070.400 2468 10 1214 16 18 20w(rad/s) phi(w)22 24 26 28

32、 30432100102030405060w(rad/s)7080901002.2.1R=input('R=');MagnitudePhaseC=input('C=');b=1;a=R*C 1;H,w=freqs(b,a); subplot(211); plot(w,abs(H); xlabel('omega(rad/s)'); ylabel('Magnitude'); title('|H(jomega)|'); grid onsubplot(212); plot(w,angle(H); xlabel('o

33、mega(rad/s)'); ylabel('Phase');title('phi(omega)'); grid on；RC=1 时|H(jw)|10.500123456w(rad/s) phi(w)789100-0.5-1-1.50123456w(rad/s)78910RC=0.1 时|H(jw)|10.500102030405060w(rad/s) phi(w)7080901000-0.5-1-1.50102030405060w(rad/s)708090100PhasePhaseMagnitudeMagnitudeRC=0.01 时|H(jw)|10

34、.500100200300400500600w(rad/s) phi(w)70080090010000-0.5-1-1.50100200300400500600w(rad/s)7008009001000RC=0.001 时|H(jw)|10.500100020003000 4000 5000 6000w(rad/s) phi(w)700080009000100000-0.5-1-1.50100020003000 4000 5000 6000w(rad/s)700080009000100002.2运用 RC=0.005 时的滤波系统t=0:0.001:0.2; b=1;a=0.001 1;sys

35、=tf(b,a); x=cos(100*t)+cos(3000*t); lsim(sys,x,t);PhaseMagnitudePhaseMagnitudeLinear Simulation Results21.510.50-0.5-1-1.5-200.020.040.060.080.10.120.140.160.180.2Time (seconds)3.3.1xn+ xn-1+ xn-2+ xn-3+ xn-4+ xn-5+ xn-6+xn-7+ xn-8=yn系统函数：H(z)=1+z-1+z-2+z-3+z-4+z-5+z-6+z-7+z-83.2b=1 1 1 1 1 1 1 1 1;

36、 a=1;impz(b,a);3.3.1b=1 1 1 1 1 1 1 1 1; a=1;H,w=freqz(b,a); plot(w/pi,abs(H)幅频特性曲线Amplitudeb=1 1 1 1 1 1 1 1 1;a=1;H,w=freqz(b,a); plot(w/pi,angle(H);相频特性b=1 1 1 1 1 1 1 1 1;a=1;H,w=freqz(b,a); zplane(b,a);零级点图4.n=-10:10;x=cos(0.3*pi*n)+0.5*cos(0.8*pi*n); stem(n,x,'filled');xlabel('n

37、9;);title('x(n)');x(n)图像n=-10:10;y=2*cos(0.3*pi*n); stem(n,y,'filled');xlabel('n');title('y(n)');stem(n,y,'filled');y(n)图像四、1.2.五、1.2.实验中所遇到的问题及解决方法在实验 3 中，错误的以为输入只有 8 项，导致图像相似但是确。在实验 4 中，对题目理解有一些问题，和心得体会交流后解决问题。第一次发现原来还可以画零极点图，对其有了更深的认识。在使用的时候，由于其所保有的函数库非常大，基

38、本涵盖了我们所有的应用，所以要勤于查阅，多使用 HELP 等帮助性功能，有助于我们学习。信号与系统实验报告实验名称：连续时间系统的复频域分析课程名称：信号与系统实验实验时间：10 月 30 日任课教师：王群实验地点：4-427实验教师：何冰松, 范哲意实验类型： 原理验证 综合设计 创新学生尹迪学号/班级：1120131225/05111361组号：学院：信息与电子学院同组搭档：专业：信息工程成绩：一、实验目的1.掌握拉斯变换及其反变换的定义，并掌握实现方法2. 学习和掌握连续时间系统系统函数的定义及复频域分析方法3. 掌握系统零极点的定义，加深理解系统零极点分布于系统特性的关系二、实验原理及

39、方法1. 利用2. 利用3. 利用实现拉斯变换实现连续时间系统的系统函数表示实现连续时间系统的零极点分析三、相应程序代码及实验结果1.h=sym('heaviside(t)-heaviside(t-2)'); h1=laplace(h); x=sym('heaviside(t)');h2=laplace(x); Y=h2.*h1; y=ilaplace(Y)ezplot(y,0,5);2.1.b=1;a=1 2 2 1;sys=tf(b,a); pzmap(sys); grid on;b=1;a=1 2 2 1;H,w=freqs(b,a); subplot(2

40、11); plot(w,abs(H);xlabel('omega(pi)');ylabel('幅频特性'); grid on; subplot(212); plot(w,angle(H);xlabel('omega(pi)');ylabel('相频特性'); grid on;冲激响应：b=1;a=1 2 2 1;r,p,k=residue(b,a) r =1.0000 + 0.0000i-0.5000 - 0.2887i-0.5000 + 0.2887ip =-1.0000 + 0.0000i-0.5000 + 0.8660i-0

41、.5000 - 0.8660ik =X(t)=(exp(-1)+(-0.5-0.2887i)*(exp(-0.5+0.866i)+(-0.5+0.2887i)*(exp(-0.5 - 0.8660i);2.b=1 0 1;a=1 2 -3 3 3 2;sys=tf(b,a); pzmap(sys); grid on;b=1 0 1;a=1 2 -3 3 3 2;H,w=freqs(b,a); subplot(211); plot(w,abs(H);xlabel('omega(pi)');ylabel('幅频特性'); grid on ; subplot(212)

42、; plot(w,angle(H);xlabel('omega(pi)');ylabel('相频特性'); grid on ;冲激响应b=1 0 1;a=1 2 -3 3 3 2;r,p,k=residue(b,a)r =0.0769 + 0.0000i-0.0300 - 0.0881i-0.0300 + 0.0881i-0.0085 - 0.1436i-0.0085 + 0.1436ip =-3.1704 + 0.0000i0.9669 + 0.9540i0.9669 - 0.9540i-0.3817 + 0.4430i-0.3817 - 0.4430ik =

43、X(t)=(0.0769)*exp(-3.1704)+(-0.03-0.0881i)*(exp(0.9669+0.9540i)+(-0.03+0.0881i)*(exp(0.9669-0.9540i)+(-0.0085 - 0.1436i)*(exp(-0.3817+0.4430i)+(-0.0085+0.1436i)*(exp(-0.3817-0.4430i);3.1.t=0:0.01:10;b=1;a=1 0;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=heaviside(t);plot(t,x);2.t=0:0.01:10;

44、b=1;a=1 2;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=exp(-2*t).*heaviside(t);plot(t,x);3.t=0:0.01:10;b=1;a=1 -2;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=exp(2*t).*heaviside(t);plot(t,x);4.t=0:0.01:10;b=1;a=1 0 4;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=(sin(2*t)/

45、2).*heaviside(t);plot(t,x);5.t=0:0.01:10;b=1;a=1 2 17;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=0.25*exp(-t).*sin(4*t).*heaviside(t);plot(t,x);6.t=0:0.01:10;b=1;a=1 -2 17;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=0.25*exp(t).*sin(4*t).*heaviside(t);plot(t,x);4.1.t=0:0.01:10;b

46、=1;a=1 2 17;sys=tf(b,a); subplot(211); pzmap(sys);subplot(212);x=0.25*exp(-t).*sin(4*t).*heaviside(t);plot(t,x);2.t=0:0.01:10;b=1 8;a=1 2 17;sys=tf(b,a); subplot(211); pzmap(sys); grid on ; subplot(212);h2=exp(-t).*(cos(4*t)+(7*sin(4*t)/4).*heaviside(t); plot(t,h2);grid on;3.t=0:0.01:10;b=1 -8;a=1 2 17;sys=tf(b,a); subplot(211);pzmap(sys);subplot(212);h3=(exp(-t).*(cos(4*t