| 4 | 1/1 | 返回列表 |
| 查看: 1113 | 回復(fù): 3 | ||
lijunjiexm鐵蟲 (小有名氣)
|
[求助]
求助一簡單matlab程序的出錯(cuò)問題,謝謝
|
|
如下m文件,為什么在命令窗口運(yùn)行 出現(xiàn)錯(cuò)誤提示: >> A=mesh2d2(10) ??? Index exceeds matrix dimensions. Error in ==> mesh2d2 at 8 A = sparse (ii(keep), jj(keep), xx(keep)); >> 請(qǐng)大家指點(diǎn),謝謝! %%%%%%%%%%%%%%%%%%%%%%% function A = mesh2d2 (n) % create an n-by-n 2D mesh for the 2nd difference operator nn = 1:n^2 ; ii = [nn-n ; nn-1 ; nn ; nn+1 ; nn+n] ; jj = repmat (nn, 5, 1) ; xx = repmat ([-1 -14-1-1]',1, n^2) ; keep = find (ii >= 1 & ii <= n^2 & jj >= 1 & jj <= n^2) ; A = sparse (ii(keep), jj(keep), xx(keep)); |
至尊木蟲 (著名寫手)
驃騎將軍
|
??? Index exceeds matrix dimensions. 下標(biāo)越界錯(cuò)誤 在mesh2d2中的第八行之前,第七行之后的分號(hào)去掉 運(yùn)行一次,把keep這個(gè)下標(biāo)打印出來看,是不是沒有find到,keep是0啊 導(dǎo)致在做稀疏矩陣A的時(shí)候ii(keep)訪問矩陣的0號(hào)元素(不存在),因?yàn)閙atlab的矩陣下標(biāo)從1開始的 |

|
把最后一行“A = sparse (ii(keep), jj(keep), xx(keep)); ”轉(zhuǎn)換成 A = sparse (ii(keep)); A = sparse (jj(keep)); A = sparse (xx(keep)); 再試試看。sparse函數(shù)轉(zhuǎn)化只能一個(gè)個(gè)來吧。 另外,你xx = repmat ([-1 -14-1-1]',1, n^2) ;是想生成2*100個(gè)元素的矩陣,不過按照你的程序應(yīng)該生成5*100矩陣吧。猜測而已。[-1 -14-1-1]應(yīng)為[-1 -1 4 -1 -1]吧?? |

|
function A = mesh2d2 (n) % create an n-by-n 2D mesh for the 2nd difference operator nn = 1:n^2 ; ii = [nn-n ; nn-1 ; nn ; nn+1 ; nn+n] ; jj = repmat (nn, 5, 1) ; xx = repmat ([-1 -14-1-1]',1, n^2) ; [keep1,keep2]= find (ii >= 1 & ii <= n^2 & jj >= 1 & jj <= n^2) ; keep=[keep1 keep2]; % A = sparse (ii(keep), jj(keep), xx(keep)); A1 = sparse (ii(keep)); A2 = sparse (jj(keep)); A3 = sparse (xx(keep)); |

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +8 | SUICHILD 2026-03-12 | 8/400 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 314求調(diào)劑 +8 | 無懈可擊的巨人 2026-03-12 | 8/400 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 308 085701 四六級(jí)已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|