【正文】 ,(t0t1)/dt+1)=1/dt。 axis([t1,t2,0,1/dt])。時間 (t)39。ylabel(39。)。單位沖激信號 39。 時域波形 ： ② 單位脈沖序列 ( ) ( )x n n?? 的實現(xiàn)程序： k1=3。 k=k1:k2。%單位脈沖出現(xiàn)的位置 f=[(kn)==0]。filled39。title(39。) xlabel(39。)。幅值 f(k)39。 時域波形 ： （ 3）單邊衰減指數(shù)信號 ( ) ( )tf t e t? ??? 和實指數(shù)序列 ()nxn a? 的實現(xiàn) ① 單邊衰減指數(shù)信號 ( ) ( )tf t e t? ??? 的實現(xiàn)程序： clear t1=1。dt=。 A1=1。 n=A1*exp(a1*t)。 axis([t1,t2,0,1])。時間 (t)39。ylabel(39。)。單邊衰減指數(shù)信號 39。 時域波形 ： ② 實指數(shù)序列 ()nxn a? 的 實現(xiàn)程序： clf。k2=10。 a=。 f=A*a.^k。filled39。 title(39。)。(k)39。ylabel(39。)。 程序 clear all。 b=3。 u=stepfun(t,t0)。 for i=1:n u(i)=b*u(i)*(t(i)t0)。 %產(chǎn)生一個正弦信號 f=y+u。 xlabel(39。)。幅值 f(t)39。title(39。)。 程序 clear all。 b=3。 u=stepfun(t,t0)。 for i=1:n u(i)=b*u(i)*(t(i)t0)。 f=y.*u。xlabel(39。)。幅值 f(t)39。title(39。) 時域波形 ： ③ 實現(xiàn)兩個連續(xù)信號的移位，即 00( ) ( )f t t f t t??， 或 者 ，常數(shù) 0 0t? 程序： clear all。 y=sin(2*pi*(t))。 plot(t,y,39。,t,yl,39。)。幅值 f(t)39。xlabel(39。)。信號的位移 39。 時域波形 ： （ 2）序列的卷積運算 實現(xiàn)兩個連續(xù)信號的卷積運算， 12( ) ( ) * ( )f t f t f t? 程序： %連續(xù)函數(shù)卷積計算 a=1000。 f1=stepfun(t1,0)。 subplot(231)。axis([5,5,0,])。時間 (t)39。 ylabel(39。)。單位階躍函數(shù) 39。 subplot(232)。ylabel(39。)。單位沖激函數(shù) 39。 y=conv(f1,f2)。t=10:1/a:10。plot(t,y)。 title(39。)。y(t)39。 f11=conv(f1,f1)。 subplot(234)。title(39。)。f11(t)39。axis([5,5,0,5000])。plot(t,f22)。f2 與 f2 的卷積 39。ylabel(39。)。 i 沖激響應 MATLAB程序： clear all。b=[3,9]。 t=0::10。 plot(t,y)。時間 (t)39。 ylabel(39。)。單位沖激響應 39。 系統(tǒng)的沖激響應圖： ii 階躍響應 MATLAB程序： clear all。b=[3,9]。 t=0::10。 plot(t,y)。時間 (t)39。 ylabel(39。)。單位階躍響應 39。 系統(tǒng)的階躍響應圖 ： ② 求系統(tǒng) y(2)(t)+y(t)=costu(t)， y(0+)=y(1)(0+)=0的全響應。a=[1 0 1]。 t=0::10。 y=lsim(sys,x,t)。 xlabel(39。)。y(t)39。title(39。)。a=[1 0 1]。 sys=ss(A,B,C,D)。 x=cos(t)。 y=lsim(sys,x,t,zi)。 xlabel(39。)。y(t)39。title(39。)。 MATLAB程序： p=。 f1=sin(k1)。 f2=1*exp(2*k1)。 f=f*p。 k3=length(f1)+length(f2)2。 subplot(2,2,1)。 title(39。)。t39。f1(t)39。 subplot(2,2,2)。 title(39。)。t39。f2(t)39。 subplot(2,2,3)。 h=get(gca,39。)。 set(gca,39。,h)。f(t)=f1(t)*f2(t)39。 xlabel(39。),ylabel(39。)。 H(s)= ssssssssss 23456 234 ????? ???? MATLAB程序 ： a=[1 463 12131 0]。 sys=tf(b,a)。 x=exp(t)。 plot(t,y)。時間 (t)39。ylabel(39。)。零狀態(tài)響應 39。 i 已知連續(xù)時間信號 f(t)=e2|t|，通過程序完成信號 f(t)的傅里葉變換 MATLAB程序 ： syms t x=exp(2*abs(t))。 subplot(2,1,1)。 subplot(2,1,2) 。 信號 f(t)的傅里葉變換如圖所示： ii 試畫出信號 f(t)= 32 ()3 tet?? 的波形及其幅頻特性曲線。 f=2/3*exp(3*t)*sym(39。) 。 subplot(2,1,1) 。 subplot(2,1,2) 。 信號 f(t)= 32 ()3 tet?? 的波形及其幅頻特性曲線如圖所示 II 傅里葉逆變換。 MATLAB程序 R=。 syms t w ifourier(1/(1+w^2),t) 運行結(jié)果 ： ans =1/2*exp(t)*heaviside(t)+1/2*exp(t)*heaviside(t) III 傅里葉變換數(shù)值計算。 MATLAB程序 ： R= 。 f=stepfun(t,1)stepfun(t,1) 。 %頻率寬度 N=500 。 W=k*W1/N 。*W)*R 。 %求 F（ jw）的實部 W=[fliplr(W),W(2 :501)] 。 %頻譜圖幅值左右對稱 subplot(2,1,1) 。 xlabel(39。) 。f(t)39。axis([2,2,2]) 。f(t)= ε(t+1)ε(t1)39。 subplot(2,1,2) 。 title(39。) 。 i f(t)= 21 ()2 tet??的頻譜。 t=5:r:5。 W=2*pi。 w=k*W/N。 F=r*f1*exp(j*t39。 F1=abs(F)。 subplot(3,1,1)。grid xlabel(39。)。f(t)39。 title(39。)。 plot(w,F1)。w39