Sensors Allocation and Observer Design for Discrete Bilateral Teleoperation Systems with Multi-Rate Sampling

This study addresses sensor allocation by analyzing exponential stability for discrete-time teleoperation systems. Previous studies mostly concentrate on the continuous-time teleoperation systems and neglect the management of significant practical phenomena, such as data-swap, the effect of sampling rates of samplers, and refresh rates of actuators on the system’s stability. A multi-rate sampling approach is proposed in this study, given the isolation of the master and slave robots in teleoperation systems which may have different hardware restrictions. This architecture collects data through numerous sensors with various sampling rates, assuming that a continuous-time controller stabilizes a linear teleoperation system. The aim is to assign each position and velocity signals to sensors with different sampling rates and divide the state vector between sensors to guarantee the stability of the resulting multi-rate sampled-data teleoperation system. Sufficient Krasovskii-based conditions will be provided to preserve the exponential stability of the system. This problem will be transformed into a mixed-integer program with LMIs (linear matrix inequalities). These conditions are also used to design the observers for the multi-rate teleoperation systems whose estimation errors converge exponentially to the origin. The results are validated by numerical simulations which are useful in designing sensor networks for teleoperation systems.


Introduction
Teleoperation systems have been extensively deployed in remote and hazardous environments during the last decades. A teleoperation system consists of master and slave mechanical parts where the master manipulates the slave via a human operator. A detailed survey of common definitions, applications, and developments of teleoperation systems is provided in [1,2]. Time delays mainly caused by data transfer across communication channels are considered, representing one of the most significant challenges in teleoperation systems since even a minor time delay can destabilize the entire teleoperation system. Motivated to solve the instability caused by time delays, a wide range of approaches have been proposed even in the recent literature for a continuous-time model of the bilateral systems. In [3], a modified wave variable compensation has been utilized to passivate the delayed communication channel. In [4], a robust controller based on a state observer is implemented on a haptic device and an asymmetric teleoperation system. The results demonstrate better tracking performance for force signals. A modified robust adaptive controller with a novel adaptive torque observer for nonlinear teleoperation systems has been proposed in [5] to overcome time-varying delays, external disturbances and uncertainties in the system dynamics. This structure is relaxed from acceleration measurement and does not require force sensors. In [6], the admittance control is applied to study the phase transition between the constraint and unconstraint motions in the teleoperator system with haptic feedback. Its stability proof is provided by hybrid systems theory and Lyapunov-Krasovskii approach. A modified motion prediction method is proposed in [7], based on a state observer with a cascade architecture. Under an adequate set of assumptions for the existence of uncertainties, the prediction errors remain bounded. These results, properly extended, have been applied to a scaled four-channel scheme under time-varying delays in [8]. Also, a novel adaptive-observer-based scaled controller for four-channel systems is presented in [9] to overcome asymmetric time-varying delays.
Besides these model-based methods, data-driven approaches have also been utilized in teleoperation systems. In [10], an adaptive neural network-based controller is designed in the existence of external disturbances, constant time delay, and internal frictions. Also, in [11], a truncated quantile critics reinforcement learning-based integrated structure is presented, leading to the quantifiable training feedback in a teleoperation system. A radial basis function neural network (RBFNN) is applied in [12] to estimate dynamic uncertainties in the existence of the actuator backlash-like hysteresis and time-varying delays.
Obviously, all of the studies above mainly have focused on the continuous-time teleoperation systems, neglecting the handling of some practical parameters such as dataswap, the sampling rate effect, and update rates of actuators on the system's stability. However, this investigation can be instructive in experimental studies such as [13], which proposes an open-source scheme of a wearable 7-DOF wireless operator arm motiontracking system for teleoperator systems or in [14], where a force tactile feedback approach for finger tracking systems is presented.
Although continuous-time teleoperator systems have received much attention, there are far fewer studies for sampled-data counterparts. However, some discrete-nature controllers have been reported for teleoperation systems in recent years. A review of discretetime control methods for sampled-data teleoperator systems has been presented in [15]. Also, a detailed comparison between the behavior of analog and sampled-data teleoperation systems is investigated in [16]. In [17], the data transfer from multiple slave robots to the master robot is regulated using the round-robin (RR) scheduling algorithm for a class of sampled-data teleoperation systems. In [18], an original method to teleoperate a discrete hyper-redundant robots is presented by utilizing mixed reality to carry out complex tasks in remote environments with more accuracy. In [19], a modified discrete sliding mode controller is proposed via a new reaching law to overcome the chattering phenomena besides disturbance rejection by utilizing an adaptive extended state observer.
A significant deficiency in previous literature on discrete-time teleoperations systems is that all transmitted signals are sampled with the same sampling rate. However, it is not a sensible assumption in the teleoperation systems due to the isolation of the master and slave robots, which may have different hardware restrictions; thus, a multi-rate approach has been proposed in this paper. In the multi-rate architecture, data are collected through numerous sensors with various sampling rates. The necessity of multi-rate sampling is not just due to different phenomena in the system. Sometimes utilizing different methods to sense the same phenomena provides unequal sampling rates in the system. As an instance in our presented study, the position signal of the master robot can be measured by either encoder, potentiometer or even a camera via image processing. Even if we relinquish applying different sensors, variable delays and data drop out in communication channels lead to different data arriving times for the controllers.
Stability analysis of the multi-rate systems and designing sampled-data controllers have received wide attention in recent years [20][21][22][23]. In an extension to our previous study in [24], which proposes Krasovskii-based conditions in the form of LMIs (Linear Matrix Inequalities) to stabilize the linear multi-rate teleoperation systems, this study focuses on sensor allocation with analyzing exponential stability. It is presumed that a continuous-time controller stabilizes the teleoperator system. The goal is to assign each position and velocity signals to sensors with different sampling rates and divide the state vector across sensors so that the multi-rate system remains exponentially stable. Different approaches for dividing the state vector across sensors can lead to stable or unstable systems. Thus, Krasovskiibased conditions satisfy exponential stability. Also, these conditions are employed to design observers on both sides of the teleoperator system. This problem will be transformed into a mixed-integer program with LMIs. The results of this paper are useful in designing sensor networks for teleoperation systems.
The remainder of this study is organized as follows. In Section 2, preliminaries and sampled-data model of the teleoperator system with multi-rate sampling approach are presented. The exponential stability analysis for sensor partitioning and observer design are provided in Sections 3 and 4, whereas in Section 5, these results are validated by numerical simulations. In Sections 6 and 7, the results are reported and concluding remarks stated.

