MATLAB实验报告1

《MATLAB实验报告1》由会员分享，可在线阅读，更多相关《MATLAB实验报告1（19页珍藏版）》请在装配图网上搜索。

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);待添加的隐藏文字内容3grid on;n=0:30;figure(2),impz(b,a,30);grid on; 以下是程序代码a1=1

10、 -0.8;a2=1 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）若

11、激励信号f1(t)如图所示，求系统的零状态响应；（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:

12、0.02:4;figure(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);y2=abs(H2)*cos(w2.*t+angle(H);y=y0+y1+y2;figure(5),plot(t,y);grid on;（2）绘图结果实验心得 通过实验，我了解了matlab软件的部分功能，并通过它完成信号波形图，相当于是对所学理论知识的一种实践和加深吧，从理论上和实际上都加深了我对这门课程的理解。而且这个软件它的强大功能，对于数学计算等类型的问题有很好的辅助作用。