matlab最佳平方逼近.docxVIP

精选文档 精选文档 可编辑 可编辑 数值分析第一次试验 最佳平方逼近试验 任兵(200820302025 ) 、问题叙述 三次最佳平方逼近多项式。求函数 f(x)=exp(x)在［-1 , 1］上的二、 三次最佳平方逼近多项式。 、问题分析 使得 由教材定义6.5有：对于给定的函数 一* ，、 一 ， ，、 S (x) Span{ 0(x), f(x) C[a, b] 1(X),L , n(X)} (x) f (x) S (x) dx min (x) f (x) a a x b a 则称 S*(x)是 f (x)在集合 Span{ 0(x),1(x), L , 显然，求最佳平万谩近函数S (x) n * aj 0 如果存在 2 . s(x) dx n(X)｝中的最佳平方逼近函数。 j(x)的问题可归结为求它的系数 f(x)2 naj j(x) dxj 0an)是关于 f(x) 2 n aj j(x) dx j 0 an)是关于 ao, a1,…，an的二次函数, 利用多元函数取得极值的必要条件, —0 ak (k = 0, 1,2,…，n) Iak b 2 (x) a f(x) n aj j(x) k (x) dx 0 j 0 a；, a*, , an ,使多元函数 b I(a0,a1, , an) a (x) a 取得极小值，也即点(a0,a［， ，2口)是19。，…，an)的极点。由于I (a。，a1，…, 得方程组 naj .j 0 b n aj . j 0 b (x) k(x) j(x)dx (x)f(x) k(x)dx, (k 0,1, 2, , n) 如采用函数内积记号 (k,j)(x) (k, j) (x) k(x) j(x)dx, (f, k) (x)f(x) k(x)dx, 那么，方程组可以简写为 n(k, j 0j)aj(f, n (k, j 0 j)aj (f, k) (k 0,1, 2,L ,n) ...(1) 这是一个包含n + 1个未知元 a。，a1, …,an的n 阶线性代数方程组，写成 矩阵形式为 0, 00,1,1)1)1,n)n)L L L L(n 0, 0 0, 1, 1) 1) 1, n) n) L L L L (n, 0) 1) (n, n) an (f, (f, M (f, ))) n) (2) 此方程组叫做求aj (j = 0, 1,2, …,n)的法方程组。 显然，其系数行列式就是克莱姆行列式 Gn = Gn ( 0, 1 ,…，n)o由于 0, 1,…，n线性无关，故Gn 0,于是上述方程组存在唯一解 ak ak(k 0,1, , n)。从而月 je了函数 f (x)在 Span{ 0(x),1(x),L , n(x)}中如 果存在最佳平方逼近函数，则必是 n _ * , 、 * ,. 一、 S (x) aj j(x) ：...??(3) j 0 三、实验程序 1、最佳平方逼近算法 （1）输入被逼近函数f（x）和对应的逼近区间［a,b］并选择逼近函数系｛a 和权函数； 解方程组（ 1 ）或（ 2 ） ，其中方程组的系数矩阵和右端的项由式（ 得到； 由式（ 3 ）得到函数的最佳平方逼近。 2、 将上述算法编写成 MATLAB 程序共需三个程序： x) } 3） 1 ）第一个程序（函数名： squar_approx.m ） 计算最佳逼近函数的系数，源代码如下： function S=squar_approx(a,b,n) %定义逼近函数 global i;global j; if nargin3 n=1;end Phi2=zeros(n+1); for i=0:n for j=0:n; ％全局变量 ％判断 生成一个 n+1*n+1 大小的全 0矩阵数组 Phi2(i+1,j+1)=quad(@rho_phi,a,b); %求rho_phi 积分 end end PhiF=zeros(n+1,1); ％生成一个 n+1*n 大小的全 0矩阵数组 for i=0:n PhiF(i+1)=quad(@fun_phi,a,b); %求fun_phi 积分 end s=Phi2\PhiF; ( 2 )第二个程序(函数名： rho_phi.m ) 代码如下： function y=rho_phi(x) global i;global j; y=(rho(x).*phi_k(x,i)).*phi_k(x,j); ( 3 )第三个程序(函数名： fun_phi.m ) function y=fun_phi(x) global i; y=(rho(x).*phi-k(x,i)).*obj(x); 四、试验结果 f(x)=exp(x) 在[-1 ， 1] 上的二、三次最佳平方逼近多项式。 (