[其他] 想請問一下多筆資料polyfit的寫法

作者: perrywu (無尾熊寶寶)   2015-04-29 19:12:14
各位先進好
小弟是matlab新手
現在有多筆資料需要做polyfit
例:現在有三個matrix(稱A,B,C)大小一樣 分別存放Data
A=[a11 a12...a1n; B=[b11 b12...b1n; C=[c11 c12...c1n;
a21 a22...a2n; b21 b22...b2n; c21 c22...c2n;
.............; .............; ..............;
am1 am2...amn]; bm1 bm2...bmn]; cm1 cm2...cmn];
現在要求出各自matrix 點對點的斜率
所以小弟的寫法是這樣
for x=1:m,
for y=1:n,
a=polyfit([1 5 10],[A(x,y) B(x,y) C(x,y)],1);
Data(x,y)=a(1);
end
end
因為實際上m和n的大小有點大 所以在這樣跑回圈算結果的時候
需要花真的非常久的時間
所以小弟想請問各位先進
有沒有比較不一樣的寫法可以較快速地完成小弟想做的事
真的很謝謝各位大大!!
作者: sunev (Veritas)   2015-04-29 22:13:00
三個點做linear fitting應該有簡單公式解Data=(13*A+B-14*C)/122;
作者: perrywu (無尾熊寶寶)   2015-04-29 22:22:00
想請問大大上面那個式子是怎麼推導的!謝謝
作者: sunev (Veritas)   2015-04-30 01:41:00
作者: perrywu (無尾熊寶寶)   2015-04-30 21:27:00
很感謝大大

Links booklink

Contact Us: admin [ a t ] ucptt.com