| 查看: 501 | 回复: 0 | ||
[求助]
Fortran转化 C++
|
|
1、试将下列描述大气对流的洛伦兹微分方程组用C 语言实现出来，要求编写出完整的程序代码。(时间步长取0.001) 其中参数值分别为： ， ===================================================================== 2、罗斯勒设计了一个新的吸引子方程组，如下： 其中参数值分别为： ， 时间步长取0.001。试将该微分方程组用C 语言实现出来，要求编写出完整的程序代码。 ============================================================================== 3.将杜芬方程组用C 语言实现出来，要求编写出完整的程序代码。时间步长取0.001。 下面是用Fortran 语言编写的程序代码，拜托转化为C++语言代码。 program lorenz implicit none integer N real x, y, z, dx, dy, dz, dt Real r, sigma, b integer i,j sigma=10.0;b=8/3; N=1E6;dt=0.001; r=28.0 open (unit=1,file='1.dat',status='unknown') x=0.10; y=0.30; z=0.2 do i=1,100000 dx=dt*( sigma*(y-x) ) dy=dt*( r*x-y-x*z ) dz=dt*( x*y - b*z ) x=x+dx y=y+dy z=z+dz if ( MOD( i, 20 )==0) then write (1,'(1x, 3f12.5)') x, y, z end if end do close (1) end program lorenz ================================================================!!! program rossler implicit none integer N integer, parameter:: NN=7 real x, y, z, dx, dy, dz, dt Real a, b, c, x_c(NN),c1,c2 integer i,j,k ! The values of parameters about the sy(i)s character aa*1, bb*1, cc*2 !------------------------------------------------------------------------------------! a=0.2; b=0.2; N=1E6; dt=0.01; !c=2.6 !,3.5, 4.1, 4.18, 4.21, 4.6, 5.7 x_c(1:7) =(/2.6, 3.5, 4.1, 4.18, 4.21, 4.6, 5.7/) !------------------------------------------------------------------------------------! x=0.10; y=0.30; z=0.2 do k=1,NN c=x_c(k) c1=(c*100)/100 c2=c*100.0-100*int(c) write(aa, '(i1.1)') int(k) write(bb, '(i1.1)') int(c1) write(cc, '(i2.2)') int(c2*1.01) open (unit=1,file='Data\rossler_c'//aa//'='//bb//'.'//cc//'.dat',status='unknown') do i=1,40000 dx=dt*( -y-z ) dy=dt*( x+a*y ) dz=dt*( b+(x-c)*z ) x=x+dx y=y+dy z=z+dz if( MOD( i, 40 )==0 ) then write (1,'(1x, 3f12.5)') x, y, z end if end do close (1) end do end program rossler integer i, j, ii character aa*1, bb*1, cc*2 ! The values of parameters about the sy(i)s pi=3.14; N=1E6 x_k(1:8)=(/0.5, 0.7, 0.9, 1.0, 1.2, 1.3, 1.4, 1.5/) do ii=1, NN k=x_k(ii) k1=(k*10.0)/10 k2=k*10.0-10*int(k) write(aa, '(i1.1)') ii write(bb, '(i1.1)') int(k1) write(cc, '(i1.1)') int(k2*1.01) open (unit=1,file='Data\bzys_k'//aa//'='//bb//'.'//cc//'.dat',status='unknown') do j=1, 20 x_0(j)=2*pi*j/20 y_0(j)=1.0*pi do i=1,100000 x(1)=x_0(j) y(1)=y_0(j) if( x(i)>=2*pi ) then x(i)=x(i)-2*pi else if( x(i)<=0.0 ) then x(i)=x(i)+2*pi end if if( y(i)>=2*pi ) then y(i)=y(i)-2*pi else if( y(i)<=0.0 ) then y(i)=y(i)+2*pi end if dx=dt* y dy=dt*( k*x-x**3 ) x=x+dx y=y+dy if ( MOD( i, 40 )==0) then write (1,'(1x, 2f12.5)') x, y end if end do end do close (1) end do end program duff |
回复此楼» 猜你喜欢
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
-
重要消息，中午系统在维护
已经有12人回复
找到一些相关的精华帖子，希望有用哦~
fortran文件读入与生成输出
已经有3人回复
- 请问下面这段c++小程序为什么不能运行
已经有6人回复
- 回帖就送金币！！请大家对CFD的编程实现给点意见！
已经有8人回复
- 相场法fortran源代码
已经有167人回复
- 关于C++中结果输出的问题
已经有8人回复
- Fortran输出二进制数据，用matlab读出，结果出错，求助问题何在？
已经有4人回复
- 编一个fortran程序
已经有17人回复
- 求助：并行计算的首选语言是C,C++还是Fortran？
已经有11人回复
- 关于fortran子程序调用内部函数
已经有7人回复
- fortran里一个子程序不运行是怎么回事
已经有13人回复
- C/C++开发专家：Visual C++代码参考与技巧大全
已经有89人回复
- C++程序怎么实现时间的转换呢？
已经有9人回复
- 如何将多模光耦合到单模光纤里，或者将其转换为单模光？
已经有28人回复
- 物理模拟是用c语言好还是fortran好？
已经有30人回复
- fortran转C或C++有什么好办法呢？要转换得比较彻底的且容易修改的。
已经有4人回复
- 如何改变fortran的浮点精度
已经有28人回复
- C语言与Fortran语言之间的转化
已经有9人回复
- 【讨论】想学lammps是否应该先学习C++
已经有28人回复
- 【求助】C++程序所需内存远远大于物理内存怎么办？
已经有8人回复
- 【分享】C/C++常用算法设计方法,fortran其他语言编程也可以参考下
已经有8人回复
- 【求助】请问一下 c++有没有按位读取文件的方法或者语句？ STL里有没有这样的功能？
已经有6人回复
- 【求助】fortran计算时出现NaN
已经有11人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
40