一
.
任务:
用
MATLAB
语言编写连续函数最佳平方逼近的算法程序(函数式
M
文件)
。并用此
程序进行数值试验,写出计算实习报告。
二
.
程序功能要求:
在书中
Page355
的程序
leastp.m
基础上进行修改,使其更加完善。要求算法程序可以适应
不同的具体函数,具有一定的通用性。所编程序具有以下功能:
1.
用
Lengendre
多项式做基,并适合于构造任意次数的最佳平方逼近多项式。
可利用递推关系
0
1
1
2
(
)
1,
(
)
(
)
(2
1)
(
)
(
1)
(
)
/
2,
3,.....
n
n
n
P
x
P
x
x
P
x
n
xP
x
n
P
x
n
n
2.
被逼近函数
f(x)
不用内联函数构造,而改用
M
文件建立数学函数。这样,此程序可通过修
改建立数学函数的
M
文以适用不同的被逼近函数(要学会用函数句柄)
。
3.
要考虑一般的情况
]
1
,
1
[
]
,
[
)
(
b
a
x
f
。因此,程序中要有变量代换的功能。
4.
计算组合系数时,计算函数的积分采用
5
点
Gauss-Legendre
求积公式。
5.
程序中应包括帮助文本和必要的注释语句。另外,程序中也要有必要的反馈信息。
6.
程序输入:
(
1
)待求的被逼近函数值的数据点
0
x
(
可以是一个数值或向量
)
(
2
)区间端点:
a,b
。
7.
程序输出:
(
1
)拟合系数:
0
1
2
,
,
,...,
n
c
c
c
c
(
2
)待求的被逼近函数值
0
0
0
0
1
1
0
2
2
0
0
(
)
(
)
(
)
(
)
(
)
n
n
s
x
c
P
x
c
P
x
c
P
x
c
P
x