应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 编一个fortran程序
181/2返回列表12下一页
查看: 2679  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

宇天旅梦

铜虫 (小有名气)

[求助] 编一个fortran程序

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

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-08-24 18:46:41
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 宇天旅梦 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