计算方法上机题文档格式.docx
- 文档编号:18723171
- 上传时间:2022-12-31
- 格式:DOCX
- 页数:15
- 大小:856.26KB
计算方法上机题文档格式.docx
《计算方法上机题文档格式.docx》由会员分享,可在线阅读,更多相关《计算方法上机题文档格式.docx(15页珍藏版)》请在冰豆网上搜索。
(2)选取距点x最近的三个结点
,
进行二次插值,即取f(x)≈
(3)n=5,
三,计算结果
(1)
0.15
0.31
0.47
0.394039
0.380067
0.356927
(2)
0.394482
0.380225
0.361399
(3)
0.394473
0.380219
0.357222
四,结果分析
编程基本按照课本套公式,方法从
(1)至(3),计算结果越来越精确。
五,源程序
#include<
stdio.h>
intmain(void)
{
inti,j;
doublea[6],b[6],x,y;
printf("
Enterx:
"
);
for(i=0;
i<
6;
i++)
scanf("
%lf"
&
a[i]);
Enterf(x):
scanf("
b[i]);
for(j=0;
j<
3;
j++){
x);
x>
a[i];
i++);
y=b[i-1]*(x-a[i])/(a[i-1]-a[i])+b[i]*(x-a[i-1])/(a[i]-a[i-1]);
x=%lf,"
x);
y=%lf\n"
y);
}
return0;
}
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]);
}
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]);
inti,j,k;
doublea[6],b[6],x,y,p;
j++)
y=0;
i++){
p=1;
for(k=0;
k<
k++){
if(k!
=i)
p=p*(x-a[k])/(a[i]-a[k]);
y=y+p*b[i];
六,附图
第九题:
已知直升机旋转机翼外形曲线轮廓线上的某些型值点及端点处的一阶导数值
y`(x0)=1.86548,y`(x18)=-0.046115,
试计算该曲线上横坐标为2,4,6,12,16,30,60,110,180,280,400,515处点的纵坐标(要求该曲线具有二阶光滑度)。
k
1
2
3
4
5
6
0.52
3.1
8.0
17.95
28.65
39.62
50.65
5.28794
9.4
13.84
20.2
24.9
28.44
31.1
13
14
15
16
17
18
364.4
416.3
468
494
507
520
20.9
14.8
7.8
3.7
1.5
0.2
7
8
9
10
11
12
78
104.6
156.6
208.6
260.7
312.5
35
36.5
36.6
34.6
31.0
26.34
根据课本图4-9流程图及追赶法解方程。
30
y
7.834331
10.476880
12.358485
16.579653
19.093173
25.386597
60
110
180
280
400
515
32.804369
36.647886
35.917141
29.368427
16.799141
0.524459
利用追赶法求M[i],比高斯消去法步骤更简单,虽然结果有微小偏差,尤其在s[515]处与课本结果偏差较大。
但在精度要求不高的情况下,追赶法可以减少工作量,提高工作效率。
doublea[19],b[19],h[19],s[18],t[18],g[19],x,y,p,q,f[19],m[20],n[20];
输入左端点一阶导:
p);
输入右端点一阶导:
q);
19;
h[i]=a[i]-a[i-1];
for(i=1;
18;
s[i]=h[i]/(h[i]+h[i+1]);
t[i]=1-s[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]=0.5;
for(i=2;
f[i]=t[i-1]/(2-s[i-1]*f[i-1]);
n[1]=g[0]/2;
20;
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];
输出M:
\n"
printf("
%lf"
m[i]);
12;
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];
f(x):
第十二题
某类疾病发生率
‰和年龄段x(每五年为一段,例如0-5岁为第一段,6-10岁为第二段)之间有形如
的关系。
试根据观测得到的如下数据表,用最小二乘法确定式中的参数a和b,并计算相应的均方误差与最大偏差。
0.898
2.38
3.07
1.84
2.02
1.94
2.22
2.77
4.02
19
4.76
5.46
6.53
10.9
16.5
22.5
35.7
50.6
61.6
81.8
(1)将非线性方程化为线性方程,对
等号两边取对数得到
(2)按照最小二乘法的原始定义,不对方程进行变换,而均方误差的方程式,通过求解非线性方程组来解决问题,此时要用到牛顿迭代法。
(1)a=0.681361,b=0.230620,e=77.307455,m=18.328522
(2)第一次:
a=0.721560,b=0.249856,e=17.925350,m=8.1115274
第二次:
a=0.239368,b=0.310747,e=5.430872,m=3.467987.
由以上结果对比,第二种方法结果更加精确,但是采用迭代法计算时,又必须多次选择初值才有可能取得较好的结果。
math.h>
inti;
doublea,b,e,m,a0,a1,b0,b1,y1,y2,x[19],y[19];
x[i]=i+1;
x:
"
x[i]);
\ny:
y[i]);
y[i]=log(y[i]);
a0=a1=b0=b1=y1=y2=0;
a1=a1+x[i];
y1=y1+y[i];
b0=b0+x[i];
b1=b1+pow(x[i],2);
y2=y2+x[i]*y[i];
a0=19;
b=(y2-y1*(b0/a0))/(b1-a1*(b0/a0));
a=(y1-b*a1)/a0;
a=exp(a);
y[i]=exp(y[i]);
e=m=0;
e=e+pow((a*exp((i+1)*b)-y[i]),2);
if(pow((a*exp((i+1)*b)-y[i]),2)>
m)
m=pow((a*exp((i+1)*b)-y[i]),2);
else;
a=%lfb=%lfe=%lfm=%lf\n"
a,b,e/19,sqrt(m));
doublea,b,e,m,f,f1,f2,g,g1,g2,c,d,x[19],y[19];
Entera0:
a);
Enterb0:
b);
f=f1=f2=g=g1=g2=0;
f=f+2*(a*exp((i+1)*b)-y[i])*exp((i+1)*b);
f1=f1+2*exp(2*(i+1)*b);
f2=f2+4*(i+1)*a*exp(2*(i+1)*b)-2*(i+1)*y[i]*exp((i+1)*b);
g=g+2*(a*exp((i+1)*b)-y[i])*(i+1)*a*exp((i+1)*b);
g1=g1+4*(i+1)*a*exp(2*(i+1)*b)-2*(i+1)*y[i]*exp((i+1)*b);
g2=g2+4*pow((i+1)*a,2)*exp(2*(i+1)*b)-2*pow((i+1),2)*a*y[i]*exp((i+1)*b);
c=f1*a+f2*b-f;
d=g1*a+g2*b-g;
b=(d-c*(g1/f1))/(g2-f2*(g1/f1));
a=(c-f2*b)/f1;
e=0;
m=0;
m=pow((a*exp((i+1)*b)-y[i]),2);
a=%lfb=%lf均方误差e=%lf最大偏差m=%lf\n"
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算方法 上机