| 查看: 1177 | 回复: 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 |
回复此楼» 猜你喜欢
filecode
已经有11人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有33人回复
-
FileCode能看出啥？
已经有19人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
用fortran编的用共轭梯度法解决线性方程组的问题。
已经有5人回复
- 关于 共轭梯度法，类牛顿法中 梯度符号判断的疑问
已经有4人回复
- 应用在分子弛豫时，共轭梯度法（CG）如何编程实现？
已经有7人回复
- 并行共轭梯度算法的C语言实现
已经有9人回复
- 关于 共轭梯度法 求助～～～好心人帮忙 有金币的哦
已经有11人回复
- 【讨论】预调件共轭梯度法（PCG）
已经有50人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
