当前位置:首页 >> 数学 >>

矩阵理论作业6:两种算法求三次最佳平方逼近多项式


两种算法求 f ( x) 的三次最佳平方逼近多项式
摘 要 对于一个较复杂的函数,往往需要求一个简单多项式来逼近。本文选取两种基函数,用两种算 法计算一个函数 f ( x)= exp( x) ? sin( x) 在 x ? [0, ? ] 上的三次最佳平方逼近多项式及其逼近误差, 然后应 2 用 matlab 进行计算作图并对比两种方法的逼近结果。 关键字:三次 最佳平方 逼近 两种算法 引言 多项式的一个重要应用就是可以用来逼近 一个区间上的连续函数,往往许多复杂的函数需 要用各种方法来进行多项式逼近。本文参考矩阵 根据 C[a, b] 上内积定义 理论讲义
[1]

(1)

x? in 对 函 数 f ( x) = e x p ( ) s x 在( )

x ? [0, ] 上进行三次最佳平方逼近,求其逼近 2
误差, 并在 matlab 中编程计算和绘图以验证和比 较逼近结果的准确性。 求最佳平方逼近的多项式

?

( f ( x), g ( x)) ? ? ? ( x) f ( x) g ( x)dx
a

b

(2)

其中权函数

? ? ( x ) =1, 在 x ? [0, ] 上 计 算 得
2

? ? /2 ?2 /8 ? 3 / 24 ? 4 / 64 ? ? a0 ? ? 2.9052 ? ? 2 ?? ? ? ? 3 4 5 ? ? / 8 ? / 24 ? / 64 ? /160 ? ? a1 ? ? ? 3.2781?(3) ?? 3 / 24 ? 4 / 64 ? 5 /160 ? 6 / 384 ? ? a2 ? ? 4.0294 ? ? 4 ?? ? ? ? 5 6 7 ?? / 64 ? /160 ? / 384 ? / 896 ? ? a3 ? ? 5.2035 ?

? f ( x)= exp( x) ? sin( x) ,x ? [0, ] 求三次的 2
最佳平方逼近多项式 (内积中的权函数 ? ( x)=1 ) 。
2 3 用两种算法实现,一是设 ? =span 1, x, x , x ,

解 得 待 定 系 数 a0 ? 0.0201 , a1 ? 0.7658 ,

a2 ? 1.5765 , a2 ? ?0.0708 。即
? ( x) ? 0.0201 ? 0.7658 x ? 1.5765 x 2 ? 0.0708 x3 (4)
误差 ? ( x) = (? ( x), ? ( x))= f ( x) 2 ? (?, f ) ? 0.0100 (5) 2 2 第二种算法:

?

?

二是设为 ? = ? L0 ( x), L1 ( x), L2 ( x), L3 ( x)? ,其中

Li ( x), i ? 0,1, 2,3 是勒让德多项式。
第一种算法:

? = ?L0 ( x), L1 ( x), L2 ( x), L3 ( x)? , 其 中

? =span ?1, x, x , x ? ,由矩阵形式
2 3

Li ( x) ,? i

0 , 1是勒让德多项式, , 2, 3 其表达式为:

Li ( x) ?

1 di [( x 2 ? 1)i ], i ? 0,1, 2,?, n (6) i i 2 i ! dx

[ pi^6/384,

pi^5/160, pi^4/64, pi^3/24, pi^2/8,

pi^3/24] pi^2/8 ] pi/2 ]

