隨機(jī)信號(hào)分析實(shí)驗(yàn)報(bào)告(編輯修改稿)
2024-11-16 23:45 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 LAＢ提供了多個(gè)設(shè)計(jì)濾波器得函數(shù),可以很方便地設(shè)計(jì)低通、帶通、高通、多帶通、帶阻濾波器。((二）實(shí)驗(yàn) 目得系統(tǒng)仿真就是信號(hào)仿真處理得一個(gè)重要部分，通過(guò)該實(shí)驗(yàn)要求學(xué)生掌握系統(tǒng)仿真得基本概念,并學(xué)會(huì)系統(tǒng)得仿真方法。((三))實(shí)驗(yàn) 結(jié)果低通濾波器帶通濾波器高通濾波器 多帶通濾波器帶阻濾波器附：源程序 １、X（n)N＝2000；fｓ=４00。Nn=ｒandom（”normａl＇,0,1,1，N)； t=（0：Ｎ—１)/fs。fi=ｒaｎｄｏm(’unｉf’,0,１,１,２)＊2＊pｉ。xn＝sin(2*ｐi*５０*ｔ+ｆｉ(1))+Nn。Ｒx=ｘcorｒ(xn，“biaｓed’)； ｍ=—N+1:N－1。Sx=abs（ｆft（xn）、^2）/Ｎ； f=(—N／2:N/2１）＊fs/N。sｕbｐlｏt（211),ｐｌot(m，Rｘ)； xlabｅl(’m’）ｙｌabel（”Rx(m）’)titlｅ（’xn 得自相關(guān)函數(shù)“)。subploｔ(212）,plｏt（f,ffｔshift（10＊lｏg10（Sx(1:N）)）)。xlabeｌ（’f/Hz”)ｙlabｅl(“Sx/dB”)tｉtle(’xn 得功率譜密度’）。低通濾波器 ｈ=ｆir1（１00,0、４)。H=ｆft(h,2＊N)。HW=aｂs(Ｈ)、＾２。Rx=xcoｒr（xn，’ｂｉasｅｄ＇)。Sx=abs（ｆfｔｓhift（fft(xn，2*N））、＾2)/（2＊Ｎ)； Ｓy＝Sx、*ＨＷ。Ｒy＝fftｓｈift（ifft（Sｙ)）。f=（Ｎ：N—1)＊ｆs/(2＊N)； m=(—N:Ｎ－1)。subｐｌot(311)。pｌot（(N:Ｎ—1)/N，fｆtshiｆt(aｂs(HW(1:2＊N))）)。tiｔlｅ(＇低通濾波器“）。ｓubplｏt(31２)，ｐｌoｔ（ｍ,Ry)。xlabｅl(”m“)ylａbeｌ（”Ry（ｍ)＇)title(’xn 經(jīng)低通濾波器得自相關(guān)函數(shù)’)； ｓubｐloｔ(31３),plot（ｆ，ffｔshift(10*log10(Ｓy（１:2*N)）))； ａxis(［—200 20０ —20 20])。xlabeｌ(“f/Hz’)ｙlabｅl(＇Sｙ/ｄＢ”)titｌe(＇xn 經(jīng)低通濾波器得功率譜密度“)； 帶通濾波器 ｈ=fｉr1(10０，［０、１ ０、5])。H=ｆft(h，2*Ｎ)。ＨW＝abs(H)、^2； Ｒx=ｘcorｒ(xn，”biａsed“）； Sｘ=aｂs（fｆtshiｆt(fｆt(ｘn,2*N））、＾2)／（2*N）； Sy=Ｓx、*ＨW； Rｙ=fｆtsｈift(ｉfft（Sｙ)）； f=(Ｎ:N1）*fs/(2＊Ｎ)。m=（N:N—1）。subｐlot(311)。ｐlot((—N:N1）/N,fftｓｈift（abs（HW（1:2＊Ｎ)))）； ｔitle(’帶通濾波器”）； ｓubｐlｏｔ（312)，ploｔ（ｍ,Ry)。xlabeｌ(’m“)yｌabel(’Ｒy(m)’)title(”xn 經(jīng)帶通通濾波器得自相關(guān)函數(shù)“)； ｓubpｌoｔ(31３)，plｏt(f,fftshift(10＊lｏg10（Sy(１：2＊N))））； ａxiｓ（［—２０0 20０ －２0 20])。xｌabｅl（’f/Hz”)yｌabｅl(“Ｓｙ/dＢ’)tｉtle（’xn 經(jīng)帶通濾波器得功率譜密度’）。高通濾波器 h=ｆｉr1（１0０,0、6,’high’）； Ｈ=fft(h,2＊N）； ＨW=abｓ（H)、^２。Ｒx=xcｏｒr(xｎ,”bｉａｓed“）。Ｓx=abs(fftｓｈift（fｆt(xn,2＊N))、^２)/（2*N)； Sｙ=Ｓx、*HW。Ｒy=ffｔshiｆt(iｆfｔ（Sｙ))。f=(Ｎ:N1）*ｆs/（2*N)。m=(—Ｎ:N—1）；sｕｂｐlｏt(31１)；plot(（N：Ｎ—1）/N,ｆftshiｆt(aｂs(HＷ（1:２*N)）))。tiｔｌe(＇高通濾波器”）。subpｌot(312），ploｔ(m,Rｙ)。xlabｅl(“ｍ’)ylaｂeｌ(’Rｙ(m)”)tiｔlｅ(＇xn 經(jīng)高通通濾波器得自相關(guān)函數(shù)’)。subplｏt（3１3）,ｐｌot(ｆ,fftshｉft(10*ｌog10(Sy（1:2＊Ｎ)）））。axis（［－２０0 2０0 —2０ 2０]）； xlａbeｌ（“ｆ/Ｈz’)ｙｌａbｅl（”Sy/ｄB“)tｉtle(＇xn 經(jīng)高通濾波器得功率譜密度＇)。多帶通濾波器 h=fｉr１(100,[０、１,0、3，0、5,0、7］)； H＝ｆfｔ(h，2*N)； HW＝ａｂｓ（Ｈ)、^2； Ｒx=xcorr（xn,＇bｉasｅd’)。Sx=abs（fｆtshift(fft(xn,2*N））、^2)/（２＊Ｎ）； Sy＝Ｓx、*HW。Ｒy=ffｔshiｆt（iｆft(Sy))。f＝（—Ｎ:N—1)＊fｓ/(2*N)。m＝（—N：Ｎ1)。sｕbｐloｔ(３１1）。ploｔ(（—N:N—１）/Ｎ,fftｓｈift（ａbs(HW（1:２＊N）)))； tiｔle(’多帶通濾波器’)； ｓｕbplot(312),pｌｏt(m,Ry）； ｘlabel（＇m’)ｙｌabｅl（”Ry(m）“）ｔitle（”xn 經(jīng)多帶通通濾波器得自相關(guān)函數(shù)“)。subpｌot（31３),plot(ｆ,fｆｔshｉｆt（1０＊ｌog10(Sｙ（1：2＊Ｎ)）))。axis（［－2０0 2０0 —20 20])。xlaｂeｌ（’ｆ／Hz”）ｙlabｅl(“Sy/dB’）ｔitle(’xn 經(jīng)多帶通濾波器得功率譜密度”）； 帶阻濾波器 h=fiｒ１(100,[０、1,0、4],’sｔoｐ’)。H=fft（h，2＊N)。HＷ=abｓ(H)、^2。Rx=ｘｃｏrr（xｎ,’ｂiaseｄ“）。Ｓx＝ａbｓ（ffｔshift(fft（ｘｎ，２*N）)、＾2)／(2*N)。Sy=Sx、＊ＨW； Ｒy＝fftshift（ｉfft(Sy））。f=(—N:N1）*ｆｓ/（2*N)。m=(N：N—1)； ｓｕｂｐlot(311）；pｌoｔ(（—N：N1）/Ｎ,fftshift(ａｂs（ＨＷ(1:２*Ｎ))))。ｔｉtｌｅ(”帶阻濾波器“）； ｓubpｌot（31２),ｐlot（m，Ｒy)； xlａbeｌ(’m’）ｙlａbel（”Ｒy(m)’)ｔitle（’ｘn 經(jīng)帶阻濾波器得自相關(guān)函數(shù)＇)； sｕbplｏｔ(3１3)，ｐlot(f,fftshift（10＊log10(Sy（１:2*N)）））。aｘis([－200 20０20 20］)。xlaｂｅl（＇f／Ｈz“)ylabｅl（”Sy/ｄB“)titlｅ（”xn 經(jīng)帶阻濾波器得功率譜密度）。第三篇：隨機(jī)信號(hào)分析基礎(chǔ)讀書(shū)報(bào)告讀書(shū)報(bào)告——隨機(jī)信號(hào)分析基礎(chǔ)本讀書(shū)報(bào)告主要分為三部分：一、自學(xué)計(jì)劃。二、理論原理知識(shí)。三、個(gè)人總結(jié)及心得體會(huì)。一、自學(xué)計(jì)劃。在研究生第一學(xué)期，開(kāi)設(shè)了隨機(jī)信號(hào)分析基礎(chǔ)課，這門(mén)課程是在信號(hào)分析基礎(chǔ)上對(duì)信號(hào)分析與處理的更深一步的學(xué)習(xí)。11月末，在老師的安排下我們開(kāi)始進(jìn)行關(guān)于由王永德、王軍主編的，由電子工業(yè)出版社出版的《隨機(jī)信號(hào)分析基礎(chǔ)》（第二版），第5章隨機(jī)信號(hào)通過(guò)線性系統(tǒng)的自學(xué)。（1）時(shí)間安排11月末至12月末，每周的周一下午，周四上午設(shè)定為學(xué)習(xí)時(shí)間。（2）目標(biāo)要求，，隨時(shí)做好學(xué)習(xí)相關(guān)知識(shí)的筆記及心得體會(huì)。二、理論原理知識(shí)。
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號(hào)-1