【正文】 0＊pi 0 2])與得頻譜比較,得頻譜函數(shù) F 1(jω)最大值就是其得１／2．（２)sｙmｓ ｔ w。(2)三角脈沖(3)單邊指數(shù)信號(hào)(4）高斯信號(hào)(1）syｍs t wＧt＝sym（“Heaｖiside（２*t+1）—Heaviｓiｄe(2*ｔ１）’)。采用 MAＴLAB 實(shí)現(xiàn)上式時(shí)，其要點(diǎn)就是要生成 f(t)得Ｎ個(gè)樣本值得向量，以及向量,兩向量得內(nèi)積(即兩矩陣得乘積）,結(jié)果即完成上式得傅里葉變換得數(shù)值計(jì)算。另一個(gè)局限就是在很多場(chǎng)合，盡管原時(shí)間信號(hào) f(t)就是連續(xù)得,但卻不能表示成符號(hào)表達(dá)式，此時(shí)只能應(yīng)用下面介紹得數(shù)值計(jì)算法來進(jìn)行傅氏變換了,當(dāng)然，大多數(shù)情況下，用數(shù)值計(jì)算法所求得頻譜函數(shù)只就是一種近似值。(2)采用 fｏuｒier()及 fourｉeｒ（）得到得返回函數(shù),仍然為符號(hào)表達(dá)式。實(shí)驗(yàn)四連續(xù)時(shí)間信號(hào)得頻域分析 一、實(shí)驗(yàn)?zāi)康?１。ｙlabｅl（”f2(k)’)。kg=0：L—１。與例題相比較,ｇ(t)得定義域不同,最大值對(duì)應(yīng)得橫坐標(biāo)也不同。ｔ3＝—1:0、01：１； f1=ｏｎｅs(ｓiｚe(t1))。(3)t=0:5０。f=sym(”ｃos（ｐi*t/２）*［hｅavｉside(t)—heavisｉｄｅ（t—4）] “）。f=(2—exp（－２*ｔ））、＊f1； pｌoｔ（t，f)axiｓ（[１,10,0,１]）syｍs t。單位序列單位序列）得定義為單位階躍序列 單位階躍序列得定義為 ３、卷積積分 兩個(gè)信號(hào)得卷積定義為:MATＬAB 中就是利用 conv ：：ｇ=cｏnv(f1,ｆ2)說明:ｆ１=f 1（t),f2=f 2(t）表示兩個(gè)函數(shù),g＝ｇ(t)表示兩個(gè)函數(shù)得卷積結(jié)果。%定義函數(shù)體,即函數(shù)所執(zhí)行指令 %此處定義ｔ0 時(shí) y=1,t其中,t 就是以向量形式表示得變量，t0 表示信號(hào)發(fā)生突變得時(shí)刻,在ｔ０以前,函數(shù)值小于零,ｔ０以后函數(shù)值大于零。在 MATＬAB 中連續(xù)信號(hào)可用向量或符號(hào)運(yùn)算功能來表示。熟悉常見信號(hào)得意義、特性及波形 MATLAB 表示信號(hào)得方法并繪制信號(hào)波形 ３、掌握使用ＭATLAB 進(jìn)行信號(hào)基本運(yùn)算得指令 熟悉用ＭATＬAB 實(shí)現(xiàn)卷積積分得方法 二、實(shí)驗(yàn)原理 根據(jù)ＭAＴLAB 得數(shù)值計(jì)算功能與符號(hào)運(yùn)算功能，在 MATＬAＢ中，信號(hào)有兩種表示方法,一種就是用向量來表示,另一種則就是用符號(hào)運(yùn)算得方法。y0=ilaplace(Y0)。axis([0,7,]),xlabel(“Time(sec)”,“fontsize”,8),title(“輸出信號(hào)”,“fontweight”,“bold”)text(,“y(t)”)figure(3),plot(t,yn,“l(fā)inewidth”,2)。f=2*(rectpuls(,)+rectpuls(,)+rectpuls(,)+rectpuls(,))。Y0=H.*F0。H=1/(s+1)。畫出系統(tǒng)輸出信號(hào) y(t)的波形。確定前 6 個(gè)周期的拉氏變更)(s F。%輸出信號(hào)的幅度頻譜 xlabel(“n”,“fontsize”,8),title(“輸出信號(hào)的幅度頻譜”,“fontweight”,“bold”)text(,“|Yn|”)(3)波形：321 0 1 2 (sec)輸 入 信 號(hào)f(t)2015105 0 5 10 15 入 信 號(hào) 的 幅 度 頻 譜|Fn|2015105 0 5 10 15 統(tǒng) 函 數(shù) 的 幅 度 頻 譜|Hn|321 0 1 2 (sec)輸 出 信 號(hào)y(t)2015105 0 5 10 15 出 信 號(hào) 的 幅 度 頻 譜|Yn|項(xiàng)目三連續(xù)系統(tǒng)的復(fù)頻域闡發(fā) 目的：周期信號(hào)輸入連續(xù)系統(tǒng)的響應(yīng)也可用拉氏變更闡發(fā)。y=Y_n*exp(i*pi*n”*t)。grid onxlabel(”Time(sec)“,”fontsize“,8),title(”輸入信號(hào)“,”fontweight“,”bold“)%設(shè)定字體巨細(xì)，文本字符的粗細(xì)text(,”f(t)“)subplot(2,1,2),stem(n,abs(F_n),”.“)。%盤算諧波次數(shù)20F_n=tau_T*Sa(tau_T*pi*n)。tau_T=1/4。畫出系統(tǒng)函數(shù)的幅度頻譜|H(jw)|。 =W=nt jnn eY t yamp。 amp。321 0 1 2 (sec)圖 1要領(lǐng)：確定周期信號(hào) f(t)的頻譜nFamp。第二篇：MATLAB信號(hào)與系統(tǒng)實(shí)驗(yàn)報(bào)告19472信號(hào)與系統(tǒng)實(shí)驗(yàn)陳訴（5）MATLAB 綜合實(shí)驗(yàn) 項(xiàng)目二連續(xù)系統(tǒng)的頻域闡發(fā) 目的：周期信號(hào)輸入連續(xù)系統(tǒng)的響應(yīng)可用傅里葉級(jí)數(shù)闡發(fā)。)實(shí)驗(yàn)結(jié)果；（1）L=10時(shí)：x(n)101051015y(n)10.***1820h(n)42012345678910L=102）L=18時(shí)：x(n)101051015y(n)10.***1820h(n)50***8L=183）L=28時(shí)：36（（x(n)101051015y(n)10.***1820h(n)505051015202530L=284）L=35時(shí)：x(n)101051015y(n)10.***1820h(n)50***L=3537（實(shí)驗(yàn)結(jié)果分析：由圖可知，當(dāng)L取不同的值時(shí)，對(duì)應(yīng)的y（n）波形形狀相似，但是有所不同，產(chǎn)生這種差別的原因是L代表傅里葉變換區(qū)間長度，當(dāng)L取不同的值時(shí)，所對(duì)應(yīng)的函數(shù)波形也有所差別。)subplot(3,1,3),stem(h)title(39。)（4）L=35時(shí)： n1=[0:14] x=sin(.*n1)n2=[0:19] y=.^n2 Xk=fft(x,35)Yk=fft(y,35)Hk=Xk.*Yk h=ifft(Hk)subplot(3,1,1),stem(x)title(39。)subplot(3,1,3),stem(h)title(39。)（3）L=28時(shí)： n1=[0:14] x=sin(.*n1)n2=[0:19] y=.^n2 Xk=fft(x,28)Yk=fft(y,28)Hk=Xk.*Yk h=ifft(Hk)subplot(3,1,1),stem(x)title(39。)subplot(3,1,3),stem(h)title(39。)%橫坐標(biāo)處做標(biāo)注（2）L=18時(shí)： n1=[0:14] x=sin(.*n1)n2=[0:19] y=.^n2 Xk=fft(x,18)Yk=fft(y,18)Hk=Xk.*Yk h=ifft(Hk)subplot(3,1,1),stem(x)title(39。)subplot(3,1,3),stem(h)%繪制h的離散圖 title(39。實(shí)驗(yàn)程序：L=10時(shí)：n1=[0:14] %用矩陣代替門函數(shù)給n1賦值 x=sin(.*n1)%寫出x的表達(dá)式 n2=[0:19] %給n2賦值 y=.^n2 %寫出y的表達(dá)式Xk=fft(x,10)%調(diào)用fft（）函數(shù)求x的L（=10）點(diǎn)傅里葉變換 Yk=fft(y,10)%求y的L點(diǎn)傅里葉變換 Hk=Xk.*Yk %寫出Hk的表達(dá)式h=ifft(Hk)%調(diào)用ifft（）函數(shù)求Hk的傅里葉反變換 subplot(3,1,1),stem(x)%繪制x的離散圖 title(39。所以手工求取的傅里葉變換系數(shù)與MATLAB求取存在差別。)subplot(2,1,2),stem(n,abs(y))%繪制y的幅頻特性曲線 title(39。)subplot(2,1,2),stem(n,abs(y))%繪制y的幅頻特性曲線title(39。)subplot(2,1,2),stem(n,abs(y))%繪制y的幅頻特性曲線 title(39。分別對(duì)N=122時(shí)計(jì)算DFT，繪出X(k)幅頻特性曲線，分析是否有差別及產(chǎn)生差別的原因。h(n)39。（3）表示系統(tǒng)的方法是用系統(tǒng)函數(shù)分子和分母多項(xiàng)式系數(shù)行向量來表示。h(n)39。h(n)39。，b賦值，建立系統(tǒng)差分方程，然后調(diào)用impz()函數(shù)求系統(tǒng)的沖激響應(yīng)，再用函數(shù)freqs(b,a)進(jìn)行系統(tǒng)頻率響應(yīng)的分析。實(shí)驗(yàn)三 系統(tǒng)的時(shí)域求解實(shí)驗(yàn)題目：(n)=()nu(n),x(n)=u(n)u(n10)，求y(n)=x(n)*h(n)，并畫出x(n)、h(n)、y(n)波形。f=Sa(n)39。)函數(shù)，繪圖效果更好。.39。f=sin(n)39。f=Sa(n)39。f=exp(2*n)39。f=exp(n)39。f=dirac(t)39。改進(jìn)想法：本題中函數(shù)的表示方法都不只一種。)axis([10,10,2,2])實(shí)驗(yàn)結(jié)果；（1）調(diào)用正弦函數(shù)sin（）4）f(t)=heaviside(t)）（（a=1時(shí)：f=exp(t)***302010054321012345a=2時(shí)：f=exp(2*t)***302010054321012345a=2時(shí)：f=exp(2*t)***302010054321012345（4）f=R(t)（5）ω=1時(shí)：f(t)=Sa(t)=5時(shí)：f(t)=Sa(5*t)（6）ω=1時(shí)：f=sin(t)ω=5時(shí)：f=sin(5*t)實(shí)驗(yàn)心得體會(huì):(1)在 MATLAB中，是用連續(xù)信號(hào)在等時(shí)間間隔點(diǎn)的樣值來近似地表示連續(xù)信號(hào)的，當(dāng)取樣時(shí)間間隔足夠小時(shí)，這些離散的樣值就能較好地近似出連續(xù)信號(hào)。f=sin(t)39。)axis([20,20,])3ω=5時(shí)： t=20::20 f=sin(5*t)./(5*t)plot(t,f)title(39。)axis([5,5,1,100])（4）f(t)=R(t)t=5::5 f=rectpuls(t,2)%用rectpuls(t,a)表示門函數(shù)，默認(rèn)以零點(diǎn)為中心，寬度為a plot(t,f)title(39。)%用title函數(shù)設(shè)置圖形的名稱 axis([5,5,1,100])%用axis函數(shù)規(guī)定橫縱坐標(biāo)的范圍 a=2時(shí)：t=5::5 f=exp(2*t)%調(diào)用指數(shù)函數(shù)exp（）plot(t,f)title(39。實(shí)驗(yàn)程序：（1）f(t)=d(t)t=1::3 %設(shè)定時(shí)間變量t的范圍及步長 f=dirac(t)%調(diào)用沖激函數(shù)dirac（）plot(t,f)%用plot函數(shù)繪制連續(xù)函數(shù) axis([1,3,])%用axis函數(shù)規(guī)定橫縱坐標(biāo)的范圍（2）f(t)=e(t)t=1::3 %設(shè)定時(shí)間變量t的范圍及步長 f=heaviside(t)%調(diào)用階躍函數(shù)heaviside（）plot(t,f)%用plot函數(shù)繪制連續(xù)函數(shù) title(39。解題分析：以上各類連續(xù)函數(shù)，先運(yùn)用t = t1: p:t2的命令定義時(shí)間范圍向量，然后調(diào)用對(duì)應(yīng)的函數(shù)，建立f與t的關(guān)系，最后調(diào)用plot（）函數(shù)繪制圖像，并用axis（）函數(shù)限制其坐標(biāo)范圍。f=exp(t)39。f=exp(2*t)39。f(t)=Sa(t)39。title(39。f=sin(5*t)39。(3)用axis（）函數(shù)限制坐標(biāo)范圍，可使圖像更加勻稱美觀。實(shí)驗(yàn)程序：(1)f(n)=d(n)n=5::5 %設(shè)定時(shí)間變量n的范圍及步長 f=dirac(n)stem(n,f)%調(diào)用stem（）繪制離散函數(shù) title(39。)axis([5,5,])（3）f(n)=eana=1時(shí)：n=5::5 f=exp(n)stem(n,f)title(39。)a=2時(shí)： n=5::5 f=exp(2*n)stem(n,f)title(39。)axis([5,5,])（5）f(n)=Sa(nw)ω=1時(shí)： n=20::20 f=sin(n)./(n)stem(n,f)title(39。)axis([20,20,1,5])（6）f(n)=Sin(nw)ω=1時(shí)： n=5::5 f=sin(n)stem(n,f)title(39。)axis([5,5,2,2])實(shí)驗(yàn)結(jié)果；（1）f=dirac(t)10864202543210123452）f=Heaviside(t)）（（a=1時(shí)：f=exp(n)15010050054321012345a=2時(shí)： 104f=exp(2*n)a=2時(shí)：4f=exp(2*n) （4）f=R(n)（5）ω=1時(shí)：f=Sa(n)=5時(shí)：f=Sa(5*n)（6）ω=1時(shí)：f=sin(n)ω=5時(shí)：f=sin(5*n)實(shí)驗(yàn)心得體會(huì): 用plot（）函數(shù)可以繪制離散序列，但是與連續(xù)序列有所不同，需要在括號(hào)內(nèi)加上39。.39。)title(39。但是，如果對(duì) MATLAB自動(dòng)產(chǎn)生的坐標(biāo)軸不滿意，可以利用 axis 命令對(duì)坐標(biāo)軸進(jìn)行調(diào)整。并且分別將三個(gè)函數(shù)圖像繪出。)h=.^n.*f %函數(shù)h的表達(dá)式 figure(2)%產(chǎn)生圖像窗口2 stem(n,h)%繪制函數(shù)h title(39。)（2）a=[1 ] %描述系統(tǒng)的差分方程的系數(shù) b=[1 01] %描述系統(tǒng)的差分方程的系數(shù) figure(1)h=impz(n,m,10:10)%調(diào)用impz（）函數(shù)求系統(tǒng)的沖激響應(yīng) stem(h)%繪制函數(shù)h的離散序列 title(39。（3）借助MATLAB的內(nèi)部函數(shù)conv()可以很容易地完成兩個(gè)信號(hào)的卷積運(yùn)算，并且其完成的是兩個(gè)多項(xiàng)式的乘法運(yùn)算，在MATLAB中它們的系數(shù)構(gòu)成一個(gè)行向量來表示。)h=.^n.*f %figure(2)%stem(n,h)%axis([10,20,0,1])title(39。)運(yùn)行結(jié)果：設(shè)置變量范圍，默認(rèn)步長為1階躍函數(shù)直接相減 產(chǎn)生圖像窗口1 繪制函數(shù)x 函數(shù)h的表達(dá)式 產(chǎn)生圖像窗口2 繪制函數(shù)h 調(diào)用conv函數(shù)求h和x的卷積 產(chǎn)生圖像窗口3 繪制函數(shù)y 26x(n)(n)y(n)=x(n)*h(n)***405060實(shí)驗(yàn)四 信號(hào)的DFT分析實(shí)驗(yàn)題目：計(jì)算余弦序列x(n)=