数值分析及MATLAB代码实现

实验二（1） 线性方程组的直接解法

实验任务：高斯顺序消元法、列选主元消元法、三角分解法

实验目的： 掌握高斯顺序消元法、列选主元消元法、三角分解法原理和应用条

件，利用 MATLAB 程序实现算法

实验课时：4 学时

预备学习：

1、 P18-20 程序 2.1 和程序 2.2 的应用举例；

2、 P27 程序 2.3 的应用举例.

n=length(b);a=[a,b];

for k=1:(n-1)

a((k+1):n,(k+1):(n+1))=a((k+1):n,(k+1):(n+1))-a((k+1):n,k)/a(k,k)*a(k,(k+1):

(n+1));

a((k+1):n,k)=zeros(n-k,1);

if flag==0,a,end

end

x=zeros(n,1);

x(n)=a(n,n+1)/a(n,n);

for k=n-1:-1:1