回覆列表
  • 1 # 高大美女的

    function [x,XA]=GaussXQByOrder(A,b)

    %高斯順序消元法

    N = size(A);

    n = N(1);

    for i=1:(n-1)

    for j=(i+1):n

    if(A(i,i)==0)

    disp("對角元素為0!"); %防止對角元素為0

    return;

    end

    l = A(j,i);

    m = A(i,i);

    A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m; %消元方程

    b(j)=b(j)-l*b(i)/m;

    end

    end

    x=SolveUpTriangle(A,b); %通用的求上三角係數矩陣線性方程組的函式

    XA = A; %消元后的係數矩陣

    function x=SolveUpTriangle(A,b)

    N=size(A);

    n=N(1);

    for i=n:-1:1

    if(i

  • 中秋節和大豐收的關聯?
  • 請問露可以組什麼詞?