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

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;

Originally posted by fspdlh at 2009-4-9 10:48:
哦，忘了说了，u,v都是列向量
N=6;
u=[1 3 3 4 5 6]';
v=[1.5 2.5 3.5 4.5 5.5 6.5]';
y=solveuv(u,v,N)

Originally posted by 飞扬2282 at 2009-4-9 11:19:


非常感谢！

其中%[A B]*[u v]'=C 是否应为

%[A B]*[v -u]'=C

？

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;