| 查看: 1139 | 回复: 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 |
回复此楼» 猜你喜欢
265求调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有21人回复
-
295材料求调剂，一志愿武汉理工085601专硕
已经有4人回复
-
0856调剂，是学校就去
已经有9人回复
-
298-一志愿中国农业大学-求调剂
已经有9人回复
-
求调剂
已经有4人回复
-
学校已经提交到NSFC，还能修改吗？
已经有7人回复
-
材料学硕297已过四六级求调剂推荐
已经有6人回复
-
295复试调剂
已经有7人回复
-
本人考085602 化学工程 专硕
已经有20人回复
找到一些相关的精华帖子，希望有用哦~
用fortran编的用共轭梯度法解决线性方程组的问题。
已经有5人回复
- 关于 共轭梯度法，类牛顿法中 梯度符号判断的疑问
已经有4人回复
- 应用在分子弛豫时，共轭梯度法（CG）如何编程实现？
已经有7人回复
- 并行共轭梯度算法的C语言实现
已经有9人回复
- 关于 共轭梯度法 求助～～～好心人帮忙 有金币的哦
已经有11人回复
- 【讨论】预调件共轭梯度法（PCG）
已经有50人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
40