1stopt里面的最小一乘法拟合曲线

qingqing5638

有谁能用1stopt里面的最小一乘法帮我拟合一个曲线啊？我想看看效果然后建议老师买这个软件！谢谢

dingd

4.0版有专门的最小一乘拟合命令。
拟合优化方面，1stOpt应该是当今最强最好用的。

re-us

最小一乘？是不是要最小绝对值，不是最小平方？如果是的话，我有个matlab的例子

qingqing5638

回复 3 # re-us 的帖子

呵呵能不能程序借我看一下，我的是多元的

re-us

这个是有一次我的homework,应该是对的，irls.m 是老师给的，大概就是迭代的算法

1. %homework for Bob's 2.1
   
2. clear;
   
3. clc
   
4. 
   
5. %precomputed data
   
6. x=1:10;
   
7. x=x';
   
8. 
   
9. y=10:10:100;
   
10. resi=[1.6734,-2.8040,-1.5749,-1.6261,3.5752,1.7820,-2.3950,-4.4920,2.0556,-2.2447];
    
11. y=y+resi;
    
12. y=y';
    
13. %outlier
    
14. y(4)=y(4)-50;
    
15. N=length(y);
    
16. 
    
17. %build the parabolic system matrix
    
18. G = [ ones(N,1) , x];
    
19. %least-squares solution
    
20. disp('Least-squares solution')
    
21. m2 = G\y
    
22. 
    
23. %1-norm solution
    
24. disp('1-norm solution')
    
25. m1=irls(G,y,1.0e-5,1.0e-5,1,25)
    
26. %Plot from the two models
    
27. mm1=G*m1;mm2=G*m2;
    
28. figure(1)
    
29. bookfonts;
    
30. plot(x,mm1,x,mm2,'k');
    
31. hold on
    
32. plot(x,y,'ok');
    
33. legend('L_1 Fit','L_2 Fit','Location','SouthEast');
    

复制代码

1. function x=irls(A,b,tolr,tolx,p,maxiter)
   
2. % x=irls(A,b,tolr,tolx,p,maxiter)
   
3. %
   
4. % Uses the iteratively reweight least squares strategy to find an
   
5. % approximate L_p solution to Ax=b.
   
6. 
   
7. [m,n]=size(A);
   
8. 
   
9. R=eye(m);
   
10. x=A\b;
    
11. 
    
12. iter=1;
    
13. while (iter <= maxiter)
    
14. iter=iter+1;
    
15. r=A*x-b;
    
16. for i=1:m
    
17. if (abs(r(i)) < tolr)
    
18. r(i)=abs(tolr)^(p-2);
    
19. else
    
20. r(i)=abs(r(i))^(p-2);
    
21. end
    
22. end
    
23. R=diag(r);
    
24. newx=(A'*R*A)\(A'*R*b);
    
25. if (norm(newx-x)/(1+norm(x)) < tolx)
    
26. return;
    
27. else
    
28. x=newx;
    
29. end
    
30. end
    
31. 
    

复制代码

re-us

大概的意思就是：如果有outliers的话，用最小1乘比较好，希望有帮助