mathematica求定积分
mathematica吧
全部回复
仅看楼主
level 1
ShellyGary 楼主
有没有大神能帮我看看呀~
最终目标是想算一个分段函数的积分,但是我写的算不出来。现在把这个分段函数分开计算,就出现了以下问题,不知道怎么解决[泪][泪][泪]具体的代码已经贴在主楼了~先谢过各位!
2022年08月11日 03点08分 1
level 1
ShellyGary 楼主
fj1 = Function[r, ArcCos[(r^2 + l2^2 - r2^2)/(2*r*l2)]];
fj2 = Function[r, ArcCos[(r^2 + l2^2 - r2^2)/(2*r*l2)]];
fj3 = Function[r,
ArcCos[(r^2 + (r1 + h1 + r2)^2 - r3^2)/(2*r*(r1 + h1 + r2))]];
*以上是三个角度的函数*
fA1 = Function[r, 2*Pi*r*kc];
fai1 = Integrate[Q/fA1[r], {r, r1, R}] // TraditionalForm
*第一段分段函数的积分,出现的结果会出现“if condition”,不知道是什么意思*
fA21 = Function[r, 2*r*(Pi - fj1[r])];
fA22 = Function[r, 0];
fA2 = Function[r, kc*fA21[r] + n*kc*fA22[r]];
fai2 = Integrate[Q/fA2[r], {r, l1, R}] // TraditionalForm
“第二段分段函数的积分,已经无法积出来”
fA31 = Function[r, 2*r*(fj2[r] - fj1[r])];
fA32 = Function[r, 2*r*(Pi - fj3[r])];
fA3 = Function[r, kc*fA31[r] + n*kc*fA32[r]];
fA41 = Function[r, 2*r*fj3[r]];
fA42 = Function[r, 2*r*(Pi - fj3[r])];
fA4 = Function[r, kc*fA41[r] + n*kc*fA42[r]];
fA5 = Function[r, 2*Pi*r*n*kc];
Piecewise[{{Q/fA1, r1 <= r < l1}, {Q/fA2, l1 <= r < l2}, {Q/fA3,
l2 <= r < l3}, {Q/fA4, l3 <= r < l4}, {Q/fA5, l4 <= r < l5}}]
*本来想直接写分段函数的,但是这样写运行不出来*
2022年08月11日 06点08分 2
吧务
level 15
…………总之你先查 ConditionalExpression 的帮助。
2022年09月02日 17点09分 3
1