MATLAB 實驗報告1

1、matlab 實驗報告 學(xué)院 _機電學(xué)院_ 姓名_ _ 班號_ 學(xué)號_信號與系統(tǒng)實驗（一）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. 已知連續(xù)lti系統(tǒng)的沖激響應(yīng)為，試繪出此系統(tǒng)的單位階躍響應(yīng)波形圖。

5、由題意可知階躍響應(yīng)由沖擊響應(yīng)卷積而成則該題程序：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(單位階躍響應(yīng));實驗截圖：信號與系統(tǒng)實習(xí)（二）2、已知調(diào)制信號試畫出此調(diào)制信號的時域波形圖及其頻譜。(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） 繪圖及運行結(jié)果顯示 3、已知某lti連續(xù)系統(tǒng)的系統(tǒng)函數(shù)為利用拉普拉斯變換法求系統(tǒng)的單位階躍響應(yīng)。（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)結(jié)果 zuoye3rt =2*t + 4/exp(t) - 3/exp(2*t) - 14、試用matlab命令畫出下列各系統(tǒng)的零極點分布圖以及對應(yīng)的單位取樣響應(yīng)h(n)的波形，并分析系統(tǒng)函數(shù)的極點對時域波形的影響。(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、已知系統(tǒng)的微分方程為（1）試畫出系統(tǒng)零極點分布圖；（2）試畫出系統(tǒng)的頻率特性曲線；（3）求系統(tǒng)的沖激響應(yīng)和階躍響應(yīng)；（4）若激勵信號f1(t)如圖

11、所示，求系統(tǒng)的零狀態(tài)響應(yīng)；（5）若激勵求系統(tǒng)的穩(wěn)態(tài)響應(yīng)。（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