- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数值分析实验
实验指导书四
理学院实验中心
数学专业实验室编写
实验四 常微分方程初值问题数值解法
【】【】【】Euler法,改进Euler法,Rung-Kutta方法求其数值解,诸如以下问题:
分别取h=0.1,0.2,0.4时数值解。
初值问题的精确解。
【】1、 熟悉各种初值问题的算法;
2、 明确各种算法的精度与所选步长有密切关系;
3、 通过计算更加了解各种算法的优越性。
【】1、 根据初值问题数值算法,编程计算;
2、 试分别取不同步长,考察某节点处数值解的误差变化情况;
3、 试用不同算法求解某初值问题,结果有何异常;
4、 分析各个算法的优缺点。
【Euler方法】
#includestdio.h
#includemath.h
/*y1=f(x,y)*/
float f(float x,float y)
{return(4*x / y-x*y);}
void Euler(float x0,float xn,float y0,int n)
{int i;
float x=x0,y=y0,h=(xn-x0)/n;
printf(x[0]=%f\ty[0]=%f\n,x,y);
for(i=1;i=n;i++)
{y=y+h*f(x,y);
x=x0+i*h;
printf(x[%d]=%f\ty[%d]=%f\n,i,x,i,y);
}}
main()
{int i,n;
float x0,xn,y0;
printf(nInput the begin and end of x:);
scanf(%f%f,x0,xn);
printf(Input the y value at %f:,x0);
scanf(%f,y0);
do {printf(Input n value[divide(%f,%f)]:,x0,xn);
scanf(%d,n);
}
while(n=1);
Euler(x0,xn,y0,n);
scanf(%f,yn);}
【Rung-Kutta方法】
#includestdio.h
#includemath.h
#define f(x,y)(y-2*x/y)
int main()
{
int m;
int i;
double a,b,y0;
double xn,yn,yn1;
double k1,k2,k3,k4;
double h;
printf(\nInput the begin and end of x:);
scanf(%lf%lf,a,b);
printf(nInput the y value at %f:,a);
scanf(%lf,y0);
printf(nInput m value[divide(%f,%f)]:,a,b);
scanf(%d,m);
if(m=0)
{
printf(Please input a number larger than 1.\n);
return 1;
}
h=(b-a)/m;
xn=a;yn=y0;
for(i=1;i=m;i++)
{
k1=f(xn,yn);
k2=f((xn+h/2),(yn+h*k1/2));
k3=f((xn+h/2),(yn+h*k2/2));
k4=f((xn+h),(yn+h*k3));
yn1=yn+h/6*(k1+2*k2+2*k3+k4);
xn+=h;
printf(x%d=%f,y%d=%f\n,i,xn,i,yn1);
yn=yn1;
}
scanf(%lf,yn);}
您可能关注的文档
- 安阳市加快发展现代职业教育任务分工及进度安排表安阳市加快发展现代职业教育任务分工及进度安排表.doc
- 安阳市2011年中考政治模拟试题及答案安阳市2011年中考政治模拟试题及答案.doc
- 安阳国际大酒店高支模施工方案安阳国际大酒店高支模施工方案.doc
- 宋秀岩:实施生态立省战略努力建设生态文明宋秀岩:实施生态立省战略努力建设生态文明.doc
- 安防技术_GE2000数字矩阵使用说明书安防技术_GE2000数字矩阵使用说明书.doc
- 宏志七年级_上册人教新目标英语模拟测试题宏志七年级_上册人教新目标英语模拟测试题.doc
- 安徽省合肥市2013届高三第一次教学质量检测数学理试题 Word版含答案安徽省合肥市2013届高三第一次教学质量检测数学理试题 Word版含答案.doc
- 完善手段不断创新努力提高工程质量监管水平完善手段不断创新努力提高工程质量监管水平.doc
- 宏观经济政策目标及其联系宏观经济政策目标及其联系.doc
- 宏观经济政策对保险业的影响20110420宏观经济政策对保险业的影响20110420.doc
文档评论(0)