逼近与机械求积题目选做

6. 求 f(x) = sin(x) 在 [0,π/2] 上的最佳一次逼近多项式，并估计误差.

背景知识
p(x) 是 f(x) ∈ C[a,b] 的 n 次最佳一致逼近多项式的充要条件是 p(x) 在 [a,b] 上至少有 n+2 个轮流为"正","负"的偏差点. 这样的点组称为Chebyshev交错点组.

本题中
f(x) = sin(x) 在 [0,π/2] 上存在无限高阶导数, 而且是单调递增的凹函数.

设 p(x) = ax+b, 则 δ(x) = f(x)-p(x) = sin(x)-ax-b 也是凹函数.
设Chebyshev交错点组为 p < q < r, 则 p=0, δ’(q)=0, r=π/2.
即 -b = -[sin(q)-aq-b] = 1-aπ/2-b, cos(q)-a = 0.
解得

$a=\frac{2}{\pi }$, $b=\frac{\sqrt{{\pi }^2-4}}{2\pi }-\frac{1}{\pi }\arccos \frac{2}{\pi }$.
$p(x)=\frac{2}{\pi }x+\frac{\sqrt{{\pi }^2-4}}{2\pi }-\frac{1}{\pi }\arccos \frac{2}{\pi }$
$\max |f(x)-p(x)|=|f(0)-p(0)|=|b|\approx 0.1053$

sqrt(pi.^2-4)./(2.*pi) - acos(2./pi)./pi

from math import sqrt, acos, pi
print(sqrt(pi**2-4)/(2*pi) - acos(2/pi)/pi)

17. 设 φ1=span{1,x}, φ2=span{x¹⁰⁰,¹⁰¹}, 分别在φ1,φ2上求一元素, 使其为x²∈C[0,1]的最佳平方逼近, 并比较其结果.

背景知识
对于基函数 e1(x), … , ed(x), f(x)∈C[a,b] 的最佳平方逼近 g(x) = k1 e1(x) + … + kd ed(x), 满足法方程:
(e1, f) = (e1, g) = (e1, Σ kj ej) = Σ kj (e1, ej)
  ⋮
(ei, f) = (ei, g) = (ei, Σ kj ej) = Σ kj (ei, ej)
  ⋮
(ed, f) = (ed, g) = (ed, Σ kj ej) = Σ kj (ed, ej)

φ1上

1/3 = 1 k1 + 1/2 k2
1/4 = 1/2 k1 + 1/3 k2

k1 = -1/6, k2 = 1

$g(x)=-\frac{1}{6}+x$
$\underset{0}{\overset{1}{\int }}(x^2-g(x){)}^2\mathrm{d}x=0.0056$

φ2上

1/103 = 1/201 k1 + 1/202 k2
1/104 = 1/202 k1 + 1/203 k2

k1 = 2009799/5356, k2 = -1004647/2678
k1 ≈ 375.2425, k2 ≈ -375.1482

$g(x)=(375.2425-375.1482x)x^{100}$
$\underset{0}{\overset{1}{\int }}(x^2-g(x){)}^2\mathrm{d}x=0.1641$

A=[1 1/2;1/2 1/3]; B=[1/3;1/4];
A=sym(A); B=sym(B);
A\B

A=[1/201 1/202;1/202 1/203]; B=[1/103;1/104];
A=sym(A); B=sym(B);
A\B
double(A\B)

f = @(x) x.^2;
g1 = @(x) -1/6 + x;
g2 = @(x) (375.2425 - 375.1482.*x) .* x.^100;
h1 = @(x) (f(x) - g1(x)).^2;
h2 = @(x) (f(x) - g2(x)).^2;

integral(h1,0,1)
integral(h2,0,1)

20. 将 f(x) = sin(x/2) 在 [-1,1] 上按Legendre多项式及Chebyshev多项式展开, 求三次最佳平方逼近多项式并画出误差图形, 再计算均方误差.

