计算方法实验报告.docx
- 文档编号:1327125
- 上传时间:2022-10-20
- 格式:DOCX
- 页数:16
- 大小:194.44KB
计算方法实验报告.docx
《计算方法实验报告.docx》由会员分享,可在线阅读,更多相关《计算方法实验报告.docx(16页珍藏版)》请在冰豆网上搜索。
计算方法实验报告
中北大学信息商务学院
计算方法实验报告
学生:
昊文学号:
1603042130
学院:
中北大学信息商务学院
专业:
电气工程及其自动化
指导教师:
薛晓健
2017年04月19日
实验一:
非线性方程的近似解法
1.实验目的
1.掌握二分法和牛顿迭代法的原理
2.根据实验容编写二分法和牛顿迭代法的算法实现
注:
(可以用C语言或者matlab语言)
2.实验设备
matlab
3.实验容及步骤
解方程f(x)=x5-3x3-2x2+2=0
4.实验结果及分析
二分法:
数据:
f=x^5-3*x^3-2*x^2+2
[nxaxbxcfc]
1-3302
2.0000-3.00000-1.50000.0313
3.0000-3.0000-1.5000-2.2500-31.6182
4.0000-2.2500-1.5000-1.8750-8.4301
5.0000-1.8750-1.5000-1.6875-2.9632
6.0000-1.6875-1.5000-1.5938-1.2181
7.0000-1.5938-1.5000-1.5469-0.5382
8.0000-1.5469-1.5000-1.5234-0.2405
9.0000-1.5234-1.5000-1.5117-0.1015
10.0000-1.5117-1.5000-1.5059-0.0343
11.0000-1.5059-1.5000-1.5029-0.0014
12.0000-1.5029-1.5000-1.50150.0150
13.0000-1.5029-1.5015-1.50220.0068
14.0000-1.5029-1.5022-1.50260.0027
15.0000-1.5029-1.5026-1.50270.0007
16.0000-1.5029-1.5027-1.5028-0.0003
17.0000-1.5028-1.5027-1.50280.0002
18.0000-1.5028-1.5028-1.5028-0.0001
19.0000-1.5028-1.5028-1.50280.0001
20.0000-1.5028-1.5028-1.5028-0.0000
牛顿迭代法
>symsx;
f=(x^5-3*x^3-2*x^2+2)
[x,k]=Newtondd(f,0,1e-12)
f=x^5-3*x^3-2*x^2+2
x=NaN
k=2
实验二:
解线性方程组的迭代法
1.实验目的
1.掌握雅克比迭代法和高斯-塞德尔迭代法的原理
2.根据实验容编写雅克比迭代法和高斯-塞德尔迭代法的算法实现
注:
(可以用C语言或者matlab语言)
2.实验设备
Matlab
3.实验容及步骤
1、分别用雅克比迭代法和高斯-塞德尔迭代法解方程Ax=b
其中A=[4-10-100
-14-10-10
0-14-10-1
-10-14-10
0-10-14-1
00-10-14]
b=[0;5;-2;5;-2;6]
4.实验结果及分析
(雅克比迭代法)
a=[4-10-100;-14-10-10;0-14-10-1;-10-14-10;0-10-14-1;00-10-14]
b=[0;5;-2;5;-2;6]
x=agui_jacobi(a,b)
a=4-10-100
-14-10-10
0-14-10-1
-10-14-10
0-10-14-1
00-10-14
b=05-25-26
k=1
01.2500-0.50001.2500-0.50001.5000
k=2
0.62501.00000.50001.00000.50001.2500
k=3
0.50001.65630.31251.65630.31251.7500
k=4
0.82811.53130.76561.53130.76561.6563
k=5
0.76561.83980.67971.83980.67971.8828
k=6
0.91991.78130.89061.78130.89061.8398
k=7
0.89061.92530.85061.92530.85061.9453
k=8
0.96261.89790.94901.89790.94901.9253
k=9
0.94901.96510.93031.96510.93031.9745
k=10
0.98261.95240.97621.95240.97621.9651
k=11
0.97621.98370.96751.98370.96751.9881
k=12
0.99191.97780.98891.97780.98891.9837
k=13
0.98891.99240.98481.99240.98481.9944
k=14
0.99621.98960.99481.98960.99481.9924
k=15
0.99481.99650.99291.99650.99291.9974
k=16
0.99821.99520.99761.99520.99761.9965
k=17
0.99761.99830.99671.99830.99671.9988
k=18
0.99921.99770.99891.99770.99891.9983
k=19
0.99891.99920.99851.99920.99851.9994
k=20
0.99961.99890.99951.99890.99951.9992
k=21
0.99951.99960.99931.99960.99931.9997
k=22
0.99981.99950.99981.99950.99981.9996
k=23
0.99981.99980.99971.99980.99971.9999
k=24
0.99991.99980.99991.99980.99991.9998
k=25
0.99991.99990.99981.99990.99981.9999
k=26
1.00001.99990.99991.99990.99991.9999
k=27
0.99992.00000.99992.00000.99992.0000
x=0.99992.00000.99992.00000.99992.0000
(高斯-赛德尔迭代法迭代法)
a=[4-10-100;-14-10-10;0-14-10-1;-10-14-10;0-10-14-1;00-10-14]
b=[0;5;-2;5;-2;6]
x=agui_GS(a,b)
a=4-10-100
-14-10-10
0-14-10-1
-10-14-10
0-10-14-1
00-10-14
b=05-25-26
k=1
Columns1through5
01.0000-0.00001.00000.0000
Column6
1.0000
k=2
Columns1through5
0.00001.00000.50001.37500.8438
Column6
1.5859
k=3
Columns1through5
0.34381.67190.40821.64890.7267
Column6
1.2837
k=4
Columns1through5
0.33021.11630.26221.82980.5574
Column6
1.4549
k=5
Columns1through5
0.73651.88910.54341.95940.0758
Column6
1.6548
k=6
Columns1through5
0.96211.39530.25241.32260.3432
Column6
1.6489
k=7
Columns1through5
0.67951.31880.07261.77380.9354
Column6
1.5020
k=8
Columns1through5
0.02311.25780.38341.83550.8988
Column6
1.3206
k=9
Columns1through5
0.02331.82640.74561.41690.3910
Column6
1.7842
k=10
Columns1through5
0.31081.36190.39071.02310.0423
Column6
1.1083
k=11
Columns1through5
0.34631.69480.20661.64880.8630
Column6
1.5174
k=12
Columns1through5
0.83591.47640.91061.90240.7240
Column6
1.1587
k=13
Columns1through5
0.84471.86980.73271.32540.3385
Column6
1.5178
k=14
Columns1through5
0.29881.34250.04641.92090.4453
Column6
1.8729
x=0.2988
1.3425
0.0464
1.9209
0.4453
1.8729
实验三:
插值与拟合
1.实验目的
1、掌握线性插值、抛物线插值、拉格朗日插值,三次样条插值与拟合
2、根据实验容,编写三次样条插值(一阶导数)的算法实现。
3、根据实验容,编写最小二乘法的算法实现。
2.实验设备
Matlab
3.实验容及步骤
(1)、给定的插值条件如下:
i
0
1
2
3
4
5
6
7
Xi
8
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算方法 实验 报告
![提示](https://static.bdocx.com/images/bang_tan.gif)