This research is performed based on the modeling of biological signals. We can produce Heart Rate (HR) and Heart Rate Variability (HRV) signals synthetically using the mathematical relationships which are used as input for the Integral Pulse Frequency Modulation (IPFM) model. Previous researches were proposed same methods such as one model of ECG signal synthetically based on RBF neural network, a model based on IPFM with random threshold, method was based on the estimation of produced signals which are dependent on autonomic nervous system using IPFM model with fixed threshold, a new method based on the theory of vector space that based on time-varying uses of IPMF model (TVTIPMF) and special functions, and two different methods for producing HRV signals with controlled characteristics and structure of time-frequency (TF) for using non-stationary HRV analysis. In this paper, several chaotic maps such as Logistic Map, Henon Map, Lorenz and Tent Map have been used. Also, effects of sympathetic and parasympathetic nervous system and an internal input to the SA node and their effects in HRV signals were evaluated. In the proposed method, output amount of integrator in IPFM model was compared with chaotic threshold level. Then, final output of IPFM model was characterized as the HR and HRV signal. So, from HR and HRV signals obtaining from this model, linear features such as Mean, Median, Variance, Standard Deviation, Maximum Range, Minimum Range, Mode, Amplitude Range and frequency spectrum, and non-linear features such as Lyapunov Exponent, Shanon Entropy, log Entropy, Threshold Entropy, sure Entropy and mode Entropy were extracted from artificial HRV and compared them with characteristics as extracted from natural HRV signal. Also, in this paper two patients that called high sympathetic Balance and Cardiovascular Autonomy Neuropathy (CAN) which is detected and evaluated by HRV signals were simulated. These signals by changing the values of the some coefficients of the normal simulated signal and with extracted frequency feature from these signals were simulated. For final generation of these abnormal signals, frequency features such as energy of low frequency band (EL), energy of high frequency band (HL), ratio of energy in low frequency band to the energy in high frequency band (EL/EH), ratio of energy in low frequency band to the energy in all frequency band (EL/ET) and ratio of energy in high frequency band to the energy in all frequency band (EH/ET) from abnormal signals were extracted and compared with these extracted values from normal signals. The results were closely correlated with the real data which confirm the effectiveness of the proposed model. Various signals derived from the output of this model can be used for final analysis of the HRV signals, such as arrhythmia detection and classification of ECG and HRV signals. One of the applications of the proposed model is the easy evaluation of diagnostic ECG signal processing devices. Such a model can also be used in signal compression and telemedicine application.