freepeople性欧美熟妇, 色戒完整版无删减158分钟hd, 无码精品国产vα在线观看DVD, 丰满少妇伦精品无码专区在线观看,艾栗栗与纹身男宾馆3p50分钟,国产AV片在线观看,黑人与美女高潮,18岁女RAPPERDISSSUBS,国产手机在机看影片

正文內(nèi)容

基于matlab的譜估計(jì)實(shí)現(xiàn)畢業(yè)設(shè)計(jì)論文-文庫(kù)吧資料

2025-06-28 21:20本頁(yè)面
  

【正文】 主瓣寬度比較窄,頻率泄露的情況相對(duì)不太明顯(也就是除主瓣以外的譜線的幅度較主瓣幅度相比較小) 。triang窗39。figure(8)plot(f7,10*log10(Pxx7))。kaiser窗39。figure(7)plot(f6,10*log10(Pxx6))。hann窗39。figure(6)plot(f5,10*log10(Pxx5))。hamming窗39。figure(5)plot(f4,10*log10(Pxx4))。blackman窗39。figure(4)plot(f3,10*log10(Pxx3))。bartlett窗39。figure(3)plot(f2,10*log10(Pxx2))。boxcar窗39。figure(2)plot(f1,10*log10(Pxx1))。)。title(39。subplot(337)。kaiser窗39。plot(f6,10*log10(Pxx6))。)。title(39。subplot(335)。hamming窗39。plot(f4,10*log10(Pxx4))。)。title(39。基于 MATLAB 的譜估計(jì)實(shí)現(xiàn)15subplot(333)。bartlett窗39。plot(f2,10*log10(Pxx2))。)。title(39。figure(1)subplot(331)。[Pxx6,f6]=periodogram(xn,window6,nfft,Fs)。[Pxx4,f4]=periodogram(xn,window4,nfft,Fs)。%直接法[Pxx2,f2]=periodogram(xn,window2,nfft,Fs)。window7=triang(length(xn))。window5=hann(length(xn))。window3=blackman(length(xn))。window1=boxcar(length(xn))。%產(chǎn)生含有噪聲的序列xn=cos(2*pi*40*n)+cos(2*pi*41*n)+3*cos(2*pi*90*n)+*randn(size(n))。 加窗周期圖的仿真采用加窗周期圖方法進(jìn)行功率譜估計(jì)的仿真程序如下:Fs=600。figure(1)plot(k,plot_Pxx) 數(shù)據(jù)重疊方法的仿真圖如圖 35 所示。k=index*Fs/nfft。[Pxx,f]=pwelch(xn,window,noverlap,nfft)。window=boxcar(100)。xn=cos(2*pi*40*n)+5*cos(2*pi*90*n)+*randn(size(n))。程序如下:Fs=600。end。 sx=sx+p。for k=1:K [p,f]=periodogram(xn(n1:(n1+L1)))。f=0。 function [sx,f]=bart(xn,K)L=floor(length(xn)/K)。%k相當(dāng)于頻率,256是PERIODOGRAM(在bart函數(shù)中調(diào)用的)函數(shù)中默認(rèn)的nfft 的值plot_Pxx=10*log10(Pxx(index+1))。index=0:round(256/2 1)。K=10。 %采樣頻率n=0:1/Fs:1。如圖 41,真實(shí)的譜線應(yīng)該是在頻率為 40Hz 和 90Hz 處為兩個(gè)沖激,而周期圖法估計(jì)出的功率譜在這兩個(gè)基于 MATLAB 的譜估計(jì)實(shí)現(xiàn)13頻率出的主瓣非常寬,失真較明顯。0 50 100 15060504030201001020圖33由前述討論可知,在周期圖方法中,無(wú)論怎樣增大數(shù)據(jù)采樣點(diǎn)數(shù) ,功率譜估N計(jì)的方差一直不為零。[Pxx,f]=periodogram(xn,window,nfft,Fs)。 xn=cos(2*pi*40*n)+5*cos(2*pi*90*n)+*randn(size(n))。仿真程序如下:Fs=600。采用加窗周期圖方法的數(shù)學(xué)表達(dá)式為 10[]NjwnnXe?????? ??其中, 為窗函數(shù)。幾種常用窗函數(shù)的數(shù)學(xué)表達(dá)式如下:矩形窗: ()1,0nN????三角窗(bartlett): 2/(1),0()/2nnN??????漢寧窗: 1()cos,21nN?????????????????哈明窗: 2(),0nN???????布萊克曼窗: 24().,011nnN?????????????????????基于 MATLAB 的譜估計(jì)實(shí)現(xiàn)11在MATLAB中還給出了幾種常用的窗函數(shù),他們的時(shí)域和頻域特性分別如下:50 0 5001boxcar域域域域域域50 0 5001bartlett域域域域域域50 0 5001blackman域域域域域域50 0 5001hamming域域域域域域50 0 5001hann域域域域域域50 0 5001kaiser域域域域域域50 0 5001triang域域域域域域圖31500 0 5000204060 boxcar域域域域域域500 0 5000102030 bartlett域域域域域域500 0 5000102030blackman域域域域域域500 0 5000102030hamming域域域域域域500 0 5000102030 hann域域域域域域500 0 50002040 kaiser域域域域域域500 0 5000102030 triang域域域域域域圖32從這幾種窗函數(shù)的頻域特性可以看出,矩形窗的主瓣寬度最小,因此運(yùn)用矩形窗函數(shù)以后,譜分辨率是最高的,但旁瓣幅度較大,功率泄露情況也比較嚴(yán)重。這就是說(shuō),若信號(hào)中包含兩個(gè)以上正弦信號(hào)且頻率非常接近時(shí),兩個(gè)分量的頻譜疊加的結(jié)果,兩個(gè)正弦信號(hào)不可能分辨,使信號(hào)模糊,降低了分辨率。而 函數(shù)的主瓣不是無(wú)限窄,主瓣寬度反比于數(shù)據(jù)記錄長(zhǎng)度 ,而在實(shí)際中又不可Sa N能獲得很長(zhǎng)的數(shù)據(jù)長(zhǎng)度。時(shí)域上乘以一矩形窗相當(dāng)于頻域與一個(gè) 函數(shù)進(jìn)行卷積。為了在平均周期圖的基礎(chǔ)上更進(jìn)一步降低用周期圖法進(jìn)行功率譜估計(jì)的方差,可以采用的方法是使每一個(gè)分段的數(shù)據(jù)有一些重疊,使得觀測(cè)的數(shù)據(jù)段的數(shù)目增多,以此來(lái)達(dá)到進(jìn)一步減小方差的目的。 平均周期圖由于周期圖估計(jì)的方差總是不趨于零,所以需要對(duì)其進(jìn)行改進(jìn),其中一種方法就是平均周期圖法。而且,無(wú)論 取何值,??4r[()]0XI???N總與 在同一數(shù)量級(jí)。,。1??根據(jù)式(32)可得 122122[()]()()jjNNNEIEXe?????????? 2[][]NNnkpqRnkpRq?基于 MATLAB 的譜估計(jì)實(shí)現(xiàn)9 (35) 12()()[][]jnkjpqEXnkpqe????其中 為矩形窗函數(shù)。按方差的定義應(yīng)有2X? 22var[()][()][()]NNNIEII????由于假設(shè) 零均值,方差為 的實(shí)高斯白噪聲序列,根據(jù)式(34)可知,當(dāng)[]nX?時(shí), 。下面討論譜估計(jì)的方差。 211[()][()][()]j jjNNNEIXeEXe?????1100[()()]jnjknk????? [][[]]jnjknkEdXedXeN??????????? ()1[[][]]jknnk ??????令 代入上式得mkn??1[()][][]]jmNmnEIdmEXne???? ???? ???? ?? ?? (33)[]jNXmbRe??????式中 1[][]]Nnbdm???? ,10N?????其 他由于 為兩個(gè)矩形窗函數(shù)的卷積,因此它是一個(gè)三角窗函數(shù),其傅里葉變[]Nbm換為基于 MATLAB 的譜估計(jì)實(shí)現(xiàn)8 22sin11()[]jnNn NBde??????????????????????式(33)為兩個(gè)函數(shù)乘積的傅里葉變換,等于他們各自傅里葉變換的卷積,即 [()]()NXEIBG???由上式可知,除非 是一個(gè)沖激函數(shù),否則 將不等于 ,故B[()]NEI()XG?周期圖作為功率譜的估計(jì)是有偏差的。由于序列的傅里葉變化可利用 計(jì)算而提高運(yùn)算效率,這是周期圖FT法的主要優(yōu)點(diǎn)。當(dāng) 為實(shí)()jNXe?[]Xn()jN[]n序列時(shí), 也為實(shí)序列,根據(jù)傅里葉變換的性質(zhì), 的傅里葉變換為 。N?[]定義矩形窗函數(shù)如下 1,01[]nNd??????其 他則式(31)可寫成 1?[][][][]XnRmdXnmN??????令 ,上式可簡(jiǎn)寫為[][]dnXY? 11?[][]XnRYY???????即 可看成是 與 的卷積。下面逐一闡述。(1)k?Xk相合性是對(duì)一個(gè)估計(jì)量的基本要求,若估計(jì)量不具有相合性,那么不論將樣本容量 取得多么大,都不能將 估計(jì)得足夠準(zhǔn)確,這樣的估計(jì)量是不可取的。12(,)n ?即,若對(duì)于任意 都滿足:對(duì)于任意 ,有???0???lim{}1nP??????則稱 是 的相合估計(jì)量。這樣,對(duì)估計(jì)量又有下述相合性的要求。?1??2 相合性無(wú)偏性和有效性都是在樣本容量 固定的前提下提出的。??可以證明,樣本均值(也就是總體的一階矩)是總體均值的無(wú)偏估計(jì)量,樣本方差是總體方差的無(wú)偏估計(jì)量。 是包含在總體 的分布中的待估參12,nX…X???X基于 MATLAB 的譜估計(jì)實(shí)現(xiàn)5數(shù),這里 是 的取值范圍。它給出了平穩(wěn)隨機(jī)過程的時(shí)域特性與頻域特性之間的聯(lián)系,它是分析隨機(jī)信號(hào)的一個(gè)最重要,最基本的公式。?4,功率譜密度可積。 2,功率譜密度是 的實(shí)函數(shù)。)XnXREXxfndx????? 離散隨機(jī)過程的功率譜密度 離散隨機(jī)過程功率譜密度定義設(shè) 為寬平穩(wěn)離散時(shí)間隨機(jī)過程,或簡(jiǎn)稱為寬平穩(wěn)隨機(jī)序列,具有零均值,[]Xn其自相關(guān)函數(shù)為基于 MATLAB 的譜估計(jì)實(shí)現(xiàn)4(){[]]}XRmEnTXm??或簡(jiǎn)寫為 (21)()[]]X當(dāng) 滿足 的條件時(shí),我們定義 的功率譜密度為 的()XRm()XR??????[]Xn()XRm離散傅里葉變換,并記為 ()G?()jmTXXmRe??????式中, 是隨機(jī)序列相鄰各值的時(shí)間間隔。[][]nmnmEXYE?nXmY離散時(shí)間隨機(jī)過程的均方值定義為 222[](。[][]nmnmXYY?(2) ,即 乘以一個(gè)常數(shù) 的均值等于 的均值乘以此常Ea?XanX基于 MATLAB 的譜估計(jì)實(shí)現(xiàn)3數(shù)。)nxnXmExfnd?????若 是單值函數(shù),則 構(gòu)成一個(gè)新的離散時(shí)間隨機(jī)過程,其均值可定義()gA{}g為 []()。本篇論文主要涉及的信號(hào)類型屬于離散隨機(jī)信號(hào),它一般是對(duì)一個(gè)連續(xù)信號(hào)進(jìn)行抽樣,對(duì)抽樣所得的序列進(jìn)行功率譜估計(jì),故它屬于連續(xù)隨機(jī)序列。0,12,kX…, ks(4) 離散型隨機(jī)序列:相應(yīng)于時(shí)間和狀態(tài)都是離散的情況。(3) 連續(xù)隨機(jī)序列,隨機(jī)過程 在任一離散時(shí)刻的狀態(tài)是連續(xù)型隨機(jī)變()t量,也就是時(shí)間離散,狀態(tài)連續(xù)的情況?;?MATLAB 的譜估計(jì)實(shí)現(xiàn)22 隨機(jī)過程及估計(jì)理論基本知識(shí) 隨機(jī)過程的基本知識(shí) 隨機(jī)過程的分類按照時(shí)間和狀態(tài)(一般稱隨機(jī)過程 在 的可能取值為它的狀態(tài))是連續(xù)()iXtit?的還是離散的來(lái)分類,可以分成以下四類:(1) 連續(xù)型隨機(jī)過程: 對(duì)于任意的 , 都是連續(xù)型隨機(jī)變量,()ttT?()Xt也就是時(shí)間和狀態(tài)都是連續(xù)的情況。在現(xiàn)代譜估計(jì)的方法中,不再認(rèn)為在觀測(cè)到的 個(gè)數(shù)據(jù)以外的數(shù)據(jù)全為零,因此克服了經(jīng)典法的缺點(diǎn),并且提高的譜估計(jì)的分辨率,能給出較好的估計(jì)。另一種是直接法,通過對(duì)觀測(cè)到的數(shù)據(jù)樣本進(jìn)行傅里葉變換,然后將所得到結(jié)果的幅值平方后得到功率譜估計(jì),這種方法又稱為周期圖法。經(jīng)典譜估計(jì)法通常又分為兩種,一種是間接法,它是由布萊克曼()和圖基()提出,又稱為 BT 法。一類是非參數(shù)化方法,另一類是參數(shù)化方法。例如,根據(jù)信號(hào)、干擾與噪聲的功率譜,可以設(shè)計(jì)適當(dāng)?shù)臑V波器,以盡量不失真地重現(xiàn)信號(hào),而最大幅度地抑制干擾與噪聲。在電子戰(zhàn)中,譜分析可用來(lái)對(duì)目標(biāo)進(jìn)行分類識(shí)別等。IIABSTRACTPower Spectral Estimation is a fundamental topic in the signal processing field. Actually
點(diǎn)擊復(fù)制文檔內(nèi)容
醫(yī)療健康相關(guān)推薦
文庫(kù)吧 www.dybbs8.com
備案圖鄂ICP備17016276號(hào)-1