matlab - 絶対値で微分方程式を解く

Question

この形式の方程式を解きたい: x' = -Ax + B.|sin(100*pi*t)| そして、私は次のようにode45を使用します:

function find_x
t = 0:0.001:10;
x0 = 0;
R1 = 90000;
R2 = 1000;
C = 0.001;
[t,x]=ode45(@rhs, t , x0);

plot(t,x);

    function dxdt = rhs(t,x)
        dxdt = -(C/R1 + C/R2)*x + C/R1*abs(sin(100*pi*t)) ;
        %It's form is dx/dt = -A.x + B.U(t)
    end
end

しかし、それは私に間違った答えを与えると思います。実際には、次のような「ダイオードブリッジとコンデンサの後の出力電圧の形を見つける」という問題からこの式を取得します。

誰かがそれを解決する別の方法を私に提案できますか? ありがとう。