我有一个依赖于 r 的符号函数，称为 u(r)。我从 u 中的微分方程得到这个函数。
否则，我有另一个符号函数，称为 sigma_r，它是 u 和 u 导数的函数。
当我显示 sigma 时，它显示为 u(r) 和 D(u)(r) 的函数。
u(r) 是 r 和 2 积分常数 C1 和 C2 的函数。
我想将 sigma_r 表示为 r、C1 和 C2 的函数。
我尝试 sigma_r = subs(sigma_r,u(r)) 但似乎 Matlab 无法替换它并且无法计算 u 导数。
这是脚本
u(r) = dsolve(diff(sigma_r) + (sigma_r - sigma_theta)/r + rho*w^2*r ==0,'IgnoreAnalyticConstraints', true) % In this differential equation sigma_r, sigma_theta are function of u, and Matlab replace their expression as function of u and solve for u. it works well.
sigma_r(r) = subs(sigma_r,u(r))
以下是 Matlab 在命令窗口中显示的内容：
u(r) =
C1/r^2.2107342132367193698883056640625 + 0.97655737574677914381027221679688*r^3 + C2*r^6.50236464850604534149169921875
sigma_r(r) =
(9671406556917033397649408*((33290247625219093223312621495090757213591748336283774753774513840*u(**0.97655737574677914381027221679688*(C7/r^2.2107342132367193698883056640625 + 0.97655737574677914381027221679688*r^3 + C8*r^6.50236464850604534149169921875)^3 +
如您所见，sigma_r 的表达式仍然包含 u。
谢谢您的帮助，
尼达尔·科赫拉德