Proposed Sampled-Data Model with Multi-Rate Architecture for Teleoperation System
In a bilateral teleoperation system, when a human manipulates the master robot, the generated motion signals are relayed to the slave robot via a wireless network, and the slave robot subsequently performs the same action. The general schematic of the teleoperation system is depicted in Figure 1.
focuses on sensor allocation with analyzing exponential stability. It is presumed that a continuous-time controller stabilizes the teleoperator system. The goal is to assign each position and velocity signals to sensors with different sampling rates and divide the state vector across sensors so that the multi-rate system remains exponentially stable. Different approaches for dividing the state vector across sensors can lead to stable or unstable systems. Thus, Krasovskii-based conditions satisfy exponential stability. Also, these conditions are employed to design observers on both sides of the teleoperator system. This problem will be transformed into a mixed-integer program with LMIs. The results of this paper are useful in designing sensor networks for teleoperation systems.
The remainder of this study is organized as follows. In Section 2, preliminaries and sampled-data model of the teleoperator system with multi-rate sampling approach are presented. The exponential stability analysis for sensor partitioning and observer design are provided in Sections 3 and 4, whereas in Section 5, these results are validated by numerical simulations. In Sections 6 and 7, the results are reported and concluding remarks stated.

Proposed Sampled-Data Model with Multi-Rate Architecture for Teleoperation System
In a bilateral teleoperation system, when a human manipulates the master robot, the generated motion signals are relayed to the slave robot via a wireless network, and the slave robot subsequently performs the same action. The general schematic of the teleoperation system is depicted in Figure 1.