[ pi^5/160, pi^4/64, [ pi^4/64, B'= pi^3/24,

勒让德正交多项式有如下递推关系

L0 ( x) ? 1 , L1 ( x) ? x ,
Li ?1 ( x) ? 2i ? 1 i xLi ( x) ? Li ?1 ( x), i ? 1, 2,?, n ? 1? 7 ? i ?1 i ?1

(exp(pi/2)*(pi^3 - 12*pi + 24))/16 (exp(pi/2)*(pi/2 - 1)*(pi/2 + 1))/2 + 1/2 (pi*exp(pi/2))/4 - 1/2 exp(pi/2)/2 + 1/2

可得 L2 ( x) ?

3 2 1 5 3 x ? , L3 ( x) ? x3 ? x 2 2 2 2

求得结果,误差为 0.3020。 算例分析 在 matlab 中编程计算(程序见附录) ,第一 种方法的结果为 >>f=@(x)exp(x).*sin(x); >> [P,error]=polyappro(f,0,pi/2) 得到结果: P= -0.0708 error = 0.0100 1.5765 0.7658 0.0201

第二种方法,勒让德多项式逼近 >> f=inline('exp(x).*sin(x)'); >> [Y,error]=LegendreApproximation(f,4) 得到结果: error = 0.3020

结论 两种方法逼近的效果比较,第一种基函数逼 近的整体效果好,误差小,但是从图上可以看出 局部细节不是很准确;而第二种勒让德多项式逼 近的结果在前半部分非常好,但后边误差较大, 总体误差大。 可以看到 A= [ pi^7/896, pi^6/384, pi^5/160, pi^4/64] 参考文献
[1]矩阵理论讲义(矩阵论第二章),19-20.

附录

方法一:求内积 >> f=inline('sin(x).*exp(x)','x'); >> fy0=quad(f,0,pi/2) fy0 = 2.9052 >> f=inline('sin(x).*exp(x).*x','x'); >> fy1=quad(f,0,pi/2) fy1 = 3.2781 >> f=inline('x.^2.*sin(x).*exp(x)','x'); >> fy2=quad(f,0,pi/2) fy2 = 4.0294 >> f=inline('x.^3.*sin(x).*exp(x)','x'); >> fy3=quad(f,0,pi/2) fy3 = 5.2035 最佳平方逼近函数: function [poly,error]=polyappro(f,a,b) syms x P=[1,x,x.^2,x.^3]; %P=[1,x,1.5*x.^2-0.5,2.5*x.^3-1.5]; m=length(P); for i=1:m for j=1:m A(i,j)=int(P(m-i+1).*P(m-j+1),a,b); end B(i)=int(f*P(m-i+1),a,b); end poly=A\B'; poly=sym2poly(poly'); xx=a:0.01:b; yy=feval(f,xx); ypoly=polyval(poly,xx); polysym=poly2sym(poly); delta=f-polysym; error=sqrt(int(delta.*delta,a,b)); error=sym2poly(error); plot(xx,yy,xx,ypoly,'r-') end

fval=feval(f,xx); for i=1:n+1 w(i)=(2*(i-1)+1)/2*quad(@(x)f(x).* mfun('p',i-1,x),-1,1); YY(i,:)=w(i).*mfun('p',i-1,xx); Y=sum(YY); error=norm(fval-Y,'inf'); plot(xx,fval,xx,Y,'r'); end end

方法二:勒让德多项式 Function [Y,error]=LegendreApproximation(f,n) xx=0:0.01:0.5*pi;


相关文章:
矩阵理论作业6:两种算法求三次最佳平方逼近多项式.doc
矩阵理论作业6:两种算法求三次最佳平方逼近多项式_数学_自然科学_专业资料。两种
计算方法各章习题及答案.doc
? x p ( x) -2 3 -1 1 0 1 1 6 2 15...x2 的三次最佳平方逼近 多项式 S3 ( x) ,并...选主元的高斯-若当消去法求矩阵 A 的逆矩阵,其中...
数值分析第6章习题.doc
数值分析第章整合版(黑组) 一、填空题 1 ? x? ? x , 1、已知 P0 ?...x3 求[-1,1]上关于 ? (x)=1 的最佳平方逼近次多项式 2 解: Pk ( ...
中科院数值分析第一次作业及参考答案1 (2014版).doc
中科院数值分析第一次作业及参考答案1 (2014版)_工...10?6 得 h ? 0.006. 9 3 3. 求 f ( x)...e 在[-1,1]上的三次最佳平方逼近多项式。 (参考...
数值分析第一次作业(2).doc
数值分析第一次作业(2)_理学_高等教育_教育专区。...ex 在[-1,1]上的三次最佳平方逼近多项式。 (...6.1 米,为控制丈量的准确性,又测量 AC ? x1 ?...
最佳平方逼近多项式.doc
最佳平方逼近的算法,能够根据给定的函数值表达求出二、三次最佳平 方逼近多项式...矩阵理论作业6:两种算法... 3页 2下载券 喜欢此文档的还喜欢 2.4...
《数值逼近》蒋尔雄第二版课后习题答案.pdf
再用 6 位数字计算出 x 与 y,已知正确解为练习...14.假设有一种算法,求 数字? 解:因为 可得到 6...: 18.求 在[0,1]上的一次最佳平方逼近多项式。 ...
数值计算方法第七章习题 2013.doc
(x) = 1 的最佳平方逼近次多项式。 20.求函数...3 0.5 1.117 4 0.6 1.223 5 0.75 1.422 xi...2011数值分析第一次作业... 13页 2下载券 ...
16春天大《数值计算方法》在线作业二.doc
16春天大《数值计算方法》在线作业二_远程、网络教育...数值稳定的算法是指舍入误差对计算结果影响不大的...[-1,1] 上的三次最佳平方逼近多项式为 0.9963+...
计算方法 第五章第二节最佳平方逼近_图文.ppt
由性质5.2.3,该方程组 的系数矩阵非奇异,因而...0.6644389 x. 例5.2.2 求函数 f ( x) ? sin ? x 在[0,1] 上的最佳平方逼近次多项式...
综合习题.pdf
7、利用Remez算法,计算函数 ,在区间[0,1] 上的...) 12、求在[-1,1]上的三次最佳平方逼近多项式。...(参考讲义与参考书) 习题二 1、用矩阵的直接三角...
Matlab第二、三次上机作业.doc
Matlab第二、三次上机作业 - 第二次上机作业 一. 任务: 用 MATLAB 语言编写连续函数最佳平方逼近算法程序(函数式 M 文件) 。并用此程序进 行数值试验,写出...
第三章-2-最佳平方逼近_图文.ppt
函数逼近与计算 3.3 最佳平方逼近 及正交多项式 ?...--Gram 矩阵 则 G 非奇异当且仅当 u1, u2, ?...(6) Tn(x) 的首项系数为 2n-1,且 |Tn(x)|...
MATLAB上机实验最佳平方逼近.doc
中国石油大学(北京)数值分析课程上机作业 一. 任务:...6. 程序输入: (1)待求的被逼近函数值的数据点 ...多项式作基的函数最佳平方逼近算法程序 LEASTP.m (...
2012数值分析第一次作业及参考答案.doc
0 .4 5 1 5 6 5( y ? 2 .2 3)( y ? 1 .1 0 ) ? 0 .2 5...e 在[-1,1]上的三次最佳平方逼近多项式。 (参考讲义与参考书,利用 x ...
最佳平方逼近算例.pdf
最佳平方逼近算例 - 例:求函数 f ( x) = e x 在[0,1]上的二次最佳平方逼近多项式,并估计平 ,小数点后保留 5 位. 方逼近误差 δ 2 2 解: (解法 ...
数值分析复习题(上).doc
三次插值多项式 P3(x)的 x3 的系数是 6,试确...x 4 } 的最佳平方逼近多项式 23.在科学计算中, ...试解释这种做法的理论依据。 ? x1 ? x 2 ? 1 ...
数值分析习题.doc
多项式, 并近似计算 cos 日二次插值) 6 已知函数...主要考察点:最小二乘法,最佳平方逼近,正交多项式的...迭代法和一般迭代法的收敛性讨论) 6 给出矩阵 A ...
计算方法复习题.doc
6.记 h? ,Newton 形式的二次插值多项式为 N 2 ...正定矩阵,则求解方程组 Ax = b 的最速下降算法为...( x)d x ,则函数 x 2 的最佳平方逼近一次...
数值分析第2次作业.doc
(-pi,pi) %计算数值 AA3=subs(A3,'t',x); ...(coff,6); 勒让德多项式 n 次最佳平方逼近 ...(f,6); 根据数值分析 P90 改进的 FFT 算法 ...
更多相关标签: