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

Originally posted by sci_chris at 2009-4-9 22:35:


按照楼主的意思，N应该是偶数了？如果是奇数，那N/2就不是一个整数，会导致u不能分成相同长度的向量，这个在matlab中是不能做的。

其实不用编成m文件进行调用，因为最小二乘可以表示成矩阵计算的，就是我在 ...

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;