符号微分方程的求解.doc

  1. 1、本文档共3页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
符号微分方程的求解

6.6 符号微分方程的求解 6.6.3 微分方程符号解示例 【 * 例 6.6.3 -1 】求 的解。 S=dsolve(Dx=y,Dy=-x); disp([blanks(12),x,blanks(21),y]),disp([S.x,S.y]) x y [ cos(t)*C1+sin(t)*C2, -sin(t)*C1+cos(t)*C2] 【 * 例 6.6.3 -2 】图示微分方程 的通解和奇解的关系。 y=dsolve(y=x*Dy-(Dy)^2,x) % 求微分方程解 clf,hold on,ezplot(y(2),[-6,6,-4,8],1) % 画奇解 cc=get(gca,Children); % 取奇解曲线的图柄 3 set(cc,Color,r,LineWidth,5) % 把奇解画成粗红线 4 for k=-2:0.5:2;ezplot(subs(y(1),C1,k),[-6,6,-4,8],1);end% 画通解 hold off,title(\fontname{ 隶书 }\fontsize{16} 通解和奇解 ) y = [ x*C1-C1^2] [ 1/4*x^2] 图 6.6.3 -2 通解和奇解曲线 ? 【 * 例 6.6.3 -3 】求解两点边值问题: 。(注意:相应的数值解法比较复杂)。 y=dsolve(x*D2y-3*Dy=x^2,y(1)=0,y(5)=0,x) y = -1/3*x^3+125/468+31/468*x^4 【 * 例 6.6.3 -4 】求边值问题 的解。(注意:相应的数值解法比较复杂)。 S=dsolve(Df=3*f+4*g,Dg=-4*f+3*g,f(0)=0,f(3)=1) S.f,S.g S = f: [1x1 sym] g: [1x1 sym] ans = exp(3*t)*sin(4*t)/sin(12)/(cosh(9)+sinh(9)) ans = exp(3*t)*cos(4*t)/sin(12)/(cosh(9)+sinh(9))

文档评论(0)

ds2fdsx + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档