数值分析计算实习题第三章第二次作业:题一:
x1021y1125x2f1polyfitxy3fpoly2symf1y1polyvalf1xx2li
space1110y2i
terp1xyx2plotxyrxy1bholdo
plotx2y2klege
d数据点3次拟合曲线3次多项式插值xlabelXylabelY
输出:f100000057520000004841
f
4591875547102675x3811296384146066816957890051440643305x25746
1469057404776431x20282409603651670423947251286016
43606096623006139007199254740992
y1
009110116002771027710116000911
03921
04611
04841
04611
03921
fx2
10000077780555603333011110777810000
01111
03333
05556
y2
00385006340063400385
01222
03000
07222
07222
03000
01222
题二:
X00010203050810Y10041050061091202246p1polyfitXY3p2polyfitXY4Y1polyvalp1XY2polyvalp2X
fplotXYrXY1bXY2gp3polyfitXY2Y3polyvalp3Xf1poly2symp1f2poly2symp2f3poly2symp3plotXYrXY1bXY2gXY3mlege
d数据点3次多项式拟合4次多项式拟合2次多项式拟合xlabelX轴ylabelY轴
输出:p1
662211281474659109266
p2288531233481627475298709427
Y109266058220454405034097302010324602
Y209427056350439905082100051986024692
p3313161240007356
fY307356064290612806454089841747726271
f17455778416425075x311258999068426241803512222945435x214073748835532840981580032809x879609302220883459537843990119007199254740992f21624271450198125x45629499534213123471944732519173x32814749767106564580931990070659x22814749767106561491459232922115x28147497671065610614094330812931125899906842624f318733x2598274179x598207333799700
题三:建立三角插值函数的m文件
ffu
ctio
ABY1Rmsa
jiaobiji
XYX1mAB分别是m阶三角多项式Tm(x)的系数aj,bj(j12m)的系数矩阵,Y1是Tm(x)在X1处的值,XY数据点Rm为均方误差
le
gthX1max1fix
12ifmmax1
mmax1e
dAzeros1m1Bzeros1m1YmY1Y
12Y1YmY
1YmA12sumY
fori1m
Bi1si
iXYAi1cosiXYe
dA2A
B2B
A1A12Y1A1fork1mY1Y1Ak1coskX1Bk1si
kX1TmA1Ak1coskXBk1si
kXkk1e
dYTmRmsumYTm2
输出:Xpi2pi33pi
fYX2cosX
ABY1Rmsa
jiaobiji
XYX116
输出:
A
1至12列
01397
4400228326
12355
0698206773066350654106474
09128
13至17列
0642606393063700635506348
0791407319
B
10e15
1至12列
0001940015000044003000082100599001330021100297
00105
00627
13至17列
001780096201049即可得16插值多项式的r