全球旧事资料 分类
数调用是不能在一个m文件中再定义的,需要再建立一个m文件,运行时用符号@标识,并为变量赋值。基本上已经熟悉了MATLAB的常用功能。

f实验五
一、实验题目
高斯列主元素消去法解方程组
解方程组:06428x103475x208468x30412703475x118423x204759x31732108468x104759x212147x308621
二、程序功能
用高斯列主元素消去法解方程组
三、算法
选主元即从增广矩阵第一列中选择绝对值最大的数所在的行与第一行互换,并将新矩阵的第二行和第三行的第一列数化为0,形成新的矩阵;步2:从新矩阵的第二行和第三行的第二列中选择绝对值最大的换到第二行,第三行第二列数化为0,形成新矩阵;步3:回代,先求出x3回代求出x2进而可得x1步1:
四、重要标识符说明A,b:系数矩阵N:b的长度X:一个长度为
,值为0的列向量C:长度为
,值为0的横向量d1:数0Aii:矩阵A中对角线的数Aji:矩阵A中第j行第i列的元素
zeros
m:生成m
阶0矩阵;abs:求绝对值;
10
f五、程序运行实例
在matlab65环境中,运行程序,结果如图:
六、源程序fu
ctio
xGaussAbA064280347508468034751842304759084680475912147b041271732108621
le
gthbXzeros
1czeros1
d10fori1
1maxabsAiimiforji1
ifmaxabsAjimaxabsAjimje
de
difmiforki
ckAik
11
fAikAmkAmkcke
dd1bibibmbmd1e
dforki1
forji1
AkjAkjAijAkiAiie
dbkbkbiAkiAiiAki0e
de
d回代求解回代求解回代求解x
b
A
fori
111sum0forji1
sumsumAijxje
dxibisumAiie
d七、个人实验总结
实验代码很长,其实到最后我也不是很明白,尤其是那些嵌套的循环语句,还好老师给了试验源代码,我们只是在此基础上稍作了些修改,最后得出了正确答案。经过五个实验的运行,已经对MATLAB环境有了进一步的了解,以后若有用到的地方也不会无从下手。
12
fr
好听全球资料 返回顶部