文档介绍：该【计算方法上机题doc 】是由【Alone-丁丁】上传分享，文档一共【14】页，该文档可以免费在线阅读，需要了解更多关于【计算方法上机题doc 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。计算方法上机题doc
第八题:
一,题目
从函数表
x
f(x)
出发,用下列方法计算f(),f(),f()的近似值:
分段线性插值;(2)分段二次插值;(3)全区间上拉格朗日插值
(要求一次性输入整张函数表,并利用计算机选择在插值计算中所需的结点)。
二,基本方法
先选取两个结点与使[,],然后在区间[,]上作线性插值,即得f(x)≈。
选取距点x最近的三个结点,,进行二次插值,即取f(x)≈。
n=5,
三,计算结果
(1)
x
f(x)
(2)
x
f(x)
printf("Enterx:");
for(i=0;i<6;i++)
scanf("%lf",&a[i]);
printf("Enterf(x):");
for(i=0;i<6;i++)
scanf("%lf",&b[i]);
for(j=0;j<3;j++){
printf("Enterx:");
scanf("%lf",&x);
for(i=0;x>a[i];i++);
if((a[i+1]-x)<(x-a[i-2]))
{
y=b[i-1]*(x-a[i])/(a[i-1]-a[i])*(x-a[i+1])/(a[i-1]-a[i+1]);
y=y+b[i]*(x-a[i-1])/(a[i]-a[i-1])*(x-a[i+1])/(a[i]-a[i+1]);
y=y+b[i+1]*(x-a[i-1])/(a[i+1]-a[i-1])*(x-a[i])/(a[i+1]-a[i]);
printf("x=%lf,",x);printf("y=%lf\n",y);
}
else{
y=b[i-2]*(x-a[i-1])/(a[i-2]-a[i-1])*(x-a[i])/(a[i-2]-a[i]);
y=y+b[i-1]*(x-a[i-2])/(a[i-1]-a[i-2])*(x-a[i])/(a[i-1]-a[i]);
y=y+b[i]*(x-a[i-2])/(a[i]-a[i-2])*(x-a[i-1])/(a[i]-a[i-1]);
printf("x=%lf,",x);printf("y=%lf\n",y);
}
}
return0;
}
(3)
#include<>
intmain(void)
{
inti,j,k;
doublea[6],b[6],x,y,p;
printf("Enterx:");
for(i=0;i<6;i++)
scanf("%lf",&a[i]);
printf("Enterf(x):");
for(i=0;i<6;i++)
scanf("%lf",&b[i]);
for(j=0;j<3;j++)
{
y=0;
printf("Enterx:");
scanf("%lf",&x);
for(i=0;i<6;i++){
p=1;
for(k=0;k<6;k++){
if(k!=i)
p=p*(x-a[k])/(a[i]-a[k]);}
y=y+p*b[i];
}
printf("x=%lf,",x);printf("y=%lf\n",y);
}
return0;
}
六,附图
(1)
(2)
(3)
第九题:
一,题目
已知直升机旋转机翼外形曲线轮廓线上的某些型值点及端点处的一阶导数值
y`(x0)=,y`(x18)=-,
试计算该曲线上横坐标为2,4,6,12,16,30,60,110,180,280,400,515处点的纵坐标(要求该曲线具有二阶光滑度)。
k
0
1
2
3
4
5
6
k
13
14
15
16
17
18
468
494
507
520
k
7
8
9
10
11
12
78
35
二,基本方法
根据课本图4-9流程图及追赶法解方程。
三,计算结果
x
2
4
6
12
16
30
y
x
60
110
180
280
400
515
y
四,结果分析
利用追赶法求M[i],比高斯消去法步骤更简单,虽然结果有微小偏差,尤其在s[515]处与课本结果偏差较大。但在精度要求不高的情况下,追赶法可以减少工作量,提高工作效率。
五,源程序
#include<>
intmain(void)
{
inti,j;
doublea[19],b[19],h[19],s[18],t[18],g[19],x,y,p,q,f[19],m[20],n[20];
printf("输入左端点一阶导:");
scanf("%lf",&p);
printf("输入右端点一阶导:");
scanf("%lf",&q);
printf("Enterx:");
for(i=0;i<19;i++)
scanf("%lf",&a[i]);
printf("Enterf(x):");
for(i=0;i<19;i++)
scanf("%lf",&b[i]);
for(i=0;i<19;i++)
h[i]=a[i]-a[i-1];
for(i=1;i<18;i++)
s[i]=h[i]/(h[i]+h[i+1]);
for(i=1;i<18;i++)
t[i]=1-s[i];
for(i=1;i<18;i++)
g[i]=6/(h[i]+h[i+1])*((b[i+1]-b[i])/h[i+1]-(b[i]-b[i-1])/h[i]);
g[0]=6/h[1]*((b[1]-b[0])/h[1]-p);g[18]=6/h[18]*(q-(b[18]-b[17])/h[18]);
f[1]=;
for(i=2;i<19;i++)
f[i]=t[i-1]/(2-s[i-1]*f[i-1]);
n[1]=g[0]/2;
for(i=2;i<20;i++)
n[i]=(g[i-1]-s[i-1]*n[i-1])/(2-s[i-1]*f[i-1]);
m[19]=n[19];
for(i=18;i>0;i--)
m[i]=n[i]-f[i]*m[i+1];
printf("输出M:\n");
for(i=1;i<20;i++)
printf("%lf",m[i]);
for(j=0;j<12;j++){
printf("\n");
printf("Enterx:");
scanf("%lf",&x);
for(i=0;x>a[i];i++);
y=m[i]*(a[i]-x)*(a[i]-x)*(a[i]-x)/h[i]/6;
y=y+m[i+1]*(x-a[i-1])*(x-a[i-1])*(x-a[i-1])/h[i]/6;
y=y+(b[i-1]-m[i]*h[i]*h[i]/6)*(a[i]-x)/h[i];
y=y+(b[i]-m[i+1]*h[i]*h[i]/6)*(x-a[i-1])/h[i];
printf("f(x):%lf",y);
}
return0;
}
六,附图
第十二题
一,题目
某类疾病发生率‰和年龄段x(每五年为一段,例如0-5岁为第一段,6-10岁为第二段)之间有形如
的关系。试根据观测得到的如下数据表,用最小二乘法确定式中的参数a和b,并计算相应的均方误差与最大偏差。
x
1
2
3
4
5
6
7
8
9
y
x
10
11
12
13
14
15
16
17
18
19
y
二,基本方法
(1)将非线性方程化为线性方程,对等号两边取对数得到。
(2)按照最小二乘法的原始定义,不对方程进行变换,而均方误差的方程式,通过求解非线性方程组来解决问题,此时要用到牛顿迭代法。
三,计算结果
(1)a=,b=,e=,m=
(2)第一次:a=,b=,e=,m=
第二次:a=,b=,e=,m=.
四,结果分析
由以上结果对比,第二种方法结果更加精确,但是采用迭代法计算时,又必须多次选择初值才有可能取得较好的结果。