如何进行多元拟合？

feiyuzhen

看maple里的帮助没有这方面的例子，只是提了一下，
请大家给个例子，谢谢，两个变量就行，用statistics包里的Fit命令

mathfan2008

用Maple时有两招必用，一是看自带的帮助系统，还有一个比较好的资源系统：
http://www.maplesoft.com/applications/
输入fit，可以找到很多曲线拟合的例子。

feiyuzhen

谢谢楼上的回答，在帮助里没有找到多元拟合，最多只有一元的，
资源系统里面也没有的

TBE_Legend

多元的不多见，1stopt应该可以。

feiyuzhen

这个已经搞定了，
多谢斑竹

qiuqiu

请问你是如何做的啊？能详细介绍一下吗？ :D

feiyuzhen

http://www.cca-es.com/bbs/dispbb ... ID=21664&page=3

TBE_Legend

Note that a model function can be nonlinear in the independent variables but linear in the model parameters.

alljoyland

用matlab做拟合,更快一些

feiyuzhen

1. restart;
   
2. with(stats);
   
3. with(fit);
   
4. with(statplots);
   
5. with(plots);
   
6. datax := [seq((rand(100))(), i = 1 .. 50)];
   
7. datay := map(proc (x) options operator, arrow; evalf(sin((1/100)*x/Pi)+(1/100)*(rand(30))(), 4) end proc, datax);
   
8. p1 := scatterplot(datax, datay);
   
9. display(p1);
   
10. leastsquare[[x, y], y = a0+a1*x+a2*x^2, {a0, a1, a2}]([datax, datay]);
    
11. f := rhs(%);
    
12. p2 := plot(f, x = 0 .. 100);
    
13. display(p1, p2);
    
14. datax1 := [1, 2, 3];
    
15. datax2 := [2, 3, 4];
    
16. datay1 := [3, 7, 8];
    
17. leastsquare[[x1, x2, y], y = a1*x1+a2*x2, {a1, a2}]([datax1, datax2, datay1]);
    
18. 
    
19. with(Statistics);
    
20. f1 := Fit(a+b*x+c*x^2, datax, datay, x);
    
21. f := Fit(a0+a1*sin((1/10)*x)+a2*sin((1/20)*x)+a3*sin((1/40)*x), datax, datay, x);
    
22. p3 := plot(f, x = 0 .. 100);
    
23. display(p1, p3);
    
24. with(linalg);
    
25. datax1 := [seq((1/100)*(rand(100))(), i = 1 .. 50)];
    
26. datax2 := [seq((1/100)*(rand(100))(), i = 1 .. 50)];
    
27. datax12 := zip(proc (x1, x2) options operator, arrow; [x1, x2] end proc, datax1, datax2);
    
28. datay := zip(proc (x1, x2) options operator, arrow; x1+2*x2+cos(x1*x2)+(1/1000)*(rand(100))() end proc, datax1, datax2);
    
29. whattype(datay);
    
30. nops(datay);
    
31. f4 := Fit(a1*x1+a2*x2+a3*cos(x1*x2)+a4*x1*x2, Matrix(datax12), datay, [x1, x2]);

复制代码