数值逼近期末复习笔记

免责声明

仅可作为期末突击复习使用,以下过程可能存在不严谨不严格的部分,但是本科期末考试一般会给分

拉格朗日插值

基本公式

会给定一系列点 (x0,f(x0)),(x1,f(x1)),,(xn,f(xn)) ,通过下面公式
$$
L_i(x) = \frac{(x-x_0)(x-x1)\cdots(x-x{i-1})(x-x_{i+1})\cdots(x-x_n)}{(x_i-x_0)(x_i-x_1)\cdots(xi-x{i-1})(xi-x{i+1})\cdots(x_i-x_n)}
$$

Pn(x)=y0L0(x)+y1L1(x)++ynLn(x)
Pn(x) 即为所求的拉格朗日多项式

例题

已知 f(1)=2,f(1)=1,f(2)=1,求 f(x) 的 Lagrange 插值多项式。
L0(x)=(x1)(x2)23=16(x23x+2)

L1(x)=(x+1)(x2)21=12(x2x2)

L2(x)=(x+1)(x1)31=13(x21)

Pn(x)=2x16(x23x+2)12(x2x2)+13(x21)

=13x212x+2312x+13x216

=18x22x+43

牛顿插值

基本公式

会给定一系列点(x0,f(x0)),(x1,f(x1)),,(xn,f(xn)),通过下面公式
差商计算公式

f[xi,xj]=f(xi)f(xj)xixj

$$
f[x_0,x_1,\ldots,x_j] = \frac{f[x_0, x1, \ldots, x{k-1}] - f[x_0, \ldots, x_k]}{x_0 - x_k}
$$

\begin{array}{|c|c|c|} \hline x_0 & f(x_0) & & \ \hline x_1 & f(x_1) & f[x_0, x_1] & \ \hline x_2 & f(x_2) & f[x_1, x_2] & f[x_0, x_1,x_2] \ \hline x_3 & f(x_3) & f[x_2, x_3] & f[x_1, x_2,x_3] & \ldots \ \hline \ldots & \ldots & \ldots & \ldots \ \hline \end{array}
Pn(x)=f(x0)+(xx0)f[x0,x1]+(xx0)(xx1)f[x0,x1,x2]+$$+ (x-x_0)(x-x1) \cdots (x-x{n-1})f[x_0, x_1, \ldots, x_n]
$$
Pn(x)为所求的插值多项式

例题

已知 f(1)=2,f(1)=1,f(2)=1,求 f(x) 的 Lagrange 插值多项式。
给出差商表

\begin{array}{|c|c|c|} \hline -1 & 2 \ \hline 1 & 1 & -\frac{1}{2} \ \hline 2 & 1 & 0 & \frac{1}{6} \ \hline \end{array}

Pn(x)=2+(x+1)(12)+(x+1)(x1)18

=212x12+(x1)18

=212x1218x218

=18x212x+43

Hermite插值

两点三次插值公式

H3(x)=(12xx0x0x1)(xx1x0x1)2y0+(12xx1x1x0)(xx0x1x0)2y1+(xx0)(xx1x0x1)2y0+(xx1)(xx0x1x0)2y1

例题

求满足条件 H(0)=1, H(0)=0.5, H(1)=2, H(1)=2 的插值多项式。
解:
直接带入公式

  • x0=0, y0=1, y0=0.5
  • x1=1, y1=2, y1=2
    将这些值代入 Hermite 插值公式:

H3(x)=(12x001)(x101)2(1)+(12x110)(x010)2(2)+(x0)(x101)2(0.5)+(x1)(x010)2(2)
最终整理得到:

H3(x)=2x33x2+2x+1

三次样条

这个学的不好
挖坑

一次最佳一致逼近多项式

基本方法

f(x) 得出 f(x)
a,b 为给定区间两端
然后 f(b)f(a)ba=a1=f(x3) ,求解出 x3 的值,代入
a0=f(a)+f(x3)2f(b)f(a)baa+x32
求解的最佳一致逼近多项式为 P(x)=a0+a1x

例题

f(x)=x2 , x[0,1] 上的一次最佳一致逼近多项式。
解:
计算 f(x):
f(x)=2x
在这个例子中,我们取 a=0b=1
计算 f(a), f(b), 以及 a1:
f(0)=02=0,f(1)=12=1
代入公式
a1=f(x3)=f(b)f(a)ba=1010=1
解方程 f(x3)=1:
2x3=1x3=12
代入公式计算 a0
a0=f(a)+f(x3)2f(b)f(a)baa+x32
接下来代入各项:
a0=f(0)+f(12)2f(1)f(0)100+122
计算各部分:
a0=0+1421122
a0=14214
a0=1814=1828=18

现在我们得到了a0a1:
a0=18,a1=1

综上所述, f(x)=x2 在区间[0,1]上的一次最佳一致逼近多项式为:
P(x)=x18

最佳平方逼近

选取 φ(x)=1, ϕ0(x)=1, ϕ1(x)=|x|xspan1,x

{(ϕ0,ϕ0)c0+(ϕ0,ϕ1)c1=(f,ϕ0) (ϕ1,ϕ0)c0+(ϕ1,ϕ1)c1=(f,ϕ1)

从而
{c0= c1=

其中 (ϕ0,ϕ1)=abϕ0(x)ϕ1(x)dx

φ(x)=c0+c1x

版权声明:除特殊说明,博客文章均为栋dong原创,依据CC BY-SA 4.0许可证进行授权,转载请附上出处链接及本声明。
如有需要,请在留言板留言,或者添加我的QQ或者微信
我只是一个学生,如有错误或者侵权,请联系我,谢!
暂无评论

发送评论 编辑评论


				
|´・ω・)ノ
ヾ(≧∇≦*)ゝ
(☆ω☆)
(╯‵□′)╯︵┴─┴
 ̄﹃ ̄
(/ω\)
∠( ᐛ 」∠)_
(๑•̀ㅁ•́ฅ)
→_→
୧(๑•̀⌄•́๑)૭
٩(ˊᗜˋ*)و
(ノ°ο°)ノ
(´இ皿இ`)
⌇●﹏●⌇
(ฅ´ω`ฅ)
(╯°A°)╯︵○○○
φ( ̄∇ ̄o)
ヾ(´・ ・`。)ノ"
( ง ᵒ̌皿ᵒ̌)ง⁼³₌₃
(ó﹏ò。)
Σ(っ °Д °;)っ
( ,,´・ω・)ノ"(´っω・`。)
╮(╯▽╰)╭
o(*////▽////*)q
>﹏<
( ๑´•ω•) "(ㆆᴗㆆ)
😂
😀
😅
😊
🙂
🙃
😌
😍
😘
😜
😝
😏
😒
🙄
😳
😡
😔
😫
😱
😭
💩
👻
🙌
🖕
👍
👫
👬
👭
🌚
🌝
🙈
💊
😶
🙏
🍦
🍉
😣
Source: github.com/k4yt3x/flowerhd
颜文字
Emoji
小恐龙
花!
上一篇
下一篇