24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2809  |  回复: 17

宇天旅梦

铜虫 (小有名气)

[求助] 编一个fortran程序

需要一个fortran程序,比较简单的那种,例如f=((5x+6z)*e^(2i) )/(2x-z),其中z=(3y-1)*(x+2)

要求输出的结果是三列数据,分别对应x,y,f .并且每列有20个数。

这三组数据是要导入orgin中进行图像模拟的。
初学fortran,很多东西还不连贯。
请大师们对上面给的式子编个程序示范一下。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

phy_yo

铁虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
open(10,file='data.dat',status='unknown')
real*8 x,y,z,x(20),y(20)
do 1 i=1,20
input (*,*)  x,y
x(i)=x
y(i)=y
1 continue
do 2 1,20
z=(3y(i)-1)*(x(i)+2)
f=((5x+6z)*e^(2i) )/(2x-z)
write(*,*) x(i),y(i),f
write(10,*) x(i),y(i),f
2 continue
pause
end
网络错误 无法获取数据
2楼2013-08-24 20:27:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宇天旅梦

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by phy_yo at 2013-08-24 20:27:45
open(10,file='data.dat',status='unknown')
real*8 x,y,z,x(20),y(20)
do 1 i=1,20
input (*,*)  x,y
x(i)=x
y(i)=y
1 continue
do 2 1,20
z=(3y(i)-1)*(x(i)+2)
f=((5x+6z)*e^(2i) )/(2x-z)
write(*,* ...

式子中的 i 是虚数
3楼2013-08-24 22:35:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phy_yo

铁虫 (小有名气)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 宇天旅梦 at 2013-08-24 22:35:22
式子中的 i 是虚数...

哦 那就把real*8改成complex*16
然后再在前面加一行parameter(i=(0,1))

[ 发自手机版 http://muchong.com/3g ]
网络错误 无法获取数据
4楼2013-08-24 23:14:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宇天旅梦

铜虫 (小有名气)

引用回帖:
4楼: Originally posted by phy_yo at 2013-08-24 23:14:52
哦 那就把real*8改成complex*16
然后再在前面加一行parameter(i=(0,1))
...

运行不了啊,提示有11个错。我用的是fortran90.
下面这些是显示的错误:
--------------------Configuration: Textxcx - Win32 Debug--------------------
Compiling Fortran...
C:\Users\Administrator\Desktop\Textxcx.for
C:\Users\Administrator\Desktop\Textxcx.for(6) : Error: Syntax error, found IDENTIFIER 'X' when expecting one of: ( % . = =>
      input (*,*)  x,y
-------------------^
C:\Users\Administrator\Desktop\Textxcx.for(10) : Error: Syntax error, found INTEGER_CONSTANT '20' when expecting one of: <END-OF-STATEMENT> ; <IDENTIFIER>
      do 2 1,20
-------------^
C:\Users\Administrator\Desktop\Textxcx.for(11) : Error: Syntax error, found IDENTIFIER 'Y' when expecting one of: * :: , <END-OF-STATEMENT> ; : ) (/ [ + - . ] /) ** / > // .LT. < ...
      z=(3y(i)-1)*(x(i)+2)
----------^
C:\Users\Administrator\Desktop\Textxcx.for(12) : Error: Unrecognized token '^' skipped
      f=((5x+6z)*e^(2i) )/(2x-z)
------------------^
C:\Users\Administrator\Desktop\Textxcx.for(12) : Error: Syntax error, found IDENTIFIER 'X' when expecting one of: * :: , <END-OF-STATEMENT> ; : ) (/ [ + - . ] /) ** / > // .LT. < ...
      f=((5x+6z)*e^(2i) )/(2x-z)
-----------^
C:\Users\Administrator\Desktop\Textxcx.for(3) : Error: A specification statement cannot appear in the executable section.
      parameter(i=(0,1))
------^
C:\Users\Administrator\Desktop\Textxcx.for(4) : Error: A specification statement cannot appear in the executable section.
complex*16 x,y,z,x(20),y(20)
--------^
C:\Users\Administrator\Desktop\Textxcx.for(7) : Error: This name has not been declared as an array.   [X]
      x(i)=x
------^
C:\Users\Administrator\Desktop\Textxcx.for(8) : Error: This name has not been declared as an array.   [Y]
      y(i)=y
------^
C:\Users\Administrator\Desktop\Textxcx.for(10) : Error: An unterminated block exists.
      do 2 1,20
^
C:\Users\Administrator\Desktop\Textxcx.for(10) : Error: This label is undefined.   [21]
      do 2 1,20
---------^
Error executing df.exe.

Textxcx.exe - 11 error(s), 0 warning(s)
5楼2013-08-24 23:27:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phy_yo

铁虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
宇天旅梦: 金币+15, ★★★★★最佳答案, 很帮助 2013-08-25 22:31:16
抱歉 之前写得太着急了 所以献丑啦 哈哈
        parameter(n=20)
        integer k
        real*8 x1,x2,y1,y2
        complex*16 z,x(n),y(n),i,f
        open(10,file='data.dat',status='unknown')
        i=(0,1)
        do 1 k=1,20
        print*, "insert xr,","xi,","yr,","yi"
        read*, x1,x2,y1,y2
        x(k)=cmplx(x1,x2)
        y(k)=cmplx(y1,y2)
1        continue
        do 2 k=1,20
        z=(3*y(k)-1)*(x(k)+2)
        f=((5*x(k)+6*z)*zexp(2*i))/(2*x(k)-z)
        write(*,*) x(k),y(k),f
        write(10,*) x(k),y(k),f
2         continue
        pause
        end
以上 fortran77测试后可用  x,y的输入是键盘输入
网络错误 无法获取数据
6楼2013-08-25 09:21:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宇天旅梦

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by phy_yo at 2013-08-25 09:21:29
抱歉 之前写得太着急了 所以献丑啦 哈哈
        parameter(n=20)
        integer k
        real*8 x1,x2,y1,y2
        complex*16 z,x(n),y(n),i,f
        open(10,file='data.dat',status='unknown')
        i=(0,1)
        do 1 k=1,20
        print*, &q ...

嗯,多谢啦。

程序运行时为什么要输入4个数,式子只要x,y啊,
print*, "insert xr,","xi,","yr,","yi"
        read*, x1,x2,y1,y2
        x(k)=cmplx(x1,x2)
        y(k)=cmplx(y1,y2)
这一段为何要这样?

如果x按照从0开始,每运行出一个结果后自动加0.1,
y从0开始,每次加0.01,这样自动输入数据,不用手动,这怎么弄啊。
7楼2013-08-25 10:50:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phy_yo

铁虫 (小有名气)

引用回帖:
7楼: Originally posted by 宇天旅梦 at 2013-08-25 10:50:28
嗯,多谢啦。

程序运行时为什么要输入4个数,式子只要x,y啊,
print*, "insert xr,","xi,","yr,","yi"
        read*, x1,x2,y1,y2
        x(k)=cmplx(x1,x2)
   ...

既然运算式里有虚数符号  x和y也就用虚数形式表示了  当然  如果他是实数  虚部为零就可以了
至于设定均匀步长的自变量  你仔细看看fortran的循环体格式  自己就能写出来
如要x从0到1.9,步长为0.1  则:
        x(1)=0
        do 10 i=1,19
        x(i+1)=x(i)+0.1
10        continue
网络错误 无法获取数据
8楼2013-08-25 14:29:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宇天旅梦

铜虫 (小有名气)

引用回帖:
8楼: Originally posted by phy_yo at 2013-08-25 14:29:45
既然运算式里有虚数符号  x和y也就用虚数形式表示了  当然  如果他是实数  虚部为零就可以了
至于设定均匀步长的自变量  你仔细看看fortran的循环体格式  自己就能写出来
如要x从0到1.9,步长为0.1  则:
        x(1) ...

这个循环现在差不多了。
还有几个小问题:
1、还是不明白为什么要输入4个数呢?
2、real*8 x1,x2,y1,y2
        complex*16 z,x(n),y(n),i,f 中,这里的8和16是什么意思?规定这样吗?
3、open语句中为何要有status='unknown',这又表示什么?

再次感谢你的回答
9楼2013-08-25 17:02:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phy_yo

铁虫 (小有名气)

引用回帖:
9楼: Originally posted by 宇天旅梦 at 2013-08-25 17:02:44
这个循环现在差不多了。
还有几个小问题:
1、还是不明白为什么要输入4个数呢?
2、real*8 x1,x2,y1,y2
        complex*16 z,x(n),y(n),i,f 中,这里的8和16是什么意思?规定这样吗?
3、open语句中为何要有 ...

虚数可以看成两个实数组成 实部一个 虚部一个
real*8 x1,x2,y1,y2 就是定义四个精度为8位的实数 这样就可以组成两个虚数
complex*16 z,x(n),y(n),i,f   就好理解了 就是定义精度为8位的虚数(实质是两个精度为8的实数)
status='unknown' 是规定   新建的文本  没有特别要求  就以这个做声明
网络错误 无法获取数据
10楼2013-08-25 19:05:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 宇天旅梦 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕322分 +8 哈哈哈吼吼吼哈 2026-04-04 8/400 2026-04-05 14:13 by brian0809
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 323分(计算机视觉和大模型项目)能直接上手 +3 chaoxiicy 2026-04-01 3/150 2026-04-05 00:50 by chongya
[考研] 333求调剂 +12 wfh030413@ 2026-04-03 13/650 2026-04-04 21:02 by jj987
[考研] 0703求调剂383分 +8 W55j 2026-04-03 8/400 2026-04-04 20:09 by xhai2011
[考研] 320求调剂 +5 振—TZ 2026-04-02 5/250 2026-04-03 14:42 by fxue1114
[考研] 土木水利328分求调剂 +6 疾风知劲草666 2026-04-02 6/300 2026-04-03 11:38 by znian
[考研] 319求调剂 +18 太容易1018 2026-04-01 18/900 2026-04-03 11:18 by linyelide
[考研] 08工科275分求调剂 +14 AaAa7420 2026-03-31 14/700 2026-04-03 11:13 by cocolv
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-04-02 6/300 2026-04-03 10:19 by 蓝云思雨
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 一志愿北京科技材料科学与工程288分,求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[考研] 346求调剂 +5 郑诚乐 2026-04-02 5/250 2026-04-02 16:38 by SZW_UJN
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 311求调剂 +14 蓝月亮亮 2026-03-30 14/700 2026-04-02 12:18 by 1753564080
[考研] 0710生物学求调剂! +6 叙述文 2026-03-31 6/300 2026-04-01 09:39 by JourneyLucky
[考研] 土木304求调剂 +5 顶级擦擦 2026-03-31 5/250 2026-04-01 08:15 by fdcxdystjk¥
[考研] 269求调剂 +4 我想读研11 2026-03-31 4/200 2026-03-31 10:04 by cal0306
[考研] 吉大生物学326分求调剂 +3 sunnyupup 2026-03-31 3/150 2026-03-31 09:28 by longlotian
信息提示
请填处理意见