Precise mathematical modeling of the primary visual cortex (V1) is still a challenging problem. Due to the high similarity of visual system of cat and human, in this paper, we present a hybrid model to track the electrical responses of neurons that are measured by a multi-electrode array implanted in cat V1. The proposed model combines a stochastic phenomenological model with a multilayer leaky integrate-and-fire (LIF) model to predict V1 responses. Since all the existing visual cortex models do not capture the stochastic properties of synaptic changes, the proposed phenomenological model provides input currents for V1 by utilizing continuous chaotic neural equations with a quantization rule. Then a multilayer LIF model is presented to mimic the functions of lateral geniculate nucleus (LGN) and V1 neurons by their corresponding differential equations. The input current in these models is from the presynaptic neurons, which are computed using the LIF model. The LGN-V1 neuronal connections are adopted from previous studies, where the receptive fields (RFs) of LGN neurons converge onto elongated spatial structures that denote RFs of V1 neurons. The main purpose of this paper is to develop a short-term plasticity model that is more consistent with the nature of the LGN and V1 responses compared to state-of-the-art models. Previous studies have not incorporated the stochastic and quantized behaviors of neurons that in the recorded data of implemented electrodes. The experimental results show the ability of the proposed model to accurately predict spikes recorded experimentally, indicating the model outperforms the state-of-the-art method.