Continuous-Time Model
In the proposed model for robotic manipulators as the master and slave subsystems, position and velocity signals have been considered as states of the system with this limitation that there is no gravitational force and friction. Also, it is presumed that the force from the human operator applied to the master robot can be modeled as a passive mass-spring-damper system. The exogenous force input exerted by the operator is presumed to be zero.
In the absence of frictions and external disturbances and assuming two degrees of freedom robots as master and slave subsystems, the linearized dynamics can be derived as:

Continuous-Time Model
In the proposed model for robotic manipulators as the master and slave subsystems, position and velocity signals have been considered as states of the system with this limitation that there is no gravitational force and friction. Also, it is presumed that the force from the human operator applied to the master robot can be modeled as a passive massspring-damper system. The exogenous force input exerted by the operator is presumed to be zero.
In the absence of frictions and external disturbances and assuming two degrees of freedom robots as master and slave subsystems, the linearized dynamics can be derived as: where subscripts m and s are notifications for master and slave robots. M, D and, S are used for mass, damping, and stiffness of robots, respectively. Position and velocity signals are mentioned by q and . q. F h and F e are imposed forces to the master and slave robots by operator and environment, respectively. F m and F s denote control signals.
The slave robot can be in a free motion or in contact with the environment. Let define the constant parameter κ = 1 for contact case and κ = 0 for free motion. The state-space equation can be stated as: where and where subscripts e and h are used for environment, and human operator, respectively.

Sampled-Data Model
Assuming the constant sampling period t s for all transmitted signals, the sampled-data counterparts of the aforementioned teleoperation system can be derived as: However, due to the reasons provided in the introduction section, a multi-rate approach is utilized in this study. Without losing the generality of the problem, the sampling rates are assumed to be the same for position signals and differ with velocity sensors in both master and slave robots sides. For this purpose, let us define the state vector as: and assume that both sides are sampled by varying sampling rates, represented by T 1 χ and T 2 χ .(χ = m for the master robot and χ = s for the slave robot). Thus, the sampling periods can be named by t 1 χk and t 2 χk . Also, for zero-order holds (ZOHs), which are indicated by Z κ , the refresh rates are donated by z χk . The schematic architecture for the proposed sampled-data bilateral system with multi-rate sampling, and 2 DOF robotic manipulators in slave and master robots sides is depicted in Figure 2. These assumptions mention that sampling and updating operations cannot occur twice at the same period. The time that has passed since the last sampling point is called the elapsed time and calculated as: The same definitions can be considered for actuators:

Assumption 2.
For actuators ε z > 0 exists which satisfies ε z < z 1 χ k+1 − z 1 χ k . These assumptions mention that sampling and updating operations cannot occur twice at the same period. The time that has passed since the last sampling point is called the elapsed time and calculated as: The same definitions can be considered for actuators: This last shows the elapsed time for ZOHs. The maximum sampling and updating periods that satisfy the stability are defined by: The duration of data transmission time from sensors to ZOHs can be calculated as: Remark 1. Unlike η χ (t) and η z (t), the function η χ sz (t) can be assigned by non-zero amounts at its discontinuities. Thus: (11) where σ χz indicates the upper bound for data transmission time from sensors to ZOHs.
In the literature, proportional-derivative (PD) continuous-time controllers are proposed to provide a stable and transparent teleoperation system. The sampled-data counterpart of these controllers are: where K p and K v are constant positive elements.
Applying the input delay approach [24], control laws are rewritten as: Remark 2. Due to the symmetric dynamics of the master and slave robots, controller gains are assumed the same in the bilateral structure.

