指定公式的MATLAB拟合

haobang008

有感于MATLAB的强大，此处记下多元回归的指定公式拟合代码，分为两个：

文件1：

1. 
   
2. <P> clc;clear all;
   
3. hc = load('HC_position.txt');   
   
4. x = hc(:,3:2:19);
   
5. y = load('mean_tide.txt');  </P>
   
6. <P>% 拟合选项一
   
7. % 根据指定的公式，进行系数的拟合
   
8. c0(9) = 0;
   
9. c = lsqcurvefit('myfun',c0,x,y)</P>
   
10. <P>% 拟合选项二
    
11. % 直接进行线性拟合</P>
    
12. <P> </P>

复制代码

文件二：

1. function val=myfun(c,x)
   
2. val = c(1)*x(:,1)+c(2)*x(:,2).^2+c(3)*x(3).^2   ...,
   
3.      +c(4)*x(:,4).^2+c(5)*x(:,5).^2+c(6)*x(:,6).^2   ...,
   
4.      +c(7)*x(:,7).^2+c(8)*x(:,8).^2+c(9)*x(:,9).^2;
   
5. end
   
6. 
   

复制代码

haobang008

lsqcurvefit 这个事最关键的万能函数

做梦的蚊子

额，没太看明白楼主这个函数是怎么拟合的。myfun是指x,y之间的函数关系式？得出结果是c1~c9这么9个值？那么文件一中的c0(9)=0的设定是用来干嘛的，进行拟合时用的c0是指的什么？能请楼主给解答下吗？谢谢啦~~~

navy_yehaijun

真是好东西，再次琢磨

haobang008

做梦的蚊子 发表于 2014-4-24 20:30
额，没太看明白楼主这个函数是怎么拟合的。myfun是指x,y之间的函数关系式？得出结果是c1~c9这么9个值？那么 ...

myfun是自己写的函数的形式，的确是指x,y之间的函数关系式，得出的结果就是需要拟合的系数c1~c9，c0(9)=0是用来设置初始的迭代值，因为matlab在进行拟合的时候是通过迭代进行计算的，需要设置初始值，一般设置为0，如果你有比较好的备选参数，也可以换掉。
另外，matlab拟合的基本原理因方法而已，不要拘泥于某种思路，因地制宜。

做梦的蚊子

haobang008 发表于 2014-5-5 07:48
myfun是自己写的函数的形式，的确是指x,y之间的函数关系式，得出的结果就是需要拟合的系数c1~c9，c0(9)=0 ...

好的，谢谢楼主~~~

haobang008

做梦的蚊子 发表于 2014-5-7 21:47
好的，谢谢楼主~~~

客气了，没啥的

daweilee1989

好东西 谢谢分享  

daweilee1989

好东西 谢谢分享  

longyongjing

谢谢楼主，学习了