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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 编一个fortran程序
北京石油化工学院2026年研究生招生接收调剂公告
182/2返回列表上一页12
查看: 2809  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

宇天旅梦

铜虫 (小有名气)
引用回帖:
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语句中为何要有 ...

我改写了一下，把x,y的虚数形式改没了。

运行没检查出错，可是运行结果什么都没有，不知道怎么回事，
代码是这样的：
      program xcx
        real*8 x(5),y(5),n
        integer k
        complex*16 i,f
        open(10,file='xcx.dat',status='unknown')
        i=(0,1)
      do 1 k=1,5
         x(1)=0
           y(1)=0
           do 3 n=0.1,0.5
              x(n+1)=x(n)+0.1
            y(n+1)=y(n)+0.1
3           continue
           read*,x,y
1     continue
      do 2 k=1,5
        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(*,*)x(k),y(k),f
2     continue
      stop
        end
赞一下 回复此楼
11楼2013-08-25 19:12:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宇天旅梦

铜虫 (小有名气)
引用回帖:
10楼: Originally posted by phy_yo at 2013-08-25 19:05:52
虚数可以看成两个实数组成 实部一个 虚部一个
real*8 x1,x2,y1,y2 就是定义四个精度为8位的实数 这样就可以组成两个虚数
complex*16 z,x(n),y(n),i,f   就好理解了 就是定义精度为8位的虚数（实质是两个精度为8的 ...

我把x,y的虚数形式换了，改了点东西，运行时没显示错误，可是运行结果没有，弹出来的就只有黑框，不知道怎么回事，代码如下：
      program xcx
        real*8 x(5),y(5),n
        integer k
        complex*16 i,f
        open(10,file='xcx.dat',status='unknown')
        i=(0,1)
      do 1 k=1,5
         do 3 n=0.1,0.5
              x(k+1)=x(k)+0.1
            y(k+1)=y(k)+0.1
3           continue
           read*,x,y
1     continue
      do 2 k=1,5
        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(*,*)x(k),y(k),f
2     continue
      stop
        end
赞一下 回复此楼
12楼2013-08-25 19:16:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phy_yo

