h***o
发帖数: 539 | 1 BBS水木清华站∶精华区
发信人: FangQ (奥萨马·本·拉登), 信区: MathTools
标 题: Mathematica函数及使用方法
发信站: BBS 水木清华站 (Wed Nov 18 21:35:55 1998)
Mathematica函数及使用方法
—————————————————————————————————————
四、解方程
Solve[eqns, vars] 从方程组eqns中解出vars
Solve[eqns, vars, elims] 从方程组eqns中削去变量elims,解出vars
DSolve[eqn, y, x] 解微分方程，其中y是x的函数
DSolve[{eqn1,eqn2,...},{y1,y2...},x]解微分方程组，其中yi是x的函数
DSolve[eqn, y, {x1,x2...}] 解偏微分方程
Eli |
|
b*****y
发帖数: 163 | 2 The following is the full turn-out in matlab 6.0:
>> y=dsolve('D2y=a*exp(4*y)+b')
Warning: Explicit solution could not be found; implicit solution returned.
> In /usr/local/matlab6/toolbox/symbolic/dsolve.m at line 292
y =
[ 2*Int(1/(2*a*exp(4*a)+8*b*a+4*C1)^(1/2),a=``..y)-t-C2=0,
-2*Int(1/(2*a*exp(4*a)+8*b*a+4*C1)^(1/2),a=``..y)-t-C2=0]
>> pretty(y)
[ y
[ /
[ | 1
[2 | -------------------------------- da - t - C2 = 0 ,
[ | |
|
s*******g
发帖数: 18 | 3 碰到了一个问题, 6个second order ODE方程组,6个未知数其中有时间的函数,p( t ), q(
t ), 方程中的未知数有diff(q( t ),t),和diff(q( t ),t,t),diff(p( t ),t,t) (二
阶导数),外加另外4个未知量a,b,c,d.
我求解的语句如下
dsolve({L1,L2,L3,L4,q(0)=0,p(0)=0},{p( t ),q( t ),a,b,c,d},numeric);
运行报错 Error, (in dsolve/numeric/process_input) invalid argument: {p( t ),q(
t ),a,b,c,d}
不知道为什么, 有人能给个hint吗? 弄了2天了,郁闷死了
附件是 |
|
l********t
发帖数: 878 | 4 mathematic 5.0 -_-b
In[4]:=
DSolve[x'[t]==-a*x[t]+a*\[ExponentialE]^(b*t),x,t]
Out[4]=
a*e^{-at+(a+b)t}/(a+b) + c1*e^{-at} |
|
b*****y
发帖数: 163 | 5 try matlab symbolic toolbox:
y=dsolve('D2y=a*exp(4*y)+b'); |
|
h*******a
发帖数: 41 | 6 比如用matlab的ode45
或maple的 dsolve
只想加边条件 y(0)=y(1),Dy(0)=Dy(1)
好像不work，大虾请指教！ |
|
s***p
发帖数: 189 | 7 请教各位数学高手
对于方程
ax''+bx'+c=0
c=f(x)
我就想推导x'和x,a,b,c的关系。
我如下用MATLAB求解的过程对吗？
let y=x'
ay'+by+c=0
y=dsolve('a*Dy'+b*y+c=0','y(0)=0')
==>
y=-c/b+exp(-b*t/a)*c/b
so
x'=-c/b+exp(-b*t/a)*c/b
=(c/b)*(exp(-b*t/a)-1)
thanks |
|
s***p
发帖数: 189 | 8 请教各位数学高手
对于方程
ax''+bx'+c=0
c=f(x)
我就想推导x'和x,a,b,c的关系。
我如下用MATLAB求解的过程对吗？
let y=x'
ay'+by+c=0
y=dsolve('a*Dy'+b*y+c=0','y(0)=0')
==>
y=-c/b+exp(-b*t/a)*c/b
so
x'=-c/b+exp(-b*t/a)*c/b
=(c/b)*(exp(-b*t/a)-1)
thanks |
|
c*******v
发帖数: 2599 | 9 你首先应该明白一个基本事实:
是你首先讽刺了我们学数学的不回答你的"问题"。
然后我才被迫反击的.
靠，这种烂人
你不能回答就别回答，反而讥刺别人问题不对
脑子白痴到底
回答个问题就那么难？
因为含有隐函数，没有去耦合，所以不行
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~```~~
你不知道dsolve能解哪种方程么？
会不会去看看帮助再来问问题?
还隐函数呢......
就TMD这么个简单答案，到你那里就谁能回答？
真是越学越白痴 |
|
J***l
发帖数: 149 | 10 用Dsolve解一个Riccati ODE,
解里面含有类似tan(sqrt(something)),
如果我想在解方程时设定“something"<0应该怎么设？
就是说我想得到something<0时的解
求助，谢谢 |
|
x********g
发帖数: 595 | 11 用mathematica,
先
DSolve[{m1'[t] == - q m1[t] - c1 m1[t] + p1 c2 m2[t],
m2'[t] == -q m2[t] - c2 m2[t] + c1 m1[t] + c3 m3[t],
m3'[t] == -q m3[t] - c3 m3[t] + p2 m2[t]}, {m1[t], m2[t], m3[t]},
t]
然后用FindFit函数. |
|
s***e
发帖数: 911 | 12 Try:
tmp=DSolve[{Dt[f[x], x] == -k1 f[x] + k1 g[x], Dt[g[x], x] == -k2 g[x] +k1 f[x],
f[0] == a, g[0] == b}, {f[x], g[x]}, x];
Simplify[tmp]. |
|
