懒惰是很奇怪的东西，它使你以为那是安逸，是休息，是福气；但实际上它所给你的是无聊，是倦怠，是消沉它剥夺你对前途的希望，割断你和别人之间的友情，使你心胸日渐狭窄，对人生也越来越怀疑。罗兰
计算方法上机实验报告1实验目的：（1）解线性方程组左除命令‘’的学习（2）自学与矩阵计算有关的一些常用函数（3）xAb与xi
vAb在耗时方面的区别（4）矩阵分解命令的学习（5）病态矩阵试验（6）SOR迭代法收敛速度受松驰因子的影响试验2实验内容：【实验一】自学与矩阵计算有关的一些常用函数（1）detAi
vAra
kArrefA求方阵A的行列式求方阵的逆求矩阵A的数值秩化矩阵A为最简阶梯形

ormAp求矩阵A的p范数co
dAp求矩阵A的p范数的条件数eigA求方阵A的特征值
【实验二】xAb与xi
vAb在耗时方面的区别解方程组Axb时尽量不要使用xi
vAb而要使用xAb下面是二者耗时方面的区别实验把下面程序拷贝为新的M文件有些命令可不关心其含义
fra
dstate0Agalleryra
dsvd2002e132产生条件数为2e13的200阶的随机矩阵xo
es2001bAxformatlo
g用15位小数显示设精确解为111
求逆法ticx1i
vAbtime1toc关闭计时器并显示耗时启动计时器
error1
ormxx1i
f计算最大误差
左除法ticx2Abtime2tocerror2
ormxx2i
f★【实验一】首先完成下面Gauss列选主元的消去法程序单独存为gaussm文件注意一定要与函数名相同
f可参考P43图32然后找一个例子调用此程序验证是否正确调用方法同Matlab内部函数调用完全一样
fu
ctio
xgaussAb
sizeAxzeros
1AugAbfork1
1行rifrk
temp1Aktemp2rkAkArrkrrArtemp1rrtemp2
增广矩阵
pivrmaxabsAugk
k
找列主元所在子矩阵的
rrk1
列主元所在大矩阵的行
对Aug实施行交换一行命令就可以了怎么写e
d
f
ifAugkk0error对角元出现0e
d程序遇到error
会中断执行并显示其中的提示内容e
d解上三角方程组AAug1
bAug
1x
fork
111e
d【实验二】编下面程序1追赶法P452Cholesky分解法P51xke
d把增广矩阵消元成为上三角forpk1
multAugpkAugkkAugpk
1消元乘子
下面是追赶法的程序你可以参考参见P45注意向量a的下标与书上不同
ffu
ctio
dtridiagabcddtridiagabcd求解三对角线性方程组的追赶法r