数值分析 插值法 计算实习题求插值
来源:学生作业帮 编辑:大师作文网作业帮 分类:数学作业 时间:2024/11/20 00:58:28
数值分析 插值法 计算实习题求插值
下列数据点的插值
x 0 1 4 9 16 25 36 49 64
y 0 1 2 3 4 5 6 7 8
可以得到平方根函数的近似,在区间[0,64]上作图.
(1)用这九个点作8次多项式插值L8(x)
(2)用三次样条(第一边界条件)程序求S(x)
从得到结果看在[0,64]上,哪个插值更精确;在区间[0,1]上,哪个插值更精确?
下列数据点的插值
x 0 1 4 9 16 25 36 49 64
y 0 1 2 3 4 5 6 7 8
可以得到平方根函数的近似,在区间[0,64]上作图.
(1)用这九个点作8次多项式插值L8(x)
(2)用三次样条(第一边界条件)程序求S(x)
从得到结果看在[0,64]上,哪个插值更精确;在区间[0,1]上,哪个插值更精确?
clear
>> x=[0 1 4 9 16 25 36 49 64];
>> y=0:8;
>> f=interp1(x,y,y);
>> F=interp1(x,y,y,'spline');
>> f,F
f =
0 1.0000 1.3333 1.6667 2.0000 2.2000 2.4000
2.6000 2.8000
F =
0 1.0000 1.5601 1.8402 2.0000 2.1706 2.3682
2.5806 2.7953
>> clear
>> x=[0 1 4 9 16 25 36 49 64];
\x0c
>> y=0:8;
>> t=0:0.1:70;
>> f=interp1(x,y,t,'spline');
>> plot(x,y,'k+',t,f,x,y,'r')
>>xlabel('x')
>>ylabel('y')
>>titel('y=x^1/2')
>>grid on
>> x=[0 1 4 9 16 25 36 49 64];
>> y=0:8;
>> f=interp1(x,y,y);
>> F=interp1(x,y,y,'spline');
>> f,F
f =
0 1.0000 1.3333 1.6667 2.0000 2.2000 2.4000
2.6000 2.8000
F =
0 1.0000 1.5601 1.8402 2.0000 2.1706 2.3682
2.5806 2.7953
>> clear
>> x=[0 1 4 9 16 25 36 49 64];
\x0c
>> y=0:8;
>> t=0:0.1:70;
>> f=interp1(x,y,t,'spline');
>> plot(x,y,'k+',t,f,x,y,'r')
>>xlabel('x')
>>ylabel('y')
>>titel('y=x^1/2')
>>grid on