soluline.m

已知两个点要求用插值的方法绘制拟和的曲线 此文件是拟和飞机下降曲线

%用来解线性方程组   index为指标变量，index=0表示计算失败，index=1表示计算成功
%A为方程组的系数矩阵 b为方程组的右端项 x为方程组的解
%de为系数矩阵A的行列式的值
function [x,de,index]=soluline(A,b)
[n,m]=size(A);nb= length(b);
index=1;de=1;x=zeros(n,1);
for k = 1:n-1
    a_max=0;
    %选主元
    for i=k:n
        if abs(A(i,k))>a_max
            a_max=abs(A(i,k));r=i;
        end
    end
   if a_max < 1e-10
       index=0;return;
   end
   %交换俩行
   if r > k
       for j=k:n
           z=A(k,j);A(k,j)=A(r,j);A(r,j)=z;
       end
       z=b(k);b(k)=b(r);b(r)=z;de=-de;
   end
   %消元过程
   for i=k+1:n
       m=A(i,k)/A(k,k);
       for j=k+1:n
           A(i,j)=A(i,j)-m*A(k,j);
       end
       b(i)=b(i)-m*b(k);
   end
   det=de*A(k,k);
end
det=de*A(n,n);
if abs(A(n,n)) < 1e-10
    index=0;return;
end
for k=n:-1:1
    for j=k+1:n
        b(k)=b(k)-A(k,j)*x(j);
    end
    x(k)=b(k)/A(k,k);
end