这篇文章长见识了啊：matlab的二次函数拟合

mytomorrow

& w0 X" I2 @$ N0 P3 l% N6 V世界如此复杂，任一元素受到太多因素的影响，因而要挑出合适的元素。例如，对于无人机价格Price，仅考虑体积V和速度S对价格的影响，构建一个函数P(V,S)。人为选定采用二次函数的方法拟合。

+ r$ j6 p, d  w: m
/ I4 f& t/ f5 _3 r: ?
9 V! u# z1 J6 l1 l收集的数据如表：

' B; x  x4 R7 B

那么就不难列出如下一个方程：
8 T0 Z" ^: M9 n( P; _7 n: u( r
价格=系数阵*数据阵   代码： c=data_new\price

6 W8 o) y1 S  z; C2 e1 \6 C6 G- W, b  q其中系数阵就是不同下标的C，其中数据阵就是【V^2,S^2，V，S，SV，1】
" D' q/ p' n! b% ]+ d( m
欲求系数阵，只需用六或以上方程（正定或超定方程），六个未知数，继而：
. ]/ h+ r, Q: o- E7 ~
5 {5 }4 _& p8 }3 ?: m6 A系数阵=价格*数据阵的逆

" x& v. V: [! W$ g: M8 C' m代码如下：

) Z7 [, f# G2 N

• %输入原始数据
• data=[1,1169,124992,16;2,948,55200,22;3,1051,89280,15;4,825,10200,8;5,921,36400,12;6,873,8500,10;7,1100,80080,19];
• %分别计算价格阵和数据阵
• price=data(:,2);
• data_new=[];
• for n=1:7
•     s=data(n,4);
•     v=data(n,3);
•     temp=[v^2,s^2,v,s,s*v,1];
•     data_new=[data_new;temp];
• end
• %计算结果
• c=data_new\price
• %验算结果
• dis=data_new*c-price
  

   
. T0 d4 \1 u  H2 }6 B3 Z
但其实，你完全可以在工具栏输入cftool唤起拟合工具栏，进行如此的设置：
9 u) W2 J& O+ H/ g
* g" }: |9 {- @; U, ?
* `! Y* y: z* y4 k$ T& }
- V* {( r7 k9 `1 i' |; h一步，就完成了全部。而且他还帮你算了各种我看不懂的系数。多么好啊？
! d! _, I0 i- c5 j- G4 P( M7 k

kinidrily

谢谢分享

dongyijie

好东东，拿来看看 好东东，拿来看看

chencol

谢谢分享