原创力文档- 1、本文档共7页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
计算物理第二次课程作业张楚珩121120173姓名-张楚珩的仓库
计算物理第⼆次课程作业 张楚珩 121120173
姓名：张楚珩
⽼师：李⽂⻜
计算物理
2014年4 ⽉17⽇
Lorenz Model
⼀、算法流程与计算公式：
利⽤了四阶Runge-Kutta⽅法，根据现已推演到的值x ,y ,z 和已知的微分⽅程，利
n n n
⽤公式
来进⾏计算。
其计算的流程图如下：
⻚码： 1/ 7
计算物理第⼆次课程作业 张楚珩 121120173
⼆、程序源代码（C语⾔）：
//
// main.c
// ComputationalPhysics
//
// Created by ZHANG CH on 14-4-16.
// Copyright (c) 2014年 NJU. All rights reserved.
//
//PROBLEM STATING:
//dx/dt = 10(y-x)
//dy/dt = -xz + rx - y
//dz/dt = xy - (8/3)z
// PRO 1:
// r = 25 Plot z-x plane
// PRO 2:
// change the val of r, plot z-x plane
//
// 1. Draw flow diagram
// 2. Make program
// 3. Give out computation results and discuss the results
//Formula implemented:
//{x_n+1} = {x_n} + (DELTA_T / 6) * (k1 + 2 k2 + 2 k3 + k4)
//k1 = f({x_n})
//k2 = f({x_n + k1*(DELTA_T/2)})
//k2 = f({x_n + k2*(DELTA_T/2)})
//k4 = f({x_n + k3*(DELTA_T)})
//input data: 25 25 0 100 \\ 12 2 9
#include stdio.h
#define DELTA_T (0.001)
#define NUM_OF_POINTS (30010)
#define SIGMA (10.0)
#define B (8.0/3.0)
double x[NUM_OF_POINTS] = {0};
double y[NUM_OF_POINTS] = {0};
double z[NUM_OF_POINTS] = {0};
int numOfLoop;
double r_min, r_max, r_step, s;
double x_ini,y_ini,z_ini;
FILE * fp;
void in
文档评论（0）