| 查看: 643 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】请教差分得到的代数方程组如何C语言实现
|
|||
|
rt，有没有这样的模板？ 我的问题具体是这样的，如下图一个电势和电流密度的差分后的形式（由电势泊松方程及电流连续性方程推导得到），怎么样写成C描述呢？ 说明：计算区域的网格划分方式，ax、ay分别为x轴和y轴方向上的网格步长，i、j为节点坐标序号。这个网格步长在C程序中如何表示？i和j又如何处理呢？ 谢谢！ [ Last edited by 夏之篮 on 2009-12-29 at 16:56 ] |
回复此楼» 猜你喜欢
实验室接单子
已经有3人回复
-
要不要辞职读博？
已经有5人回复
-
不自信的我
已经有10人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
hakekill
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1709.3
- 帖子: 131
- 在线: 19小时
- 虫号: 246017
- 注册: 2006-04-29
- 性别: GG
- 专业: 水环境与生态水利
★ ★ ★ ★ ★ ★ ★ ★
nono2009(金币+2,VIP+0):谢谢应助！ 12-29 20:36
夏之篮(金币+6,VIP+0):谢谢，基本完成我的提问，不好意思，留俩回答我接下来的问题吧，不排除增加悬赏哈~ 12-29 21:24
nono2009(金币+2,VIP+0):谢谢应助！ 12-29 20:36
夏之篮(金币+6,VIP+0):谢谢，基本完成我的提问，不好意思，留俩回答我接下来的问题吧，不排除增加悬赏哈~ 12-29 21:24
|
这个都是通过定义变量来完成的，如网格步长都取0.5 ax=0.5; ay=0.5; i，j分别是你的网格节点数，假如网格节点为4x4的 rows = 4; cols = 4; float u[3,3], rho[3, 3]; max_iter = 10; for(k=1, k for(j = 0, j < cols, j++){ ax2 = ax * ax; ay2 = ay * ay; epsilon = ??; u[i,j]=1/2*(ax2+ay2)*(ax2*(u[i+1,j]+u[i-1,j])+ay2*(u[i,j+1]+u[i,j-1])+ax2*ay2*rho[i,j]/epsilon); .. } } } 大致情况是这样的 [ Last edited by hakekill on 2009-12-29 at 20:07 ] |
4楼2009-12-29 20:05:31
hakekill
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1709.3
- 帖子: 131
- 在线: 19小时
- 虫号: 246017
- 注册: 2006-04-29
- 性别: GG
- 专业: 水环境与生态水利
2楼2009-12-29 16:08:16
3楼2009-12-29 19:48:18
5楼2009-12-29 21:26:24