Fortran | 普通差分法解非线性常微分方程组（洛伦兹对流方程组）

0.写在前面

今天上机，老师直接让解方程组，我们上节课只教了龙格-库塔解法，我直接一脸懵逼......好在想到了，就是用最一般的差分解。

1.待求方程组

2.初始条件和要求

σ=10，ρ=28，β=8.0/3

t0=0，x0=0，y0=1，z0=0

Δt=0.01 N=3000

3.代码

program exp1implicit nonereal x0,x(1:3000),y0,y(1:3000),z0,z(1:3000),t0,t(1:3000),ts,a,b,cinteger N,its=0.01n=3000x0=0.0y0=1.0z0=0.0t0=0.001a=10.0b=28.0c=8.0/3x(1)=x0+ts*a*(y0-x0)y(1)=y0+ts*x0*(b-z0)-y0*tsz(1)=z0+ts*(x0*y0-c*z0)t(1)=t0+tsdo i=2,nx(i)=x(i-1)+ts*a*(y(i-1)-x(i-1))y(i)=y(i-1)+ts*x(i-1)*(b-z(i-1))-y(i-1)*ts z(i)=z(i-1)+ts*(x(i-1)*y(i-1)-c*z(i-1))t(i)=t(i-1)+tsenddoopen(1,file='./output_exp1.csv')write(1,'(3(f13.8,a),f13.8)') t0,',',x0,',',y0,',',z0do i=1,nwrite(1,'(3(f13.8,a),f13.8)') t(i),',',x(i),',',y(i),',',z(i)enddoclose(1)
end program exp1

#因为可视化的代码是老师给的，不好放出来，大家看个效果就行

#exp2的初值z0改为0.001，其他条件不变

Fortran | 普通差分法解非线性常微分方程组（洛伦兹对流方程组）相关推荐

用matlab画多普勒加宽线性函数,洛伦兹线性函数
对于 CEST MRI 图像中的任一个像素,利用伪佛克脱线型(PVP)代替洛伦兹作为拟合函数, 使拟合方法适合在更大的饱和功率和组织中磁化转移含量更大的情况下拟合,...... 3.3谱线加宽和线 ...
c语言模拟洛伦兹吸引子,洛伦兹吸引子相图的简易实现.pdf
洛伦兹吸引子相图的简易实现.pdf 洛伦兹吸引子相图的简易实现晏志武 ( 临沧市临翔区第一中学云南临沧 ) ( 收稿日期 ) 摘要用平面动画软件绘制出高质量的洛伦兹吸引子相图关键词洛伦兹 ...
matlab求洛伦兹方程的解,[转载]用Matlab求解洛伦兹方程
1. 洛伦兹方程求解本文说明用Matlab工具箱求解洛伦兹方程的过程,并给出吸引子的三维动态图象.洛伦兹方程如下: (1)这是一个自洽的方程组,求解过程如下: (1) 建立自定义函数 functio ...
【Scipy高级计算】(2) 常微分方程、洛伦兹吸引子，附python完整代码
大家好,在上一篇博文中,我介绍了如何使用Scipy库计算定积分和二重积分,感兴趣的可以看一下:https://blog.csdn.net/dgvv4/article/details/124226759 ...
高斯分布+柯西-洛伦兹分布+三种光谱线型函数（洛伦兹线型函数+多普勒[高斯]线型函数+vogit 线型函数）
更新时间:2020/4/26 文章目录一.两种分布 1.1 高斯分布(正态分布) 1.2 柯西分布(洛伦兹分布) 二.光谱线型函数 2.1 光谱线型函数的几种形式 2.2 洛伦兹线型函数与多普勒(高 ...
泊松回归、gamma回归、Tweedie回归等广义线性回归模型GLM的评估指标：校准曲线、洛伦兹曲线、卡方检验、AIC、BIC、偏差（Deviance）指标
泊松回归.gamma回归.Tweedie回归等广义线性回归模型GLM的评估指标:校准曲线(Calibration curve). 洛伦兹曲线(Lorenz Curve).卡方检验.AIC.BIC.偏差 ...
MATLAB 数据处理（二）非线性拟合——洛伦兹拟合（Lorentz fit）
文章目录一.拟合示例二.单峰洛伦兹 2.1 洛伦兹线型函数表达式与物理含义 2.2 lsqcurvefit非线性拟合 2.3 代码实现三.双峰洛伦兹 3.1 洛伦兹线型函数表达式与物理含义 3. ...
相对论【3】洛伦兹速度变换和婓索实验
前言洛伦兹变换有一个洛伦兹因子γ\gammaγ,通过推导得到相对论条件下,有这样一个压缩因子,之后推导一是可以从头推,二是可以在伽利略变换的方程基础上加上这个因子 . 洛伦兹速度变换之前进行洛伦兹 ...
R语言建模收入不平等：分布函数拟合及洛伦兹曲线(Lorenz curve)
最近我们被客户要求撰写关于洛伦兹曲线的研究报告,包括一些图形和统计输出. 洛伦兹曲线来源于经济学,用于描述社会收入不均衡的现象.将收入降序排列,分别计算收入和人口的累积比例. 本文,我们研究收入和不平 ...

Fortran | 普通差分法解非线性常微分方程组（洛伦兹对流方程组）

0.写在前面

1.待求方程组

2.初始条件和要求

3.代码

Fortran | 普通差分法解非线性常微分方程组（洛伦兹对流方程组）相关推荐

最新文章

热门文章