1024手机基地看电影,午夜福利视频导航,国产精品福利在线一区,亚洲欧美日韩另类成人,在线观看午夜日本理论片,成年超爽免费网站,国产精品成人免费,精品动作一级毛片,成人免费观看网站,97精品伊人久久大香蕉

設(shè)為首頁(yè)收藏本站

開(kāi)啟輔助訪問(wèn) 切換到寬版

快捷導(dǎo)航

»論壇 › 電子技術(shù)分類(lèi)討論區(qū) › 資料共享 › 高斯消元法

返回列表

查看: 1249|回復(fù): 0

高斯消元法

[復(fù)制鏈接]

電梯直達(dá)

樓主

ID:334306 發(fā)表于 2018-5-20 17:03 | 只看該作者回帖獎(jiǎng)勵(lì)

|倒序?yàn)g覽 |閱讀模式

function x=gauss_elim(a,b)
% Gauss 消去法
[n,m]=size(a);
nb=length(b);
det=1; %存儲(chǔ)行列式值
x=zeros(n,1);
for k=1:n-1
for i=k+1:n
      if a(k,k)==0
      return
      end
      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=det*a(k,k); %計(jì)算行列式
end
det=det*a(n,n);
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