数值分析上机(四五章).doc

数值分析上机题目

学院：土木工程学院

姓名：黎亮（B）

学号：140973

第四章

一、题目

（1）编制求第一型3次样条插值函数的通用程序；

（2）已知汽车门曲线型值点的数据如下：

i

0

1

2

3

4

5

xi

0

1

2

3

4

5

yi

2.51

3.30

4.04

4.70

5.22

5.54

i

6

7

8

9

10

xi

6

7

8

9

10

yi

5.78

5.40

5.57

5.70

5.80

端点条件为，用所编程序求车门的3次样条插值函数S（x），并打印出S（i+0.5）（i=0,1,……,9）。

二、通用程序（本题我是用matlab编的，如下所示：）

clc

%%输入相关参数

n=input(请输入X的小标i的最大值n=);

n=n+1;

xn=zeros(1,n);

yn=zeros(1,n);

xn(1,:)=input(请以矩阵的形式输入X值的大小：);

yn(1,:)=input(请以矩阵的形式输入与X值对应的Y：);

daoshu0=input(请输入起点的1阶导数);%输入边界条件

daoshun=input(请输入终点的1阶导数);

%%求dj的值

d=zeros(n,1);

h=zeros(1,n-1);

chashang1=zeros(1,n-1);

chashang2=zeros(1,n-2);

fori=1:n-1

h(i)=xn(i+1)-xn(i);

chashang1(i)=(yn(i+1)-yn(i))/h(i);%求出一阶差商

end

fori=2:n-1

chashang2(i)=(chashang1(i)-chashang1(i-1))/(xn(i+1)-xn(i-1));%求二阶差商

d(i)=6*chashang2(i);

end

d(1)=6*(chashang1(1)-daoshu0)/h(1);

d(n)=6*(daoshun-chashang1(n-1))/h(n-1);

%%求M值

B=zeros(n);

miu=zeros(1,n-2);

lmd=zeros(1,n-2);

fori=1:n-2

miu(i)=h(i)/(h(i)+h(i+1));%求μ和λ

lmd(i)=1-miu(i);

end

%%给矩阵赋值

B(1,2)=1;

B(n,n-1)=1;

fori=1:n

B(i,i)=2;

end

fori=2:n-1

B(i,i-1)=miu(i-1);

B(i,i+1)=lmd(i-1);

end

M=B\d;

%%求插值函数并输出结果

disp(S(x)=);

symsx;

fori=1:n-1

Sx(1,i)=(yn(i)+(chashang1(i)-(M(i)/3+M(i+1)/6)*h(i))*(x-xn(i))...

+0.5*M(i)*(x-xn(i))^2+(M(i+1)-M(i))/(6*h(i))*(x-xn(i))^3);

fprintf(%s(%d,%d)\n,char(vpa(Sx(1,i),5)),xn(i),xn(i+1))

disp(================================================================);

end

%%输出S(i+0.5)的值

fori=1:n-1

x=xn(i)+0.5;

Y(i)=yn(i)+(chashang1(i)-(M(i)/3+M(i+1)/6)*h(i))*(x-xn(i))...

+0.5*M(i)*(x-xn(i))^2+(M(i+1)-M(i))/(6*h(i))*(x-xn(i))^3;

end

disp(S(i+0.5));

disp(ix(i+0.5)S(i+0.5));

fori=1:n-1

fprintf(%d%.5f%.5f\n,i-1,xn(i)+0.5,Y(i));

End

三、求解结果

请输入X的小标i的最大值n=10

请以矩阵的形式输入X值的大小：[012345678910]

请以矩阵的形式输入与X值对应的Y：[2.513.304.044.705.225.545.785.405.575.705.80]

请输入起点的1阶导数0.8

请输入终点的1阶导数0.2

S(x)=

-0.008514*x^3-0.001486*x^2+0.8*x+2.51(0,1)

=============================================================