| 查看: 1201 | 回复: 0 | ||
[求助]
关于共轭梯度算法（CGM)的一些困惑
|
|
关于图片中的那一条算法。该如何快速的实现呢。 假设一个M*N的网格区域。共有M*N个点。即产生的矩阵qx，qy，p都为M*N size。这些节点共有两种状态，黏着状态stick和滑移状态slip。分别用集合Ist和Isl表示。Tst和Isl由0和1组成。处在黏着区的节点(i,j)满足Ist（i，j）=1，不满足为0. Isl也同理。 我对图片中的操作是这样的： q=abs(sqrt(qx.^2+qy.^2)); Ist(q>mu*p)=0;%q>mu*p时，黏着状态转为滑移状态 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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有13人回复
找到一些相关的精华帖子，希望有用哦~
用fortran编的用共轭梯度法解决线性方程组的问题。
已经有5人回复
- 关于 共轭梯度法，类牛顿法中 梯度符号判断的疑问
已经有4人回复
- 应用在分子弛豫时，共轭梯度法（CG）如何编程实现？
已经有7人回复
- 并行共轭梯度算法的C语言实现
已经有9人回复
- 关于 共轭梯度法 求助～～～好心人帮忙 有金币的哦
已经有11人回复
- 【讨论】预调件共轭梯度法（PCG）
已经有50人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
