 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (3346)
- >導(dǎo)師招生 (276)
- >基金申請(qǐng) (241)
- >考研 (233)
- >文獻(xiàn)求助 (127)
- >蟲(chóng)友互識(shí) (112)
- >考博 (101)
- >休閑灌水 (45)
- >招聘信息布告欄 (27)
- >論文投稿 (26)
- >論文道賀祈福 (25)
- >博后之家 (19)
- >碩博家園 (17)
- >找工作 (13)
- >教師之家 (11)
- >SciFinder/Reaxys (10)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1114 | 回復(fù): 3 | |||
lijunjiexm鐵蟲(chóng) (小有名氣)
|
[求助]
求助一簡(jiǎn)單matlab程序的出錯(cuò)問(wèn)題，謝謝
|
|
如下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)); |
回復(fù)此樓» 猜你喜歡
焦慮
已經(jīng)有13人回復(fù)
-
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑
已經(jīng)有8人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有9人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有5人回復(fù)
-
285化工學(xué)碩求調(diào)劑（081700）
已經(jīng)有11人回復(fù)
-
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑
已經(jīng)有9人回復(fù)
-
274求調(diào)劑
已經(jīng)有6人回復(fù)
-
求材料調(diào)劑
已經(jīng)有10人回復(fù)
-
本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑
已經(jīng)有4人回復(fù)
-
國(guó)自科面上基金字體
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求matlab求解一個(gè)方程組的問(wèn)題
已經(jīng)有6人回復(fù)
- 求助使用MATLAB計(jì)算綠色葉片的覆蓋度的問(wèn)題
已經(jīng)有24人回復(fù)
- 求助Matlab解方程問(wèn)題
已經(jīng)有6人回復(fù)
- 【求助】弱弱的問(wèn)一下Matlab數(shù)據(jù)提取的問(wèn)題
已經(jīng)有5人回復(fù)
- 【求助】初學(xué)matlab 請(qǐng)教一個(gè)很簡(jiǎn)單的問(wèn)題 非常感謝高手們
已經(jīng)有3人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題
已經(jīng)有21人回復(fù)
- 【求助】大家看看，matlab安裝過(guò)程的截圖，找出錯(cuò)
已經(jīng)有14人回復(fù)
- 【求助】MATLAB多個(gè)mat中 數(shù)據(jù)畫(huà)圖問(wèn)題
已經(jīng)有7人回復(fù)
- 【求助】問(wèn)一個(gè)關(guān)于matlab中l(wèi)egend的問(wèn)題
已經(jīng)有11人回復(fù)
libralibra
至尊木蟲(chóng) (著名寫(xiě)手)
驃騎將軍
- 程序強(qiáng)帖: 40
- 應(yīng)助: 817 (博后)
- 金幣: 12914.1
- 紅花: 64
- 帖子: 2238
- 在線: 287.3小時(shí)
- 蟲(chóng)號(hào): 696514
- 注冊(cè): 2009-02-05
- 專(zhuān)業(yè): 計(jì)算機(jī)軟件
【答案】應(yīng)助回帖
★ ★
余澤成(金幣+2): 謝謝參與應(yīng)助！ 2011-11-26 23:43:22
余澤成(金幣+2): 謝謝參與應(yīng)助！ 2011-11-26 23:43:22
|
??? Index exceeds matrix dimensions. 下標(biāo)越界錯(cuò)誤 在mesh2d2中的第八行之前,第七行之后的分號(hào)去掉 運(yùn)行一次,把keep這個(gè)下標(biāo)打印出來(lái)看,是不是沒(méi)有find到,keep是0啊 導(dǎo)致在做稀疏矩陣A的時(shí)候ii(keep)訪問(wèn)矩陣的0號(hào)元素(不存在),因?yàn)閙atlab的矩陣下標(biāo)從1開(kāi)始的 |
2樓2011-11-26 17:53:51
【答案】應(yīng)助回帖
★ ★
余澤成(金幣+2): 謝謝參與應(yīng)助！ 2011-11-26 23:43:34
余澤成(金幣+2): 謝謝參與應(yīng)助！ 2011-11-26 23:43:34
|
把最后一行“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è)來(lái)吧。 另外，你xx = repmat ([-1 -14-1-1]',1, n^2) ;是想生成2*100個(gè)元素的矩陣，不過(guò)按照你的程序應(yīng)該生成5*100矩陣吧。猜測(cè)而已。[-1 -14-1-1]應(yīng)為[-1 -1 4 -1 -1]吧？？ |
3樓2011-11-26 23:11:33
【答案】應(yīng)助回帖
★ ★
余澤成(金幣+2): 謝謝參與應(yīng)助！ 2011-11-26 23:43:44
余澤成(金幣+2): 謝謝參與應(yīng)助！ 2011-11-26 23:43:44
|
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樓2011-11-26 23:20:31
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 lijunjiexm 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
