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

两种算法求的三次最正确平方逼近多项式

摘要

对于一个较复杂的函数，往往需要求一个简单多项式来逼近。本文选取两种基函数，用两种算法计算一个函数在上的三次最正确平方逼近多项式及其逼近误差，然后应用matlab进行计算作图并比照两种方法的逼近结果。

关键字：三次最正确平方逼近两种算法

引言

多项式的一个重要应用就是可以用来逼近一个区间上的连续函数，往往许多复杂的函数需要用各种方法来进行多项式逼近。本文参考矩阵理论讲义对函数在上进行三次最正确平方逼近，求其逼近误差，并在matlab中编程计算和绘图以验证和比拟逼近结果的准确性。

求最正确平方逼近的多项式

，求三次的最正确平方逼近多项式〔内积中的权函数〕。用两种算法实现，一是设，二是设为，其中是勒让德多项式。

第一种算法：

，由矩阵形式

(1)根据上内积定义

(2)

其中权函数，在上计算得解得待定系数，，，。即

(4)

误差(5)

第二种算法：

，其中是勒让德多项式，其表达式为：

(6)

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

可得，

求得结果，误差为0.3020。

算例分析

在matlab中编程计算〔程序见附录〕，第一种方法的结果为

f=@(x)exp(x).*sin(x);

[P,error]=polyappro(f,0,pi/2)

得到结果：

P=

-0.07081.57650.76580.0201

er