Question

我想在Maxima内的微分方程中更改变量，例如：

(%i1) diff(y(x),x);
                                   d
(%o1)                              -- (y(x))
                                   dx
(%i2) subst([x=a*z],%);
                                  d
(%o2)                            ---- (y(a z))
                                 da z

但是我不知道如何继续并获得diff(y(a*z),z)/a。 Maxima手册mentions the question，但指的是未处理此主题的at命令。

Answer 1

在maxima中有几种方法可以解决这个问题，让我们使用pdiff。首先我们在最大值中解决两个颂歌

(%i1) eq:'diff(y,x) = -y;
                                   dy
(%o1)                              -- = - y
                                   dx
(%i2) ode2(eq, y, x);
                                          - x
(%o2)                            y = %c %e
(%i3) sol: ic1(%, x= 1, y= 8);
                                         1 - x
(%o3)                            y = 8 %e

maxima对复杂的颂歌没有问题：

(%i4) eqcv : 'diff(y,x) = (x*y) / (3*x^2 - y^4);
                                dy      x y
(%o4)                           -- = ---------
                                dx      2    4
                                     3 x  - y
(%i5) ode2(eqcv, y, x);
                                   4    2
                                  y  - x
(%o5)                           - ------- = %c
                                      6
                                   2 y
(%i6) sol: ic1(%, x= 2, y= 1);
                                    4    2
                                   y  - x    3
(%o6)                            - ------- = -
                                       6     2
                                    2 y

但是让我们使用pdiff

在这个等式上使用变量

y = z ^ a

(%i1) load(pdiff)$

(%i5) eqcv : diff(y(x),x) = (x*y(x)) / (3*x^2 - y(x)^4);

                                         x y(x)
(%o5)                       y   (x) = ------------
                             (1)         2    4
                                      3 x  - y (x)

(%i6) eqcv, y(x) := z(x)^a;

                                                 a
                           a - 1              x z (x)
(%o6)                a z(x)      z   (x) = --------------
                                  (1)         2       4 a
                                           3 x  - z(x)

然后我们可以简化：

(%i7) % / z(x)^(a -1);
                                          x z(x)
(%o7)                     a z   (x) = --------------
                             (1)         2       4 a
                                      3 x  - z(x)
(%i8) % / a;
                                         x z(x)
(%o8)                    z   (x) = ------------------
                          (1)            2       4 a
                                   a (3 x  - z(x)   )

现在我们可以应用a = 1/2

(%i50) %, a = 1/2;
                                        2 x z(x)
(%o50)                      z   (x) = ------------
                             (1)         2    2
                                      3 x  - z (x)

开始解决：

(%i51) eq : %;
                                        2 x z(x)
(%o51)                      z   (x) = ------------
                             (1)         2    2
                                      3 x  - z (x)
(%i52) eq;
                                        2 x z(x)
(%o52)                      z   (x) = ------------
                             (1)         2    2
                                      3 x  - z (x)


(%i56) eq : convert_to_diff(eq);
                           d             2 x z(x)
(%o56)                     -- (z(x)) = ------------
                           dx             2    2
                                       3 x  - z (x)

(%i57) depends(z,x);
(%o57)                              [z(x)]

(%i58) eq : subst(z,z(x),eq);
                                dz     2 x z
(%o58)                          -- = ---------
                                dx      2    2
                                     3 x  - z

(%i59) ode2(eq,z,x);

                                   2    2
                                  z  - x
(%o59)                          - ------- = %c
                                     3
                                    z

您也可以将此应用于偏微分方程example from a mailing list：

我正在处理下面的偏微分方程   函数T（r，t）：

diff（T，t） - T ^ 2 / r ^ 2 * diff（（r ^ 2 * T ^（sigma）* diff（T，r）+ q /   T），r）= 0

我现在想介绍一个新变量

eta = r / t ^（1/3）

和一个新功能

Psi（η）= t ^（1/3）* T ^（1 + sigma）

我希望Maxima能够改变偏微分方程   对于T（r，t），将其描述为Psi（η）的微分方程。

溶液：

(%i1) load(pdiff)$
(%i2) de : diff(f(x,y),x,2) + diff(f(x,y),y,2)$
(%i3) de, f(x,y) := g(x/y);
(%o3) g[(2)](x/y)/y^2+(2*x*g[(1)](x/y))/y^3+(x^2*g[(2)](x/y))/y^4
(%i4) y^2 * %;
(%o4) (g[(2)](x/y)/y^2+(2*x*g[(1)](x/y))/y^3+(x^2*g[(2)](x/y))/y^4)*y^2
(%i5) ratsubst(w,x/y,%);
(%o5) (w^2+1)*g[(2)](w)+2*w*g[(1)](w)
(%i6) nde : convert_to_diff(%);
(%o6) 2*w*('diff(g(w),w,1))+(w^2+1)*('diff(g(w),w,2))
(%i7) depends(g,w)$
(%i8) nde : subst(g,g(w),nde);
(%o8) ('diff(g,w,2))*(w^2+1)+2*('diff(g,w,1))*w
(%i9) ode2(nde,g,w);
(%o9) g=%k1*atan(w)+%k2

所以在你的情况下：

(%i11) load(pdiff)$

(%i12) eq : diff(y(x),x);
(%o12)                              y   (x)
                                     (1)
(%i13) eq, y(x) := t(a*x);
(%o13)                            a t   (a x)
                                     (1)

Maxima CAS：如何用衍生形式代替自变量

1 个答案: