matlab中自相關(guān)函數(shù)xcorr

1、2627matlab中自相關(guān)函數(shù)xcorr1、介紹一下MATLAB函數(shù)xcorr的使用2、自相關(guān)是什么意思3、兩個(gè)數(shù)相關(guān)為什么是兩個(gè)數(shù)的乘積xcorr(2,2)=4%?給幾個(gè)簡單例子講解一下自相關(guān)函數(shù)是描述隨機(jī)信號X(t)在任意兩個(gè)不同時(shí)刻tl,t2的取值之間的相關(guān)程度.設(shè)原函數(shù)是f(t),則自相關(guān)函數(shù)定義為R(u)=f(t)*f(-t),其中*表示卷積.給個(gè)例子:dt=.1;t=0:dt:100;x=cos(t);a,b=xcorr(x,unbiased);plot(b*dt,a)上面代碼是求自相關(guān)函數(shù)并作圖,matlab中查看幫助時(shí),helpxcorr解釋其意思是：C(m)=EA(n+m)

2、*conj(B(n)=EA(n)*conj(B(n-m)；但是,在調(diào)用xcorr函數(shù)求自相關(guān)時(shí)，有scaleopt參數(shù)r=xcorr(s,SCALEOPT)SCALEOPT有biased-scalestherawcross-correlationby1/M.unbiased-scalestherawcorrelationby1/(M-abs(lags).coeff-normalizesthesequencesothattheauto-correlationsatzerolagareidentically1.0.none-noscaling(thisisthedefault).注意觀察下面的測試

3、:s=123r=xcorr(s);r=3.00008.000014.00008.00003.0000當(dāng)用r=xcorr(s,unbiased)時(shí)就能得到r=3.00004.00004.66674.00003.00002627Fft歷程主要內(nèi)容(signal_analyzer.rar)如下:testl.txt第一列為時(shí)間t,0.001s為步長；第二列為系統(tǒng)的輸入信號；第三列為系統(tǒng)的輸出信號。signal_analyzer.m首先將testl.txt的數(shù)據(jù)讀入workspace，然后進(jìn)行FFT,分析出兩信號的頻率，輸出對于輸入的幅值衰減(db),輸出對輸入的相位差。信號信號8求付利葉變換實(shí)際頻率點(diǎn)付利葉變換1,141求能量最大點(diǎn)fz1,f1=max(abs(sinfftfz2,f2=max(abs(cosfft19求兩信號能量最大點(diǎn)的頻率得到實(shí)際頻率2223求兩信號能量最大點(diǎn)的相位28信號頻率信號頻率信號比信號幅值衰減信號比信號相位差度復(fù)制代碼