MATLAB Answers

how to detemine this residual matrix ?

15 ビュー (過去 30 日間)
Ker Nad
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
Peng Li 2020 年 8 月 10 日
You question isn't clear enough for people to provide valuable comments honestly.
Ker Nad
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
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 件のコメント

表示 1 件の古いコメント
Peng Li
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
Ker Nad 2020 年 8 月 11 日
Sincerely, I have no idea, I have never programmed the equations in Matlab.
Thank you anyway.
Peng Li
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