Sensor Allocation with Exponential Stability Analysis
This section provides a basis for sensor partitioning to guarantee the exponential stability of the sampled-data system. According to the previous section, it is assumed that four sensors with various sampling rates are available, thus n s = 4.
Data-drop out is modelled via switches in Figure 2.; let nominate sensing blocks by S i , i ∈ IS where IS = {1, . . . , n s }. We define P = {IS i |i ∈ IS } as a partition of sensors with the following properties: Each sensor in S i , i ∈ IS is dedicated to sample one of the transmitted signals on both sides of the system. To parameterize the partition matrix, we define Y ∈ {0, 1} n x ×n s ; its rows correspond to the states, as well as the columns to the n s sensors, belonging to the set of IS i , i ∈ IS. The matrix's elements are defined in the Boolean set {0, 1}. When Y(i, j) = 1; it means that the j-th state is sampled by i-th sensor or, in other words, j ∈ IS i ; in the case of Y(i, j) = 0, we have j / ∈ IS i . The maximum admissible sampling periods are known and named by σ i where i ∈ IS. To imply this partitioning to the multi-rate controller the diagonal matrices C i ∈ {0, 1} n x ×n x , i ∈ IS is considered as: It is presumed that sensors have a time-driven structure, and ZOHs are event-driven and update themselves immediately after receiving new data.
The following proposed theorem provides LMIs to assign state vectors to each sensor and assures that the stability does not jeopardize. Theorem 1. Consider the sampled-data dynamic of the teleoperation system in (6) and (13). For positive scalars α and β and positive definite matrices E W o , N, Q, Q, and W i (i ∈ IS), there exists a sensor network parametrized by a matrix Y that guarantees the exponential stability with a rate greater than σ/2, if it satisfies the following LMIs: Proof. To verify the theorem above, we present an LKF function and show that the LMIs in (16) and (17) satisfy the stability preliminaries. Let us define the Lyapunov function as: where x(s))ds We demonstrate that LMIs in Theorem 1 are sufficient conditions for . V + αV ≤ 0 in the period among two successive sampling rates. The derivative of V 1 , V 2 and V 4 are straightforward and can be calculated as: . .
Considering t − η z = t n for an arbitrary vector k(t) . V 2 can be simplified as: .
A similar inequality can be assigned for the derivative of V 3 as follows: .
where W = ∑ i∈IS W i and x(t) = ∑ i∈IS C i x(t − η i (t)).
Let define the augmented parameters as: .
Based on the theorem above, the problem of state assignments for sensors that satisfies the stability can be restated as finding Y ∈ {0, 1} 4 * 4 subject to the (16) and (17).

Observer Design for the Multi-Rate Teleoperator System
In this section, observer design for the multi-rate teleoperator system has been addressed. For the calculated maximum admissible sampling periods, the contribution of this section is providing Krasovskii-based conditions to design linear observers on both sides of the teleoperator system. The estimation error of the observer should converge to the origin exponentially. The extension schematic of the sampled-data system with multi-rate observers has been depicted in Figure 3. Based on the theorem above, the problem of state assignments for sensors that satisfies the stability can be restated as finding (16) and (17).

Observer Design for the Multi-Rate Teleoperator System
In this section, observer design for the multi-rate teleoperator system has been addressed. For the calculated maximum admissible sampling periods, the contribution of this section is providing Krasovskii-based conditions to design linear observers on both sides of the teleoperator system. The estimation error of the observer should converge to the origin exponentially. The extension schematic of the sampled-data system with multi-rate observers has been depicted in Figure 3. In this scheme, l K is the gain of Luenberger observer. 11 22 1 , , A A B and, 2 B have been defined in (3) and (4), respectively. The estimation error of the observer is notified by Thus, the error dynamics can be defined as: In this scheme, K l is the gain of Luenberger observer. A 11 , A 22 , B 1 and, B 2 have been defined in (3) and (4), respectively. The estimation error of the observer is notified by e(t) = q(t) −q(t). Thus, the error dynamics can be defined as: where: Thus, Theorem 1 can be extended to calculate the observer gain via provided LMIs.

