mathematica 解微分方程组_mathematica吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：10,788贴子：
mathematica 解微分方程组收藏
求高手用软件或是手解这个微分方程组，其中B,C为已知的参数，本人思考有4,5天，始终解不出，坐等高手！！！
处女星号邮轮由上海出发前往大阪,畅享大阪自然美景和饕餮美食
两个方程四个未知数，怎么解？
M1(x)和M2(x)也是已知函数，要求的是q1（x）和q2（x）吗？
好像还真可以解。先把两个方程加起来，可得p1+p2满足形式y'[x] + f[x] y[x] == 0DSolve[y'[x] + f[x] y[x] == 0, y[x], x]可以给出通解，再分别带入——好像手算都行~~
用DSolve[{p12'[x] + (2 c + \[Mu]1[x] + \[Mu]2[x]) p12[x] ==
2 (p12[x]) b}, {p12[x]}, x]得p12=p1+p2再用DSolve[{Derivative[1][p1][x] + (c + \[Mu]1[x])*p1[x] ==
E^Integrate[2*b - 2*c - \[Mu]1[K[1]] - \[Mu]2[K[1]], {K[1], 1, x}]*C[1]*b},
{p1[x]}, x]
其实上面的是原本式子，我化简了下。带*是laplace表达式，能帮我把最后的式子算出来吗，最后我自己可以化简的，非常感谢！
请问一下问题出在哪里？
看来我二了下面换一种思路。做代换p1[x]-&p10 Exp[k[x]],p2[x]-&p20Exp[k[x]],其中p10,p20为常数。带入原方程组equ = {D[p1[x], x] + (c + u1[x]) p1[x] == (p1[x] + p2[x]) b,
D[p2[x], x] + (c + u2[x]) p2[x] == (p1[x] + p2[x]) b}equ2 = (equ /. {p1 -& (p10 Exp[k[#]] &), p2 -& (p20 Exp[k[#]] &)}) //
Simplify可得一关于p10,p20的二元齐次方程，令其行列式为0，得一关于k'[x]的二次方程。解之。带回，可得p1[x],p2[x]的两组线性无关解。
一组解是{p1 -& Function[x, E^(
1/2 \[Integral](2 b - 2 c - u1[x] - u2[x] - Sqrt[
4 b^2 + u1[x]^2 - 2 u1[x] u2[x] +
u2[x]^2]) \[DifferentialD]x)],
p2 -& Function[x, (1/b)
E^(1/2 \[Integral](2 b - 2 c - u1[x] - u2[x] - Sqrt[
4 b^2 + u1[x]^2 - 2 u1[x] u2[x] +
u2[x]^2]) \[DifferentialD]x) (-b + c + u1[x] +
1/2 (2 b - 2 c - u1[x] - u2[x] - Sqrt[
4 b^2 + u1[x]^2 - 2 u1[x] u2[x] + u2[x]^2]))]}
大家看一下哪里有问题：Clear["Global`*"];equ = {D[p1[x], x] + (c + u1[x]) p1[x] == (p1[x] + p2[x]) b,
D[p2[x], x] + (c + u2[x]) p2[x] == (p1[x] + p2[x]) b};k0 = k' /. Solve[(c + u1[x] - b + k') (k' + c + u2[x] - b) == b^2, k'];ks0 = {Integrate[k0[[1]], x], Integrate[k0[[2]], x]};solution11 = {p1 -& Function[x, Exp[ks0[[1]]] // Evaluate],
p2 -& Function[
x, (c + u1[x] - b + k0[[1]]) Exp[ks0[[1]]]/b // Evaluate]};solution22 = {p1 -& Function[x, Exp[ks0[[2]]] // Evaluate],
p2 -& Function[
x, (c + u1[x] - b + k0[[2]]) Exp[ks0[[2]]]/b // Evaluate]};(equ /. solution11) // Simplify(equ /. solution22) // Simplify
登录百度帐号推荐应用
为兴趣而生，贴吧更懂你。或}