| 查看: 1103 | 回复: 7 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】c++编程，请高手帮我看看这个程序应该怎么写
|
|||
|
程序中的公式为： p[I][j]=(p[j]+p[i-1][j])/2; s[I][j]=-2*miu*h*p[I][j]*p[I][j]*(1-p[I][j])/(4*dx*4*dx); fw[i-1][j]=((k[j]+k[i-1][j])/2*v[I][j]+(k[j]+k[i-2][j])/2*v[I-1][j])/2; fe[j]=((k[i+1][j]+k[j])/2*v[I+1][j]+(k[j]+k[i-1][j])/2*v[I][j])/2; aw[I-1][j]=dw[i-1][j]+fw[i-1][j]/2; ae[I+1][j]=de[j]-fe[j]/2; 在这个程序中包含i，j，I,J。我可以像下边这样写程序吗？ for(I=3;I<=m;I++) for(J=3;J<=n;J++) { for(i=2;i<=m;i++) for(j=2;j<=n;j++) { p[I][j]=(p[j]+p[i-1][j])/2; s[I][j]=-2*miu*h*p[I][j]*p[I][j]*(1-p[I][j])/(4*dx*4*dx); fw[i-1][j]=((k[j]+k[i-1][j])/2*v[I][j]+(k[j]+k[i-2][j])/2*v[I-1][j])/2; fe[j]=((k[i+1][j]+k[j])/2*v[I+1][j]+(k[j]+k[i-1][j])/2*v[I][j])/2; aw[I-1][j]=dw[i-1][j]+fw[i-1][j]/2; ae[I+1][j]=de[j]-fe[j]/2; } } [ Last edited by 飘雪139 on 2008-3-14 at 20:47 ] |
回复此楼» 猜你喜欢
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
时间戳又变了8-15
已经有15人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有9人回复
2楼2008-03-15 08:53:11
3楼2008-03-15 17:46:06
|
for(I=3;I<=m;I++) for(J=3;J<=n;J++) { for(i=2;i<=m;i++) for(j=2;j<=n;j++) { p[I][j]=(p[j]+p[i-1][j])/2; s[I][j]=-2*miu*h*p[I][j]*p[I][j]*(1-p[I][j])/(4*dx*4*dx); dw[i-1][j]=miu*k[i-1][j]/dx; de[j]=miu*k[j]/dx; fw[i-1][j]=((k[j]+k[i-1][j])/2*v[I][j]+(k[j]+k[i-2][j])/2*v[I-1][j])/2； fe[j]=((k[i+1][j]+k[j])/2*v[I+1][j]+(k[j]+k[i-1][j])/2*v[I][j])/2; aw[I-1][j]=dw[i-1][j]+fw[i-1][j]/2; ae[I+1][j]=de[j]-fe[j]/2; au[I][j]=(k[j]+k[i-1][j])/2; } } |
4楼2008-03-15 17:54:57
zsuzwr
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 贵宾: 0.15
- 金币: 2811.7
- 散金: 10
- 帖子: 554
- 在线: 275小时
- 虫号: 209571
- 注册: 2006-03-06
- 专业: 凝聚态物性I:结构、力学和
5楼2008-03-15 18:02:51
6楼2008-03-15 18:09:06
sun20082002
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2392.3
- 帖子: 50
- 在线: 73.6小时
- 虫号: 218450
- 注册: 2006-03-13
- 性别: GG
- 专业: 金属材料的凝固与结晶学
|
网络不好，发重了，希望版主给删除一个。 先看一点点 for(I=3;I<=m;I++) for(J=3;J<=n;J++) { for(i=2;i<=m;i++) for(j=2;j<=n;j++) { p[I][j]=(p[j]+p[i-1][j])/2; …… } p[I][j] 的值最后都是 (p[m][n]+p[m-1][n])/2, 应该不对吧？ 而且在第二次 ij 循环中计算多次。 p[I][j] 与 p[j] 的本质没变，还是同一个数组，你用 I 并不代表它会改变。 对于这样的循环 要不重新找一个数组存放数值，要不 只用 一个 i j 的循环，在内部处理好I 与i 的关系。 [ Last edited by sun20082002 on 2008-5-15 at 10:25 ] |
7楼2008-05-15 10:17:32
sun20082002
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2392.3
- 帖子: 50
- 在线: 73.6小时
- 虫号: 218450
- 注册: 2006-03-13
- 性别: GG
- 专业: 金属材料的凝固与结晶学
★ ★ ★
lei0736(金币+3,VIP+0):谢谢分析
lei0736(金币+3,VIP+0):谢谢分析
|
先看一点点 for(I=3;I<=m;I++) for(J=3;J<=n;J++) { for(i=2;i<=m;i++) for(j=2;j<=n;j++) { p[I][j]=(p[j]+p[i-1][j])/2; …… } p[I][j] 的值最后都是 (p[m][n]+p[m-1][n])/2, 应该不对吧？ 而且在第二次 ij 循环中计算多次。 p[I][j] 与 p[j] 的本质没变，还是同一个数组，你用 I 并不代表它会改变。 对于这样的循环 要不重新找一个数组存放数值，要不 只用 一个 i j 的循环，在内部处理好I 与i 的关系。 |
8楼2008-05-15 10:21:51
