[問題] 請問我程式哪裡出問題已回收

看板MATLAB作者 (e^(ix)=COS(x)+iSIN(x))時間16年前 (2009/04/02 19:55), 編輯推噓1(103)
留言4則, 2人參與, 最新討論串1/1
這是作業 是要我們寫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
我剛剛試了一下沒問題阿 該不會在command window打成
04/04 12:56, 1F

04/04 12:57, , 2F
[L U] = ul(A) 吧...@@
04/04 12:57, 2F

04/04 12:59, , 3F
還是因為truncate error?
04/04 12:59, 3F

04/05 18:06, , 4F
我下週上機再試試看 謝謝 ^^
04/05 18:06, 4F
文章代碼(AID): #19rAUyv6 (MATLAB)
文章代碼(AID): #19rAUyv6 (MATLAB)