微分方程计算方法：四阶龙格-库塔法求解Lorenz方程及敏感性分析.pdfVIP

1、用四阶龙格-库塔法解著名的Lorenz方程，

x′=−sx+sy





y′=−xz+rx−y

′

z=xy−bz



其中b=8/3，普朗特数s=10，雷诺数r=28，初始值取[10,10,10]，

给出x-z平面上的轨线分布图；

分别调整初始值及参数r、s、b，讨论Lorenz方程的解对初始变量和参数的敏

感性。

一、流程图

1、用四駖龙格‑库塔法解著名的洛伦兹方程，

x′=−sx+sy





y′=−xz+rx−y

′

z=xy−bz



份位b=8/3,普朗特数s=10,雷诺数r=28,值取[10,10,10],

x‑z处理院线图像图;

和参数r、s、b，讨论洛伦兹方程对初变量的解和参数的重力。

一、流程图

二、源代码

#includestdio.h

#includestdlib.h

#includemath.h

#defineb8./3

#defines10.

#definer28.

#defineh0.01

doublef(doublexx,doubleyy,doublezz,inti)

{

doublex[3]={xx,yy,zz};

if(i0)return-s*x[0]+s*x[1];

if(i1)return-x[0]*x[2]+r*x[0]-x[1];

if(i2)returnx[0]*x[1]-b*x[2];

}

intmain()

{

doublex[3]={10.,10.,10.};

doublek[5][3],t[5]={0,0,h/2.,h/2.,h};

inti,j;

intm;

FILE*out;

out=fopen(rl.txt,w);

for(m=1;m=1000;m++)

{

for(i=1;i5;i++)//calk

{

for(j=0;j