[其他] 想请问一下多笔资料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