MATLAB##高斯消元##G-S迭代法##二分法##牛顿法

时间： 2021年5月26日21:48:57

上机练习02：利用MATLAB程序完成线性方程组数值求解（高斯消元法、G-S迭代法）、非线性方程求解（二分法、牛顿法），共4种计算方法，需完成4个程序，具体题目可从教材中选取。

2. 计算内容

高斯消元法：

function x=colGauss(A,b)
n=length(b);%取数组b的长度
k=1 ;
A (1:n, n+1)=b ;%将数组b放入到数组A的最后一列
while k<=n-1%while循环开始
b=abs(A(k:n,k)) ;%取A数组中第一列取绝对值赋给数组b
[cmax,ik ]=max(b);%取数组b中最大值赋给cmax，并将最大值所在行索引号赋给ik
if ik>k
for i=k:n+1%将数组A中每一列中绝对值最大的数作为主元
t=A(k,i) ;
A(k, i)=A(ik, i) ;
A(ik, i)=t ;
end
end
for i=k+1:n%消元过程
t=A(i,k)/A(k,k);
A(i,k:n+1)=A(i,k:n+1)-t*A(k, k:n+1);
end
k=k+1 ;
end %while循环结束
x (n)=A(n, n+1)/A(n, n) ;%求解最后结果
for k=n-1:-1 : 1