Question

我想增加确定性输入的噪音。

例如

# Input stimulus 
def stim(t): 
    if 0.0 < t < 1.0: 
        return 50.0 
    return 0.0

此函数在第一毫秒内返回一个脉冲-但我想添加一个随机噪声信号

当我这样做的时候

平均值= 0 标准= 1

# Input stimulus 
def stim(t): 
    if 0.0 < t < 1.0: 
        return 50.0 + np.random.normal(mean,std)
    return 0.0

那我真的在整个轴上都产生了白噪声？？但是我只添加一个随机数，而不是沿所有t进行处理。

目标是估计白噪声输入对总膜电压V（t）的影响，这是在ODE系统中计算的。

def compute_derivatives(y, t0): 
    dy = np.zeros((4,)) 

    Vm = y[0] 
    n = y[1] 
    m = y[2] 
    h = y[3] 


     # dVm/dt 
    GK = (gK / Cm) * np.power(n, 4.0) 
    GNa = (gNa / Cm) * np.power(m, 3.0) * h 
    GL = gL / Cm 

    dy[0] = (stim(t0) / Cm) - (GK * (Vm - VK)) - (GNa * (Vm - VNa)) - (GL * (Vm - Vl)) 

     # dn/dt 
    dy[1] = (alpha_n(Vm) * (1.0 - n)) - (beta_n(Vm) * n) 

    # dm/dt 
    y[2] = (alpha_m(Vm) * (1.0 - m)) - (beta_m(Vm) * m) 

     # dh/dt 
    dy[3] = (alpha_h(Vm) * (1.0 - h)) - (beta_h(Vm) * h) 

    return dy 

 # State (Vm, n, m, h) 
Y = np.array([0.0, n_inf(), m_inf(), h_inf()]) 


 # Solve ODE system 
 # Vy = (Vm[t0:tmax], n[t0:tmax], m[t0:tmax], h[t0:tmax]) 
Vy = odeint(compute_derivatives, Y, T)

如何添加噪声信号以及如何估计影响？也许是沿着t的RMS振幅？

白噪声输入

0 个答案: