弹簧问题(拟合的妙用).docVIP

内容：进一步学习最小二乘原理及其Matlab实现 一、复习 二、问题 X x1 x2 x3 x4 x5 x6 x7 x8 x9 Y y1 y2 y3 y4 y5 y6 y7 y8 y9 x,y分别表示弹簧弹力试验中的拉伸长度和弹力，由物理知识知在弹性限度内，y=kx，在弹性限度外，y=ax^2+bx+c，图中显示前五个点大致呈直线，后五个点大致呈抛物线。 要求拟合出这两段曲线的函数表达式，满足 1、在x5处连接（即连续） 2、在x5处光滑（理解为左导数等于右导数即可导） 二、思路： 1、难点：若仅按刚讲过的原理和方法，可以拟合出两段曲线所代表的函数，即解超定方程组，求出最小二乘解： 对前五个点方程组为 进而，X=k=A\b 对于后五个点方程组为 进而，X=A\b求出abc 但是，无论从结果还是从原理分析，以上方法拟合的结果都无法保证曲线经过某个点，只能保证总体上距离平方和最小，这样就不可能满足1和2 2、办法：怎么办呢？拿到这个问题，我首先想到，可以用条件极值（或叫到约束的优化）来做，即在满足条件1和2的前提下，使J(a,b,c)最小。但是，由于这一部分至今没有讲，暂时不能用。其次考虑插值，两点原因注定不能用（1）测量点不准确，不适合用插值（2）插值不能产生满足要求的函数表达式。我们能用的方法目前只有再次考虑最小二乘拟合 原理已经复习过了，还是使得该超定方程组的两端差距