MATLAB 实验报告1

1、matlab 实验报告 学院 _机电学院_ 姓名_ _ 班号_ 学号_信号与系统实验（一）1. 试绘出下列各信号的波形图：（1）实验程序：n=-3:15;m=n-1;a=2;x1=a.m;stem(n,x1,fill),grid on,xlabel(m),title(x(m)=2m)实验截图（2）试验程序：n=-10:30;f1=sin(pi/5*n);f2=cos(pi/3*n);x2=f1+f2;stem(n,x2,fill),grid on,xlabel(n),title(f1+f2) 实验截图：（3）实验程序：t=linspace(-2,2,400); x1=sinc(t/pi); x

2、2=cos(10*pi*t); f2=x1.*x2; plot(t,f2,t,x1/2,r:,t,-x1/2,r:); axis(-2,2,-2,2); title(f2)实验截图：（4）试验程序：t=-2:0.01:4;a=1;a=-1;xt1=a*exp(a*t);xt2=sin(10*pi*t);a=1;a=-1/2;xt3=a*exp(a*t);xt4=sin(9*pi*t);f1=xt1.*xt2+xt3.*xt4;plot(t,f1,t,xt1.*xt2/2,r:,t,-xt1.*xt2/2,r:); axis(-2,2,-2,2); title(f1)实验截图：2. 试绘出右图所

3、示信号 f(t) 的波形图，并绘出 f(-2t+2) 及其偶分量 evf(t) 和奇分量odf (t)的波形图。 试验程序：t=linspace(-5,5,300);n=-2.*t+2;m=-t;f1=t.*heaviside(t)+(1-t).*heaviside(t-1)-heaviside(t-2);f2=n.*heaviside(n)+(1-n).*heaviside(n-1)-heaviside(n-2);f=m.*heaviside(m)+(1-m).*heaviside(m-1)-heaviside(m-2);f3=(f1+f)./2;f4=(f1-f)./2;subplot(2

4、,2,1);plot(t,f1);xlabel(t);ylabel(f(t);title(f(t);subplot(2,2,2);plot(t,f2);xlabel(t);ylabel(f(-2.*t+2;);title(f(t)=f(-2.*t+2);subplot(2,2,3);plot(t,f3);axis(-5,5,-1,1);xlabel(t);ylabel(evf(t);subplot(2,2,4);plot(t,f4);axis(-5,5,-1,1);xlabel(t);ylabel(odf(t);实验截图：3. 已知连续lti系统的冲激响应为，试绘出此系统的单位阶跃响应波形图。

5、由题意可知阶跃响应由冲击响应卷积而成则该题程序：t=linspace(-1,30,300);f=(3./2)-2.*exp(-1.*t)+(1./2).*exp(-2.*t).*heaviside(t);plot(t,f);xlabel(t);ylabel(r(t);title(单位阶跃响应);实验截图：信号与系统实习（二）2、已知调制信号试画出此调制信号的时域波形图及其频谱。(1)程序代码：w0=12.*pi;a=4;tao=1./2;syms t w;f=a.*(heaviside(t+tao./2)-heaviside(t-tao./2)*cos(w0.*t);t=-1:0.01:10;

6、figure(1);ezplot(f);ft=sym(12.*(heaviside(t+0.25)-heaviside(t-0.25)*cos(12.*pi.*t);fw=simple(fourier(ft);figure(2);ezplot(abs(fw),-4.*pi,4.*pi);axis(-4.*pi,4.*pi,-0.5,0.5);grid on;（2） 绘图及运行结果显示 3、已知某lti连续系统的系统函数为利用拉普拉斯变换法求系统的单位阶跃响应。（1）程序代码%a=4 4 4;%b=1 3 2 0;%sys=tf(b,a);syms s t;h=(4*s2+4*s+4)/(s3+

7、3*s2+2*s);hs=sym(h);ht=ilaplace(hs);et=heaviside(t);es=laplace(et);%rs=hs*es;rs=h./s;rt=ilaplace(rs)t=-10:0.1:20;%figure(1),plot(t,ht);%figure(2),plot(t,rt);（2)结果 zuoye3rt =2*t + 4/exp(t) - 3/exp(2*t) - 14、试用matlab命令画出下列各系统的零极点分布图以及对应的单位取样响应h(n)的波形，并分析系统函数的极点对时域波形的影响。(1)a=1 -0.8;b=1 0;figure(1),zpla

8、ne(b,a);grid on;n=0:30;figure(2),impz(b,a,30);grid on;(2)a=1 0.8;b=1 0;figure(1),zplane(b,a);grid on;n=0:30;figure(2),impz(b,a,30);grid on;（3）a=1 -1.2 0.72;b=1 0;figure(1),zplane(b,a);grid on;n=0:30;figure(2),impz(b,a,30);grid on;（4）a=1 -1;b=1 0;figure(1),zplane(b,a);grid on;n=0:30;figure(2),impz(b,

9、a,30);grid on;（5）a=1 -1.6 1;b=1 0;figure(1),zplane(b,a);grid on;n=0:30;figure(2),impz(b,a,30);grid on;（6）a=1 -1.2;b=1 0;figure(1),zplane(b,a);grid on;n=0:30;figure(2),impz(b,a,30);grid on;（7）a=1 -2 1.36;b=1 0;figure(1),zplane(b,a);grid on;n=0:30;figure(2),impz(b,a,30);grid on; 以下是程序代码a1=1 -0.8;a2=1

10、0.8;a3=1 -1.2 0.72;a4=1 -1;a5=1 -1.6 1;a6=1 -1.2;a7=1 -2 1.36;a=1 -0.8 1 0.8;%a=a1;a2;a3;a4;a5;a6;a7;b=1 0;m=0;for k=1:7; b=1 0;a=a(k); m=m+1; figure(m),zplane(b,a); grid on; m=m+1; n=0:30; figure(m),impz(b,a,30); grid on;end5、已知系统的微分方程为（1）试画出系统零极点分布图；（2）试画出系统的频率特性曲线；（3）求系统的冲激响应和阶跃响应；（4）若激励信号f1(t)如图

11、所示，求系统的零状态响应；（5）若激励求系统的稳态响应。（1）程序代码：a=1 6 5;b=2 6;zs=roots(b);ps=roots(a);plot(real(zs),imag(zs),blacko,real(ps),imag(ps),blackx,markersize,12);grid on;axis(-6,1,-1,1);w=-3*pi:0.01:3*pi;h=freqs(b,a,w);subplot(211),plot(w,abs(h);grid on;subplot(212),plot(w,angle(h);grid on;sys=tf(b,a);t=0:0.02:4;figu

12、re(2),impulse(sys,t);figure(3),step(sys,t);syms a t;f1=heaviside(t)-heaviside(t-1);ysz=lsim(sys,f1,t);figure(4),lsim(sys,f1,t);f2=3+cos(2.*t)+cos(5.*t);w0=0;w1=2;w2=5;h0=(2.*1i.*w0)+6)./(1i.*w0).2+6.*1i.*w0+5);h1=(2.*1i.*w1)+6)./(1i.*w1).2+6.*1i.*w1+5);h2=(2.*1i.*w2)+6)./(1i.*w2).2+6.*1i.*w2+5);y0=abs(h0)*cos(w0.*t+angle(h);y1=abs(h1)*cos(w1.*t+angle(h);y