背景知识
Legendre多项式: 当区间为[-1,1], 权函数ρ(x)=1时, 由{1,x,x², … ,xⁿ}正交化得到的多项式$L_n(x)=\frac{1}{2^{n}n!}\frac{{\mathrm{d}}^n}{\mathrm{d}{x}^n}(x^2-1{)}^n$.

$L_0(x)=1$
$L_1(x)=x$
$L_2(x)=\frac{3x^2-1}{2}$
$L_3(x)=\frac{5x^3-3x}{2}$

$k_0=\frac{(f,L_0)}{(L_0,L_0)}=0$
$k_1=\frac{(f,L_1)}{(L_1,L_1)}=12\sin \left(\frac{1}{2}\right)-6\cos \left(\frac{1}{2}\right)\approx 0.4876$
$k_2=\frac{(f,L_2)}{(L_2,L_2)}=0$
$k_3=\frac{(f,L_3)}{(L_3,L_3)}=826\cos \left(\frac{1}{2}\right)-1512\sin \left(\frac{1}{2}\right)\approx -0.0082$

$g(x)=\sum _0^3{k}_i{L}_i(x)=\left(2065\cos \left(\frac{1}{2}\right)-3780\sin \left(\frac{1}{2}\right)\right)x^3+\left(2280\sin \left(\frac{1}{2}\right)-1245\cos \left(\frac{1}{2}\right)\right)x\approx 0.5x-0.02x^3$

$\Vert \delta \Vert =\sqrt{\underset{-1}{\overset{1}{\int }}(f(x)-g(x){)}^2\mathrm{d}x}=\sqrt{(f,f)-(f,g)}=\sqrt{229160\cos \left(1\right)+356879\sin \left(1\right)-424119}\approx 1.3966\times 10^{-5}$

x = sym('x');
f = sin(x/2);

L0 = sym_legendre(0);
L1 = sym_legendre(1);
L2 = sym_legendre(2);
L3 = sym_legendre(3);

L0 = collect(L0,x)
L1 = collect(L1,x)
L2 = collect(L2,x)
L3 = collect(L3,x)

k0 = int(L0*f,-1,1) / int(L0*L0,-1,1)
k1 = int(L1*f,-1,1) / int(L1*L1,-1,1)
k2 = int(L2*f,-1,1) / int(L2*L2,-1,1)
k3 = int(L3*f,-1,1) / int(L3*L3,-1,1)

double(k1)
double(k3)

g = k0.*L0 + k1.*L1 + k2.*L2 + k3.*L3;
g = collect(g, x)
vpa(g)

fplot(f-g,[-1,1])
title('mean square error')

err2 = sqrt( int((f-g).^2,-1,1) )
double(err2)

err2 = sqrt( int(f.^2,-1,1) - int(f.*g,-1,1) )
double(err2)

function l = sym_legendre(n)
    x = sym('x');
    f = @(x) (x.^2-1).^n;
    l = diff(f(x), n);
    l = l / (factorial(n) .* 2.^n);
end

22. 用最小二乘法求一个形如 y=a+bx² 的经验公式, 使它数据相拟合, 并求均方误差.

19	25	31	38	44
19.0	32.3	49.0	73.3	97.8

观察发现曲线具有凸性, 因此考虑二次多项式和指数函数. 按照题目要求考虑 y=a+bx².

φ=span{1,x²}

(f,e1) = (e1,e1) k1 + (e1,e2) k2
(f,e2) = (e1,e2) k1 + (e2,e2) k2

271.4 = 5 k1 + 5327 k2
369321.5 = 5327 k1 + 7277699 k2

y ≈ g(x) = 0.9726 + 0.0500 x²

$\Vert \delta {\Vert }_2=\sqrt{\sum _i(Y_i-g(X_i){)}^2}=\sqrt{0.0150}=0.1226$

X = [ 19   25   31   38   44 ];
Y = [19.0 32.3