铁虫 (小有名气)
引用回帖:
12楼: Originally posted by 宇天旅梦 at 2013-08-25 19:16:10
我把x,y的虚数形式换了，改了点东西，运行时没显示错误，可是运行结果没有，弹出来的就只有黑框，不知道怎么回事，代码如下：
      program xcx
        real*8 x(5),y(5),n
        integer k
        complex*16 i,f
        open(10,f ...

我把x,y的虚数形式换了，改了点东西，运行时没显示错误，可是运行结果没有，弹出来的就只有黑框，不知道怎么回事，代码如下：
      program xcx
        real*8 x(5),y(5)
        integer k，n
        complex*16 i,f
        open(10,file='xcx.dat',status='unknown')
        i=(0,1)
        x(1)=0
        y(2)=0
         do 3 n=1,4               !/数组变量必须是整数才行 你可以想一下为什么这里是4/
            x(n+1)=x(n)+0.1
            y(n+1)=y(n)+0.1
3        continue

      do 2 k=1,5
        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     !/这样数据才会写入xcx.dat/
2     continue
       pause
        end
赞一下 回复此楼
网络错误 无法获取数据
13楼2013-08-25 19:46:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宇天旅梦

铜虫 (小有名气)
引用回帖:
13楼: Originally posted by phy_yo at 2013-08-25 19:46:30
我把x,y的虚数形式换了，改了点东西，运行时没显示错误，可是运行结果没有，弹出来的就只有黑框，不知道怎么回事，代码如下：
      program xcx
        real*8 x(5),y(5)
        integer k，n
        comp ...

显示啦
n=1,4是因为前面已经给出了x(1)=0， y(2)=0，而每个只容纳5个数。
这里改成了y(2)=0？  原来的y(1)=0不对吗？这是任意的吗？

write(10,*)x(k),y(k),f  这里的10 是什么含义？我见过    write(*,10)x,y,f
                                                                       10  format(......)

现在的输出结果怎么才可以弄成x一列，y一列，f一列呢？

DAT文件的打开方式设置为“origin",却没有一点反应，这又怎么导入？
赞一下 回复此楼
14楼2013-08-25 20:12:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phy_yo

铁虫 (小有名气)
引用回帖:
14楼: Originally posted by 宇天旅梦 at 2013-08-25 20:12:50
显示啦
n=1,4是因为前面已经给出了x(1)=0， y(2)=0，而每个只容纳5个数。
这里改成了y(2)=0？  原来的y(1)=0不对吗？这是任意的吗？

write(10,*)x(k),y(k),f  这里的10 是什么含义？我见过    write ...

y(2)=0 是我打错了 应该是y(1)=0 不过在这里不影响  这里你即便不写 默认的x(1),y(1)都是0  改用其他数据的的话 要记得定义数组的首项
format的用法我就不细说了 这里用不到  你回头自己看书吧
dat文件可以用记事本打开 就像txt一样 复制 然后粘贴到oringin就行了
赞一下 回复此楼
网络错误 无法获取数据
15楼2013-08-25 20:24:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宇天旅梦

铜虫 (小有名气)
引用回帖:
15楼: Originally posted by phy_yo at 2013-08-25 20:24:55
y(2)=0 是我打错了 应该是y(1)=0 不过在这里不影响  这里你即便不写 默认的x(1),y(1)都是0  改用其他数据的的话 要记得定义数组的首项
format的用法我就不细说了 这里用不到  你回头自己看书吧
dat文件可以用记事 ...

这一部分明白了，可是自己写了一个和这大同小异的程序，出错了。

部分程序如下：
do 2 k=1,100
        m=(KY*50+i*SQRT(KX*KX+KY*KY))
     $/(SQRT(KX*KX+KY*KY+KY*KY*50*50)-50*50/(6.582122*10**28))
      
        n=(KY*50-i*SQRT(KX*KX+KY*KY))/
     $(SQRT(KX*KX+KY*KY+KY*KY*50*50)-50*50/(6.582122*10**28))
     
      w=(KY*50+i*SQRT(1+KY*KY))/
     $(SQRT(KX*KX+KY*KY+KY*KY*50*50))

      z=(KY*50-i*SQRT(1+KY*KY))/
     $(SQRT(KX*KX+KY*KY+KY*KY*50*50))
        f=((m-n)*(w-z)*EXP(-i*SQRT(1+KY*KY))/
     $((m-w)*(n-z)*e**(2*i*SQRT(KX*KX+KY*KY)-(m-z)(n-z))
其中关系是 m(KX,KY)    n(KX,KY)     w(KX,KY)  ,  z(KX,KY)

f(m(KX,KY), n(KX,KY), w(KX,KY), z(KX,KY))

给出的错误是：
--------------------Configuration: Text1 - Win32 Debug--------------------
Compiling Fortran...
C:\Program Files\DevStudio\MyProjects\Text1.for
C:\Program Files\DevStudio\MyProjects\Text1.for(28) : Error: Syntax error, found '(' when expecting one of: * :: , <END-OF-STATEMENT> ; : ) (/ [ + - . ] /) ** / > // .LT. < ...
     $(((m-w)*(n-z)*exp(2*i*SQRT(KX*KX+KY*KY))-(m-z)(n-z))
----------------------------------------------------^
C:\Program Files\DevStudio\MyProjects\Text1.for(28) : Error: Syntax error, found END-OF-STATEMENT when expecting one of: , )
     $(((m-w)*(n-z)*exp(2*i*SQRT(KX*KX+KY*KY))-(m-z)(n-z))
----------------------------------------------------------^
C:\Program Files\DevStudio\MyProjects\Text1.for(15) : Error: The shapes of the array expressions do not conform.   [M]
m=(KY*50+i*SQRT(KX*KX+KY*KY))/
--------^
C:\Program Files\DevStudio\MyProjects\Text1.for(18) : Error: The shapes of the array expressions do not conform.   [N]
n=(KY*50-i*SQRT(KX*KX+KY*KY))/
--------^
C:\Program Files\DevStudio\MyProjects\Text1.for(21) : Error: The shapes of the array expressions do not conform.   [W]
      w=(KY*50+i*SQRT(1+KY*KY))/
------^
C:\Program Files\DevStudio\MyProjects\Text1.for(24) : Error: The shapes of the array expressions do not conform.   [Z]
      z=(KY*50-i*SQRT(1+KY*KY))/
------^
Error executing df.exe.

Text1.exe - 6 error(s), 0 warning(s)
赞一下 回复此楼
16楼2013-08-25 21:45:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phy_yo

铁虫 (小有名气)
引用回帖:
16楼: Originally posted by 宇天旅梦 at 2013-08-25 21:45:24
这一部分明白了，可是自己写了一个和这大同小异的程序，出错了。

部分程序如下：
do 2 k=1,100
        m=(KY*50+i*SQRT(KX*KX+KY*KY))
     $/(SQRT(KX*KX+KY*KY+KY*KY*50*50)-50*50/(6.582122*10**28))
      
...

你这部分有点部分过头了吧  循环变量k在循环体内根本没看到  你应该不是故意删的……所以你的这些表达式不会是直接复制粘贴的吧……还全是大写
KX,KY是是自变量 要做成数组形式 如KX（20），KY（20）之类的  再套入循环变量就行了  你自己慢慢琢磨吧  我感觉我前面写的那些代码够你做参考了~或者你可以考虑用matlab   作这种计算他比fortran方便
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

宇天旅梦
网络错误 无法获取数据
17楼2013-08-25 22:20:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宇天旅梦

铜虫 (小有名气)
送红花一朵
引用回帖:
17楼: Originally posted by phy_yo at 2013-08-25 22:20:40
你这部分有点部分过头了吧  循环变量k在循环体内根本没看到  你应该不是故意删的……所以你的这些表达式不会是直接复制粘贴的吧……还全是大写
KX,KY是是自变量 要做成数组形式 如KX（20），KY（20）之类的 ...

不是粘的啊，我一个一个敲进去的。。。。。敲了几次，大小写混了。。。

我琢磨琢磨吧。

非常非常感谢你的帮助啦，送一朵红花！
赞一下 回复此楼
18楼2013-08-25 22:26:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 宇天旅梦 的主题更新
182/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 找调剂 +10 楚乔乔 2026-04-01 10/500 2026-04-05 22:19 by syh9288
[考研] 一志愿同济大学323分（080500）求调剂 +8 yikeniu 2026-04-01 8/400 2026-04-05 18:15 by cql1109
[考研] 285求调剂 +11 哦呦呼o 2026-04-04 11/550 2026-04-05 17:59 by 猪会飞
[考研] 367求调剂 +3 谢28 2026-03-30 3/150 2026-04-05 13:27 by huangmoli
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 085410人工智能 初试316分 求调剂 +6 残星拂曙 2026-03-31 6/300 2026-04-05 12:15 by rainbow11
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 一志愿南农090401，268，求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
[考研] 求调剂 +3 农业工程与信息� 2026-04-04 3/150 2026-04-04 12:19 by 舍而后得
[考研] 材料调剂 +11 吴棂颖！ 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 求材料调剂，一志愿郑州大学289分 +15 硕星赴 2026-04-03 15/750 2026-04-04 01:01 by userper
[考研] 305求调剂 +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
[考研] 求调剂不挑专业 +3 xrh030412 2026-04-01 3/150 2026-04-03 14:40 by 氮气气气
[考研] 260求调剂 +6 朱芷琳 2026-04-02 6/300 2026-04-02 20:27 by 6781022
[考博] 材料工程专业硕士申博 +3 麟正宇 2026-03-30 3/150 2026-04-02 15:04 by greychen00
[考研] 0856初试324分求调剂 +6 想上学求调 2026-04-01 6/300 2026-04-02 11:42 by 星空星月
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 358求调剂 +3 王向阳花 2026-03-31 3/150 2026-04-01 09:56 by zzchen2000
[考研] 本2一志愿C9-333分，材料科学与工程，求调剂 +9 升升不降 2026-03-31 9/450 2026-03-31 18:01 by 无际的草原
信息提示
关闭
请填处理意见
关闭