原创力文档- 1、本文档共3页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
微分方程符号解和数值解
问题1：求 满足初始条件 的符号解。
在文本框中输入：
s=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x)
运行得：s =
3*exp(-2*x)*sin(5*x)
问题2：求解微分方程 .先求符号解,再求数值解, 并作图进行比较.
在文本框中输入：
s=dsolve(Dy=(y^2-t-2)/(4*(t+1)),y(0)=2)
运行得：s =1+(t+1)^(1/2)
所以求得符号解s =1+(t+1)^(1/2)
为了求数值解,先编写M-文件fun.m
在文本框中输入：
function y=fun(t,y)
y=(y^2-t-2)/(4*(t+1));
保存，再运行新命令，在文本框中输入：
clear;close;t=0:10;
y=1+(t+1).^(1/2);
plot(t,y) %画符号解的图形
hold on %保留已画好的图形,
[t,y]=ode45(fun,[0,10],2);
plot(t,y,ro); %画数值解图形,用红色小圆圈
xlabel(t),ylabel(y)
运行结果见下图, 可见符号解和数值解吻合得很好.
问题3：求解微分方程 .先求符号解,再求数值解, 并进行比较.
在文本框中输入：
s=dsolve((1+x^2)*D2y=2*x*Dy,y(0)=1,Dy(0)=3,x)
运行得： s =1+x^3+3*x
所以求得符号解s =1+x^3+3*x
为了求数值解,先编写M-文件fun.m
在文本框中输入：
function z=ill(x,y)
z=[2*x*y(1)/(1+x^2),y(1)];
保存，再运行新命令，在文本框中输入：
clear
x=0:0.5:10;
y=1+x.^3+3*x;
plot(x,y)
hold on
xs=[0,10];
y0=[3,1];
[x,y]=ode45(ill,xs,y0);
plot(x,y(:,2),ro)
运行结果见下图, 可见符号解和数值解吻合得很好.
文档评论（0）