The analysis of the uterine contraction have become a general practice in an effort to improve the clinical management of uterine contractions during pregnancy and labour in human beings. The fluctuations in uterine activity may occur without affecting progress of gestation, however the painful and fashion contractions may be the first threat of misscarage. While pigs were considered as an referential preclinical model, the computational modelling of spontaneous myoelectrical activity of complex systems of porcine myometrium in peri-fertilization period has been proposed. The higher order statistic, multivariate cumulants and Joint Skewness Band Selection method, have been applied to study the dependence structure of electromyographic (EMG) signal with an effective EMG feature. Than the model of recognition of multivariate, myoelectricaly changes according to crucial stages for successful fertilization and early pregnancy maintenance has been estimated. We found that considering together time and frequency features of EMG signal was extremely non-Gaussian distributed and the higher order multivariate statistics such as cumulants, have to be used to determine the pattern of myoelectrical activity in reproductive tract. We confirmed the expectance that the probabilistic model changes on a daily base. We demonstrated the changes in proposed model at the crucial time points of in peri-fertilization period. We speculate the activity of the middle of uterine horn and the power (minimum and maximum) and pauses between myoelectrical burst features are essential for the functional role of uterine contractility in peri-fertilization period.