NewsomeTask.py 文件源码

def get_sqe_above(self,high_p,rich_return=False):
        p = dict()
        p['g_exc'] = high_p[0] # nS
        p['sen2dec_exc'] = high_p[1] # pA/Hz
        if len(high_p)>2:
            p['frac_proj'] = high_p[2]
        else:
            p['frac_proj'] = 0.1

        win1_all = np.zeros((6,6))
        for i_m in xrange(len(self.m_plot)):
            for i_c in xrange(len(self.c_plot)):
                m = self.m_plot[i_m]
                c = self.c_plot[i_c]
                win1_all[i_m,i_c] = self.quick_Newsome_model(p,m,c)

        sqe = 0
        sqe += sum((self.p_m_m-win1_all.mean(axis=1))**2)
        sqe += sum((self.p_m_c-win1_all.mean(axis=0))**2)

        if rich_return:
            return p
        else:
            return sqe