找回密码
 注册
Simdroid-非首页
查看: 146|回复: 7

[符号计算] 求解微分方程

[复制链接]
发表于 2008-7-8 14:40:48 | 显示全部楼层 |阅读模式 来自 上海
用mathematica5.0求解一常微分方程,方程形式见附件,但求解结果总是把输入的信息再重新显示一遍。是不是求解不出来呀?还有别的命令能求解这个方程吗?
谢谢!
file:///C:/DOCUME%7E1/ADMINI%7E1/LOCALS%7E1/Temp/moz-screenshot.jpgfile:///C:/DOCUME%7E1/ADMINI%7E1/LOCALS%7E1/Temp/moz-screenshot-1.jpgfile:///C:/DOCUME%7E1/ADMINI%7E1/LOCALS%7E1/Temp/moz-screenshot-2.jpg

本帖子中包含更多资源

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

×
发表于 2008-7-8 15:08:33 | 显示全部楼层 来自 美国
Simdroid开发平台
Mathematica里面变量不能用下标。 你把R_0改成R0,以及alpha_0改成alpha0就好了。
不过得到的解很长,没什么意义。不利于理解结果的意义。
不如先简化以及给定一些初始值,再解方程。
回复 不支持

使用道具 举报

 楼主| 发表于 2008-7-10 09:32:08 | 显示全部楼层 来自 上海
试了一下还是不行。
可能是因为微分方程的系数是变量,并不是常系数的微分方程,所以无法求解。
简化的方法也找了好多,都凑不成一个常系数的微分方程的形式,
郁闷至极!
回复 不支持

使用道具 举报

发表于 2008-7-10 12:40:28 | 显示全部楼层 来自 美国
根据我说的,能算出来的,大概要3分钟的样子(Timing给的是90.609秒)。不用我在给你打一遍吧。

  1. eq = Cos[\[Alpha]]^2/R - Cos[\[Alpha]0]^2/R0 -
  2.   1/R (Cos[\[Alpha]]^2 ((
  3.        2 Cos[\[Alpha]]^2 Cos[\[Theta]] \[Rho]'[\[Theta]])/R + (
  4.        Sin[\[Theta]] \[Rho]''[\[Theta]])/R)) +
  5.   P Sin[\[Alpha]] R/
  6.    B (1 + (Sin[\[Theta]] \[Rho][\[Theta]])/R - (
  7.      Cos[\[Alpha]]^2 Cos[\[Theta]] \[Rho]'[\[Theta]])/R)
  8. DSolve[eq == 0, \[Rho][\[Theta]], \[Theta]] // Timing
复制代码

评分

1

查看全部评分

回复 不支持

使用道具 举报

发表于 2008-7-10 16:41:33 | 显示全部楼层 来自 江苏无锡
I think so big result is not useful at all.

It's two D equation,  " Sin[\[Theta]] \[Rho]''[\[Theta]])/R))
Check Dsolve help and guide/EquationSolving TUTORIALS

Better plot it first for understanding than hard CPU solve.
回复 不支持

使用道具 举报

 楼主| 发表于 2008-7-11 10:10:38 | 显示全部楼层 来自 上海
smarten,的确是解不出看来,我把结果贴上来,计算时间不到一秒钟,0.375second.是不是设置有问题,材鸟,学习中!

本帖子中包含更多资源

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

×
回复 不支持

使用道具 举报

发表于 2008-7-11 11:12:30 | 显示全部楼层 来自 美国
你的Mathematica是什么版本,我算出来了,不过就是没什么意义,除非你先给一些参数赋值。
回复 不支持

使用道具 举报

 楼主| 发表于 2008-7-11 20:53:43 | 显示全部楼层 来自 上海
我的版本是5.0的,pojieban。该问题属于边值问题,求解的其实是一个挠度曲线,有四个边界条件,可以确定相关参数。为什么我的版本不能求解呢?我同学装的也不能求解,很奇怪。 smarten 你的版本是多少?正版?有没有需要设置的地方?
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-27 19:40 , Processed in 0.048025 second(s), 17 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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