zoukankan      html  css  js  c++  java
  • 二次函数三点式

    拉格朗日插值法:

    已知二次函数过三个点 $(x1,y1)$,$(x2,y2)$,$(x3,y3)$
    求函数的解析式。

    你还在 高斯消元 ?!

    二次函数的三点式表示法,让你远离$EPS$蒙不对的烦恼。

     $f(x)$ $=$ $frac{(x-x_{2})(x-x_{3})}{(x_{1}-x_{2})(x_{1}-x_{3})}*y_{1}$ $+$ $frac{(x-x_{1})(x-x_{3})}{(x_{2}-x_{1})(x_{2}-x_{3})}*y_{2}$ $+$ $frac{(x-x_{1})(x-x_{2})}{(x_{3}-x_{1})(x_{3}-x_{2})}*y_{3}$

     $f(x)$ $=$ $frac{x^{2}-(x_{2}+x_{3})x+x_{2}x_{3}}{(x_{1}-x_{2})(x_{1}-x_{3})}*y_{1}$ $+$ $frac{x^{2}-(x_{1}+x_{3})x+x_{1}x_{3}}{(x_{2}-x_{1})(x_{2}-x_{3})}*y_{2}$ $+$ $frac{x^{2}-(x1+x2)x+x_{1}x_{2}}{(x_{3}-x_{1})(x_{3}-x_{2})}*y_{3}$

    我们发现,分母通分的话,拆分后在以$x^{n}$为关键字合并同类项可得。

    ps:懒得打LATEX,所以部分过程省略。

     $(x_{1}-x_{2})(x_{1}-x_{3})(x_{2}-x_{3})f(x)$ $=$

    $[x_{1}(y_{3}-y_{2})+x_{2}(y_{1}-y_{3})+x_{3}(y_{2}-y{1})]*x^{2}$ $+$

    $[x_{1}^{2}(y_{2}-y_{3})+x_{2}^{2}(y_{3}-y_{1})+x_{3}^{2}(y_{1}-y{2})]*x$ $+$

    $x_{2}x_{3}y_{1}*(x_{2}-x_{3})+x_{1}x_{3}y_{2}*(x_{3}-x_{1})+x_{1}x_{2}y_{3}*(x_{1}-x_{2})$

    这样的话,各个系数就直接出来了嘛。。。如果我们让 $x_{1}>=x_{2}>=x_{3}$ 的话就更好了。

    二次项系数:
     $frac{[x_{1}(y_{3}-y_{2})+x_{2}(y_{1}-y_{3})+x_{3}(y_{2}-y{1})]}{(x_{1}-x_{2})(x_{1}-x_{3})(x_{2}-x_{3})}$

    一次项系数:
     $frac{[x_{1}^{2}(y_{2}-y_{3})+x_{2}^{2}(y_{3}-y_{1})+x_{3}^{2}(y_{1}-y{2})]}{(x_{1}-x_{2})(x_{1}-x_{3})(x_{2}-x_{3})}$

    常数项:
     $frac{x_{2}x_{3}y_{1}*(x_{2}-x_{3})+x_{1}x_{3}y_{2}*(x_{3}-x_{1})+x_{1}x_{2}y_{3}*(x_{1}-x_{2})}{(x_{1}-x_{2})(x_{1}-x_{3})(x_{2}-x_{3})}$

    好了。就解到这里吧。剩下的东西就交给$DEVC++$来做了。

  • 相关阅读:
    声明
    Random类——获取随机数
    JDK_API的使用方法
    Ajax?有谁开始学习了吗?
    用xslt循环xml同一节点的不同子节点
    在Repeater控件中嵌套一个Repeater控件
    html+css的一些技巧.收集中...
    记录一下: onbeforeunload()方法, 避免未保存而关闭页面.
    简单的C# Socket编程
    不实用的UriBuilder类
  • 原文地址:https://www.cnblogs.com/arcturus/p/9180349.html
Copyright © 2011-2022 走看看