原创力文档- 1、本文档共48页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
一、经典四阶龙格-库达法解一阶微分方程组
1．摘要：已知点近似精度求解微分方程组
2. 关键词：龙格-库达法(R-K法)
3. 正文
3.1 算法原理：
从点（t0，y0）开始迭代， yk+1=yk+h(f1+2f2+2f3+f4)/6，
其中：
f1=f(tk,yk)
f2=f(tk+h/2,yk+hf1/2)
f3=f(tk+h/2,yk+hf2/2)
f4=f(tk+h,yk+hf3)
3．2 程序流程图：
3.3 程序代码：
#include iostream
#includeiomanip
using namespace std;
float G(float pointT,float pointX,float pointY);
float F(float pointT,float pointX,float pointY);
void function(float * t,float * x,float * y,float h);
struct aa
{
float t;
float x;
float y;
struct aa*next;
};
void statement(void)
{
cout******************endl;
cout*计算微分方程组：*endl;
cout*dx/dt=x+2y *endl;
cout*dy/dt=3x+2y *endl;
cout*的龙格-库塔解 *endl;
cout******************endl;
}
float G(float pointT,float pointX,float pointY)
{
return(pointX+2*pointY);
}
float F(float pointT,float pointX,float pointY)
{
return(3*pointX+2*pointY);
}
void function(float * t,float * x,float * y,float h)
{
float g[4],f[4];
g[0]=G(*t,*x,*y);
f[0]=F(*t,*x,*y);
g[1]=G(*t+h/2,*x+h*g[0]/2,*y+h*f[0]/2);
f[1]=F(*t+h/2,*x+h*g[0]/2,*y+h*f[0]/2);
g[2]=G(*t+h/2,*x+h*g[1]/2,*y+h*f[1]/2);
f[2]=F(*t+h/2,*x+h*g[1]/2,*y+h*f[1]/2);
g[3]=G(*t+h,*x+h*g[2],*y+h*f[2]);
f[3]=F(*t+h,*x+h*g[2],*y+h*f[2]);
*t=*t+h;
*x=*x+h*(g[0]+2*g[1]+2*g[2]+g[3])/6;
*y=*y+h*(f[0]+2*f[1]+2*f[2]+f[3])/6;
}
void main()
{
int step=0,count=0;
float a=0,b=0,h=0;
struct aa * head=NULL,* point=NULL,*nextpoint=NULL;
float * pointT=NULL,* pointX=NULL,* pointY=NULL;
char flag=y;
while (flag==y||flag==Y)
{
point=(struct aa *)malloc(sizeof(struct aa));
point-next=NULL;
head=point;
nextpoint=head;
statement();
cout请输入微分区间[a,b]端点:;
cinab;
point-t=a;
cout请输入函数在端点a=%d处的初值(x,y):;
cinpoint-xpoint-y;
cout请输入微分区间的区间段个数step:;
cinstep;
h=(b-a)/step;
pointT=(float *)malloc(sizeof(float));
pointX=(float *)malloc(sizeof(float));
pointY=(float *)malloc(sizeof(float));
* pointT=point-t;
* pointX=point-x;
* pointY=point-y;
for(count=0;countstep;count++)
{
point=(struct aa *)malloc(sizeof(struct
文档评论（0）