| 查看: 696 | 回复: 14 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hym5701木虫 (小有名气)
|
[交流]
[求助]再求解下面的线性方程组fortran程序
|
||
|
已知u，v，w，s，d，e怎么求a，b sigma(an*un+bn*vn)=dn sigma(an*wn+bn*sn)=en 不知道怎么粘贴公式，就只能像上面那样写了，sigma是求和符号，n是下标，从0到无穷大。测试通过15个金币送上，谢谢。 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
728728
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 4940.1
- 散金: 239
- 帖子: 519
- 在线: 397.7小时
- 虫号: 92818
- 注册: 2005-09-11
- 性别: GG
- 专业: 可再生与替代能源利用中的
★ ★ ★ ★ ★ ★ ★ ★
hym5701(金币+5,VIP+0):谢谢你的解释，现在我假设右边dn和en是常数的话，能提供解这样的程序吗，我对fortran不是很熟悉，以前学过一点点，这个问题是我同事给我的，他说一定要用fortran解，没办法。你说可以改下，可以帮我改下吗。谢谢
lei0736(金币+3,VIP+0):谢谢
hym5701(金币+5,VIP+0):谢谢你的解释，现在我假设右边dn和en是常数的话，能提供解这样的程序吗，我对fortran不是很熟悉，以前学过一点点，这个问题是我同事给我的，他说一定要用fortran解，没办法。你说可以改下，可以帮我改下吗。谢谢
lei0736(金币+3,VIP+0):谢谢
|
我没有现成的程序。 如果n是一个较大的数的话，假设它的上限是N，那么方程租右侧应该是N而不是n，这样就是一个二元线性方程组。如果我这么理解对的话，那直接先算出a，b的系数（un，vn之类的直接求和），直接用消去法就可以了。 如果右侧确实是n的话，那么就是2N维线性方程组，但自变量只有2个，成了超定方程，这个不能用解方程的方法求解了，只能用最小二乘法或者曲线拟合的方法了。这些算法在Fortran算法集里面都有，稍加改动就可以用，你在论坛上就能找到算法集。 这是我的理解，不知道理解对了没有 |
4楼2008-12-25 15:50:12
728728
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 4940.1
- 散金: 239
- 帖子: 519
- 在线: 397.7小时
- 虫号: 92818
- 注册: 2005-09-11
- 性别: GG
- 专业: 可再生与替代能源利用中的
2楼2008-12-24 17:23:32
hym5701
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4204.3
- 散金: 3
- 帖子: 174
- 在线: 75.8小时
- 虫号: 574617
- 注册: 2008-06-16
- 专业: 凝聚态物性 II ：电子结构
3楼2008-12-24 19:51:04
728728
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 4940.1
- 散金: 239
- 帖子: 519
- 在线: 397.7小时
- 虫号: 92818
- 注册: 2005-09-11
- 性别: GG
- 专业: 可再生与替代能源利用中的
★ ★ ★ ★ ★ ★ ★
coldwind042(金币+2,VIP+0):非常感谢！
lei0736(金币+5,VIP+0):谢谢
coldwind042(金币+2,VIP+0):非常感谢！
lei0736(金币+5,VIP+0):谢谢
|
这么说的话就是第一种理解？那就简单了，程序如下d*MAXN就是dn，常数，以此类推。MaxN是N的上限。 Integer N,MaxN Real U,V,W,S,D,E,a,b Real un,vn,wn,sn Real uu,vv,ww,ss,NN,deno Read (*,*) MAXN,U,V,W,S,D,E un = 0 vn = 0 wn = 0 sn = 0 Do N = 1,MAXN NN = N*N un = un + U*NN vn = vn + V*NN wn = wn + W*NN sn = sn + S*NN End do deno = un*sn - wn*vn If (abs(deno) >= 0.000001 ) then a = (sn*D*MAXN - vn*E*MAXN) / deno b = (D*MAXN - un*a) / vn write(*,*) "a = ",a," b = ",b Else write(*,*) "please check the input data!" End if End |
5楼2008-12-25 18:36:02