原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Mat lab微分方程的求解实验报告
《数学实验》报告
实验名称?Mat lab常微分方程的求解?
学 院 计算机与通信工程学院
专业班级? 计1103????????
姓??名????????
学??号???????
??
2013年6 月
【实验目的】
通过练习，熟悉Mat lab的求解常微分方程，函数文件的创建等。了解Mat lab的命令窗口及其基本操作和常用命令。
通过练习，熟悉Mat lab的一些基本操作，掌握符号解法和数值解法，以及其中常用的方法。
【实验任务】
求解微分方程y=xsin（x）/cos（y）。
用数值方法求解下列微分方程，用不同颜色和线形将y和y画在同一个图形窗口里：y+ty-y=1-2t，初始时间：t0=0；终止时间：tf=π；初始条件：y|t=0=0.1，y|t=0=0.2。
【实验程序】
题一：
y= dsolve(Dy=x*sin(x)/cos(y),x)
题二:
令：a=y，b=y=a，b=y则 ：
a=0*a+1*b+0;b=1*a-t*b+(1-2*t)；
[a;b]=[0 1;1 -t][a;b]+[0;1](1-2*t)；
故化为一阶微分方程有：
x=[0 1;1 -t]x+[0;1]u；其中：x=[a;b],u=(1-2*t)。
初始条件：当t=0时，a=0.1；b=0.2。
function xd = mainfun( t,x )
%UNTITLED Summary of this function goes here
% Detailed explanation goes here
u=1-2.*t;
xd=[0 1;1 -t]*x+[0;1]*u;
end
clf;
t0=0;
tf=pi;
x0t=[0.1;0.2];%初始条件
[t,x]=ode23(mainfun,[t0,tf],x0t)
y=x(:,1);
Dy=x(:,2);
plot(t,y,g--,t,Dy,k-.)
xlabel(t轴),ylabel(Y轴)
legend(原始函数y,一阶导数Dy)
【实验结果】
题一：
y =
asin(C3 + sin(x) - x*cos(x))
题二：
fun
t =
0
0.0145
0.0873
0.2014
0.3259
0.4621
0.6121
0.7778
0.9621
1.1482
1.2767
1.4053
1.5188
1.6706
1.8601
2.0891
2.3569
2.6546
2.9687
3.1416
x =
0.1000 0.2000
0.1030 0.2158
0.1214 0.2883
0.1598 0.3798
0.2116 0.4479
0.2756 0.4847
0.3485 0.4813
0.4245 0.4282
0.4939 0.3168
0.5388 0.1599
0.5513 0.0320
0.5466 -0.1069
0.5272 -0.2358
0.4780 -0.4127
0.3788 -0.6347
0.2037 -0.8952
-0.0742 -1.1797
-0.4677 -1.4646
-0.9691 -1.7290
-1.2793 -1.8586
【实验总结】
通过本次实验，熟悉了Mat lab的一些基本操作，通过练习，会求解常微分方程，并熟练掌握符号解法和数值解法，以及其中常用的方法，包括欧拉方法、梯形公式、龙格库塔方法等。
通过练习，熟悉Mat la
integrated energy, chemicals and textile Yibin city, are the three core pillars of the industry. In 2014, the wuliangye brand value to 73.58 billion yuan, t
文档评论（0）