多项式插值的振荡现象
实验内容及要求:
设区间[-1,1]上函数 f(x)=1/(1+25x2)。
考虑区间[-1,1]的一个等距划分,分点为xi= -1 + 2i/n,i=0,1,2,…,n,
则拉格朗日插值多项式为
.
其中,li(x),i=0,1,2,…,n是n次Lagrange插值基函数。
function t_charpt2
% 数值实验二:含“实验2.1:多项式插值的震荡现象”和“实验2.2:样条插值的收敛”
% 输入:实验选择,函数式选择,插值结点数
% 输出:拟合函数及原函数的图形
result = inputdlg({
'请选择实验,若选2.1,请输入1,否则输入2:'},'charpt_2', 1,{
'1'});
Nb = str2num(char(result));
if(Nb ~= 1) & (Nb ~=2) errordlg('实验选择错误!');return;end
promps = {
'请选择实验函数,若选f(x),请输入f,若选h(x),请输入h,若选g(x),请输入g:'};
titles = 'charpt_2';
result = inputdlg(promps,'charpt 2',1,{
'f'});
Nb_f = ch