Wolfram mathematica 在mathematica中求解方程组中的问题
我已经为我的问题写了这组动力学方程 聚变中的磁化目标。我只是想知道为什么Mathematica 无法解决它。 我以前做过这个“NDSolve”，我得到了很好的答案。但是当我改变我的公式时，它不能解决什么问题。我可以把全部代码发给你Wolfram mathematica 在mathematica中求解方程组中的问题,wolfram-mathematica,Wolfram Mathematica,我已经为我的问题写了这组动力学方程 聚变中的磁化目标。我只是想知道为什么Mathematica 无法解决它。 我以前做过这个“NDSolve”，我得到了很好的答案。但是当我改变我的公式时，它不能解决什么问题。我可以把全部代码发给你 `bal = {(3/2)*ne[t]*k*Te'[t] == \[Eta]d*wd + wie - wb + f\[Alpha]*\[Eta]f*wf - whe, (3/2)*ni[t]*k* Ti'[t] == (1 - \[Eta]d)*
`bal = {(3/2)*ne[t]*k*Te'[t] == \[Eta]d*wd + wie - wb +
f\[Alpha]*\[Eta]f*wf - whe, (3/2)*ni[t]*k*
Ti'[t] == (1 - \[Eta]d)*wd - wie + f\[Alpha]*(1 - \[Eta]f)*wf -
whi, nd'[t] = -nd[t]*nT[t]*\[Sigma],
nT'[t] = -nd[t]*nT[t]*\[Sigma],
n\[Alpha]'[t] = nd[t]*nT[t]*\[Sigma], Te[0] = 1, Ti[0] = 1,
nd[0] == nT[0] == \!\(TraditionalForm\`
\*FractionBox[\(1.4447999999999998`*^26\), \(2\)]\), n\[Alpha][0] = 0}
sol = NDSolve[bal, {Te, Ti, nd, nT, n\[Alpha]}, t]
NDSolve::deqn: Equation or list of equations expected instead of -
((2.85474*10^-12 E^(-19.983 ((1<<1>>Plus[<<3>>]<<1>>Power<<1>>
<<1>>])/Ti[t])^(1/3)) nd[t] nT[t])/(Ti[t]^(2/3) (1-(15.136 Ti[t]+4.6064
Ti[<<1>>]^2-0.10675 Ti[<<1>>]^3)/(1000+75.189 Ti[<<1>>]+13.5
Power[<<2>>]+0.01366 Power[<<2>>]))^(5/6))) in the first argument
{4.8*10^-9 nd[t] (Te^\[Prime])[t]==900000000000000000-8.70051*10^-25
(nd[t]+nT[t])^2 Sqrt[Te[t]]-(5.2266*10^46 <<1>>^<<1>> (11.92
+1.69505*10^-9 <<1>>^3))/((nd[t]+nT[t]) (3.77 +<<21>>
<<1>>+1.32084*10^-19 Power[<<2>>]))+(8.7331*10^17 (24-Log[Times[<<4>>]])
nd[t]^2 (-Te[t]+Ti[t]))/(1.09626*10^24 Te[<<1>>]+5.97059*10^20
<<1>>)^(3/2)+(5.152*10^-16 E^(-19.983 Times[<<2>>]^(1/3)) nd[t] nT[t]
(8/3 (4.32916*10^-7+Times[<<3>>])+64/9 Plus[<<2>>]^2))/((1+104/27
Plus[<<2>>]+64/9 Power[<<2>>]) (32+Te[t]) Ti[t]^(2/3) (1-Plus[<<3>>]
Power[<<2>>])^(5/6)),<<7>>,0}.`enter code here`
NDSolve:：deqn:预期的方程式或方程式列表，而不是-
（（2.85474*10^-12 E^（-19.983）（（1个像素[]功率
])/Ti[t]^（1/3））nd[t]nT[t]/（Ti[t]^（2/3）（1-（15.136 Ti[t]+4.6064
Ti[]^2-0.10675 Ti[]^3）/（1000+75.189 Ti[]+13.5
第一个参数中的幂[]+0.01366幂[]）^（5/6）））
{4.8*10^-9第[t]（Te^\[Prime]）[t]==90000000000000000-8.70051*10^-25
(nd[t]+nT[t]^2平方米[Te[t]]-(5.2266*10^46)(11.92)
+1.69505*10^-9^3）/（nd[t]+nT[t]）（3.77+
+1.32084*10^-19幂[]）+（8.7331*10^17（24对数[次[]））
nd[t]^2（-Te[t]+Ti[t]）/（1.09626*10^24 Te[]+5.97059*10^20
)^（3/2）+（5.152*10^-16 E^（-19.983倍[]^（1/3））nd[t]nT[t]
（8/3（4.32916*10^-7+倍[]）+64/9加[]^2））/（1+104/27
加【】+64/9次方【】）（32+Te[t]）Ti[t]^（2/3）（1加【】
Power[]）^（5/6）），0}。`在此处输入代码`
==传统形式的k=1;\[Eta]d=1;wd=1;wie=1;wb=1;f\[Alpha]=1;\[Eta]f=1;wf=1;whe=1;whi=1;\[Sigma]=1;
ne[t_]:=2t+1;ni[t_]:=3t+2;
bal = {(3/2)*ne[t]*k*Te'[t] == \[Eta]d*wd + wie - wb + f\[Alpha]*\[Eta]f*wf - whe,
(3/2)*ni[t]*k*Ti'[t] == (1 - \[Eta]d)*wd - wie + f\[Alpha]*(1 - \[Eta]f)*wf - whi,
nd'[t] == -nd[t]*nT[t]*\[Sigma],
nT'[t] == -nd[t]*nT[t]*\[Sigma],
n\[Alpha]'[t] == nd[t]*nT[t]*\[Sigma],
Te[0] == 1,
Ti[0] == 1,
nd[0] == nT[0] == 1.4447999999999998`*^26/2,
n\[Alpha][0] == 0};
sol = NDSolve[bal, {Te,Ti,nd,nT,n\[Alpha]}, {t,0,1}];
Plot[{Te[t],Ti[t]}/.sol[[1]],{t,0,1}]
现在用你的实际值替换所有这些变量，用你的实际函数替换ne[t]和ni[t]，看看你得到了什么。我把一些
==传统的形式k=1;\[Eta]d=1;wd=1;wie=1;wb=1;f\[Alpha]=1;\[Eta]f=1;wf=1;whe=1;whi=1;\[Sigma]=1;
ne[t_]:=2t+1;ni[t_]:=3t+2;
bal = {(3/2)*ne[t]*k*Te'[t] == \[Eta]d*wd + wie - wb + f\[Alpha]*\[Eta]f*wf - whe,
(3/2)*ni[t]*k*Ti'[t] == (1 - \[Eta]d)*wd - wie + f\[Alpha]*(1 - \[Eta]f)*wf - whi,
nd'[t] == -nd[t]*nT[t]*\[Sigma],
nT'[t] == -nd[t]*nT[t]*\[Sigma],
n\[Alpha]'[t] == nd[t]*nT[t]*\[Sigma],
Te[0] == 1,
Ti[0] == 1,
nd[0] == nT[0] == 1.4447999999999998`*^26/2,
n\[Alpha][0] == 0};
sol = NDSolve[bal, {Te,Ti,nd,nT,n\[Alpha]}, {t,0,1}];
Plot[{Te[t],Ti[t]}/.sol[[1]],{t,0,1}]