找回密码
 欢迎注册
查看: 14291|回复: 4

[提问] 如何用mathematica求解如下的常微分方程

[复制链接]
发表于 2010-12-28 09:38:38 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
如何用mathematica求解如下的常微分方程

DSolve[{y''[t] + a*y'[t] + b*y[t] == 0}, y[t], t]
其中a^2-4*b<0
要求表达成三角函数正弦和余弦的线性组合的形式,
我不是不会求解这个方程,我只是用mathematica办到这事。
请给出命令吧
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2011-1-22 00:43:34 | 显示全部楼层
1# mathematica
  1. Refine[DSolve[{y''[t] + a*y'[t] + b*y[t] == 0}, y[t], t] // ComplexExpand, a^2 - 4 b < 0]
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2011-1-22 10:52:53 | 显示全部楼层
1# mathematica Refine[DSolve[{y''[t] + a*y'[t] + b*y[t] == 0}, y[t], t] // ComplexExpand, a^2 - 4 b < 0]
wayne 发表于 2011-1-22 00:43


wayne果然是大神呀!
我当时不会求解,然后就把a换成一个符号常数,比如E(2.71828.....),
b换成另外一个符号常数,比如Pi(3.1415926...)这样的话就会满足
a^2-4b<0,同时由于符号常数在进行运算的时候,不会被加减乘除运算“吸收”掉,
等运算完,我再把E和Pi换成我自己想要的数。
这是我的策略。
不过看起来还是wayne牛呀,真的很佩服,不知道wayne都是从什么地方学来的?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2011-1-23 00:11:05 | 显示全部楼层
4# mathematica
  1. Refine[y[t] /. DSolve[{y''[t] + a*y'[t] + b*y[t] == 0}, y[t], t] // ComplexExpand, a^2 - 4 b < 0] /. {C[1] -> (c1 + c2 I)/2, C[2] -> (c1 - c2 I)/2} // TrigReduce
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-5-18 18:44 , Processed in 0.056464 second(s), 18 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表