最佳平方逼近方法.docVIP

2016-2017（1）专业课程实践论文 用最佳平方逼近法求逼近函数 肖夏， 29，R数学12-1班 一、 算法理论 设函数组φ0,φ1,…，φ H=Span{ 则H中的任意一个元素为 p 对空间C[a,b]的任意两个函数f，g，定义内积 f,g 对于给定的函数f(x)∈C[a,b]， f- 则称p*x为子空间H中对于f 特别地，若φjx=xj，j=0,1,…m则称满足条件的p*x 设f(x)∈C[a,b]，p*x∈H是子空间H中对于f(x 求最佳平方逼近元素p*x= 因 f- 得 k=0 得 φ 求出ck*，带入p 二、算法框图 开始 开始 定义权函数ω( 定义权函数ω(x) 输入a 输入a，b，m 否m≥0 否 m≥0 是 for k=1:n+1 for k=1:n+1 for j=1:n+1结束 for j=1:n+1 结束 lk,j= l 结束结束 结束 结束 yk y 结束 c(i) c(i) 结束p 结束 p 三、算法程序 function S=abc(n,a,b) //创建一个函数，里面填入次数，和区间范围 base=inline(x^(j-1),x,j);///定义 quan=inline(1,x); for k=1:(n+1) for j=1:(n+1) syms x l(k,j)=int(base(x,k)*base(x,j)*quan(x),x,a,b); end y(k)=int(base(x,k)*(sqrt(x^2+1)),x,a,b);//红色字体是f end l; y; c=vpa(inv(l)*y,3) p=0; for i=1:(n+1) p=p+c(i)*base(x,i); end p 四、算法实现 例1.求fx=x2 解： f, f, 由方程组 1 c0=0.934，c 第一题的解： 例2. 求fx=sinx在 解： f, f, 由方程组 1 c0=0.036，c 第二题的解： 例3. 求fx=arctan?x在0 解： f, f, f, 由方程组 1 c0=-0.005，c 第三题的解： 例4.求fx=e2x在0 解： f, f, 由方程组 1 c0=0.195，c 第三题的解：