matlab最佳平方逼近[精].docVIP

数值分析第一次试验 最佳平方逼近试验 任 兵（200820302025） 问题叙述 求函数f(x)=exp(x)在[-1，1]上的二、三次最佳平方逼近多项式。 二、问题分析 由教材定义6.5有：对于给定的函数，如果存在 使得 则称S*(x)是f (x)在集合中的最佳平方逼近函数。 显然，求最佳平方逼近函数的问题可归结为求它的系数，使多元函数 取得极小值，也即点()是I (a0, …，an)的极点。由于I (a0, a1, …，an)是关于a0, a1, …，an的二次函数，利用多元函数取得极值的必要条件， (k = 0, 1, 2, …, n) 即 得方程组 如采用函数内积记号 那么，方程组可以简写为 ……………...（1） 这是一个包含n + 1个未知元a0, a1, …, an的n + 1阶线性代数方程组，写成矩阵形式为 …………（2） 此方程组叫做求aj (j = 0, 1, 2, …, n)的法方程组。 显然，其系数行列式就是克莱姆行列式Gn = Gn ((0, (1, …, (n)。由于(0, (1, …, (n线性无关，故Gn ( 0，于是上述方程组存在唯一解。从而肯定了函数f (x)在中如果存在最佳平方逼近函数，则必是 ……………………………...（3） 三、实验程序 最佳平方逼近算法 输入被逼近函数f(x)和对应的逼近区间[a,b]并选择逼近函数系{∮（x）}和权函数； 解方程组（1）或（2），其中方程组的系数矩阵和右端的项由式（3）得到； 由式（3）得到函数的最佳平方逼近。 将上述算法编写成MATLAB程序共需三个程序： （1）第一个程序（函数名：squar_approx.m） 计算最佳逼近函数的系数，源代码如下： function S=squar_approx(a,b,n) ％定义逼近函数 global i;global j; ％全局变量 if nargin3 n=1;end ％判断 Phi2=zeros(n+1); ％生成一个*n+1大小的全0矩阵数组for i=0:n for j=0:n; Phi2(i+1,j+1)=quad(@rho_phi,a,b); ％求rho_phi积分 end end PhiF=zeros(n+1,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]上的二、三次最佳平方逼近多项式。 （1）编写下面三个函数：外部函数；多项式函数；被逼近函数； function y=rho(x) ％外部函数 y=1; function y=phi_k(x,k) ％多项式函数 if k==0 y=ones(size(x)); else y=x.^k; end function y=obj(x) ％被逼近函数 y=exp(x) （2）当求的是二次逼近时得到如下结果 clear S=squar_approx(-1,1,2) S = 0.9963 1.1036 0.5367 绘制两者的图形： fun=exp(x); fplot(fun,[-1,1]) hold on xi=-1:0.1:1; yi=polyval(S,xi); plot(xi,yi,r:) 得到如下结果： （3）当求的是三次逼近时得到如下结果 S=squar_approx(-1,1,3) S = 0.9963 0.9980 0.5367