Theorem 2.
Given σ i x i ∈ (1, 2, 3, 4) and α > 0, there exists an observer gain K l = E −1 Y which the estimation error of the observer is exponentially stable for that if there exist positive definite matrices E and W 0 and matrices Y, Q and Q, with proper dimensions, satisfying the following LMIs: where Γ 1 , Γ 2 and, Γ 3 were defined in (19) and (20) and F = diag(exp(ασ i )I) i ∈ (1, 2, 3, 4) Proof. The proof is straightforward and similar to the proof of Theorem 1. It is enough to replace the Lyapunov function in (21) with the following terms and proceed with the same steps in Theorem 1.
The theorem above provides sufficient conditions for exponential convergence of estimation error to the origin for the multi-rate observer. It is assumed that the upper bound of sampling intervals for all sensors is available. The aim is to design observers whose gains are calculated via an optimization process stated in the forms of LMIs. Indeed, we first find Y subject to E > 0 and W 0 > 0, and then compute the observer gain as K l = E −1 Y.

Numerical Simulation
This section illustrates the circumstance of states assignment between sensors for the teleoperator system via numerical simulations. Two DOF robotic manipulators are considered as master and slave subsystems with the following dynamic parameters: The exerted force is modelled by mass-spring dynamics with 15 N/m for spring gain and 2 Ns/m for the damper gain; this force is depicted in Figure 4. It is assumed that the slave robot's end-effecter has touched the environment at 2 rad, as well as the controller gains are set to K p = 15 and K v = 0.3. The tuning process is based on the method proposed in [25]. In this approach, the integral of the absolute error is considered a cost function for the optimization algorithm to find optimal gain values for the PD controller to achieve the best tracking performance. According to Theorem 1, the sensor partitioning problem that preserves the stability can be transformed into a mixed-integer program with LMIs. This last is in contrast with our previous study, where the maximum admissible sampling periods and refresh rates had been calculated.  Two scenarios are defined to evaluate the sensor partitioning problem. In the first scenario, the master robot position signal for the first state and second state is presumed to be acquired via encoders at sampling periods up to 0.4 s and 0.6 s, and the rate of the velocity states on both sides will be measured by inertial measurement unit is 0.8 s. Thus, with master robot position signals and master-slave robots velocity signals assigned to the sensors, the sensor network matrix Y is defined as: where the rows of this matrix indicate position signals for two links. The first column is related to the encoder with a 0.4 s sampling rate and the second column with a 0.6 s sampling rate. We solve the LMIs in Theorem 1 to find which encoder should be utilized to sample the position signals in the slave side such that the stability is preserved. Solving (16) and (17)   Two scenarios are defined to evaluate the sensor partitioning problem. In the first scenario, the master robot position signal for the first state and second state is presumed to be acquired via encoders at sampling periods up to 0.4 s and 0.6 s, and the rate of the velocity states on both sides will be measured by inertial measurement unit is 0.8 s. Thus, with master robot position signals and master-slave robots velocity signals assigned to the sensors, the sensor network matrix Y is defined as: where the rows of this matrix indicate position signals for two links. The first column is related to the encoder with a 0.4 s sampling rate and the second column with a 0.6 s sampling rate. We solve the LMIs in Theorem 1 to find which encoder should be utilized to sample the position signals in the slave side such that the stability is preserved. Solving (16) and (17) with the MATLAB toolbox YALMIP [26] yields: y 1 = 1, y 2 = 0, y 3 = 0, and y 4 = 1.
This means that the first and second position states of the slave robot should be assigned to the first encoder and the second encoders, respectively, to provide exponential stability. To validate this concept, the position and velocity trajectories are depicted in Figures 5-8.   Figures 5 and 6 show the instability of the teleoperator system when the second encoder is utilized for the first position state of the slave robot side. In Figures 7 and 8, both teleoperator systems' stability and transparency are obtained considering the proper sensor partitioning on the slave robot side.   Figures 5 and 6 show the instability of the teleoperator system when the second encoder is utilized for the first position state of the slave robot side. In Figures 7 and 8, both teleoperator systems' stability and transparency are obtained considering the proper sensor partitioning on the slave robot side.  The root mean squared error (RMSE) is the square root of the quadratic mean of all errors. The calculated RMSE error related to Figure 8 is equal to 0.0003416.
To illustrate the efficiency of multi-rate sampling in comparison with the single-rate sampling, the maximum admissible sampling period based on the LMIs in Theorem 1 is calculated with the YALMIP software [26] and demonstrated in the following Table 1.    The root mean squared error (RMSE) is the square root of the quadratic mean of all errors. The calculated RMSE error related to Figure 8 is equal to 0.0003416.

