通过Matlab ode45设计求解可激励电路

Question

我正在研究Excitable Dynamical系统的非线性电路。简单地说，我们可以使用初始条件通过Matlab“Ode45”求解其非线性方程。但我需要通过直接输入（外部刺激;不是初始条件）来刺激它。如何在第一个等式（DV）中将方波脉冲定义为电流源（具有优选的时间开始）？

clear;close all;clc 
gfi = 0.2308;      gso = 0.769;        Efi = 17.333;    Eso = 0;  TAu=3; 

Km = 0.842;         v12m = 5.85;       Kn = 4.4;      v12n = 9.667;  
C =0.01;  a=6;

DV  =  @(v,no,ni)   (-1/C)*((gfi*(v-Efi))/(1+exp((v12m-v)/Km))+gso*(no*ni)*(v-Eso)); 
Dno =  @(v,no,TAu) (1/TAu)*((1/(1+exp((v12n-v)/Kn)))-no); 
Dni =  @(v,ni,TAu) (1/a*TAu)*((1/(2+exp((v12n-v)/-Kn)))-ni); 

EDS= @(v,no,ni,TAu)[DV(v,no,ni);Dno(v,no,TAu);Dni(v,ni,TAu)]; 

[T,Y] = ode45(@(t, x) EDS(x(1), x(2) ,x(3),TAu),[0 25], [11,0.1,0.90]);