:round(nFFT/2-1);k=t*Fs/nFFT;p實(shí)現(xiàn)特征向量法的功率譜估計(jì)調(diào)用方式（1）s=peig(x,psin(2*pi*100*n)+4*sin(2*pi*500x1);xlabel('NW=2');subplot(3,1x:round(nFFT/2-1);k=t*Fs/nFFT;p實(shí)現(xiàn)特征向量法的功率譜估計(jì)調(diào)用方式（1）s=peig(x,psin(2*pi*100*n)+4*sin(2*pi*500x1);xlabel('NW=2');subplot(3,1x信號(hào)運(yùn)算MATLA現(xiàn):x=x1+x22、信號(hào)延遲3、信號(hào)乘4、信號(hào)變化幅度n2nn1n2nMATLA現(xiàn)：Ex=sum(abs(x)^2)n2))信號(hào)采樣積數(shù)學(xué)描述：y信號(hào)能量nn1n2nn1x(nx1);xlabel('NW=2');subplot(3,1應(yīng)非線性組合算法，默認(rèn)值Unity:相同加權(quán)的線性組合n2))信號(hào)采樣積數(shù)學(xué)描述：y信號(hào)能量nn1n2nn1x(nx1);xlabel('NW=2');subplot(3,1應(yīng)非線性組合算法，默認(rèn)值Unity:相同加權(quán)的線性組合EigA);EV=fliplr(EV(:,I));%MUSICfo1N1NMATLAB矩形窗w=boxcar(n)布萊克曼窗w=blackman(n)漢寧窗w=hanning(n)e(iang);a=exp(-j*twpi*d*sin(phSD估計(jì)，結(jié)果如圖所示：randn('state');n=0);[p1,f]=pcov(x,20,[],Fs);[p2,w的范圍為[0,Pi];e(iang);a=exp(-j*twpi*d*sin(phSD估計(jì)，結(jié)果如圖所示：randn('state');n=0);[p1,f]=pcov(x,20,[],Fs);[p2,w的范圍為[0,Pi];若x為復(fù)信號(hào)，則w的范圍為[0,2*數(shù)字濾波器的特性分析2、頻率響應(yīng)（幅頻響應(yīng)和相頻響應(yīng)）e(iang);a=exp(-j*twpi*d*sin(phSD估計(jì)，結(jié)果如下圖所示：randn('state',1);x,w]=peig(……)(5)[Pxx,f]=peig(…為[0,Fs];（5）[Pxx,f]=pmtm(……，F(xiàn)s,100e(iang);a=exp(-j*twpi*d*sin(phSD估計(jì)，結(jié)果如下圖所示：randn('state',1);x,w]=peig(……)(5)[Pxx,f]=peig(…為[0,Fs];（5）[Pxx,f]=pmtm(……，F(xiàn)s,10011001MhPh10-10-50-10010-1putSNR(dB)%A=exp(-j*twpi*d.'*sep+1;%un=urn+juinurn=normrnd(0);[p1,f]=pcov(x,20,[],Fs);[p2,exp(-j*2*pi*d*sin(doa(k))/lambputSNR(dB)%A=exp(-j*twpi*d.'*sep+1;%un=urn+juinurn=normrnd(0);[p1,f]=pcov(x,20,[],Fs);[p2,exp(-j*2*pi*d*sin(doa(k))/lamb功率譜估（ＰＳＤ）1、mean函數(shù)（1）y=mean(X)：當(dāng)X為向量時(shí)，此函數(shù)結(jié)果為X的均值；當(dāng)X為矩陣時(shí)，（3）dim=1時(shí)，函數(shù)結(jié)果為一個(gè)行向量，其元素分別為矩陣每列元素的均值。則它的每一列相當(dāng)于一個(gè)變量，函數(shù)返回結(jié)果為該矩陣的列與列之間的A=cov(X)；%計(jì)算協(xié)方差度為2M-1的互相關(guān)函數(shù)序列；當(dāng)A和B長(zhǎng)度不同時(shí)，則要對(duì)長(zhǎng)度小的進(jìn)行補(bǔ)零；如果A為列向量，則C也為列向量，如果A為行向量，則C也為行向量。矩陣，該矩陣的列是由矩陣A所有列之間的互相關(guān)函數(shù)構(gòu)成。對(duì)序列進(jìn)行歸一化處理計(jì)算序列的非歸一化相關(guān)果為X的方差；當(dāng)X為矩陣時(shí)，則它的每一列相當(dāng)于一個(gè)變量，函數(shù)…,Fs)(6)[s,f]=peig(……，果為X的方差；當(dāng)X為矩陣時(shí)，則它的每一列相當(dāng)于一個(gè)變量，函數(shù)…,Fs)(6)[s,f]=peig(……，NW，noverartlett(n)三角窗w=triang(n)布萊克曼窗wmehod):Method有：adapt:Thomson自適Periodogram函數(shù)：功率譜估計(jì)Pxx=periodogram(x,window):參數(shù)window用來(lái)指定所采用的窗函數(shù)若x為復(fù)信號(hào)，則Pxx的長(zhǎng)度為NFFT；同時(shí)返回和估計(jì)PSD的位置一一對(duì)應(yīng)的線性頻率f,參數(shù)Fs為采樣頻率。Fs=2000;NFFT=1024;t(f,Pxx1,'r:',f,Pxx2,'g--');ylNn0MATLA現(xiàn)：Px=sum((abs(x)^2)/NM全部功率譜值，需要使用另一個(gè)輸入?yún)?shù)whole.（3）s=pt(f,Pxx1,'r:',f,Pxx2,'g--');ylNn0MATLA現(xiàn)：Px=sum((abs(x)^2)/NM全部功率譜值，需要使用另一個(gè)輸入?yún)?shù)whole.（3）s=p;學(xué)習(xí)好資料歡迎下載endPmusicf(1,i)=10*loPo00例2、利用FFT直接法計(jì)算上面噪聲信號(hào)的功率譜Fs=2000;nFFT=1024;-10-20-30-40”功率譜，若x為復(fù)數(shù)，則返回結(jié)果為“雙邊”功率譜。（2”功率譜，若x為復(fù)數(shù)，則返回結(jié)果為“雙邊”功率譜。（2）Pxx,w]=peig(……)(5)[Pxx,f]=peig(…實(shí)現(xiàn)特征向量法的功率譜估計(jì)調(diào)用方式（1）s=peig(x,p返回結(jié)果為該矩陣的列與列之間的協(xié)方差矩陣，diag(cov(Fs=2000;nFFT=1024;自回歸功率譜估計(jì)的改進(jìn)的協(xié)方差方法Pmcov函數(shù)應(yīng)用實(shí)例：比較兩種方法在噪聲信號(hào)的功率譜估計(jì)中的效果，發(fā)現(xiàn)兩種方法基Fs=1000;x=periodogram(x):返回向量x的功率譜估計(jì)向量函數(shù)調(diào)用方式（1x=periodogram(x):返回向量x的功率譜估計(jì)向量函數(shù)調(diào)用方式（1）y=cov(X):當(dāng)X為向量時(shí)，函數(shù)返回結(jié)orr’)(3)s=peig(x,p,nFFT)(4)[Pxf]=pmtm(x,2,nFFT,Fs);[P2,f]=pm-20協(xié)方差方法-30-40-50-60-70-80-90-1000時(shí)間-帶寬的乘積NW，窗口數(shù)=2*NW-1（2）Pxx=pmtm(x,NW):（3）Pxx=pmtm(x,NW,NFFT):參數(shù)NFFT用來(lái)指定FFT運(yùn)算所采用的點(diǎn)數(shù)。若x為復(fù)信號(hào)，則Pxx的長(zhǎng)度為NFFT；NFFT默認(rèn)值為256果為X的方差；當(dāng)X果為X的方差；當(dāng)X為矩陣時(shí)，則它的每一列相當(dāng)于一個(gè)變量，函數(shù)(0,delta/2,1,sample);un=urn+m*):用特征向量法對(duì)離散時(shí)間信號(hào)x進(jìn)行功率譜估計(jì)。p是信號(hào)x中波峰變寬，這說(shuō)明頻譜泄露增大了。MUSIC法估計(jì)Pmusic同時(shí)返回和估計(jì)PSD的位置一一對(duì)應(yīng)的線性頻率f,參數(shù)Fs為采樣頻率。應(yīng)用說(shuō)明：利用Multitaper進(jìn)行PSD估計(jì)，并比較NW取不同數(shù)值時(shí)的結(jié)果。Fs=1000;nFFT=1024;料歡迎下載4035302520151050010020030旁瓣的高度。n料歡迎下載4035302520151050010020030旁瓣的高度。n一定時(shí)，beta越大，其頻譜的旁瓣越小，但主瓣SD的位置一一對(duì)應(yīng)的線性頻率f,參數(shù)Fs為采樣頻率。若x為實(shí)離散時(shí)間信號(hào)x進(jìn)行功率譜估計(jì)，若x為實(shí)數(shù)，則返回結(jié)果為“單邊-20-40-20-40-20-40NW=2NW=4NW=102、MUSIC法估計(jì)Pmusic函數(shù)：實(shí)現(xiàn)MUSIC法的功率譜計(jì)算列都進(jìn)行功率譜估計(jì)。注意：為了返回實(shí)信號(hào)的全部功率譜值，需要使用另ambda/2;%陣元間距snr=15;%%信噪比%%%ambda/2;%陣元間距snr=15;%%信噪比%%%導(dǎo)向(……)(8)peig(……)應(yīng)用說(shuō)明例：用特征向量法進(jìn)行Pheangleofsi字濾波器的特性分析脈沖響應(yīng)：impz函數(shù)調(diào)用方式：（1）[hNW默認(rèn)值為2*p,參數(shù)noverlap的默認(rèn)值為NW-1一對(duì)應(yīng)的特征值所組成的向量；而e為相關(guān)矩陣的特征值向量。例：用MUSIC法進(jìn)行PSD估計(jì)，結(jié)果如圖所示：ww-10-20iaMUSIC方法，它主要適用于混有在噪聲的正弦信號(hào)的功率譜估計(jì)）Peig函數(shù)：實(shí)現(xiàn)特征向量法的功率譜估計(jì)它主要適用于混有在噪聲的正弦信號(hào)的功率譜估計(jì)）Peig函數(shù)：度。例：[b,a]=butter(4,0.05);impz(而e為相關(guān)矩陣的特征值向量。應(yīng)用說(shuō)明例：用MUSIC它主要適用于混有在噪聲的正弦信號(hào)的功率譜估計(jì)）Peig函數(shù)：度。例：[b,a]=butter(4,0.05);impz(而e為相關(guān)矩陣的特征值向量。應(yīng)用說(shuō)明例：用MUSIC法進(jìn)行P以看出：NW數(shù)值越大，曲線越平滑，說(shuō)明方差比較小。但同時(shí)看到例：用特征向量法進(jìn)行PSD估計(jì)，結(jié)果如下圖所示：P0-10-20-300MUSIC算法);a(d)=1000;endifP>1disp('Thea相頻y=abs(x):計(jì)算);a(d)=1000;endifP>1disp('Thea相頻y=abs(x):計(jì)算x各元素的絕對(duì)值。當(dāng)x為一個(gè)復(fù)數(shù)時(shí)da*[0:M-1]);endB=B';%%%導(dǎo)向向量xx=好資料歡迎下載0-10-20-30-40-50angle(d%zn%xn%*********************************************************數(shù)：協(xié)方差函數(shù)估計(jì)（3）相關(guān)系數(shù)估計(jì)計(jì)算學(xué)習(xí)好資料歡迎下載C);subplot(3,1,3);plot(f,Pxx3);數(shù)：協(xié)方差函數(shù)估計(jì)（3）相關(guān)系數(shù)估計(jì)計(jì)算學(xué)習(xí)好資料歡迎下載C);subplot(3,1,3);plot(f,Pxx3);A);EV=fliplr(EV(:,I));%MUSICfo噪聲，%%%平均后結(jié)果更好[pp,ppp]=size(x);50-5-10-15-20X))是該矩陣每一個(gè)列向量的方差。（2）y=cov(X,Y)0;nFFT=1024;t=0:1/Fs:1;x=sin(2end%xnforn=1:samplexn(n)=2*cosX))是該矩陣每一個(gè)列向量的方差。（2）y=cov(X,Y)0;nFFT=1024;t=0:1/Fs:1;x=sin(2end%xnforn=1:samplexn(n)=2*cosmehod):Method有：adapt:Thomson自適0-10-20-30-40-50%deg->rad%陣列數(shù)量%space%%numberofDOA%角度%inputSNR(dB)%X=A*S;%%%%[EV,D]=eig(Rxx);%%%%%MUSIC10,nFFT,Fs);Pxx1=10*log10(P1);包含的復(fù)數(shù)正弦波信號(hào)的個(gè)數(shù)。(2)s=peig(r,p,’10,nFFT,Fs);Pxx1=10*log10(P1);包含的復(fù)數(shù)正弦波信號(hào)的個(gè)數(shù)。(2)s=peig(r,p,’cf,fs)(2)模擬濾波器頻率響應(yīng)：freqs函數(shù)調(diào)用方式：);subplot(3,1,3);plot(f,Pxx3);end%SP=abs(SP);SPmax=max(SP);SP=10*log10(SP/SPmax);2、用ESPRIT方法求DOAN=200;%%快拍數(shù)M=8;%%陣元數(shù)snr=15;%%信噪比%%%導(dǎo)向向量%%%導(dǎo)向向量%%%%噪聲平均因?yàn)閙atlab產(chǎn)生的噪聲不太好，不是嚴(yán)格意義上的白噪聲