示例 4 步，在同一个会话中按顺序运行。每一步都在引擎上执行过，下面的输出就是它实际产生的。

第 1 步

vi[t_] = Sin[t];
vr[t_] = Abs[vi[t]]//PiecewiseExpand[#, t∈Reals]&;
Plot[{vi[t], vr[t]}, {t, 0, 16}]

输出

-Graphics-

第 2 步

system = {vo'[t] == charging[t]vr'[t] - (1 - charging[t])vo[t] / (r c), vo[0] == 0};
control = {charging[0] == 1, WhenEvent[vr'[t] < 0, charging[t] -> 0], WhenEvent[And[vo[t] < vr[t], charging[t] == 0, vr'[t] > 0], charging[t] -> 1]};

无输出——这一步是在为下一步做准备。

第 3 步

pars = {r -> 0.1, c -> 1000};
sol = NDSolve[{system, control} /. pars, vo, {t, 0, 16}, DiscreteVariables -> charging∈{0, 1}];

无输出——这一步是在为下一步做准备。

第 4 步

Plot[{vo[t] /. sol, vr[t]}, {t, 0, 16}]

输出

-Graphics-

在 Qurak 上运行全部操作示例

AC-to-DC full-wave rectifier

用到的函数

相关操作示例