# how to detemine this residual matrix ?

15 ビュー (過去 30 日間)
Ker Nad 2020 年 8 月 10 日
コメント済み: Peng Li 2020 年 8 月 12 日
I calculate the cumulative sum of a matrix :Gm,n(i,j)
Now, the next step is '' we adopt the simplest function of a plane-fitting (i.e.Gm,n(i,j)=ai+bj+c ) to fit the trending for each surface Gm,n and determine the residual matrix ym,n(i,j).''
Please, can someone give me the code how to detemine this residual matrix, please ?

#### 2 件のコメント

Peng Li 2020 年 8 月 10 日
You question isn't clear enough for people to provide valuable comments honestly.
Ker Nad 2020 年 8 月 10 日
It is explained like this in the paper.
To estimate the eq 2, I need to determine ''the residual matrix''.
Sincerely, I myself didn't understand, especially the fact that I'm not a mathematician, but I need to program this method.

サインインしてコメントする。

### 回答 (1 件)

Peng Li 2020 年 8 月 10 日
Interesting. This is an application of the detrended fluctuation analysis (DFA) to a 2D image. Based on what your screenshot shows, it implements the algorithm similarly like being implemented to a time series -- cut into segments based on a time scale s (or here a time-spatial scale), integration (cumulative sum), linear fitting to get residual, and finally there should be a log-log fit between s and F(s).
To answer the specific question you asked, the residual is nothing but the point difference between G_{m,n} and \tilde{G}_{m,n}. This means, the y variable y_{m,n} = G_{m,n} - \tide{G}_{m,n}. F is simply the average of sum of squared y.

#### 4 件のコメント

Peng Li 2020 年 8 月 10 日
\tilde{G}_{m,n} is the linear 2D fitting ai+bj+c. Basically c is the intercept, a and b being the slope along i and j respectively. Could you try coding it by yourself first? It should be straightforward now if you understand the algorithm. It's a bit wierd I guess if I do that for you lol
Ker Nad 2020 年 8 月 11 日
Sincerely, I have no idea, I have never programmed the equations in Matlab.
Thank you anyway.
Peng Li 2020 年 8 月 12 日
understand and I would just like to encourage you to try it out and ask back here if you get any issues. That way you probably have higher chance to get a solution. I will see what I can do but honestly I have short of time this week and next to meet multiple deadlines, conferences ...

サインインしてコメントする。

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by