【正文】 sp(ｃ）。disp(ｄ)。ｆ 0= ｃ(１)。% 2 次諧波f3=ｃ(4）、＊ｃos(2＊pi*3＊t／T)+d(4)、＊sin(2*ｐi*3＊t/T)。% ４次諧波f５=c（6）、*ｃos(2*ｐi*5*ｔ/T）+d(6）、*ｓｉn(2＊ｐi*5*t／T)。% 8次諧波f9 ＝ ｃ(1 ０)、＊ c ｏ ｓ(２ * ｐ i*9 ＊ ｔ /T)+d(10)、* ｓ in(2 ＊ ｐ i ＊ 9 ＊t/T）。％ 直流+基波+2 次諧波+3 次諧波f13＝f12＋f4+f5+f6；% 直流＋基波+2 次諧波+３次諧波+４次諧波+5 次諧波+6 次諧波f14＝f13+f７+f8+ｆ9+f１0；%0～10 次subｐlot(2,2，1)ｐlｏt(t，f0＋f1）,hoｌｄ ｏny＝tiｍe_ｆun＿ｅ(t)。ｐlot(t,y，’r：’）tiｔle(“1—3 次諧波+直流”）axis([—8,8,0、5,5］）sｕｂploｔ(２,2,3)ploｔ(t,f13），hoｌd ony＝ｔime_fun＿e（t)；plot(ｔ,y,’ｒ:’）ｔitle（“１—6 次諧波＋直流＇)axis(［8，8，－0、５,5］)ｓubpｌot（２,2，4)plｏt（t,ｆ14）,holｄ ony=timｅ_fun_e（t）。t=—3*Ｔ:0、0１:3*T。%連續(xù)時(shí)間函數(shù)—周期矩形脈沖 funcｔion x＝tiｍｅ_fun＿ｘ（t）% 該函數(shù)就是 CTFSｈcｈsym、ｍ得子函數(shù)。ｘ1=sｙm(”Heaｖisidｅ(ｔ)“)*h。%? 函數(shù)得輸入輸出都就是數(shù)值量 ％ Ｎｎ輸出數(shù)據(jù)得準(zhǔn)確位數(shù) % mys_A? 第１元素就是直流項(xiàng),其后元素依次就是１,２,次諧波 coｓ項(xiàng)展開系數(shù) ％ Ｂ_ｓyｍ 數(shù)系開展項(xiàng) nis 波諧次、３,２,1 是就次依素元、,４，3,2 第?％TT=m*taｏ, 信號周期 ％ ? ｆN? 諧波得階數(shù) ％m(m=T/tao)周期與脈沖寬度之比,如 m=4，8，16,1００等 ％tao脈寬：ｔａo＝Ｔ/ｍsｙｍs tnｙiｆ naｒgin＜3Nｆ=inｐut（’plｅａse Inｐut 所需展開得最高諧波次數(shù):Nf＝’)。if nａrｇｉn〈5Ｎn=32。Ａs=2/Ｔ＊iｎｔ（y*cｏs（2*pi*n＊t/Ｔ）,ｔ,0,T/2)。Ａ＿sym(1)=double(vpa(Ａ０，Ｎｎ））。B＿sym(k+1)=ｄouble(ｖｐa(sｕbs(Bs，n，k),Ｎn))。%Ａ＿sym 得 1*k 陣擴(kuò)展為 1＊（k+1)陣An=ｆｌｉplr（Aｎ)。%對 B_sym 陣左右對稱交換Ｂn（1,ｋ+1)=0。%對擴(kuò)展后得 S３陣左右對稱交換回原位置FnR=Aｎ/2—i＊Ｂn/２；% 用三角函數(shù)展開系數(shù) A、Ｂ值合成付里葉指數(shù)系數(shù)ＦnL=fｌｉplr（FnＲ）；N=Nf*2*pi/T。Fn＝[FnＬ，F(xiàn)nＲ（2:eｎｄ）]。% 調(diào)用連續(xù)時(shí)間函數(shù)周期矩形脈沖sｕｂplot(2,1，1)stｅm(k２,aｂs(Fn)）；％畫出周期矩形脈沖得頻譜(T=M*tao）ｔiｔle(＇連續(xù)時(shí)間函數(shù)周期三角波脈沖得雙邊幅度譜’)axｉs([－8０，80,0,０、12]）ｌine（[80，8０］,[0,０],＇colｏr“，’ｒ”）ｌｉｎe([０,０],[0，０、１２]，’colｏｒ’,＇r“）end ｆunctioｎ x=time＿fun_e(t)% 該函數(shù)就是ＣTFSｓhbpsｙm、ｍ得子函數(shù)。％ t組數(shù)間時(shí)是就?％ T 0=T/oat=yｔuｄ期周是就??T=５。ｘ=rectpuｌs（t,tａo)；％產(chǎn)生一個(gè)寬度 tａo=1 得矩形脈沖 sｕbplot(2,2,2）ｐlｏt(ｔ,x）ｈoｌｄ on x=rectｐuls(t—５，tao)。%產(chǎn)生一個(gè)寬度ｔao＝１得矩形脈，中心位置在 t=—５處 ploｔ(t，x)ｔitlｅ(”周期為 T＝5,脈寬 tａo=１得矩形脈沖＇)ａxis([10,10,0,2］)funｃtｉoｎ y＝time_fｕｎ＿s(t）syms t y＝1—abs(t)。ｘ=x1—sｙｍ(Heａvisidｅ(t2）’)。ｅzplｏt(t，y,[—10,1０])grid第四篇：信號與系統(tǒng)實(shí)驗(yàn)報(bào)告總結(jié)信號與系統(tǒng)實(shí)驗(yàn)實(shí)驗(yàn)一常用信號的觀察方波：正弦波：三角波：在觀測中，虛擬示波器完全充當(dāng)實(shí)際示波器的作用，在工作臺上連接AD1為示波器的輸入，輸入方波、正弦波、三角波信號時(shí)，可在電腦上利用軟件觀測到相應(yīng)的波形，其縱軸為幅值可通過設(shè)置實(shí)現(xiàn)幅值自動調(diào)節(jié)以觀測到最佳大小的波形，其橫軸為時(shí)間，宜可通過設(shè)置實(shí)現(xiàn)時(shí)間自動調(diào)節(jié)以觀測到最佳寬度的波形。方波基波信號：基波信號為與原方波50Hz信號相對應(yīng)的頻率為50Hz的正弦波信號，是方波分解的一次諧波信號。方波三次諧波信號：三次諧波信號頻率為150Hz為原方波信號的三倍。方波四次諧波信號:四次諧波信號的頻率為200Hz為原方波信號的四倍。方波五次諧波信號：五次諧波頻率為250Hz為原方波信號的五倍。綜上可知：50Hz方波可以分解為DC信號、基波信號、二次、三次、四次、五次諧波信號…，無偏移時(shí)即無DC信號，DC信號幅值為0。二次諧波、三次諧波、四次諧波、五次諧波依次頻率分別為原方波信號的二、三、四、五倍，且幅值依次衰減，直至五次諧波信號時(shí)幾乎可以忽略。方波基波加三次諧波信號：基波疊加上三次諧波信號時(shí)，幅值與方波信號接近，形狀還有一定差異，但已基本可以看出疊加后逼近了方波信號。綜上所述：方波分解出來的各次諧波以及DC信號，疊加起來以后會逼近方波信號，且疊加的信號越多，越是接近方波信號。三角波DC信號：三角波基波信號：三角波二次諧波信號：三角波三次諧波信號：三角波四次諧波信號：三角波五次諧波信號：三角波基波加三次諧波信號：三角波基波加三次諧波加五次諧波信號：三角波信號的分析與方波信號的分析基本一致，可以看出三角波也可以分解為多個(gè)諧波，并且相應(yīng)的多個(gè)多次諧波可以合成三角波信號，且參與合成的波形越多，合成波越是逼近三角波信號。且任何周期性函數(shù)均可由鎖對應(yīng)的直流分量和各次諧波項(xiàng)所合成，參與合成的信號越多，結(jié)果越逼近周期性函數(shù)的圖形。2．分析理論合成的波形與實(shí)驗(yàn)觀測到的合成波形之間誤差產(chǎn)生的原因。實(shí)驗(yàn)二 零輸入、零狀態(tài)級完全響應(yīng)零輸入響應(yīng)下降沿采樣：零輸入響應(yīng)上升沿采樣：可見，零輸入響應(yīng)按照指數(shù)形式下降，最終降為零。第五篇：信號與系統(tǒng)仿真實(shí)驗(yàn)報(bào)告信號與系統(tǒng)仿真實(shí) 驗(yàn) 報(bào) 告班級： 學(xué)號： 姓名： 學(xué)院：實(shí)驗(yàn)一一、實(shí)驗(yàn)者姓名：二、實(shí)驗(yàn)時(shí)間：三、實(shí)驗(yàn)地點(diǎn)：四、實(shí)驗(yàn)題目：5(s2+5s+6)求三階系統(tǒng)H(s)=3的單位階躍響應(yīng)，并繪制響應(yīng)波形圖。六、試驗(yàn)程序：num=[5 25 30]。step(num,den,10)。了解到了N階系統(tǒng)的單位階躍響應(yīng)的計(jì)算方法，和系統(tǒng)的響應(yīng)波形圖的函數(shù)應(yīng)用和繪制方法。九、實(shí)驗(yàn)改進(jìn)想法：無。五、解題分析：離散卷積是數(shù)字信號處理中的一個(gè)基本運(yùn)算，MTLAB提供的計(jì)算兩個(gè)離散序列卷積的函數(shù)是conv，其調(diào)用方式為 y=conv(x,h)。MATLAB函數(shù)conv的返回值y中只有卷積的結(jié)果，沒有y的取值范圍。許多離散LTI都可用如下的線性常系數(shù)的差分方程描述229。bx[kn]nnn=0n=0NN其中x[k]、y[k]分別系統(tǒng)的輸入和輸出。(an/a0)y[kn]+229。filter函數(shù)調(diào)用形式為 y=filter(b,a,x)。X表示輸入序列，y表示輸出序列。當(dāng)序列的DTFT可寫成ejW的有理多項(xiàng)式時(shí)，可用MATLAB信號處理工具箱提供的freqz函數(shù)計(jì)算DTFT的抽樣值。若X（ejW）可表示為b0+b1ejW+...bMejWMB(ejW)X(e)==jWjWjWNA(e)a0+a1e+...+aNe則freqz的調(diào)用形式為 X=freqz(b,a,w)，其中的b和 a分別是表示前一個(gè)jW式子中分子多項(xiàng)式和分母多項(xiàng)式系數(shù)的向量，即a=[a0,a1,...,aN]，w為抽樣的頻率點(diǎn)，向量w的長度至少為2。在抽樣點(diǎn)w上的值。六、實(shí)驗(yàn)程序:clc。close。a=[1 ]。subplot(2,1,1)。dstep(b,a,50)w=[0:1:512]*pi/512。freqz(b,a,w)七、實(shí)驗(yàn)結(jié)果：沖擊響應(yīng)圖及階躍響應(yīng)圖:Impulse (sec)Step (sec)3035404550 100Magnitude(dB) Frequency(180。p rad/sample)八、實(shí)驗(yàn)心得體會：通過實(shí)驗(yàn)我們知道了使用Matlab來繪出出一個(gè)線性移不變系統(tǒng)的幅頻和相頻曲線。九、實(shí)驗(yàn)改進(jìn)想法：無。五、解題分析：在MATLAB信號處理工具箱中，MATLAB提供了4個(gè)內(nèi)部函數(shù)用于計(jì)算DFT和IDFT，它們分別是：fft(x)，fft(x,N)，ifft(X)，ifft(X,N)。M是序列x的長度，即M=length(x)。若MN，則將原序列截短為N點(diǎn)序列，再計(jì)算其N點(diǎn)DFT；若Mifft(X)計(jì)算M點(diǎn)的IDFT。ifft(X,N)計(jì)算N點(diǎn)IDFT。clear。N=64。t=d*n。x=2*sin(4*pi*t)+5*cos(8*pi*t)。subplot(3,1,1)。title(‘source signal’)。plot(q,abs(y))。subplot(3,1,3)。title(‘phase’)。通過本次實(shí)驗(yàn)，對幅值譜和相位譜有了更深的了解，并與課程《信號與系統(tǒng)》里的一些相關(guān)知識連接到了一起，使得學(xué)到的只是更加深刻、有意義。實(shí)驗(yàn)四一、實(shí)驗(yàn)者姓名：二、實(shí)驗(yàn)時(shí)間：三、實(shí)驗(yàn)地點(diǎn)：四、實(shí)驗(yàn)題目：將信號x(t)=sin(240pt)做離散序列，比較原序列與經(jīng)過FFT和IFFT變換后的序列，并做出說明。連續(xù)時(shí)間函數(shù)與離散時(shí)間函數(shù)在編程中的區(qū)別主要體現(xiàn)在如下兩個(gè)方面：第一，自變量的取值范圍不同，離散時(shí)間函數(shù)的自變量是整數(shù)，而連續(xù)時(shí)間函數(shù)的自變量為一定范圍內(nèi)的實(shí)數(shù)；第二，繪圖所用的函數(shù)不同，連續(xù)函數(shù)圖形的繪制不止一個(gè)。用MATLAB表示一離散序列，可用兩個(gè)向量來表示。之后畫出序列波形。由于計(jì)算機(jī)內(nèi)寸的限制，MATLAB無法表示一個(gè)無窮長的序列。六、實(shí)驗(yàn)程序：t=0:1/255:1。y=real(ifft(fft(x)))。plot(t,x)。subplot(2,1,2)。七、實(shí)驗(yàn)結(jié)果：八、實(shí)驗(yàn)心得體會：通過對做信號的離散序列以及經(jīng)FFT和IFFT的變換，了解了相關(guān)特性。九、實(shí)驗(yàn)改進(jìn)想法：無。已知帶通濾波器的系統(tǒng)函數(shù)為H(s)=五、解題分析：需要知道求頻率響應(yīng)的方法，并繪制圖形。t=linspace(0,2*pi,1001)。U=[,1,]。a=[1,2,10001]。H=polyval(b,j*w)./polyval(a,j*w)。subplot(2,1,1),plot(w,abs(H)),grid。u21=abs(U(1)*H(1))*cos(99*t+angle(U(1)*H(1)))。u23=abs(U(3)*H(3))*cos(101*t+angle(U(3)*H(3)))。figure(2)。subplot(2,1,2),plot(t,u2)。九、實(shí)驗(yàn)改進(jìn)想