は公表数字に私はMATLABで遅れて微分方程式を解くためにしようとしている

Question

を再現するためのMatlabで遅延微分方程式を解きますこの方程式をシミュレートして、さまざまなパラメータ（すなわち、異なる時間遅延、ODEなどとの比較）でどのように振る舞うかを見て回ることができます。私はコード例hereの成功を制限しています。これまで

私のコードは次のとおりです。

function General_mRNA_DDE 
    sol = dde23(@General_mRNA_DDE2,2,@input_function,[0,5]) 


    figure; 
    plot(sol.x,sol.y) 


    function dydt = General_mRNA_DDE2(t,y,z) 
    k0=1; 
    k=10; 
    mRNA0=1; %initial concentration of mRNA 
    gamma=0.1; 
    z 
    dydt= [k0 + k*z - gamma*y]; 
    end 

    function hist = input_function(t) 
     hist = 1; 
    end 


    end 

しかし、私は基本的に非常に急な指数曲線のように見えてきました。ここで私が再現しようとしているものです：

この論文のDOIから：10.15252/msb.20177554（http://msb.embopress.org/content/msb/13/5/928.full.pdf）

誰が私は正確に姿を再現するために何かアドバイスはありますか？誘導体および未知のmRNAの値として、これは時間遅延微分方程式ではない事前

Answer 1

に

おかげで、同じ時間から取られます。アクチベーター値がその前の時間のmRNA値に依存しない場合、制御機能の値が遅延時間からのものであることは重要ではない。

新しい微分方程式

(exp(gamma*t) * mRNA(t))' = exp(gamma*t) * (k0 + k*Activator(t-delta_t)) 

は、ESP、簡単な統合によって解決することができるように。あなたは積分因子exp(gamma*t)を適用することができますActivator関数が区分的に定数である場合。