数值计算偏重分析-插值和Newton-Cotes-matlab

《数值计算偏重分析》作业

1、

用三次样条插值spline

x=1:12
y=[5,8,9,15,25,29,31,30,22,25,27,24]
xx=1:0.1:12
p=spline(x,y)
yy=ppval(p,xx)

插值结果为：

5 5.624206104 6.154452386 6.598956604 6.965936515 7.263609878 7.500194451 7.683907991 7.822968257 7.925593008 8 8.054406992 8.097031743 8.136092009 8.179805549 8.236390122 8.314063485 8.421043396 8.565547614 8.755793896 9 9.305165927 9.673420644 10.10567536 10.60284129 11.16582963 11.79555161 12.49291842 13.25884129 14.09423141 15 15.9739293 17.00128568 18.06420655 19.1448293 20.22529134 21.28773008 22.31428291 23.28708724 24.18828047 25 25.70911687 26.32143663 26.84749845 27.29784151 27.68300499 28.01352808 28.29994994 28.55280977 28.78264673 29 29.21360322 29.42496781 29.63379967 29.83980466 30.04268868 30.2421576 30.43791731 30.6296737 30.81713263 31 31.17547026 31.33069212 31.45030289 31.51893984 31.52124029 31.44184151 31.2653808 30.97649545 30.55982275 30 29.29051576 28.46026371 27.54698879 26.58843596 25.62235017 24.68647636 23.81855949 23.05634451 22.43757636 22 21.7704667 21.73225305 21.85774195 22.11931632 22.48935905 22.94025305 23.44438124 23.97412653 24.50187181 25 25.44561743 25.83472408 26.1680434 26.44629878 26.67021365 26.84051143 26.95791554 27.02314939 27.0369364 27 26.9130636 26.77685061 26.59208446 26.35948857 26.07978635 25.75370122 25.3819566 24.96527592 24.50438257 24

2、

三次样条插值

x=[0,3,5,7,9,11,12,13,14,15]
y=[0,1.2,1.7,2.0,2.1,2.0,1.8,1.2,1.0,1.6]
xx=1:0.1:15
p=spline(x,y)
yy=ppval(p,xx)
xlswrite('C:\Users\86159\Desktop\data2',yy)

0.466537249 0.50936064 0.551507218 0.592981886 0.633789547 0.673935105 0.713423461 0.75225952 0.790448184 0.827994355 0.864902937 0.901178833 0.936826946 0.971852178 1.006259433 1.040053613 1.073239622 1.105822362 1.137806736 1.169197648 1.2 1.230218695 1.259858636 1.288924726 1.317421868 1.345354965 1.37272892 1.399548635 1.425819014 1.45154496 1.476731375 1.501383163 1.525505226 1.549102468 1.572179791 1.594742098 1.616794292 1.638341277 1.659387954 1.679939228 1.7 1.719574454 1.738663893 1.757268898 1.775390052 1.793027937 1.810183136 1.826856232 1.843047805 1.85875844 1.873988718 1.888739222 1.903010533 1.916803235 1.93011791 1.94295514 1.955315507 1.967199595 1.978607984 1.989541259 2 2.009983488 2.019485794 2.028499683 2.037017925 2.045033287 2.052538535 2.059526439 2.065989764 2.07192128 2.077313753 2.082159951 2.086452641 2.090184591 2.093348569 2.095937342 2.097943678 2.099360344 2.100180108 2.100395738 2.1 2.098991593 2.097392933 2.095232369 2.092538248 2.089338917 2.085662723 2.081538014 2.076993138 2.072056441 2.066756271 2.061120976 2.055178903 2.0489584 2.042487813 2.03579549 2.028909779 2.021859027 2.014671582 2.00737579 2 1.99244448 1.984097184 1.974217988 1.962066769 1.946903402 1.927987763 1.904579727 1.875939171 1.84132597 1.8 1.751566427 1.697011572 1.637667047 1.574864465 1.509935436 1.444211572 1.379024485 1.315705786 1.255587087 1.2 1.150089814 1.106256528 1.068713822 1.037675371 1.013354855 0.99596595 0.985722334 0.982837686 0.987525682 1 1.020474318 1.049162314 1.086277666 1.13203405 1.186645145 1.250324629 1.323286178 1.405743472 1.497910186 1.6

3、

x=[200,220,250,270,280]
y=[4,4.5,4.7,4.8,5.2]
xx=240
p=spline(x,y)
yy=ppval(p,xx)
%xlswrite('C:\Users\86159\Desktop\data2',yy)

结果为：4.6924

4、

并分析误差

f=@(x) log(x)

%梯形公式
a=1;b=2
int1=(b-a)/2*(f(a)+f(b))
error1=(b-a)^3/12*1

%simpson
int2=(b-a)/6*(f(a)+4*f((a+b)/2)+f(b))
error2=(b-a)/180*((b-a)/2)^4*6

%cotes
t=linspace(a,b,5)
int3=(b-a)/90*(7*f(a)+32*f(t(2))+12*f(t(3))+32*f(t(4))+7*f(b))
error3=2*(b-a)/945*((b-a)/4)^6*120

%导数用以下方法求
% syms m
% u=log(m)
% w=diff(u,6)

梯形公式求得的积分值为：0.3466，截断误差不超过：0.0833

Simpson公式求得积分值：0.3858，截断误差不超过：0.0021

Cotes公式求得积分值为：0.3863，截断误差不超过：6.2004e-05

833

Simpson公式求得积分值：0.3858，截断误差不超过：0.0021

Cotes公式求得积分值为：0.3863，截断误差不超过：6.2004e-05