文档介绍：2016-2017第一学期数值分析
上机实验报告
姓名：XXX
学号：20162
学院：上木工程学院
导师： •
联系电话： •
指导老师:
目录
第一题 1
1
1
)) ; $向里的2范数控制渓差
P=H‘ ：
if (wc<wcx)
break;
end if (k==n)
fprintf ('迭代发散'): return;
end
end
x=p：
HDZ=A*x:
%解向里 q迭代灰数
5FC
HD2
end
$误差限
%返回检验值
E] function DI2TI (n)
为b列向里的第一个元素值
-%
A=5*eye(50):
A<1,2)=-1;
A(50, 49)=-1;
Hfor(j=2:49)
A(j, [j-1, j+l])=-l；
曰id $该部分输入A中元素
□ for j=l:50
B(l, j)=j:
end
3 for
i=2:50
H
for j=l:50
B=zeros (50, 50):
%输入B中策一行元秦
B(i, j)+l;
end
中每行元寿等于上行加一
end
]for
i=l:50 for
j=l:50
j)>50)
B(i, j)-50: %蛮找B中大于50的數,减去刊，再附给该位置
end
end
end
®该部分将所有b敢入E中，输入B中元素
b=E(Il:50]?n); ykb (A, b)
备调用雅可比函数，输出i咕HD2
end
（命令行窗口显示内容）数据太多截取一部分
【= HDZ =
1. 3384
5. 0000
—
□
6. 0000
6
2. 3212
7. 0000
■
i
2. 6641
8. 0000
s
2. 9995
9. 0000
9
3. 3332
10. 0000
10
3. 6666
11. 0000
11
4. 0000
12. 0000
12
4. 3333
13. 0000
13
4. 6667
14. 0000
14
□. 0000
15. 0000
15
5. 3333
16. 0000
16
5. 6667
17. 0000
17
6. 0000
18. 0000
18
6. 3333
19
6. 6667
20. 0000
20
7. 0000
21. 0000
21
7. 3333
22. 0000
22
7. 6667
23. 0000
23
8. 0000
24. 0000
24
8. 3333
25. 0000
25
8. 6667
26. 0000
26
9. 0000
27. 0000
27
9. 3333
28. 0000
28
9. 6667
29. 0000
29
10. 0000
30. 0000
30
10. 3333 k =
wc =
31. 0000
31
10. 6667
32. 0000
32
11. 0000 18
6. 5841e-06
33. 0000
33
11. 3333
34. 0000
34
第三题
3・1题目要求
对函数/(%) = 在区间卜1,1]上取血=-1 + / (i = 0,l,2 -,10),
对函数进行多项式插值和三次样条插值，并画出插值函数及f(x)的函 数；
对函数求其三次拟合曲线并画出拟合曲线的图像，与(a)中结果进行 比较。
□I function Z] = DI3TI (x, y)
%本函数是对数堀和y，邊行线性，三次样条，3次拟合，作出囹的函数
x= input ('输入实验数据’'
X
y= input ('输入实验数据’ry ?
y
N=length(x):
X=x(l):0. 01:x(N):
峨性釦页式怔值函数值 %三次样条插值函数值 辂次拟台多顷式系数 頤次拟台函数值
f=l./(l+25*(X. ”2));
f1=interpl(x5 y, X):
f 2= spline (x9 y9 X):
a=polyf it (x, y, 3):
f3=poLyvraL(a, X):
plot (x,y/+9,X, f/k\X7fl/g\X