数值分析龙格库塔

前言随着计算机的迅速开展和广泛应用，在众多领域内，我们越来越认识到科学计算是科学研究的重要方法。数值计算方法是一种利用计算机解决数学问题的数值近似解方法，特别是无法用人工过计算器计算的数学问题。数值计算方法常用于矩阵高次代数方程矩阵特征值与特征向量的数值解法，插值法，线性方程组迭代法，函数逼近，数值积分与微分，常微分方程初值问题数值解等。作为数学与计算机之间的一条通道，数值计算的应用范围已十分广泛，作为用计算机解决实际问题的纽带，数值算法在求解线性方程组，曲线拟合、数值积分、数值微分，迭代方法、插值法、拟合法、最小二乘法等应用广泛。数值计算方法是和计算机紧密相连的，现代计算机的出现为大规模的数值计算创造了条件，集中而系统的研究适用于计算机的数值方法是十分必要的。数值计算方法是在数值计算实践和理论分析的根底上开展起来的。通过数值计算方法与实验将有助于我们理解和掌握数值计算方法根本理论和相关软件的掌握，熟练求解一些数学模和运算。并提高我们的编程能力来解决实际问题。摘要对于本次计算方法与实习的实践环节，我们采用改良欧拉〔Euler〕方法对给定的数据进行分析，改良的欧拉〔Euler〕方法是解决常微分方程初值问题常用的数值解法，本文在简要介绍改良欧拉〔Euler〕方法及四阶龙格库塔公式的根底上，通过编写C语言程序实现两种数值算法。通过本次实践环节，我们很好的了解了常微分方程数值解法的原理。出色的完成了本次课程设计。关键词：欧拉方法；四阶龙格--库塔；C语言；数值分析目录前言 -0-摘要 -1-实验设计内容 -3-一．实验目的 -3-二．实验内容 -3-三．实验算法 -3-四．实验程序 -4-⑴改良欧拉方法 -4-⑵四阶龙格库塔方法 -5-实验心得 -7-实验设计内容一．实验目的熟悉求解常微分方程初值问题的有关方法和理论，主要是改良欧拉法、四阶龙格-库塔法；会编制上述方法的计算程序，包括求解微分方程组的计算程序。二．实验内容分别用改良的欧拉方法与四阶龙格库塔公式取〔h=0.1〕求解些列微分方程初值问题：三．实验算法解一阶常微分方程初值问题将区间[a,b]作n等分，取步长。欧拉公式为梯形公式为改良欧拉法，采用公式或表常用的是四阶龙格-库塔公式取步长为h，由初值出发，可得未知函数在区间[a,b]上各节点处的近似值。四．实验程序⑴改良欧拉方法#include<stdio.h>#defineN20voidModEuler(float(*f)(float,float),floatx0,floaty0,floatxn,intn){ inti; floatyp,yc,x=x0,y=y0,h=(xn-x0)/n; printf("x[0]=%f\ty[0]=%f\n",x,y); for(i=1;i<=n;i++) { yp=y+h*(*f)(x,y); x=x0+i*h; yc=y+h*(*f)(x,yp); y=(yp+yc)/2; printf("x[%d]=%f\ty[%d]=%f\n",i,x,i,y); }}floatf(floatx,floaty){ returnx*x+y*y;}voidmain(){ floatxn=1,x0=0,y0=0; ModEuler(f,x0,y0,xn,N); getch();}2.运行结果⑵四阶龙格库塔方法1．龙格库塔程序#include<stdio.h>voidRunge_Kutta(float(*f)(floatx,floaty),floata,floatb,floaty0,intN){ floatx=a,y=y0,K1,K2,K3,K4; floath=(b-a)/N; inti; printf("x[0]=%f\ty[0]=%f\n",x,y); for(i=1;i<=N;i++) { K1=(*f)(x,y); K2=(*f)(x+h/2,y+h*K1/2); K3=(*f)(x+h/2,y+h*K2/2); K4=(*f)(x+h,y+h*K3); y=y+h*(K1+2*K2+2*K3+K4)/6; x=a+i*h; printf("x[%d]=%f\ty[%d]=%f\n",i,x,i,y); }}floatf(floatx,floaty){ returnx*x+y*y;}voidmain(){ floata=0,b=1,y0=0; Runge_Kutta(f,a,b,y0,20); getch();}2.龙格库塔结果实验心得课程设计是培养学生综合运用所学知识,发现,提出,分析和解决实际问题,锻炼实践能力的重要环节,是对学生实际工作能力的具体训练和考察过程.随着科学技术开展的日新日异，程序软件已经成为当今计算机应用中空前活泼的领域，在生活领域中的应用可以说得是无处不在。因此作为二十一世纪的大学来说掌握简单的软件的开发技术是十分重要的。回忆起此次数值分析的课程设计，至今我仍感慨颇多，确实，从理论到实践，在整整一个星期的日子里，可以说得是苦多于甜，但是可以学到很多很多的的东西，同时不仅可以稳固了以前所学过的知识，而且学到了很多在书本上所没有学到过的知识。通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正为社会效劳，从而提高自己的实际动手能力和独立思考的能力。在设计的过程中遇到问题，可以说得是困难重重，这毕竟第一次做的，难免会遇到过各种各样的问题，同时在设计的过程中发现了自己的缺乏之处，对以前所学过的知识理解得不够深刻，掌握得不够牢固，通过这次课程设计之后，一定把以前所学过的知识重新温故。我得到了很多同学的帮助。我想没有他们我可能都要放弃了，因为我本人对数值分析的算法本来就不是很熟悉，学的东西也不能很好的灵活应用，理论联系不了实际。以前的汇编语言没学好，一开始的程序这块儿就要令我抓狂了。通过这段时间的学习，我认为要学好C语言编程这门课程，不仅要认真阅读课本知识，更重要的是要通过上机实践才能增强和稳固我的知识。特别是作为大学生，更要注重实践这一环节，只有这样我们才能成为合格的计算机人材。整个过程不断的调试，在调用子程序的时候遇到了局部问题，最后幸好都解决了。通过这次课程设计也稳固了一些已经学习的数值分析的算法，同时了解了