MATLAB實(shí)驗(yàn)報(bào)告

1、matlab 實(shí)驗(yàn)報(bào)告 學(xué)院 _機(jī)電學(xué)院_ 姓名_ _ 班號_ 學(xué)號_信號與系統(tǒng)實(shí)驗(yàn)（一）1. 試?yán)L出下列各信號的波形圖：（1）實(shí)驗(yàn)程序：n=-3:15;m=n-1;a=2;x1=a.m;stem(n,x1,'fill'),grid on,xlabel('m'),title('x(m)=2m')實(shí)驗(yàn)截圖（2）試驗(yàn)程序：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+

2、f2') 實(shí)驗(yàn)截圖：（3）實(shí)驗(yàn)程序：t=linspace(-2,2,400); x1=sinc(t/pi); x2=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')實(shí)驗(yàn)截圖：（4）試驗(yàn)程序：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.

3、*xt4;plot(t,f1,t,xt1.*xt2/2,'r:',t,-xt1.*xt2/2,'r:'); axis(-2,2,-2,2); title('f1')實(shí)驗(yàn)截圖：2. 試?yán)L出右圖所示信號 f(t) 的波形圖，并繪出 f(-2t+2) 及其偶分量 evf(t) 和奇分量odf (t)的波形圖。 試驗(yàn)程序：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).*

4、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,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.*

5、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)');實(shí)驗(yàn)截圖：3. 已知連續(xù)lti系統(tǒng)的沖激響應(yīng)為，試?yán)L出此系統(tǒng)的單位階躍響應(yīng)波形圖。由題意可知階躍響應(yīng)由沖擊響應(yīng)卷積而成則該題程序：t=linspace(-1,30,300);f=(3./2)-2.*exp(-1.*t)+(1./2

6、).*exp(-2.*t).*heaviside(t);plot(t,f);xlabel('t');ylabel('r(t)');title('單位階躍響應(yīng)');實(shí)驗(yàn)截圖：信號與系統(tǒng)實(shí)習(xí)（二）2、已知調(diào)制信號試畫出此調(diào)制信號的時(shí)域波形圖及其頻譜。(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;figure(1);ezplot(f);ft=sym('12.*(heavis

7、ide(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） 繪圖及運(yùn)行結(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+3*s2+2*s);hs=sym('h');ht

8、=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)的零極點(diǎn)分布圖以及對應(yīng)的單位取樣響應(yīng)h(n)的波形，并分析系統(tǒng)函數(shù)的極點(diǎn)對時(shí)域波形的影響。(1)a=1 -0.8;b=1 0;figure(1),zplane(b,a);grid o

9、n;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,a,30);grid on;

10、（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;待添加的隱藏文字內(nèi)容1figure(2),impz(b,a,30);grid on; 以下是程序代碼a1=1 -0.8;a2=1 0.8

11、;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)零極點(diǎn)分布圖；（2）試畫出系統(tǒng)的頻率特性曲線；（3）求系統(tǒng)的沖激響應(yīng)和階躍響應(yīng)；（4）若激勵(lì)信號f1(t)如圖所示，

12、求系統(tǒng)的零狀態(tài)響應(yīng)；（5）若激勵(lì)求系統(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;

13、sys=tf(b,a);t=0: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(