[交流] 【求助】Matlab最小二乘法简单编程代码？

Originally posted by sci_chris at 2009-4-9 15:08:
用最小二乘法最简单：把

x1i = (uj-ui);
x2i = -(vj-vi);

放到一个矩阵中：X = [x1i x2i]；
而u,v可以看作是系数向量：b = [u;v]；
yi = vi*uj-vj*ui；

则可以建立方程：Xb = y；

利用最小二乘法的 ...

Originally posted by fspdlh at 2009-4-8 18:46:
function y=solveuv(u,v,N)
%Av-Bu=C
ind=floor(N/2);
A=u(ind+1:ind*2)-u(1:ind);
B=v(ind+1:ind*2)-v(1:ind);
C=v(1:ind).*u(ind+1:ind*2)-v(ind+1:ind*2).*u(1:ind);
%[A B]*[u v]'=C
y=([A B])\C;