- 1、本文档共9页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实用标准文案
文档
常微分方程的求解与定性分析实验报告
一、实验综述
1、实验目的及要求
归纳和学习求解常微分方程(组)的基本原理和方法;
掌握解析、数值解法,并学会用图形观察解的形态和进行解的定性分析;
熟悉MATLAB软件关于微分方程求解的各种命令;
通过范例学习建立微分方程方面的数学模型以及求解全过程;
通过该实验的学习,使学生掌握微分方程(组)求解方法(解析法、欧拉法、梯度法、改进欧拉法等),对常微分方程的数值解法有一个初步了解,同时学会使用MATLAB软件求解微分方程的基本命令,学会建立微分方程方面的数学模型。这对于学生深入理解微分、积分的数学概念,掌握数学的分析思维方法,熟悉处理大量的工程计算问题的方法是十分必要的。
2、实验仪器、设备或软件
电脑 、matlab7.0
二、实验过程(实验步骤、记录、数据、分析)
实验内容:
根据实验内容和步骤,完成以下实验,要求写出实验报告(实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论)
1.求微分方程的解析解,并画出它们的图形。
y = y + 2 x, y (0) = 1, 0 x 1;
m=dsolve(Dy=y+2*x,y(0)=1,x)
ezplot(m,[0 1])
m =
3*exp(x) - 2*x – 2
1.求微分方程 的数值解,要求编写求解程序。
function dy=vdp1000(t,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=-y(1)+0.1*y(1)^3;
[T,Y]=ode15s(vdp1000,[0 10],[0 0]);
plot(T,Y(:,1),-)
3.Rossler微分方程组:
当固定参数b=2,c=4时,试讨论随参数a由小到大变化(如 a∈(0,0.65))而方程解的变化情况,并且画出空间曲线图形,观察空间曲线是否形成混沌状?
function r=rossler(t,x)
global a;
global b;
global c;
r=[-x(2)-x(3);x(1)+a*x(2);b+x(3)*(x(1)-c)];
global a;
global b;
global c;
b=2;
c=4;
t0=[0,200];
for a=0:0.1:0.6
[t,x]=ode45(rossler,t0,[0,0,0]);
subplot(1,2,1);
plot(t,x(:,1),r,t,x(:,2),g,t,x(:,3),b);
title(x(红色),y(绿色),z(蓝色)随t的变化情况);xlabel(t);
subplot(1,2,2);
plot3(x(:,1),x(:,2),x(:,3))
title(相图);xlabel(x);ylabel(y);zlabel(z);
pause
end
结果显示:
a=0:
a=0.1:
a=0.2:
a=0.3:
a=0.4:
a=0.5:
结果分析:从图像可以看出,当a=0时,微分方程的解(x,y,z)收敛与(0,0.5,0.5);当a=0.1时,(x,y,z)仍收敛与(0,0.5,0.5),只是收敛速度减慢;当a=0.2时,(x,y,z)已发散,周期性变化;随着a的增大,(x,y,z)接近其极限环的速度加快,空间曲线成混沌状。
4.炮弹发射角的确定
炮弹发射视为斜抛运动,已知初始速度为200m/s,若要击中水平距离360m、在垂直距离160m的目标,当忽略空气阻力时,发射角应为多大?此时炮弹的运行轨迹如何?
要求:
(1) 建立在忽略空气阻力情况下的描述炮弹发射轨迹的数学模型;
(2) 用Matlab 软件求解方程和微分方程;
(3) 结合实际对解的合理性进行分析。
进一步思考:
如果要考虑水平方向的阻力,且设阻力与(水平方向)速度成正比,系数为0.1(1/s ),结果又如何?此时炮弹的运行轨迹如何?
解:
(1)忽略空气阻力时,设发射角为a,炮弹的飞行时间为t
水平方向:Vx=V0*cos a=200*cos(a)
竖直方向:Vy=V0*sin a=200*sin(a)
得到:x=Vx*t=200*cos(a)*t=360
y=Vy*t-1/2*g*t^2=200*sin(a)*t-1/2*9.8*t^2=160
得到炮弹的路程为Y=360*tan(a)-4.9*(360/200/cos(a)).^2-160
编程为:
function Y=fun2(a)
Y=360*tan(a)-4.9*(360/200/cos(a)).^2-160;
function Y=fun3(a0,a1,n,tol)
a(1)=a0;
a(2)=a1;
b=1;
i=2;
while(a
您可能关注的文档
- QCT413-2002汽车电气设备基本技术条件.doc
- 岩石力学研究地现状和未来.doc
- 三角形三边关系教学设计课题.doc
- 2016移动吊装作业案列分析报告题.doc
- 分层分类地战略性人力资源管理系统模式.doc
- 2、工法正文--覆土顶板沉降后浇带先封闭后浇筑施工工法(1).doc
- 2012年数学一轮复习精品试题第十、十一模块概率与统计算法初步、数系地扩充与.doc
- 八年级数学下册一次函数专题练习---2.doc
- 2015年普通高等学校招生全国统一考试语文试题(福建卷,含解析汇报).doc
- 层流车URS用户需求实用标准.doc
- 两个相干源产生波具有3b chmc.pdf
- wifi高性能模块ce fcc rohs认证-xrf report电子档.pdf
- linux pci multiport installation manual多端口安装手册.pdf
- 二章基本概念理论.pdf
- 互联网工程任务组ietfs turnerrfc5940 txt.pdf
- mtz0数学高级证明集合关系群mathematics paper 3 sets relations and groups french.pdf
- 网上书店需求规格说明.pdf
- 系统测试用例pcweb终端.pdf
- apv1 keil环境下接口兼容问题.pdf
- 2023届江西省八所重点中学高三3月联考英语试题(原卷版).docx
文档评论(0)