- restart:
- Cd := Vector([.686, .693, .702, .714, .728, .745, .764,
- .785, .594, .597, .601, .607, .614, .621, .63, .641,
- .545, .546, .549, .553, .557, .562, .566, .574, .514,
- .515, .518, .52, .523, .526, .53, .535, .494, .495,
- .496, .498, .5, .503, .504, .505, .479, .48, .481, .482,
- .484, .485, .488, .49]):
- S := Vector(10^3*[.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, .5, 1,
- 1.5, 2, 2.5, 3, 3.5, 4, .5, 1, 1.5, 2, 2.5, 3, 3.5, 4,
- .5, 1, 1.5, 2, 2.5, 3, 3.5, 4, .5, 1, 1.5, 2, 2.5, 3,
- 3.5, 4, .5, 1, 1.5, 2, 2.5, 3, 3.5, 4]):
- V := Vector([40, 40, 40, 40, 40, 40, 40, 40, 60, 60,
- 60, 60, 60, 60, 60, 60, 80, 80, 80, 80, 80, 80, 80, 80,
- 100, 100, 100, 100, 100, 100, 100, 100, 120, 120, 120,
- 120, 120, 120, 120, 120, 140, 140, 140, 140, 140, 140,
- 140, 140]):
- W := Statistics:-Fit(a+b*s+c*v+d*s^2+e*s*v+f*v^2,
- Matrix([S, V]), Cd, [s, v]);
复制代码
数据是随便取的,拟合问题可以这样试试 |