数学实验matlab解常微分方程

数学实验matlab解常微分方程

ID:28060867

大小:636.00 KB

页数:12页

时间:2018-12-07

数学实验matlab解常微分方程_第1页
数学实验matlab解常微分方程_第2页
数学实验matlab解常微分方程_第3页
数学实验matlab解常微分方程_第4页
数学实验matlab解常微分方程_第5页
资源描述:

《数学实验matlab解常微分方程》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

1、GUIZHOUUNIVERSITY实验课程名称数学实验实验项目名称MATLAB解常微分方程年级2010级专业信计(101)学生姓名成富学1007010167理学院实验时间:20012年5月9日学生实验室守则一、按教学安排准时到实验室上实验课,不得迟到、早退和旷课。二、进入实验室必须遵守实验室的各项规章制度,保持室内安静、整洁,不准在室内打闹、喧哗、吸烟、吃食物、随地吐痰、乱扔杂物,不准做与实验N容无关的事,非实验用品一律不准带进实验室。三、实验前必须做好预习(或按要求写好预习报告),未做预习者不准参加实验。四、实验必须服从教师的安排和指导,认真按规程操作,未经教师允许不得擅

2、自动用仪器设备,特别是与本实验无关的仪器设备和设施,如擅自动用或违反操作规程造成损坏,应按规定赔偿,严重者给予纪律处分。五、实验中要节约水、电、气及其它消耗材料。六、细心观察、如实记录实验现象和结果,不得抄袭或随意更改原始记录和数据,不得擅离操作岗位和干扰他人实验。七、使用易燃、易爆、腐蚀性、有毒有害物品或接触带电设备进行实验,应特别注意规范操作,注意防护;若发生意外,要保持冷静,并及时向指导教师和管理人员报告,不得自行处理。仪器设备发生故障和损坏,应立即停止实验,并主动向指导教师报告,不得自行拆卸查看和拼装。八、实验完毕,应清理好实验仪器设备并放回原位,清扫好实验现场,经

3、指导教师检杳认可并将实验记录交指导教师检杳签字后方可离去。九、无故不参加实验者,应写出检查,提出申请并缴纳相应的实验费及材料消耗费,经批准后,方可补做。十、自选实验,应事先预约,拟订出实验方案,经实验室主任同意后,在指导教师或实验技术人员的指导下进行。十一、实验室内一切物品未经允许严禁带出室外,确需带出,必须经过批准并办理手续。姓名成富学号1007010167实验组实验时间2012-5-9指导教师刘轶中成缋实验项目名称MATLAB解微积分学生所在学院:理学院专业:信计班级:101实验目的及要求:实验目的:1:学习MATLAB有关常微分方程计算指令;2:学会常微分方程(组)的

4、数値求解Euler法和刚性方程组问题:3:分析数值'汁算法求解数学问题的局限性:4:建模实验研宂常微分方程。实验要求:在MATLAB环境下熟练运用各知识作图实验(或算法〉原理:根据函数的表达式,建立M文件规则MATLAB中语句的功能法输入相应的句子。建立M文件。运用数学函数中基本函数作图。、运用plot(x,y),ezplot(x,y)输入相应的变量。MATLAB将会根据输入的变作出相应图案实验硬件及软件平台:PCMATLAB实验步骤:1打幵MATLAB,进入到窗口状态。2根据要求输入要求的题目。3运用MATLAB特殊的函数编程。4检验程序是否正确。5运用结果。6分析得到的

5、实验数裾。实验内容(包括实验具体内容、算法分析、源代码等等):1(1):>>fun=inline(,x+y,,,x,,,y,);»[t,y]=ode45(fun,[0123],1)01231.00003.436611.778136.17111(2)。»fun=inline(,[-2*y(l)-3*y(2);2*y(l)+y(2)]*/t,y,);»[t,y]=ode45(fun,[010],[-2.7;2.8]);plot(y(:,l),y(:,2))»plot(y(:,l),y(:,2))1(3):»clear;close;fun=@(t,y)[y(2);0.01*y(

6、2)A2-2*y(l)+sin(t)];[t,y]=ode45(fun,[05],[0;l]);plot(t,y(:,1))1(4):»clear;close;fun=@(t,y)[y(2);2.5*y(2)+1.5*y(l)+45*exp(2*t)];[t,y]=ode45(fun,[02],[2;1]);plot(t,y(:,1))1(5):»clear;close;fun=@(t,y,mu)[y(2);-mu*(y⑴八2-1)*y(2)-y(l)];[t,y]=ode45(fun,[020],[2;0],[],l);mu=l.plot(y(:,l),y(:,2));h

7、oldon;[t,y]=ode45(fun,[020],[2;0],[],2);mu=2.plot(y(:,l),y(:,2),’r’);holdoff;mu=1mu=21(6):»clear;close;sinit=bvpinit(l:0.5:3,[2;0])odefun=inline(’[y(2);-2/t*y(2)+(2*y(l)+10*cos(log⑴))/t/t]',’t’,’y’);bcfun=inlineC[ya(l)-l;yb(l)-3],/ya,,,yb,);sol=bvp4c(odefun,bc

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。