| 5 | 1/1 | 返回列表 |
| 查看: 840 | 回復(fù): 3 | |||
| 【懸賞金幣】回答本帖問題,作者我是陳小生將贈送您 100 個金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
我是陳小生金蟲 (著名寫手)
|
[求助]
求助:已知高度值矩陣z(x,y),如何計(jì)算其沿x、y方向的自相關(guān)長度correlation length? 已有1人參與
|
||
|
如題,已知高度值矩陣z(x,y),如何計(jì)算其沿x、y方向的自相關(guān)長度correlation length? 有大神做過嗎? @beefly @月只藍(lán) 發(fā)自小木蟲Android客戶端 |
|
自相關(guān)矩陣和互相關(guān)矩陣的matlab實(shí)現(xiàn) 一維實(shí)值x的自相關(guān)矩陣Rxx應(yīng)為實(shí)對稱的toeplitz矩陣,而一維實(shí)值信號x,y的互相關(guān)矩陣Rxy為非對稱的toeplitz陣,matlab提供的corrmtx產(chǎn)生的并非通常意義下的autocorrelation matrix 事實(shí)上,我們可以利用xcorr+toeplitz和corrmtx兩種方法實(shí)現(xiàn)自相關(guān)陣Rxx和互相關(guān)陣Rxy 一、 Rxx 1)% implementation with xcorr and toeplitz m= ;% dfine the time lag m+1, and m+1<=n; n=length(x);%location of rxx(0); rx=xcorr(x);%length of rx is 2*n-1; Rxx=toeplitz(rx(n:n+m))/n; 2)%implementation with corrmtx m= ;% dfine the time lag m+1,and m+1<=n rx=corrmtx(x,m); Rxx=rx'*rx; 二、Rxy 1)% implementation with xcorr and toeplitz m= ;% dfine the time lag m+1, and m+1<=n; n=max(length(x),length(y));location of rxy(0); rxy=xcorr(x,y);%length of rxy is 2*n-1; RR=toeplitz(rxy)/n;%RR is a (2*n-1)*(2*n-1) matrix Rxy=RR(1:m,n:n+m);%the exact location of Rxy in RR; 2)% implementation with corrmtx m= ;% dfine the time lag m+1, and m+1<=n; rx=corrmtx(x,m); ry=corrmtx(y,m); Rxy=rx'*ry; %on the other hand, Ryx=Rxy' 上面的方法實(shí)現(xiàn)了自相關(guān)和互相關(guān)的有偏矩估計(jì)(實(shí)際是用實(shí)現(xiàn)卷積的前提下做到的),也是做統(tǒng)計(jì)分析的常用手段,當(dāng)然除了有偏矩估計(jì),corrmtx還有很多可選參數(shù),以供不同目的使用。 還有Rxy=Ryx',也就是說要求Ryx,只需要計(jì)算Rxy即可。 需要注意的是,當(dāng)時延m+1接近于信號長度n的時候,xcorr后面的值(rx(m+1),rx(m), r(m-1)...)估計(jì)的并不準(zhǔn)確,這將嚴(yán)重影響器等后續(xù)工作的效果,可以確信的是,n-m-1>100的時候,Rxx是可信的。 另外,自相關(guān)陣Rxx是一個主對角線絕對占優(yōu)陣,也就是說,主對角線的值遠(yuǎn)大于其它對角線 這樣svd(Rxx)得到的奇異值和eig(Rxx)得到的特征值幾乎相同,這也是為什么有的用svd方法求主分量,有的則用evd(eigen value decomposition)求主分量了 發(fā)自小木蟲Android客戶端 |



| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +4 | mingmingry 2026-03-09 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 家人們 調(diào)劑不迷路 看這里 +8 | likeihood 2026-03-09 | 13/650 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工296分求調(diào)劑 +8 | 稻妻小編 2026-03-09 | 8/400 |
|
|
[考研]
|
簡木ChuFront 2026-03-09 | 4/200 |
|
|
[考博] 找博導(dǎo) +3 | 小呆呆熊 2026-03-07 | 3/150 |
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
[考研] 288求調(diào)劑 +5 | 好好- 2026-03-07 | 5/250 |
|
|
[考研] 22408 275分求調(diào)劑 +3 | 宇智波比 2026-03-03 | 4/200 |
|
|
[考研] 332材料求調(diào)劑 +6 | zjy101327 2026-03-05 | 7/350 |
|
|
[考研] 考研282分求調(diào)劑,接受跨專業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|
|
[考研] 一志愿西交化工專碩288專業(yè)課93求調(diào)劑 +6 | 好運(yùn)好運(yùn)接接 2026-03-04 | 6/300 |
|
|
[考研] 材料328求調(diào)劑 +10 | 一個蘿卜02 2026-03-03 | 10/500 |
|
|
[考研] 070300化學(xué) 280 一志愿太原理工 求調(diào)劑 +8 | 拾玖壹 2026-03-04 | 8/400 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-04 | 3/150 |
|
|
[考研] 085600 材料與化工 298 +14 | 小西笑嘻嘻 2026-03-03 | 14/700 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|
|
[考研] 293求調(diào)劑 +4 | 是樂渝哇 2026-03-03 | 4/200 |
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|