【正文】 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% T4_2a=[ ]。filled39。 %函數(shù)表達(dá)式為 f(t)=[u(t0)u(t1)][u(t1)u(t2)]，兩個(gè)門t2=linspace(1,2,101)。)。grid on。14（2）執(zhí)行結(jié)果：（3）分析： ，操作窗口求出了周期矩形脈沖信號(hào)的Fourier級(jí)數(shù)（復(fù)指數(shù)形式）的系數(shù)，（其中n=9，9次諧波）（復(fù)指數(shù)形式）系數(shù)的振幅譜和相位譜圖像，（其中n=9，9次諧波）。)。ylabel(39。原信號(hào)頻譜 39。 Fy=fourier(y)。 %化為頻域 ff1=subs(f1,t,t1)。 xlabel(39。t/s39。title(39。 b=[0 , 1]。fx139。19(2) 執(zhí)行結(jié)果:（3）分析：，表示出了RC=1，2,3,4,5，時(shí)的幅頻曲線和相頻曲線，可以觀察到RC越大，曲線開始下降的頻率越低，即RC與截止頻率成反比。)。)。)。for a=3::3 R11=subs(Fe,v,a)。)。)。for a=3::3 R11=subs(Fe,v,a)。 axis([3,3,4,]) subplot(2,3,5)。)。 C(i)=angle(R11)。 subplot(2,3,5)。22（2）執(zhí)行結(jié)果23（3）分析： 信號(hào)在通過(guò)線性無(wú)失真系統(tǒng)后，波形不失真；幅度譜增加了增益系數(shù)K，相位譜由常數(shù)變?yōu)橐淮魏瘮?shù)；時(shí)域上有時(shí)延，頻域上有相移。 ezplot(f,[5,5]) title(39。fc=0。end subplot(3,2,6)。 ff=fourier(f)。ezplot(Fc,[40,40])title(39。 f1=ifourier(Fc)/fs。（1）源程序28%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 對(duì)傳輸函數(shù) ，選取不同的參數(shù)設(shè)計(jì)% 低通、高通、帶通濾波器。bode(bt,at)。Rp=3。figure(5)。figure(7)。hold on。)。30(2)執(zhí)行結(jié)果31(3)分析 ，使用方法二得出的傳遞函數(shù)和題目所給的二階Butterworth 低通濾波器相同，均為，截止頻率為W0=1 rad/s。 % 利用 freqs 求出系統(tǒng)的頻率響應(yīng)，數(shù)值解，比起bode函數(shù)來(lái)利于作圖和標(biāo)注 magH1=abs(H1)。)。)。grid。[z,p,k]=buttap(2)。通過(guò)本實(shí)驗(yàn)了解低通、高通、帶通、全通濾波器的性能及特點(diǎn)。采樣信號(hào)頻譜，fs2fm(fs=)39。 G=heaviside(w4)+heaviside(w+4)。subplot(4,2,1)。)。 ezplot(fc,[40,40]) title(39。close。 plot(a,C)title(39。 ezplot(r1,[3,3])title(39。)。 subplot(2,3,6)。 subplot(2,3,4)。幅度譜39。 i=i+1。 r1=ifourier(R1,t)。幅度譜39。 end b=3::3。 %系統(tǒng)函數(shù)(1)R1=Fe*H1。 Fe=fourier(e,v)。)。xlabel(39。提示：為最大值的處對(duì)應(yīng)的頻率為通帶截止頻率，首先求取并找到和RC關(guān)系，然后根據(jù)題意選定即可確定RC值。)。 %低通后取出的信號(hào)，時(shí)域subplot(4,2,7)。 Ff1a=abs(Ff1)。 調(diào)制后的信號(hào)頻譜 39。)。ylabel(39。 %畫原信號(hào)的時(shí)域圖像plot(t1,x1)。 %符號(hào)函數(shù)可以表示成2*u(t)1FFw=abs(Fw)。ylabel(39。)。n39。圖10（1）源程序：13%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 求如圖10所示周期矩形脈沖信號(hào)%的Fourier級(jí)數(shù)表達(dá)式，畫出頻譜圖，%并用前N次諧波合成的信號(hào)近似。for i=N+1:N+length(n) y(i)= a1*y(iN:i1)39。 title(39。系統(tǒng)的階躍響應(yīng)y(t)39。title(39。title(39。) %縱坐標(biāo)定義title(39。運(yùn)用課堂上學(xué)到的理論知識(shí)，從RC、RL一階電路的響應(yīng)中正確區(qū)分零輸入響應(yīng)、零狀態(tài)響應(yīng)、自由響應(yīng)與受迫響應(yīng)。 %求卷積和k=0:length(f)1。filled39。grid on % 子圖4：作出序列圖f[k+2]subplot(2,3,4)。f39。ylabel(39。 k=6:5。stem(k,f,39。)。)。)。)。xlabel(39。 %脈寬為1f=rectpuls(t,w)。f=A*sin(omega0*t+phi)39。f=A*cos(omega0*t+phi)。alpha=。)。t39。ylabel(39。xlabel(39。圖3（1）源程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%t=20::20。grid on % 子圖2：作出f（t）subplot(2,3,2)。grid on % 子圖4：作出f（12*t）subplot(2,3,4)。grid on4 （2）執(zhí)行結(jié)果如上圖，作出了1個(gè)原函數(shù)圖像，和4個(gè)變換后的函數(shù)圖像實(shí)驗(yàn)二 離散時(shí)間信號(hào)時(shí)域分析一、 實(shí)驗(yàn)內(nèi)容用MATLAB表示離散信號(hào)：。xlabel(39。f= A.*sin(k)。stem(k,f,39。filled39。)。 title(39。)。 title(39。dt=。b=[1]。te=5。dt=。（1）源程序 11%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 假設(shè)某系統(tǒng)的單位函數(shù)響應(yīng) ，% 系統(tǒng)激勵(lì)信號(hào) ，求系統(tǒng)的零狀態(tài)響應(yīng)，% 并畫圖表示。 b=[ 0 1 0]。)。 %劃分為兩個(gè)門信號(hào)疊加n=0:9。title(39。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%用前N次諧波合成的信號(hào)近似。試用fourier()函數(shù)求下列信號(hào)的傅里葉變換,并畫出(1) (2) (1)源程序：15%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% (1)求下列信號(hào)的傅里葉變換 ,并作圖 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% syms t 。ylabel(39。x(t)39。)。 % 調(diào)制后信號(hào)的頻譜Fya=abs(Fy)。 subplot(4,2,5)。w (rad/s)39。)。低通取出的信號(hào)頻譜 39。freqs(b,a,0::2*pi)。)。RC=1/wc, 取wc =1000 rad/s,則 RC=103, 。 axis([3,3,0,2])。 axis([3,3,2])。axis([3,3,3,3])。 C(i)=angle(R11)。 axis([3,3,1,2])。axis([3,3,3,3])。 C(i)=angle(R11)。 ezplot(39。 %（4）兩個(gè)條件均不滿足時(shí)clear。 i=i+1。 ezplot(39。 ，信號(hào)通過(guò)只滿足幅值條件 的系統(tǒng)，波形失真，幅度譜不變，相位譜變化。時(shí)域信號(hào)39。for n=6:6 fc=fc+(pi/2)*fs*(heaviside(w2n*2*pi*fs)+heaviside(w+2n*2*pi*fs))。ezplot(fc,[40,40])title(39。 subplot(4,2,2)。重建信號(hào)頻譜,fs=2fm(fs=3)39。 subplot(4,2,6)。繪制波特圖。title(39。Rs=60。bode(b,a)。bode(b,a)。 % 畫出幅頻曲線end title(39。 %模擬濾波器定階 得到n=11。grid。 % 將傳遞函數(shù)由零極點(diǎn)式轉(zhuǎn)化成多項(xiàng)式形式 [H1,w1]=freqs(b,a,w)。帶阻濾波器39。Chebyshev Highpass39。)。W0=1。27實(shí)驗(yàn)八 傳輸函數(shù)與系統(tǒng)特性一、實(shí)驗(yàn)?zāi)康膶W(xué)習(xí)由系統(tǒng)函數(shù)確定系統(tǒng)頻率特性的方法。ezplot(fc,[40,40]) title(39。)。 f=sin(2*t)/(2*t)。采樣信號(hào)頻譜，采樣頻率339。 end subplot(3,2,3)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear。enda=3::3。 subplot(2,3,4)。幅度譜39。 end b=3::3。 r1=ifourier(R1,t)。 ezplot(Fe,[3,3]) title(39。 C(i)=angle(R11)。D=abs(R1)。 ezplot(Fe,[3,3])title(39。 i=i+1。 H1=2*exp(j*v*1)。 %第一幅圖像ezplot(e,[3,3]) %ezplot畫符號(hào)函數(shù)的幅度譜axis([3,3,2])。fy139。plot(t1,fx1)。(b) 信號(hào)包含了一個(gè)低頻分量和一個(gè)高頻分量，確定適當(dāng)?shù)腞C值，濾除信號(hào)中的高頻分量并畫出信號(hào)和在s范圍內(nèi)的波形。w (rad/s)39。 %通過(guò)低通濾波器，取出原信號(hào) f2 = ifourier(Ff2,t)。)。title(39。y(t)39。)。 %原信號(hào)的時(shí)域表示，數(shù)值表示subplot(4,2,1)。Fw=fourier(2*heaviside(t)1)。)。Theta39。xlabel(39。二、