版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

beckzx

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 36.6
帖子: 41
在線: 15.6小時
蟲號: 3402977
注冊: 2014-09-07
性別: GG
專業(yè): 固體力學(xué)

[求助] 關(guān)于共軛梯度算法（CGM)的一些困惑

關(guān)于圖片中的那一條算法。該如何快速的實現(xiàn)呢。
假設(shè)一個M*N的網(wǎng)格區(qū)域。共有M*N個點。即產(chǎn)生的矩陣qx，qy，p都為M*N size。這些節(jié)點共有兩種狀態(tài)，黏著狀態(tài)stick和滑移狀態(tài)slip。分別用集合Ist和Isl表示。Tst和Isl由0和1組成。處在黏著區(qū)的節(jié)點(i,j)滿足Ist（i，j）=1，不滿足為0.  Isl也同理。
我對圖片中的操作是這樣的：
   q=abs(sqrt(qx.^2+qy.^2));
      Ist(q>mu*p)=0;%q>mu*p時，黏著狀態(tài)轉(zhuǎn)為滑移狀態(tài)
      Isl=Ic-Ist;
   qx(Isl>0)=Isl.*((mu*p.*Isl)./sqrt(((qx.*Isl+1e-6).^2+(qy.*Isl+1e-6).^2)).*(qx.*Isl));
   qy(Isl>0)=Isl.*((mu*p.*Isl)./sqrt(((qx.*Isl+1e-6).^2+(qy.*Isl+1e-6).^2)).*(qy.*Isl));

可是運行，報錯：
qx(Isl>0)=Isl.*((mu*p.*Isl)./sqrt(((qx.*Isl+1e-6).^2+(qy.*Isl+1e-6).^2)).*(qx.*Isl));
Error in test (line 15)

不知道錯在哪里。
還是說有更好的辦法處理這個問題。
請教大神。

111.png

回復(fù)此樓

» 猜你喜歡

307求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑已經(jīng)有15人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有14人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有18人回復(fù)
265求調(diào)劑已經(jīng)有3人回復(fù)
296求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有13人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有15人回復(fù)
生物學(xué)調(diào)劑招人�。。� 已經(jīng)有4人回復(fù)

1樓 2015-01-25 22:41:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

用fortran編的用共軛梯度法解決線性方程組的問題。已經(jīng)有5人回復(fù)
關(guān)于共軛梯度法，類牛頓法中梯度符號判斷的疑問已經(jīng)有4人回復(fù)
應(yīng)用在分子弛豫時，共軛梯度法（CG）如何編程實現(xiàn)？已經(jīng)有7人回復(fù)
并行共軛梯度算法的C語言實現(xiàn) 已經(jīng)有9人回復(fù)
關(guān)于共軛梯度法求助～～～好心人幫忙有金幣的哦已經(jīng)有11人回復(fù)
【討論】預(yù)調(diào)件共軛梯度法（PCG）已經(jīng)有50人回復(fù)

點擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 beckzx 的主題更新