| 4 | 1/1 | 返回列表 |
| 查看: 841 | 回復: 3 | |||
| 【懸賞金幣】回答本帖問題,作者我是陳小生將贈送您 100 個金幣 | |||
我是陳小生金蟲 (著名寫手)
|
[求助]
求助:已知高度值矩陣z(x,y),如何計算其沿x、y方向的自相關長度correlation length? 已有1人參與
|
||
|
如題,已知高度值矩陣z(x,y),如何計算其沿x、y方向的自相關長度correlation length? 有大神做過嗎? @beefly @月只藍 發(fā)自小木蟲Android客戶端 |


|
自相關矩陣和互相關矩陣的matlab實現 一維實值x的自相關矩陣Rxx應為實對稱的toeplitz矩陣,而一維實值信號x,y的互相關矩陣Rxy為非對稱的toeplitz陣,matlab提供的corrmtx產生的并非通常意義下的autocorrelation matrix 事實上,我們可以利用xcorr+toeplitz和corrmtx兩種方法實現自相關陣Rxx和互相關陣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' 上面的方法實現了自相關和互相關的有偏矩估計(實際是用實現卷積的前提下做到的),也是做統(tǒng)計分析的常用手段,當然除了有偏矩估計,corrmtx還有很多可選參數,以供不同目的使用。 還有Rxy=Ryx',也就是說要求Ryx,只需要計算Rxy即可。 需要注意的是,當時延m+1接近于信號長度n的時候,xcorr后面的值(rx(m+1),rx(m), r(m-1)...)估計的并不準確,這將嚴重影響器等后續(xù)工作的效果,可以確信的是,n-m-1>100的時候,Rxx是可信的。 另外,自相關陣Rxx是一個主對角線絕對占優(yōu)陣,也就是說,主對角線的值遠大于其它對角線 這樣svd(Rxx)得到的奇異值和eig(Rxx)得到的特征值幾乎相同,這也是為什么有的用svd方法求主分量,有的則用evd(eigen value decomposition)求主分量了 發(fā)自小木蟲Android客戶端 |

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703求調劑 +5 | jtyq001 2026-03-10 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 301求調劑 +5 | Liyouyumairs 2026-03-11 | 5/250 |
|
|
[考研] 一志愿中科大材料與化工。353分能調劑到哪? +5 | MayUxw1 2026-03-08 | 7/350 |
|
|
[考研] 274環(huán)境工程求調劑 +8 | 扶柳盈江 2026-03-05 | 9/450 |
|
|
[考研] 一志愿河海大學085900土木水利專碩279求調劑不挑專業(yè) +3 | SunWwWwWw 2026-03-10 | 3/150 |
|
|
[考研] 標題:撿漏預警|08工科/09農學調劑!英語要求低,過線即有機會! +7 | 馬超放煙花 2026-03-07 | 12/600 |
|
|
[考研] 求調劑! +3 | 朔朔話 2026-03-09 | 3/150 |
|
|
[考研] 0817化學工程319求調劑 +7 | lv945 2026-03-08 | 9/450 |
|
|
[考研]
|
likeihood 2026-03-06 | 16/800 |
|
|
[考研] 294 英二數二物化 求調劑 +6 | 米飯團不好吃 2026-03-09 | 6/300 |
|
|
[考研] 數二英二309分請求調劑 +3 | dtdxzxx 2026-03-09 | 4/200 |
|
|
[考研] 337求調劑 +3 | 睡醒,。 2026-03-09 | 3/150 |
|
|
[考研] 269求調劑 +3 | 朔朔話 2026-03-08 | 4/200 |
|
|
[考研] 347求調劑 +4 | 浮云滿足 2026-03-07 | 4/200 |
|
|
[考研] 085701環(huán)境工程317分求調劑 +9 | 6汆尼9 2026-03-07 | 9/450 |
|
|
[考研] 材料化工求調劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 346分材料求調劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|
[考博] 2026申博自薦 六級440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|
|
[考研] 334求調劑 +6 | Trying] 2026-03-05 | 8/400 |
|
|
[考研] 282求調劑 +7 | 夕~日 2026-03-05 | 8/400 |
|