[問題] 請問我程式哪裡出問題已回收
這是作業
是要我們寫LU程式
我是GaussNaive的M檔去做修改
程式碼如下
function [U,L] = ul(A)
[m,n] = size(A);
L=eye(n);
if m~=n, error('Matrix A must be square'); end
nb = n;
Aug = [A];
% forward elimination
for k = 1:n-1
for i = k+1:n
factor = Aug(i,k)/Aug(k,k);
Aug(i,k:nb) = Aug(i,k:nb)- factor*Aug(k,k:nb);
L(i,k) = factor;
end
end
U = Aug
==
結果出現一些問題
LU的矩陣相乘結果不是原矩陣
秀給老師看 (但可能是他趕時間吧)
老師他也看不出哪出問題
--
哇!73級分耶!我的孩子真不是蓋的。去讀台大吧!
什麼?你想讀XX系?這樣就上不了台大了。
什麼興趣?興趣不重要啦!你老子我的面子才重要。
我生你養你讓你去補習考上大學,你就是要聽我的。
想想我跟人家說:"我孩子是讀台大的" 多驕傲啊!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 163.25.118.130
推
04/04 12:56, , 1F
04/04 12:56, 1F
→
04/04 12:57, , 2F
04/04 12:57, 2F
→
04/04 12:59, , 3F
04/04 12:59, 3F
→
04/05 18:06, , 4F
04/05 18:06, 4F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章