【正文】 y=simple(y)。)*(4*t/a1)。 %y=y1sym(39。Heaviside(ta1)39。Heaviside(t+a1)39。a=T/2。 x=sawtooth(t2*T/3,1)。) axis([0,20,0,]) end function x=sjb_timefun T=10。 title(39。T=10，脈寬為 1的周期三角波脈沖 39。t=2*T::2*T。 S7=[S6,S5(2:end)]。 N=Nf*2*pi/T。end if nargout==0 S1=fliplr(A_sym) S1(1,k+1)=A_sym(1) S2=fliplr(1/2*S1) S3=fliplr(1/2*B_sym) S3(1,k+1)=0 S4=fliplr(S3) S5=S2i*S4。 for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn))。 Bs=int(2*y*sin(2*pi*n*t/T)/T,t,0,T)。 A0=2*int(y,t,0,T)/T。Nn=32。)。Nf=input(39。 源程序五： （典型周期脈沖的頻譜分析 —— 三角波） function [A_sym,B_sym]=sjb syms t n k y 26 T=10。duty=50。 function x=fangbo_2 T=5。 y=subs(y,a1,a)。Heaviside(t+a1)39。)。)*2sym(39。 y1=sym(39。) axis([20,20,0,]) end function y=fangbo_1 syms a a1 T=5。 title(39。T=5，占空比為 50%的周期方波脈沖 39。t=2*T::2*T。 S7=[S6,S5(2:end)]。 N=Nf*2*pi/T。end if nargout==0 S1=fliplr(A_sym) S1(1,k+1)=A_sym(1) S2=fliplr(1/2*S1) S3=fliplr(1/2*B_sym) S3(1,k+1)=0 S4=fliplr(S3) S5=S2i*S4。 for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn))。 Bs=int(2*y*sin(2*pi*n*t/T)/T,t,0,T)。 A0=2*int(y,t,0,T)/T。Nn=32。)。Nf=input(39。 end subplot(3,1,1) plot(t,x) hold on axis([20,20,0,]) 源程序四： （典型周期脈沖的頻譜分析 —— 方波） function [A_sym,B_sym]=fangbo syms t n k y T=10。 for i=1:n。 t=n*T::n*T。tao=T/5。 y=simple(y)。)*A。 y=y1sym(39。Heaviside(t+a1)39。a=tao/2。 A=1。周期與脈沖寬度之比 M=39。)。 axis([150,150,4,4]) end function y=fun_in(t) syms a a1 T=input(39。 x=fun_mc subplot(3,1,2) stem(k2,abs(S7))。 k2=N:2*pi/T:N。 S6=fliplr(S5)。 B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn))。 A_sym(1)=double(vpa(A0,Nn))。 As=int(2*y*cos(2*pi*n*t/T)/T,t,0,T)。end y=fun_in(t)。 if nargin5。pleas Input 信號(hào)的周期 T=39。)。Nf=input(39。 for n=1:2:N fN=fN+cos(pi*n*t)*sinc(n/2)。 % N為輸入要達(dá)到的最高次諧波的次數(shù) c0=。N=39。 源程序二：（吉布斯效應(yīng)） t=2::2。Heaviside()39。)*h。 x1=sym(39。 y=h.*(e1e2)。 e1=1/2+1/2.*sign(t+tao/2)。 tao=*T。T=5?；?+2次諧波 +3次諧波 +6次諧波 39?；?+2次諧波 +3次諧波 39?；?+2次 諧波 39。基波 39。 f8=f7+f4+f5。 f6=f1+f2。 f4=2*(.*cos(2*pi*4*t/5)+0.*sin(2*pi*4*t/5))。 f2=2*(.*cos(2*pi*2*t/5)+0.*sin(2*pi*2*t/5))。disp(d) t=8*a::Ta。 end if nargout==0 c=A_sym。 for k=1:Nf %Nf 為最高次諧波的次數(shù) A_sym(k+1)=vpa(subs(As,n,k),Nn)。 Bs=int(2*x*sin(2*pi*n*t/T)/T,t,a,Ta)。 %求傅里葉級(jí)數(shù)展開式的系數(shù)。end x=time_fun_x(t)。end if nargin5。 if nargin4。tao=*T。對于那些只需要該參數(shù)就可以實(shí)現(xiàn)的程序也給只出主要程序。 沒有他們的幫助我不會(huì)那么快克服那些困難，也不會(huì)這么快學(xué)到這么多的知識(shí)。通過這次課設(shè)讓我明白了理論聯(lián)系實(shí)踐的重要性，書本上的理論知識(shí)學(xué)了不少，我們必須得應(yīng)用到實(shí)踐當(dāng)中，做到學(xué)以致用，這樣我們才能有不斷的創(chuàng)新。這次課程設(shè)計(jì)讓我深深的明白了自己以后該做什么，該怎么去做。在以后的學(xué)習(xí)中應(yīng)該多多的注意實(shí)踐知識(shí)的訓(xùn)練和積累。要想真正掌握它還需要繼續(xù)努力學(xué)習(xí)它。對 word 也有了進(jìn)一步的掌握。對以前的理論知識(shí)有了更進(jìn)一步的認(rèn)識(shí)和理解。以后應(yīng)該多注意，充分合理的利用網(wǎng)絡(luò)，通過網(wǎng)絡(luò)來學(xué)習(xí)東西。 在寫相關(guān)源程序的時(shí)候，我還收索了大量的網(wǎng)站，在網(wǎng)上收索了很多關(guān)于 MATLAB 的資料。設(shè)計(jì)的核心內(nèi) 容就是 利用 MATLAB 強(qiáng)大的圖形處理功能，符號(hào)運(yùn)算功能以及數(shù)值計(jì)算功能，實(shí)現(xiàn)連續(xù)時(shí)間周期信號(hào)頻域分析的仿真波形。主瓣寬度與周期和脈沖寬度的關(guān)系與方波時(shí)的規(guī)律基本一致，這里不再贅述。 在 附錄的源程序五 中， 只需要 給出不同的 T值和 width值， 就會(huì)得到 三角波的 不同的 雙邊和單邊頻譜圖。 15 周期三角波脈沖頻譜的 MATLAB 實(shí)現(xiàn) 周期三角波脈沖如圖 ，周期 T=5，其幅度為 1。 周期方波信號(hào)頻譜與周期矩形脈沖信號(hào) 具有相同的規(guī)律， 這里不再贅述。 將在下一小節(jié)中給出不同參數(shù)時(shí)的頻譜圖。 1 0 8 6 4 2 0 2 4 6 8 101 0 . 8 0 . 6 0 . 4 0 . 200 . 20 . 40 . 60 . 81周期方波脈沖 圖 周期方波脈沖 編寫 ，源程序文件見附錄程序 四 。脈沖周期 T 愈大，譜線間隔愈小，頻譜越稠密；反之，則越稀疏。如果周期無限增長（這時(shí)就成為非周期信號(hào)），那么，相鄰譜線的間隔將趨近于零，周期信號(hào)的離散頻譜就過渡到非周期信號(hào)的連續(xù)頻譜。因?yàn)棣又挡蛔?，改?T值就等于改變了 T=ατ中比例系數(shù)α的大小。x1/2)+1/2.*(x==1/2)固定為 1/2，然后分別將 a,b值變?yōu)椋?a=4,b=4 和 a=6,b=6，來研究方波周期對其頻譜的影響。T0=ba。在 10 MATLAB 代碼中 a=5?？梢?，信號(hào)的頻帶寬度與脈沖寬度 τ 成反比 [2]。 9 圖 不同τ值時(shí)的頻譜 容易看出，在 T不變的情況下，減小τ值，可以使頻譜變得更密集，增大τ值則可以使頻譜變得稀疏，因此，需要在不同的情況下選擇不 同的τ值，才能是系統(tǒng)變得更加符合實(shí)際需要。這個(gè)語句是控制τ值的，現(xiàn)在的參數(shù)是 1/2。xsqual=(x)1/2.*(x==1/2)+1.*(x1/2amp。改變脈沖寬度，就是在圖 中 T 值不變的情況下，改變的τ值的大小，同時(shí)τ必須小于 T。 單邊頻譜和雙邊頻譜的區(qū)別就是求值的范圍不同，單邊頻譜求的是頻率大于 0的情況，而雙邊頻譜求的是所有頻率的情況，即包括頻率小于 0 的情況，這個(gè)區(qū)別在上面的兩張圖中可以非常明顯地看出來。 由圖可見，周期信號(hào)的譜線只出現(xiàn)在頻率為 0, ,2 ,...?? 等原周期信號(hào)頻率的整數(shù)倍的離散頻率上，即周期信號(hào)的頻譜是離散譜。圖 相位譜為雙邊相位譜。 類似地，也可畫出各諧波初相角 n? 與角頻率的關(guān)系圖，如圖 n?與角頻率的關(guān)系圖，稱為相位頻譜，簡稱相 位譜。圖 譜（用 nA 繪制的頻譜）。圖中每條豎線代表該頻率分量的幅度，稱為譜線。 2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 200 . 51諧波次數(shù) N = 5 0 時(shí)的合成波形2 1 . 5 1 0 . 5 0 0 . 5 1 1 .