这个资源是_entirely_ MATLAB。枫有时这里讨论因为枫以前用于MATLAB符号工具箱。
我将尝试你的系统在枫今晚当我回家,但这将是几个小时。
数学的问题
4视图(30天)
显示旧的评论
好吧,事情是这样的,数学不返回dif.equations这个系统解决方案
DSolve[{侯' [t] = = 2 * p12 [t] - (1 / r) * p12 [t] ^ 2,
p12的[t] = =
(第22位[t]) - w ^ 2 * 11赛[t] - 2 * z * w * p12 [t] - (1 / r) * p12 [t] *第22位[t],
第22位' [t] = = 2 * (w ^ 2) * p12 [t] -
4 * z * w *第22位[t] - (1 / r) *第22位[t] ^ 2 + (w ^ 4) * q,侯[0]= = 0,p12 [0] = = 0,
第22位[0]= = 0},{赛[t], p12 [t],第22位[t]}, t]
谁能请试穿你的系统并告诉我如果你有同样的问题吗?
我也尝试NDSolve函数
NDSolve[{侯' [t] = = 2 * p12 [t] - (1 / r) * p12 [t] ^ 2,
p12的[t] = =
(第22位[t]) - w ^ 2 * 11赛[t] - 2 * z * w * p12 [t] - (1 / r) * p12 [t] *第22位[t],
第22位' [t] = = 2 * (w ^ 2) * p12 [t] -
4 * z * w *第22位[t] - (1 / r) *第22位[t] ^ 2 + (w ^ 4) * q,侯[0]= = 0,p12 [0] = = 0,
第22位[0]= = 0},{赛[t], p12 [t],第22位[t]}, {t 0 10}]
但我得到错误消息:NDSolve:: ndnum:遇到非数值值导数在t = = 0. '
我非常感谢你的帮助
有没有。
(我使用mathematica 8.0)
我知道这主要是关于matlab的问题,但我希望应该没有问题
2的评论
安德烈亚斯•戈斯
2012年6月14日
一则轶事:我记得当我在做技术和安装支持自己,客户叫我们其他产品有时当我们建议给供应商打电话,客户(他还使用MATLAB对一些应用程序)说:“是金宝app的,我知道,但是你的服务要好得多”:-)下载188bet金宝搏
答案(1)
沃尔特·罗伯森
2012年6月14日
在枫,该表达式
dsolve ([diff(赛(t), t) = 2 * p12 (t) -p12 (t) ^ 2 / r, diff (p12 (t), t) =第22位(t) - w ^ 2 * 11赛(t) 2 * z * w * p12 (t) -p12 (t) *第22位(t) / r, diff(第22位(t), t) = 2 * w ^ 2 * p12 (t) 4 * z * w *第22位(t)第22位(t) ^ 2 / r + w ^ 4 * q, p12(0) = 0,侯(0)= 0,第22位(0)= 0])
枫说这没有解决方案。
如果你减少p12初始条件(0)= 0,让侯(0)和第22位(0)待定,则系统有两个解决方案,金宝搏官方网站
赛(t) = (2 * z * r + (4 * r ^ 2 * z ^ 2 + w ^ 2 * q * r) ^ (1/2)) / w
p12 (t) = 0
第22位(t) = w * (2 * z * r + (4 * r ^ 2 * z ^ 2 + w ^ 2 * q * r) ^ (1/2))
和
赛(t) = - z (2 * * r + (4 * r ^ 2 * z ^ 2 + w ^ 2 * q * r) ^ (1/2)) / w
p12 (t) = 0
第22位(t) = - w * (2 * z * r + (4 * r ^ 2 * z ^ 2 + w ^ 2 * q * r) ^ (1/2))}
(两个是相似的,但有一些符号的变化。)
如果一个一个检查这些解决方案然后t的注意,这金宝搏官方网站些都是独立的,所以如果一个强加侯(0)= 0,意味着侯(t)和第22位(t)都必须是0(不是不可能考虑到公式,但它需要r = 0或需要一些特殊的w和q的关系)
还有一个解决方案没有初始值条件如果p12的方程(t) = 2 * r。
还有第三个解决方案没有初始值的方程的条件相当复杂,似乎是递归地定义的。我不懂什么是枫想说在它的输出。在任何情况下显然从我所理解,结合{侯(0)= 0,p12(0) = 0,第22位(0)= 0}作为初始条件方程组无解。
