线性代数实验：矩阵求解方法比较与应用.pdfVIP

实验九线性代数实验

班级：统计0502学号：1303050206：欧钊锋

实验目的：

掌握MATLAB求矩阵的秩、行列式和逆矩阵令

学会用逆矩阵法、左除法、克拉默法则求解方程Ax=B并通过比较知道哪种方法用时最

少，最精确

学会运用以上的方法求解给定的恰定，超定和欠定线性方程组

实验理论：

m*n线性方程组Ax=B的求解是用矩阵除来完成的，X=A\B，当m=n且A可逆时，给出唯

一解。这时矩阵除A\B相当于inv(A)*B;当nm时，矩阵除给出方程的最小二乘解；当nm

时，矩阵除给出方程的最小范数解。恰定方程组Ax＝b的求解比较简单．一般可用两种方

法：一种是利用逆矩阵求解x＝inv(A)b；另一种是用除法求解X＝A＼b．两种方法的异同

点是：算法上都采用Guass消去法，但用除法求解时，无需求A的逆，这样可以很好地保

证求解时的计算精度，还能节省大量的计算时间．当然也可以用Cramer法则求解方程组。

由于超定方程组无解，而在实际应用中，能求得其最小二乘解也是有意义的，这时，方程组

的求解仍可用除法和广义逆矩阵法，不过这样求得的解不会满足Ax＝b，而是其最小二乘意

TT

义下的解，即解x＝inv(AA)Ab．从理论上说，欠定方程组有无穷多组解。这时，我们仍

然可用广义逆矩阵法和矩阵除法求解方程组，但前者所求的解是所有解中范数最小的一个，

而后者所求的解是所有解中含零个数最多的一个，因而两种解法的结果一般会有所不同。

但它们求出的解是所有解系中的一个特解，要求其通解，则先解其对应的齐次方程的通解，

然后再加上自身的一个特解即可。值得注意的是，用除法求解欠定方程组时，会提出警告，

并给出系数矩阵的秩，但警告并不会影响方程求解。

实验内容：

（一）比较用逆矩阵法、左除法、克拉默法则求解方程Ax=b的用时和误差

程序设计：

%求解用逆矩阵法、左除法、克拉默法则求解方程Ax=b的用时和误差，

%并相互比较，可知左除法最方便

A=rand(100)*1.e2;

x=ones(100,1);

b=A*x;

tic

y=inv(A)*b;

toc

err=norm(y-x)

res=norm(A*y-b)

tic

y=A\b;

toc

Experiment9LinearAlgebraExperiment

Class:Statistics0502StudentID:1303050206Name:OuZhaofeng

Experiment

purpose：

MastertheMATLABcommandsforfindingtherank,determinantandinversematrixofamatrix

Learntousetheinversematrixmethod,leftdivisionmethod,andCramersruletosolvetheequationAx=

B,andknowwhichmethodtakestheleasttimeandisthemostaccuratethroughcomparison

Learntousetheabovemethodstosolvegivendefinite,overdeterminedand

underdeterminedlinearequations

Experimental

Theory：

Thesolutionofthem*nlinearequationsystemAxBiscompletedbymatrixdivision,XA\B,whenm=

nandAarereversible,theuniquesolutionisgiven.Atthistime,matr