Single-Rate Sampling
To illustrate the efficiency of multi-rate sampling in comparison with the single-rate sampling, the maximum admissible sampling period based on the LMIs in Theorem 1 is calculated with the YALMIP software [26] and demonstrated in the following Table 1.    Figures 5 and 6 show the instability of the teleoperator system when the second encoder is utilized for the first position state of the slave robot side. In Figures 7 and 8, both teleoperator systems' stability and transparency are obtained considering the proper sensor partitioning on the slave robot side.

Single-Rate Sampling
The root mean squared error (RMSE) is the square root of the quadratic mean of all errors. The calculated RMSE error related to Figure 8 is equal to 0.0003416.
To illustrate the efficiency of multi-rate sampling in comparison with the single-rate sampling, the maximum admissible sampling period based on the LMIs in Theorem 1 is calculated with the YALMIP software [26] and demonstrated in the following Table 1. It is deduced that by applying multi-rate samplers, there will be fewer conservative conditions for sampling instants or update rates to guarantee the system's stability. Specifically, while the maximum admissible sampling period for position sensor to preserve the stability is 150 ms in the single-rate method, multi-rate architecture can increase this upper bound to 352 ms.
The second scenario assumes that the transmitted signals are assigned to sensors with different and nonuniform sampling rates. These sensors can be categorized as a fast-rate sensor with a maximum sampling interval of 0.2 s, a medium-rate sensor with a maximum sampling interval of 0.6 s and a slow-rate sensor with a maximum sampling interval of 1 s. The aim is to assign these sensors to the transmitted signals on both sides, preserving stability. In this scenario, the sensor network matrix should have eight rows for four position and four-velocity signals and three columns for fast, medium and slow rate sensors. Thus, the number of Boolean elements stored in Y is 24. Solving LMIs in Theorem 1, the sensor network matrix Y is defined as: It is deduced that by applying multi-rate samplers, there will be fewer conservative conditions for sampling instants or update rates to guarantee the system's stability. Specifically, while the maximum admissible sampling period for position sensor to preserve the stability is 150 ms in the single-rate method, multi-rate architecture can increase this upper bound to 352 ms.
The second scenario assumes that the transmitted signals are assigned to sensors with different and nonuniform sampling rates. These sensors can be categorized as a fast-rate sensor with a maximum sampling interval of 0.2 s, a medium-rate sensor with a maximum sampling interval of 0.6 s and a slow-rate sensor with a maximum sampling interval of 1 s. The aim is to assign these sensors to the transmitted signals on both sides, preserving stability. In this scenario, the sensor network matrix should have eight rows for four position and four-velocity signals and three columns for fast, medium and slow rate sensors. Thus, the number of Boolean elements stored in Y is 24. Solving LMIs in Theorem 1, the sensor network matrix Y is defined as: In the proposed sensor allocation scheme, position signals are sampled with the fast-rate samplers according to their significant roles in stabilizing the system. Position and velocity trajectories for this case are depicted in Figures 9 and 10. The bilateral system becomes unstable if slow-rate sensors are used to sample the position signals.
The position estimation errors for the teleoperator system are depicted in Figure 11, which confirms that the observer gain in (39) preserves the exponential stability of the estimation error.  The position estimation errors for the teleoperator system are depicted in Figure 11, which confirms that the observer gain in (39) preserves the exponential stability of the estimation error.
The position estimation errors for the teleoperator system are depicted in Figure 11, which confirms that the observer gain in (39) preserves the exponential stability of the estimation error. Figure 11. Position estimation error trajectory for observer gain calculated in (39). Figure 11. Position estimation error trajectory for observer gain calculated in (39).

