24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 803  |  回复: 11
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

icecream8508

金虫 (著名写手)

[交流] 【求助】新人请教fortran编程问题![求助完成]

刚接触Fortran不久。因为自己求解的东西需要迭代,而自己对迭代的认识就是附初值,然后更新每次的值,直到收敛。大致思路也该如此,但是在程序实现的时候,总是不对的。
因此,今天就想从最简单的迭代开始算起。想先把最基本的一个迭代的例子走通,但是结果仍然不对。
也问了几个人,但是大家对于程序的编写都有自己习惯的语言,觉得问完之后,更茫然了。
在检查程序的时候,不知道错在哪里了?
不知,大家在审查程序时有什么好方法吗?

[ Last edited by wangmonk on 2009-6-7 at 20:44 ]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

具体的错误以及问题,可以帮你解决,某些问题,不可言传。
4楼2009-06-03 21:39:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

icecream8508

金虫 (著名写手)

我今天做的简单例子的编程附着如下:
!******* 练习迭代算法*******

!******** 所需要求解的方程******
!  x(1)=0.2*x(2)+0.1*x(3)+0.3
!   x(2)=0.2*x(1)+0.1*x(3)+1.5
!   x(3)=0.2*x(1)+0.1*x(2)+2
  
!******* 1.jACOBI 迭代*******
program main

real B(3)
real A(3,3),x(3),s(30)
integer N,k
integer m
real(kind=8) :: eps
real(kind=8) :: xo
integer,parameter:: NM=1000

eps=1.0e-6
!***********系数矩阵*************
B(1)=0.3
B(2)=1.5
B(3)=2

A(1,2)=0.2
A(1,3)=0.1
A(2,1)=0.2
A(2,3)=0.1
A(3,1)=0.2
A(3,2)=0.1

!***************迭代公式************
! do k=1,NM,1
!Do i=1,3
!   x(i)=0.0
!end do

20 Do k=1,NM,1

Do i=1,3
  s(i)=0.0  
  Xo=x(i)
  x(i)=0.0
   Do j=1,3
      If(i.NE.j) then
!       S(i)=S(i)+A(i,j)*Xo
       S(i)=S(i)+A(i,j)*X(j)
      End if
   End do
x(i)=B(i)-s(i)

If(abs(xo-x(i)).gt.eps) then
Go to 20
else
  exit
end if

end do
end do

!write(*,*) (x(i),i=1,3)
do j=1,3
  write(*,*) x(j)
end do

end program main
2楼2009-06-03 11:50:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
你是想知道问题出现在哪里,还是想知道怎么找错误,
这个不是在地图上找某个城市那么公式化和简单化,有一般的调试流程,但有很大的经验成份,还有对该程序的熟悉程度以及算法,程序理论依据的熟悉程度。不是你想象的那么简单。
3楼2009-06-03 21:38:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

icecream8508

金虫 (著名写手)

先谢谢楼上的。
一是先想找错误
二是,想知道怎么找错误的?自己已经检查了好几遍了。这个还是一个初步的迭代。我自己的问题要比这复杂的多,头大。

再者,对于我们刚接触的fortran来说,应该怎么做,最好哪?!
5楼2009-06-03 22:28:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 300求调剂,材料科学英一数二 +8 leaflight 2026-03-24 8/400 2026-03-29 01:31 by fmesaito
[考研] 调剂求院校招收 +5 鹤鲸鸽 2026-03-28 5/250 2026-03-29 01:10 by 我是小康
[考研] 各位老师好,我的一志愿为北京科技大学085601材料专硕 +7 Koxui 2026-03-28 7/350 2026-03-29 00:58 by 我是小康
[考研] 086000生物与医药调剂 +3 Feisty。 2026-03-28 7/350 2026-03-29 00:41 by 544594351
[考研] 295求调剂 +4 wei-5 2026-03-26 4/200 2026-03-28 23:20 by 小木虫tim
[考研] 学硕274求调剂 +9 Li李鱼 2026-03-26 9/450 2026-03-28 21:42 by bymhappy
[考研] 0703本科郑州大学求调剂 +3 nhj_ 2026-03-25 3/150 2026-03-28 13:24 by Iveryant
[考研] 材料与化工(0856)304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 331环境科学与工程求调剂 +3 熠然好运气 2026-03-27 3/150 2026-03-28 04:11 by fmesaito
[考研] 285求调剂 +4 AZMK 2026-03-27 7/350 2026-03-27 20:59 by AZMK
[考研] 324求调剂 +5 hanamiko 2026-03-26 5/250 2026-03-27 10:33 by wangjy2002
[考研] 284求调剂 +11 junqihahaha 2026-03-26 12/600 2026-03-27 04:37 by wxiongid
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 一志愿 南京邮电大学 288分 材料考研 求调剂 +3 jl0720 2026-03-26 3/150 2026-03-26 13:39 by zzll406
[考研] 化学调剂一志愿上海交通大学336分-本科上海211 +4 小鱼爱有机 2026-03-25 4/200 2026-03-26 10:19 by aa331100
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +6 jiajunX 2026-03-22 6/300 2026-03-25 23:05 by licg0208
[考研] 生物技术与工程 +3 1294608413 2026-03-25 4/200 2026-03-25 18:02 by 1294608413
[考研] 网络空间安全0839招调剂 +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
信息提示
请填处理意见