数值分析实验二 解线性方程组

一、实验目的

1.了解求线性方程组的直接法的有关理论和方法；会编制列主元消去法的程序；
2. 熟悉迭代法的有关理论和方法；会编制雅可比迭代法；
3. 掌握所用迭代方法的收敛性及其收敛速度问题。

二、实验题目

2.对于第1题的弹道轨迹抛物线方程，根据已知的三个弹道轨迹经过点（1，6）、（3，5）、（7，2）建立的线性方程组，采用迭代法求解，得到此弹道轨迹抛物线方程。
采用雅可比迭代法求其近似解，并画出解随着迭代次数变化的趋势线（允许的最大迭代次数N，近似解

三、实验原理

1、高斯列主元消去法的精度eps，由用户设定）。（50分）

四、实验内容与结果

实验题目1 (列主元高斯消元法)

1. 程序源码

function [augmentedMatrix,ans] = ColMain(A,b)
%COLMAIN 此处显示有关此函数的摘要
%   此处显示详细说明
 
%常量定义
ZERO = 10^(-15);
 
%增广矩阵
augmentedMatrix = [A b];
 
 
[row_a,col_a] = size(A);
[row_b,col_b] = size(b);
 
%解
ans = zeros(col_a,1);
 
%先判断A是不是方阵
if row_a ~=col_a
    sprintf('矩阵A不是方阵')
    return;
end
 
%判断b是否符合要求
if row_b ~= row_a
    sprintf('b的格式为列向量且b的行数和A的行数要相同')
    return;
end
 
%判断b的行和A的行是否一致
if row_a ~= row_b
    sprintf('向量b长度必须和矩阵A的行数相等')
    return;
end
 
%进行消元
for col = 1:col_a-1
    %找每一列的最大值
    maxCol = max