北航数值B大作业二Word格式.docx
- 文档编号:18366226
- 上传时间:2022-12-15
- 格式:DOCX
- 页数:20
- 大小:73.08KB
北航数值B大作业二Word格式.docx
《北航数值B大作业二Word格式.docx》由会员分享,可在线阅读,更多相关《北航数值B大作业二Word格式.docx(20页珍藏版)》请在冰豆网上搜索。
doublel[5],s;
for(m=0;
m<
c;
m++)
for(k=0,l[m]=1;
k<
k++)
{if(k!
=m)l[m]*=(b-X[a+k])/(X[a+m]-X[a+k]);
for(m=0,s=0;
s+=l[m]*Y[a+m];
returns;
voiderci()//二次插值
/prcline1.txt"
if(i<
=0.5*(X[1]+X[2]))
{y[i]=lagrange(0,i,3);
elseif(i>
0.5*(X[10]+X[11]))
{y[i]=lagrange(10,i,3);
else
{for(j=0;
j<
10;
j++)
{if(i>
0.5*(X[j]+X[j+1])&
=0.5*(X[j+1]+X[j+2]))
y[i]=lagrange(j,i,3);
}
voidsanci()//三次插值
/prcline2.txt"
=X[2])
{y[i]=lagrange(0,i,4);
X[10])
{y[i]=lagrange(9,i,4);
x=%.2lf,y=%9.6lf\n"
{for(j=1;
9;
X[j+1]&
=X[j+2])
y[i]=lagrange(j,i,4);
voidsanciyangtiao()//三次样条插值
inti,j,k,t;
doublea[N][N]={0},m[N],b[N],c,d,h[N],x[M],y[M];
/prcline3.txt"
a[0][1]=b[0]=b[N-1]=a[N-1][N-2]=0;
//确定系数矩阵
a[0][0]=a[N-1][N-1]=2;
for(i=1;
N;
h[i]=X[i]-X[i-1];
N-1;
{a[i][i]=2;
a[i][i+1]=h[i+1]/(h[i]+h[i+1]);
a[i][i-1]=1-a[i][i+1];
b[i]=6/(h[i]+h[i+1])*((Y[i+1]-Y[i])/h[i+1]-(Y[i]-Y[i-1])/h[i]);
for(k=0;
k++)//LU法解方程组
{for(j=k;
{for(t=0,c=0;
t<
k;
t++)c+=a[k][t]*a[t][j];
a[k][j]=a[k][j]-c;
if(k<
N-1)
{for(t=0,d=0;
t++)d+=a[j+1][t]*a[t][k];
a[j+1][k]=(a[j+1][k]-d)/a[k][k];
i;
t++)c+=a[i][t]*b[t];
b[i]=b[i]-c;
m[N-1]=b[N-1]/a[N-1][N-1];
for(i=N-2;
i>
=0;
i--)
{for(t=i+1,d=0;
t++)d+=a[i][t]*m[t];
m[i]=(b[i]-d)/a[i][i];
{x[i]=i;
for(j=1;
{if(i>
=X[j-1]&
=X[j])
y[i]=m[j-1]/(6*h[j])*pow((X[j]-i),3)+m[j]/(6*h[j])*pow((i-X[j-1]),3)+(Y[j-1]/h[j]-m[j-1]*h[j]/6)*(X[j]-i)+(Y[j]/h[j]-m[j]*h[j]/6)*(i-X[j-1]);
main()//主函数
xianxing();
erci();
sanci();
sanciyangtiao();
return0;
三、程序运行结果
x
线性插值
二次插值
三次插值
三次样条
1
1.122363
1.334295
1.408378
1.250028
2
2.244726
2.555257
2.651036
2.464375
3
3.367089
3.662886
3.741956
3.607359
4
4.489451
4.657184
4.695122
4.6433
5
5.471849
5.538148
5.524515
5.536721
6
6.055882
6.30578
6.244118
6.274823
7
6.639916
6.96008
6.867913
6.887646
8
7.22395
7.284531
7.409883
7.409992
9
7.807983
7.834371
7.884009
7.876665
10
8.303684
8.359434
8.355296
8.321714
11
8.711053
8.859719
8.846586
8.761856
12
9.118421
9.335226
9.313015
9.196485
13
9.525789
9.785956
9.755631
9.624236
14
9.933158
10.211907
10.175483
10.043746
15
10.340526
10.485162
10.573621
10.453654
16
10.747895
10.876094
10.951092
10.852596
17
11.155263
11.253878
11.308947
11.23921
18
11.562632
11.618513
11.648233
11.612132
19
11.97
20
12.19
12.308338
12.29998
12.31175
21
12.41
12.633527
12.616236
12.637517
22
12.63
12.945568
12.919033
12.947734
23
12.85
13.24446
13.208638
13.242833
24
13.07
13.530203
13.485316
13.52325
25
13.29
13.802798
13.749331
13.789416
26
13.51
14.062244
14.000949
14.041765
27
13.73
14.308541
14.240436
14.280731
28
13.95
14.54169
14.468057
14.506747
29
14.17
14.572632
14.684076
14.720247
30
14.39
14.783684
14.888761
14.921663
31
14.61
14.985789
15.082375
15.111429
32
14.83
15.178947
15.265184
15.289979
33
15.05
15.363158
15.437454
15.457746
34
15.27
15.538421
15.59945
15.615162
35
15.49
15.704737
15.751438
15.762663
36
15.71
15.862105
15.893681
15.900681
37
15.93
16.010526
16.026447
16.029648
38
16.15
39
16.2
16.280526
16.26723
16.262128
40
16.25
16.402105
16.375734
16.366266
41
16.3
16.514737
16.475734
16.462605
42
16.35
16.618421
16.567452
16.551336
43
16.4
16.713158
16.651108
16.632653
44
16.45
16.798947
16.726925
16.706746
45
16.5
16.875789
16.795125
16.773808
46
16.55
16.943684
16.855928
16.834031
47
16.6
17.002632
16.909557
16.887607
48
16.65
16.863158
16.956233
16.934728
49
16.7
16.908421
16.996177
16.975585
50
16.75
16.948947
17.029612
17.010371
51
16.8
16.984737
17.056759
17.039278
52
16.85
17.015789
17.077839
17.062497
53
16.9
17.042105
17.093075
17.080221
54
16.95
17.063684
17.102687
17.092641
55
17.080526
17.106898
17.099949
56
17.05
17.092632
17.105928
17.102339
57
17.1
58
17.06
17.102632
17.095983
17.093115
59
17.02
17.100526
17.087341
17.081823
60
16.98
17.093684
17.074183
17.066252
61
16.94
17.082105
17.05662
17.046531
62
17.065789
17.034765
17.022788
63
16.86
17.044737
17.008726
16.995152
64
16.82
17.018947
16.978615
16.96375
65
16.78
16.988421
16.944543
16.928711
66
16.74
16.953158
16.90662
16.890165
67
16.818421
16.864958
16.848238
68
16.66
16.775789
16.819668
16.803061
69
16.62
16.730526
16.770859
16.75476
70
16.58
16.682632
16.718643
16.703465
71
16.54
16.632105
16.66313
16.649303
72
16.578947
16.604432
16.592404
73
16.46
16.523158
16.542659
16.532896
74
16.42
16.464737
16.477922
16.470907
75
16.38
16.403684
16.410332
16.406565
76
16.34
77
16.273684
16.272022
16.271311
78
16.16
16.204737
16.20144
16.200485
79
16.07
16.133158
16.128283
16.127483
80
15.98
16.058947
16.052576
16.052262
81
15.89
15.982105
15.974349
15.974782
82
15.8
15.902632
15.893629
15.895002
83
15.820526
15.810443
15.812882
84
15.62
15.735789
15.72482
15.72838
85
15.53
15.648421
15.636787
15.641456
86
15.44
15.534737
15.546371
15.552069
87
15.35
15.442632
15.453601
15.460178
88
15.26
15.348421
15.358504
15.365742
89
15.17
15.252105
15.261108
15.26872
90
15.08
15.153684
15.16144
15.169071
91
14.99
15.053158
15.059529
15.066755
92
14.9
14.950526
14.955402
14.961731
93
14.81
14.845789
14.849086
14.853958
94
14.72
14.738947
14.740609
14.743394
95
14.63
96
14.5
14.518947
14.512299
14.513774
97
14.37
14.405789
14.392604
14.394878
98
14.24
14.290526
14.271025
14.273513
99
14.11
14.173158
14.147673
14.149879
100
13.98
14.053684
14.022659
14.024179
101
13.85
13.932105
13.896094
13.896612
102
13.72
13.808421
13.768089
13.767381
103
13.59
13.682632
13.638753
13.636687
104
13.46
13.554737
13.508199
13.50473
105
13.33
13.376537
13.371712
106
13.2
13.243878
13.237835
107
13.110332
13.103298
108
12.94
12.976011
12.968304
109
12.81
12.841025
12.833054
110
12.68
12.705485
12.697749
111
12.55
12.569501
12.56259
112
12.42
12.433186
12.427778
113
12.29
12.296648
12.293514
114
12.16
115
12.03
12.031662
12.027387
116
11.9
11.903296
11.895625
117
11.77
11.774875
11.764616
118
11.64
11.646371
11.634261
119
11.51
11.517756
11.50446
120
11.38
11.389003
11.375114
121
11.25
11.260083
11.246125
122
11.12
11.13097
11.117394
123
10.99
11.001634
10.98882
124
10.86
10.883684
10.87205
10.860306
125
10.73
10.753158
10.742188
10.731752
126
10.6
10.622105
10.612022
10.60306
127
10.47
10.490526
10.481524
10.474129
128
10.34
10.358421
10.350665
10.344861
129
10.21
10.225789
10.219418
10.215157
130
10.08
10.092632
10.087756
10.084918
131
9.95
9.958947
9.955651
9.954045
132
9.82
9.824737
9.823075
9.822439
133
9.69
134
9.55
9.554737
9.556399
9.556651
135
9.41
9.418947
9.422244
9.4224
136
9.27
9.282632
9.287507
9.287275
137
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 北航 数值 作业