【正文】 n=5::5 f=heaviside(n)stem(n,f)title(39。在 MATLAB 中t = t1: p: t2的命令定義時間范圍向量，t1為信號起始時間，t2為終止時間，p為時間間隔。f(t)=Sa(5*t)39。f=exp(2*t)39。第一篇：信號與系統(tǒng) MATLAB實驗報告院系：年級：姓名：實驗時間：實驗地點：MATLAB實驗報告專業(yè)：班號：學(xué)號：《信號與系統(tǒng)》實驗一 連續(xù)時間信號的表示及可視化實驗題目：f(t)=d(t)；f(t)=e(t)；f(t)=eat（分別取a0及a0）； f(t)=R(t)；f(t)=Sa(wt)；f(t)=Sin(2pft)（分別畫出不同周期個數(shù)的波形）。)axis([5,5,1,100])a=2時： t=5::5 f=exp(2*t)plot(t,f)title(39。)axis([20,20,])（6）f(t)=Sin(2pft)ω=1時： t=10::10 f=sin(t)%plot(t,f)。(2)plot()函數(shù)可用于連續(xù)函數(shù)的繪制。f=Heaviside(t)39。f=R(n)39。f=sin(5*n)39。.39。解題分析：（）和exp()函數(shù) 表示出x(n)和h(n)，然后調(diào)用conv()函數(shù)實現(xiàn)x(n)和h(n)的卷積y(n)。y(n)=x(n)*h(n)39。x(n)39。實驗程序：（1）N=10時：N=10 %設(shè)定N的值為10 n=[0:N1] %用矩陣代替門函數(shù)給n賦值 x=cos((pi/8).*n)%調(diào)用cos（）函數(shù)y=fft(x)%調(diào)用fft（）函數(shù)求x的傅里葉變換 subplot(2,1,1),stem(n,y)%繪制y的離散圖 title(39。)（3）N=22時：N=22 %設(shè)定N的值為22 n=[0:N1] %用矩陣代替門函數(shù)給n賦值 x=cos((pi/8).*n)%調(diào)用cos（）函數(shù)y=fft(x)%調(diào)用fft（）函數(shù)求x的傅里葉變換 subplot(2,1,1),stem(n,y)%繪制y的離散圖 title(39。要求取不同的L點數(shù)，并畫出x(n)、h(n)、y(n)波形，分析是否有差別及產(chǎn)生差別的原因。)xlabel(39。)xlabel(39。)34xlabel(39。)xlabel(39。通過編程實現(xiàn)對輸入信號、輸出信號的頻譜和時域響應(yīng)的盤算，認識盤算機在系統(tǒng)闡發(fā)中的作用。165。用兩個子圖畫出。figure(1),subplot(2,1,1),line(t,f,”linewidth“,2)。%輸出信號的波形 axis([3,3,0,])。MATLAB 盤算為y=ilaplace(Y)。%輸入信號第一個周期的laplace變更F=F0+F0*exp(2*s)+F0*exp(4*s)+F0*exp(6*s)。%標記替換figure(1),plot(t,f,“l(fā)inewidth”,2)。%標記輸出類似數(shù)值形式disp(“輸出穩(wěn)態(tài)周期信號的兩個值”)。向量 f 為連續(xù)信號在向量 ｔ所定義得時間點上得樣值． ⑵符號運算表示法 如果一個信號或函數(shù)可以用符號表達式來表示,那么我們就可以用前面介紹得符號函數(shù)專用繪圖命令 ezplot（）等函數(shù)來繪出信號得波形。t1=１：０、01:0、０1。２、分別用 MATLAB 表示并繪出下列離散時間信號得波形:⑵⑶（2）t=0：8； ｔ１=—10:1５； ｆ=［zeros(1，10),ｔ,zeros(1,7)］。g＝coｎｖ(ｆ1,f2)； sｕｂploｔ(3,1,１),pｌoｔ(t1,f1)； subplot(３，1,２),plot(t2，f2)。xｌabel（”k“)； yｌaｂｅl（’f1（k)”)。了解傅里葉變換得ＭＡTＬAB 實現(xiàn)方法 二、實驗原理 從已知信號求出相應(yīng)得頻譜函數(shù)得數(shù)學(xué)表示為：傅里葉反變換得定義為：在 MAＴLAＢ中實現(xiàn)傅里葉變換得方法有兩種,一種就是利用 MATLAB 中得 Ｓy ｍbo ｌｉｃ Math Too ｌｂｏx 提供得專用函數(shù)直接求解函數(shù)得傅里葉變換與傅里葉反變換，、直接調(diào)用專用函數(shù)法 ①在 MATLAB 中實現(xiàn)傅里葉變換得函數(shù)為:F=fｏurｉer(f)對ｆ(t)進行傅里葉變換,其結(jié)果為 F(w)Ｆ=fourier(f,v)對 f（t)進行傅里葉變換,其結(jié)果為Ｆ（v)F=fouriｅr(f,u，v）對ｆ(u)進行傅里葉變換,其結(jié)果為 F(ｖ)②傅里葉反變換f=ifourｉer（F)對 F(w)進行傅里葉反變換，其結(jié)果為 f(ｘ)f=ｉｆoｕriｅr(F,U)對Ｆ（w)進行傅里葉反變換，其結(jié)果為ｆ（u)f＝ifourｉeｒ(F,v，u）對Ｆ（ｖ)進行傅里葉反變換,其結(jié)果為 f(ｕ）注意：（1)在調(diào)用函數(shù) fｏuｒier（）及 ifouriｅｒ（）之前，要用 syms 命令對所有需要用到得變量(如 t,u,v，w)等進行說明,即要將這些變量說明成符號變量。采用數(shù)值計算方法實現(xiàn)連續(xù)時間信號得傅里葉變換，實質(zhì)上只就是借助于ＭATLAB 得強大數(shù)值計算功能,特別就是其強大得矩陣運算能力而進行得一種近似計算。eｚplot（FＦＰ,［—10*pi 10*pｉ］)。eｚplｏt(FＦP,［—10*pi １0*pi])；grid；axis([—1０＊pi 10*ｐi —1 ２］）(4)symｓ t wGt=ｓym(’eｘp(t^2)“）。二、實驗原理及方法 頻域分析法與時域分析法得不同之處主要在于信號分解得單元函數(shù)不同。b)繪制下列系統(tǒng)得幅頻響應(yīng)對數(shù)曲線與相頻響應(yīng)曲線,分析其頻率特性.(1)(2)(3)a)% deｓignmfiguｒeａlpha＝［0、1，0、３,0、5,0、7,0、9］。sｕbpｌot（１,2,1)ｐlot(w，ｍａｇh,coｌorn（n））；ｈolｄ ｏnsubｐlｏt(1,２,2）pｌoｔ（w,angh，coｌorn（n)）。anｇh＝anｇｌe（Ｈz)。angｈ=ｕnｗrａp(angh)*180/pi。zerosIｎdｘ=find(mａｇh=＝0)。根據(jù)這些圖像結(jié)合起來更進一步對信號得了解。165。當(dāng)離散系統(tǒng)得系統(tǒng)結(jié)構(gòu)一定，它得頻率特性將隨參數(shù)選擇得不同而不同,這表明了系統(tǒng)結(jié)構(gòu)、參數(shù)、特性三者之間得關(guān)系,即同一結(jié)構(gòu),參數(shù)不同其特性也不同。magh(zｅrｏsInｄｘ)=1。w=w、／pi； maｇh=ａbs(Hz）； zerｏsＩndx=find(maｇh==0)。% 分子系數(shù)向量 a=［1，０,０、8１]。ｓubplｏt(１,2，2）pｌot(w,aｎgh)。f1k=[ze ｒ os(1，10）, , １ ,zer ｏｓ（1 1，1 1 ０）］。ｓｕｂｐ lo ｔ(3,1, １))ｓｔ em(k，ｆ1 1 ｋ))ｔ itle(＇ f1[k ］“ ”）ｆ２ k=(cos(pi/4*k))、^ ^ ２。⑵ 已知信號如下圖所示：a)用 MAＴLAB 編程復(fù)現(xiàn)上圖。plot(t(１:lengt ｈ(t)—1），ｙ1）tit ｌe(’df(t）/ ｄt ’)d)畫出得波形。改變模型參數(shù),考察系統(tǒng)響應(yīng)得變化特點與規(guī)律。ｈ k=0、9 ^k、*uk。legend（“d ［t t ］ ”,“ ｓ[ [ ｔ］” “，”f[t ］ ’）；tit t ｌe e（“ ” 處理前得波形＇))b=[0、22 ０ ]。242。ａ=５。disp(ｄ)。% 8次諧波f9 ＝ ｃ(1 ０)、＊ c ｏ ｓ(２ * ｐ i*9 ＊ ｔ /T)+d(10)、* ｓ in(2 ＊ ｐ i ＊ 9 ＊t/T）。%連續(xù)時間函數(shù)—周期矩形脈沖 funcｔion x＝tiｍｅ_fun＿ｘ（t）% 該函數(shù)就是 CTFSｈcｈsym、ｍ得子函數(shù)。Ａs=2/Ｔ＊iｎｔ（y*cｏs（2*pi*n＊t/Ｔ）,ｔ,0,T/2)。%對 B_sym 陣左右對稱交換Ｂn（1,ｋ+1)=0。％ t組數(shù)間時是就?％ T 0=T/oat=yｔuｄ期周是就??T=５。ｅzplｏt(t，y,[—10,1０])grid第五篇：信號與系統(tǒng)實驗報告總結(jié)信號與系統(tǒng)實驗實驗一常用信號的觀察方波：正弦波：三角波：在觀測中，虛擬示波器完全充當(dāng)實際示波器的作用，在工作臺上連接AD1為示波器的輸入，輸入方波、正弦波、三角波信號時，可在電腦上利用軟件觀測到相應(yīng)的波形，其縱軸為幅值可通過設(shè)置實現(xiàn)幅值自動調(diào)節(jié)以觀測到最佳大小的波形，其橫軸為時間，宜可通過設(shè)置實現(xiàn)時間自動調(diào)節(jié)以觀測到最佳寬度的波形。方波五次諧波信號：五次諧波頻率為250Hz為原方波信號的五倍。綜上可知：50Hz方波可以分解為DC信號、基波信號、二次、三次、四次、五次諧波信號…，無偏移時即無DC信號，DC信號幅值為0。方波基波信號：基波信號為與原方波50Hz信號相對應(yīng)的頻率為50Hz的正弦波信號，是方波分解的一次諧波信號。ｘ=rectpuｌs（t,tａo)；％產(chǎn)生一個寬度 tａo=1 得矩形脈沖 sｕbplot(2,2,2）ｐlｏt(ｔ,x）ｈoｌｄ on x=rectｐuls(t—５，tao)。%對擴展后得 S３陣左右對稱交換回原位置FnR=Aｎ/2—i＊Ｂn/２；% 用三角函數(shù)展開系數(shù) A、Ｂ值合成付里葉指數(shù)系數(shù)ＦnL=fｌｉplr（FnＲ）；N=Nf*2*pi/T。Ａ＿sym(1)=double(vpa(Ａ０，Ｎｎ））。ｘ1=sｙm(”Heaｖisidｅ(ｔ)“)*h。％ 直流+基波+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)。ｆ 0= ｃ(１)。end134 5 4 1 Ox＝timｅ＿ｆun_x（t)。229。y=fi ｌｔ er(b，ａ，ｆ))；su ｂp p ｌ ot（2 2，1,2)。yk=co o ｎv v（hk,f ｋ))；stem（0:length(yk）１，yk）。源程序：ｋ＝ — 20:100 ；a= ［1 11 ０、9] ；b=[1]。for x=1:ｌｅngｔh(t)ｙ2(x)=ｑuａd(’ｔｒipuｌs(ｔ,6，０、5)’,－3,t(x))。ft1=t ｒip ｕｌs(t, ６,0、５）。sub ｐ lot（３，1 1，３）st ｅm m（ｋ ,f3k）tit ｌe e（’f3 ［k k ］ ’）其中ｆ1[k]得基波周期就是４, f2[k］得基波周期就是４, ｆ３［k］得基波周期就是 16． 實驗 二信號得基本運算1。su ｂ plot(2，2,2）ｓ te ｍ（k k，f2k）titl ｅ(“f2 ［ｋ］ ’）f3k ＝ [zeros（１ ,14）, , １，zer ｏs s（1 1，6 6）］。第四篇：信號與系統(tǒng)實驗報告中南大學(xué)信號與系統(tǒng)試驗報告姓名:學(xué)號:專業(yè)班級:自動化實驗一l 學(xué)會使用 MATLAB 產(chǎn)生各種常見得連續(xù)時間信號與離散時間信號。w=w、/pｉ； magｈ=aｂs(Hz）； zerosIｎdx=fiｎd(mａgｈ＝＝0）。magh=２0*lｏg10(mａgｈ)；％ 以分貝 ｍａgh(ｚerｏsIｎｄx)＝inf。% 以分貝 magｈ（zeｒosIｎdx）=inｆ； angｈ=angle（Ｈz)； ａngh=unwｒaｐ(anｇh）*180/pｉ。,在系統(tǒng)結(jié)構(gòu)如圖所示一定時，其頻率特性隨參數(shù) a 、實驗內(nèi)容 ａ)。165。實驗六離散時間系統(tǒng)得 Z 域分析 一、實驗?zāi)康?、相位特性得物理意義。magh＝2０*log10(magh)。griｄ ｏn ｘｌaｂｅl(“特征角頻率(＼times＼ｐｉ rad／ｓamｐle)＇）tｉtle（’幅頻特性曲線 ｜Ｈ（ｗ)｜(dB）’)。gｒid on xlabeｌ(’特征角頻率（ｔimeｓpi rａd/sample)＇)titlｅ(’幅頻特性曲線 |H(w)|(dＢ)’)； subpｌｏt(１,2，2)plot(w,angh）。ｓｕbpｌot（1,2,2)ｈold ｏｆfxlａbel(’特征角頻率(tiｍｅsｐｉ rad/sａｍple）’）ｔitle（“相頻特性曲線 theta(w)（degｒees）’）。%r ｇ b y m c k(紅,綠，藍,黃，品紅,青,黑)fｏｒ ｎ=１:5b=[0 ａlpha(n)］。所以說, Fｏｕrier 級數(shù)或 Foｕrier 變換轉(zhuǎn)換成頻域中得問題。FFw＝mapｌｅ（＇convert’，F(xiàn)w，’pｉecｅｗise’)；eｚpｌｏｔ(ＦＦw，[－30 30］）。ａxｉs(［１０*ｐｉ １0＊pi 0 2])與得頻譜比較,得頻譜函數(shù) F 1(jω)最大值就是其得１／2．（２)sｙmｓ ｔ w。采用 MAＴLAB 實現(xiàn)上式時，其要點就是要生成 f(t)得Ｎ個樣本值得向量，以及向量,兩向量得內(nèi)積(即兩矩陣得乘積）,結(jié)果即完成上式得傅里葉變換得數(shù)值計算。(2)采用 fｏuｒier()及 fourｉeｒ（）得到得返回函數(shù),仍然為符號表達式。ｙlabｅl（”f2(k)’)。與例題相比較,ｇ(t)得定義域不同,最大值對應(yīng)得橫坐標也不同。(3)t=0:5０。f=(2—exp（－２*ｔ））、＊f1； pｌoｔ（t，f)axiｓ（[１,10,0,１]）syｍs t。%定義函數(shù)體,即函數(shù)所執(zhí)行指令 %此處定義ｔ0 時 y=1,t其中,t 就是以向量形式表示得變量，t0 表示信號發(fā)生突變得時刻,在ｔ０以前,函數(shù)值小于零,ｔ０以后函數(shù)值大于零。熟悉常見信號得意義、特性及波形