Discussion
In this paper, the sensor allocation approach is designed for sampled-data bilateral teleoperation systems. Based on Theorem 1 and its proof, this problem is converted into a mixed-integer program with LMIs that, if satisfied, the teleoperation system's stability will be preserved. The solution is stated as a sensor network matrix, that its rows represent the transmitted states and its columns show the number of sensors with varying sampling rates. The elements of this matrix are defined in the Boolean set {0, 1}, where 1 means that the state is sampled by the related sensor.
In contrast to in-depth studies on continuous-time teleoperation models in the literature, there is less research on sampled-data teleoperation systems. However, as deduced from the results and simulations in Section 5, sampling has a significant impact on the stability and neglecting it can undermine both transparency and stability of the teleoperator systems. Compared to similar studies that take the sampling analysis into account and provide a sampled-data model for teleoperator systems such as [19,27,28], the novelty of the proposed method relies on designing a sensor allocation basis for the multi-rate structure of bilateral teleoperator systems. This proposed method covers the significant deficiency in previous literature of discrete-time teleoperations systems, where all transmitted signals are sampled with the same sampling rates. Uniform sampling is not sensible for teleoperation systems since, besides the possibility of utilizing various sensors with varying sampling rates on both sides, master and slave robots may have different hardware restrictions regarding the existence distance. Therefore, a multi-rate approach has been presented in this study to collect data through numerous sensors with various sampling rates, as an extension to our previous investigations, where the maximum sampling interval of sensors has been calculated to guarantee the stability of the teleoperator system. Indeed, in [29], we introduced a stability framework for single-rate uniform and nonuniform sampling intervals for teleoperator systems. Also, in [22], the transmission rates' upper bounds were calculated to guarantee the teleoperator system stability. However, the sensor allocation problem among different possible scenarios was an open question covered in the proposed research work. For this purpose, an additional term is added to Lyapunov-Krasovskii function in (21) based on the diagonal matrix C to consider the multi-rate partitioning of sensors. Thus, the stability concept is stated in LMIs, whose solution yields the sensor network matrix. It is concluded that stability can be guaranteed with higher admissible sampling intervals and less determinative conditions. Compared to other LMI schemes, such as input delay in [24], the LMIs presented in Theorem 1 impose less determinative conditions on the maximum permissible sampling interval regarding the defined Lyapunov-Krasovskii function.
Finally, the contributions of this study can be summarized as follows: 1. Designing a sensor allocation structure for bilateral teleoperator systems.

2.
Providing exponential stability analysis for the proposed multi-rate sampling approach to collect data through numerous sensors with various sampling rates in teleoperation systems.

3.
Designing linear observers for master and slave robots sides of sampled-data teleoperation systems relied on Krasovskii-based conditions and maximum admissible sampling periods.
These evaluations can aid in the development of guidelines for more transparent and stable teleoperator systems.

Conclusions
Unlike continuous-time teleoperation systems, there are far fewer studies for sampleddata counterparts. A significant deficiency in previous literature on discrete-time teleoperations systems is that all transmitted signals are sampled with the same sampling rate. This study has investigated the effect of sampling rates of samplers and update rates of ZOHs on the teleoperator system's stability. Variable sampling intervals were considered for transmitted signals due to the isolation of the master and slave robots and inner time delays and data dropout in the communication channel. Also, in this paper, the sensor allocation issue was addressed for preserving the exponential stability of the sampled-data teleoperator systems using multi-rate sampling. This structure leads to less conservative conditions in the sense of stability and maximum admissible sampling periods. Sufficient conditions based on the Krasovskii functions were utilized to partition the transmitted signals according to the provided LMIs. Besides, these conditions were employed to design observers for determining the maximum admissible sampling periods for each sensor. As future works, these results can be extended to four-channel teleoperator systems with nonlinear dynamics.