经典数学实验 4个 MTALAB 二次曲线拟合 迭代法解线性方程组

function gaosi %本程序解决输入任意矩阵求解其根 n=input('请输入未知数的个数:'); A=zeros(n,n+1); A=input('请输入增广矩阵的 ：'); ans=gs(A,n); disp('则方程的解为：'); if ans==0 disp('error'); else for i=1:n disp(sprintf('x( %d )= %16.16f\n',i,ans(i))); end end function x=gs(A,n) %下面采用列主元高斯消元法求解线性方程 for k=1:n-1 r=k; %找出列元中的绝对值最大值 for i=k:n if abs(A(i,k))>abs(A(r,k)) r=i; end end if abs(A(r,k))<eps disp('No unique slution!'); exit (0); end if r~=k %r,k进行换行 for j=k:n+1 c=A(k,j); A(k,j)=A(r,j); A(r,j)=c; end end %进行消元计算 for i=k+1:n c=A(i,k)/A(k,k); for j=k+1:n+1 A(i,j)=A(i,j)-c*A(k,j); end end end if abs(A(n,n))<eps disp('No slution!'); exit (0); end %进行回代求解 for k=n:-1:1 x(k)=A(k,n+1); for j=k+1:n x(k)=x(k)-A(k,j)*x(j); end x(k)=x(k)/A(k,k); end