J = 0.05;
B = 0.02;
Tload[t_] := 0.0;
R1 = 2;
\[Alpha] = 30*\[Pi]/180;
d = 0.05;
g = 0.001;
Nturns = 200;
\[Mu]0 = 4*\[Pi]*10^(-7);
r = 0.03;
lm = 0.02;
d = 0.1;
W = Sqrt[2*r^2 - 2*r^2*Cos[2*a]];
Rgmax = g/(\[Mu]0*r*d);
Rm = lm/(\[Mu]0*W*d);
square[x_, x1_, x2_] := UnitStep[x - x1]*(1 - UnitStep[x - x2])
Rg[\[Theta]_] :=
Rgmax/((2 \[Alpha] - Mod[\[Theta][t], \[Pi]])*
square[Mod[\[Theta][t], \[Pi]], -2 \[Alpha],
2 \[Alpha]] + (2 \[Alpha] - Mod[-\[Theta][t], \[Pi]])*
square[Mod[-\[Theta][t], \[Pi]], -2 \[Alpha], 2 \[Alpha]])
L[\[Theta]_] :=
Nturns^2/(
2*Rg[\[Theta]] +
Rm) *(square[Mod[\[Theta][t], \[Pi]], -2 \[Alpha], 2 \[Alpha]] +
square[Mod[-\[Theta][t], \[Pi]], -2 \[Alpha], 2 \[Alpha]])
Plot[L[\[Theta]], {\[Theta][t], -8 \[Alpha], 8 \[Alpha]},
PlotRange -> All]
我很抱歉,但我无法粘贴,因为mathematica可以显示(这里是新手)。然而,我尝试了几种解决方案来摆脱L0的1/0 inf表达,但我没有得到它。
请复制并粘贴到笔记本中并运行它。
我将Lθ乘以我的平方函数以使未定义的区域等于零(f.e 2a,Pi-2a应该为零作为其他区间)但是它不起作用。
如何正确识别此功能?
提前致谢。