文档介绍：数值分析实验
董海云
数理学院数学实验教学中心
0 Matlab介绍入门知识 3
1绪论 17
1 例题解答 17
Matlab中数值计算精度 20
2线性方程组的直接解法 22
1例题解答 22
2 Mat ，只要内存或文件中存在即可，变量名之，保存或装入全部变量.-ascii选项使文件以 ASCII格式处理，-append选项控制将变量追加到MAT文件中.
⑴向量的创建
用步长生成法：
数组=初值:步长（增量）:终值
» a= 1::3
a =

用linspace生成：
数组=linspace（初值，终值，等分点数目）
» b=linspace（l,3,5） b =

列向量用分号（;）作为分行标记：
» c=[l;2;3;4;]
C =
1
2
3
4
若不想输出结果，在每一条语句后用分号作为结束符，若留空或用逗号结束，则在执行该语句后会把结果输出来.
» a+b;
» a+b
ans =
3 4 5 6
⑵矩阵的创建
直接输入：
：将矩阵的元素用方括号括起来，按矩阵行的顺序输入各元素，同一行的各元素之间用空格或逗号分隔，不同行的元素之间用分号分隔.
» A=[l 2 3;4 5 6;2 3 5] A =
12 3
4 5 6
2 3 5
利用矩阵函数创建：
» B=magic（3）%魔方阵 B =
8 16
5 7
9 2 » C=hilb(3)%3 阶 Hilbert 矩阵
C =

Matlab中用％引导注释
其它创建矩阵函数还有：
eye(m,n):生成m行n列单位矩阵.
zeros(m,n):生成m行n列全零矩阵.
ones(m,n):生成全1矩阵.
rand(m,n):生成随机矩阵.
rand:生成一个随机数.
diag(A):取A的对角线元素.
tril(A):取A的下三角元素.
triu(A):取A的上三角元素.
hilb(n):生成 n 维 Hilbert 矩阵.
randn (n):产生均值为0,方差为1的标准正态分布随机矩阵.
vander (V):生成以向量V为基础向量的范得蒙矩阵.
invhilb(n):求n阶的希尔伯特矩阵的逆矩阵.
toeplitz(x, y):生成一个以x为第一列,y为第一行的托普利兹矩阵
compan(p):生成伴随矩阵，p是一个多项式的系数向量，高次幕系数排在前，低次幕排在后.
pascal (n):生成一个n阶帕斯卡矩阵.
compan:生成伴随矩阵
⑶矩阵运算
MATLAB的基本算术运算有:+(加)、一(减)、*(乘)、/(右除)、'(左除)、
"(乘方).
加法：
» A+B
ans =
9 3 9
7 10 13
6 12 7
减法：
» B-A
ans =
7-13
-10 1
2 6-3
乘法：
» A*B
ans =
26
38
26
71
83
71
45
62
43
除法：
» magic(3)/hilb(3)
ans =
+OO3 *

-

-

-

-
在MATLAB中，有一种特殊的运算，因为其运算符是在有关算术运算符前面加点，.*、./、.'和两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维参数相同.
» A. *B
ans =
8
2
18
12
25
42
8
27
10
MATLAB提供了6种关系运算符：〈（小于）、〈=（小于或等于）、＞（大于）、＞=（大
于或等于）、
=（等于）、
〜=（不等于）.
» A>B
ans =
0
1
0
1
0
0
0
0
1
MATLAB提供了3种逻辑运算符：& （与）、|（或）和〜（非）.
在逻辑运算中，确认非零元素为真，用1表示，零元素为假，