An Off-Body Narrowband and Ultra-Wide Band Channel Model for Body Area Networks in a Ferryboat Environment

: In the article an off-body narrowband and ultra-wide band channel model for body area networks in a ferryboat environment is described. Considering the limited number of publications there is a need to develop an off-body channel model, which will facilitate the design of radio links, both from the multimedia services provider and the security point of view, for body area networks in this atypical environment. A mobile heterogeneous measurement stand, using radio distance measurements, which consists of three types of devices: miniaturized mobile nodes, stationary reference nodes, and a data acquisition server, was developed. A detailed analysis of both radio channels’ parameters was carried out. An analysis of system loss for off-body communication, including mean system loss, large-scale fading (corresponding to body shadowing), and small-scale fading (associated with the multipath phenomenon), both for 868 MHz narrowband and for 6489 MHz ultra-wide band channels, was performed. A statistical analysis of the obtained system loss model parameters was also carried out; good ﬁt to the empirical data is observed.


Introduction
The development of mobile radio communication systems is related to constantly growing expectations and the needs of the users. The first generation of mobile telephony (1G) was launched in 1979 and was based on analog transmission. At the beginning of the 1990s, there was a significant breakthrough in the standardization of the mobile telephony. The implementation of the second-generation (2G) Global System for Mobile Communications (GSM) network offered not only voice services, but also transfer of data, text, and multimedia messages, indicated a milestone in the technological progress. Then, at the turn of the 21st century, the Universal Mobile Telecommunications System (UMTS) system was introduced, which is the third-generation (3G) mobile telephony standard. Almost a decade later, miniaturization of electronic devices and technological progress in terms of increase in computational power, as well as growing expectations for various services offered by service providers, initiated the development of the Long-Term Evolution (LTE) and LTE Advanced systems, which are recognized as the fourth-generation system (4G). A fast growth of data transmission and spectral efficiency, as well as the reduced delays of point-to-point transmission in the whole network, are some of the main goals achieved during the implementation of this system [1 -3].
The trend of introducing the next generations of systems every decade can be seen, so it is expected to launch the first fifth-generation (5G) systems around 2020. It is supposed to be a significant technological leap in relation to existing systems, not only due to a significant increase in the data rate

Related Work
The rapid technological progress in development of small size wearable devices dedicated for WBANs does not eliminate the key problem of ensuring a stable radio link. In these types of networks one has to deal with high system losses, i.e., introduced by the human body and multipath propagation in the environment, especially for applications with dynamic users. The indoor propagation issues, concerning WBANs, are widely analyzed by researchers. The majority of currently designed WBANs are used in indoor environments in which radio wave propagation is significantly disturbed compared to open space, mainly due to the stronger multipath effect. In the current state of the art, many indoor off-body channel models were proposed, but still there is a need to perform additional measurements and analytical studies for other types of environments. Therefore, it was decided to carry out measurements in a passenger ferryboat, which is an atypical and very harsh environment in terms of a radio wave propagation. Due to its metal structure, it cannot be treated as a typical indoor environment.
The obtained measurements allowed for development of an off-body channel model, which will facilitate the design of WBANs. It is very important, both from the perspective of multimedia service providers, as well as the security point of view. The practical example of WBAN operating in the ferryboat is a system for remote tracking and monitoring of people, developed at the Department of Radio Communication Systems and Networks, Gdansk University of Technology [11]. This system was developed for the Polish Maritime Regional Unit of Border Guard, and its goal was to increase the security of the officers on duty. It consists of off-body personal identification modules, which are mounted directly on the human body, while reference nodes are attached to the ferryboat's construction elements. The telemetry data about the officers are sent to the operating center.
Several previous research studies in this environment analyzed ultra-wide band (UWB) channels operating in the 6-9 GHz band and narrowband (NB) channels operating in the 2.4 GHz band. In [12], measurements of NB 2.4 GHz radio channel in a similar ferryboat environment were carried out.
Despite the environments' similarities, the main goal of the analysis was to evaluate a simple model for stationary surveillance network. In [13], measurements were carried out for an on-body UWB channel. The proposed model consists of a normally distributed random variable, which represents the received signal variations due to fading. The work in [14] is focused on obtaining a statistical model for the UWB propagation channel in a ferryboat passenger cabin. Moreover, the measurements were focused on determining and further modeling of the channel impulse response, not the system loss. In [15], the developed path loss model was evaluated for UWB off-body communication, where the fading is modeled by a Gaussian variable, but the model was not investigated in the ferryboat environment. The most complex model for WBANs in a ferryboat is proposed in [16]. It accounts for the mean system loss (MSL), and attenuation due to fading. The measurement scenarios were dynamic, but only the NB channel at 2.45 GHz is analyzed.
Considering the limited number of publications related to channel modeling in the passenger ferryboat and the growing interest in the use of WBANs in such an environment, the conducted research and proposed novel system loss model is unique. It should also be noted that the authors have used custom measurement equipment, which uses radio distance measurements (RDMs) to dynamically estimate the distance between the moving person and the static reference node (RN). This approach gives accurate data despite the different user motion velocities during the measurements.

Measurement Equipment
Conducting measurements of NB and UWB in the real environment required the development of heterogeneous nodes adapted to work around the human body. The proposed measurement equipment consists of three types of devices: mobile node (MN), stationary R,N and data acquisition server (DAS).
The MN and the RN have similar hardware design. The only difference, from the electrical diagram point of view, is the existence of the RS-232 wired interface in RN. In Figure 1, a block diagram of the MN is presented. carried out. Despite the environments' similarities, the main goal of the analysis was to evaluate a simple model for stationary surveillance network. In [13], measurements were carried out for an onbody UWB channel. The proposed model consists of a normally distributed random variable, which represents the received signal variations due to fading. The work in [14] is focused on obtaining a statistical model for the UWB propagation channel in a ferryboat passenger cabin. Moreover, the measurements were focused on determining and further modeling of the channel impulse response, not the system loss. In [15], the developed path loss model was evaluated for UWB off-body communication, where the fading is modeled by a Gaussian variable, but the model was not investigated in the ferryboat environment. The most complex model for WBANs in a ferryboat is proposed in [16]. It accounts for the mean system loss (MSL), and attenuation due to fading. The measurement scenarios were dynamic, but only the NB channel at 2.45 GHz is analyzed.
Considering the limited number of publications related to channel modeling in the passenger ferryboat and the growing interest in the use of WBANs in such an environment, the conducted research and proposed novel system loss model is unique. It should also be noted that the authors have used custom measurement equipment, which uses radio distance measurements (RDMs) to dynamically estimate the distance between the moving person and the static reference node (RN). This approach gives accurate data despite the different user motion velocities during the measurements.

Measurement Equipment
Conducting measurements of NB and UWB in the real environment required the development of heterogeneous nodes adapted to work around the human body. The proposed measurement equipment consists of three types of devices: mobile node (MN), stationary R,N and data acquisition server (DAS).
The MN and the RN have similar hardware design. The only difference, from the electrical diagram point of view, is the existence of the RS-232 wired interface in RN. In Figure 1, a block diagram of the MN is presented. Each MN and RN is heterogeneous, in the sense that it has two simultaneously working, time synchronized radio interfaces:  NB radio module CC1120 from Texas Instruments, working in the ultra-high frequency (UHF) 868 MHz band [17], Each MN and RN is heterogeneous, in the sense that it has two simultaneously working, time synchronized radio interfaces: • NB radio module CC1120 from Texas Instruments, working in the ultra-high frequency (UHF) 868 MHz band [17], • UWB radio module DWM1000 from DecaWave, working in the super-high frequency (SHF) 6489 MHz band [18].
Both radio interfaces work in accordance with the IEEE 802.15.4 standard, including the operating frequencies in the Wireless Medical Telemetry Service (WMTS) and UWB bands for Europe, which is indicated as possible for use in WBANs [19].
The selected parameters of the 868 MHz module should be mentioned: center frequency 868.3 MHz, 5 dBm of output power, 2-GFSK (Gaussian frequency shift keying) modulation with frequency deviation of 50 kHz, 50 kbps bitrate, 67 B packet data length with data whitening included. The antenna for the NB interface is Yageo's compact ceramic antenna ANT1204F007R0870A [20], with 1.67 dBi gain, resonance frequency equal to 870 MHz, and reflection loss −28 dB. The choice of antenna was motivated by its physical dimension and the possibility of using it in miniaturized devices.
The selected UWB parameters for the measurements are as follows: 6489 MHz center frequency, −41.3 dBm/MHz transmitted power spectral density and 499.2 MHz of the occupied bandwidth. The ACS5200HFAUWB ceramic omnidirectional antenna element for the UWB interface is integrated in the DW1000 module [21].
The DAS includes a GETAC F110 industrial tablet [22], together with a dedicated DAS server application. The custom developed dedicated software is a multi-threaded application, which allows for continuous monitoring of the obtained telemetry data and the status of each network node.

Measurement Scenarios
Propagation measurements were conducted on the passenger ferryboat, named Motor Ferry WAWEL (MF WAWEL), during mooring at the port quay. During the measurements, only the crew was present on the ferryboat. Structural elements of the environment (i.e., walls, ceiling, floor, doors, and handrails) are made of steel of various thicknesses.
For all considered scenarios, the stationary RN was attached to the corridor wall at a height of 1.2 m, on the seventh passenger deck. A single MN was placed on the user moving along the axis of a narrow corridor. The separation between human skin and radiating elements was about 1 cm. The measurements were carried out by one user: a man with a height of 1.72 m, a weight of 60 kg and a body mass index (BMI) of 20.3.
In Figure 2, a schematic diagram of the measurement equipment is shown. It illustrates the places where MN were mounted. From many possible mounting locations, three were selected as representative for WBANs applications, including smart watches, heart rate monitors, multimedia glasses, etc. [7,23,24]. These are: It is known that static scenarios are characterized by a smaller spread of propagation attenuation around the mean value and a smaller error in distance estimation using RDMs due to the absence of the small-scale fading component associated with the user's body motion, or lack of changes of its position. For this reason, it was decided to analyze only dynamic scenarios, which are much more complex in terms of radio wave propagation [25,27,28].
The measurements were carried out for two dynamic scenarios (SCNs): approaching to (APR), and departing from (DEP), the RN. These two scenarios correspond to a typical human behavior in indoor environments, and they are widely considered in WBAN research. They allow for the analysis of the measurement data independently of the side of the body on which the node is located, i.e., regardless of the antenna's placement on the left or right wrist. For the DEP scenario, the user started Typically, in other works [23][24][25], the distance between measurement points was calculated on the basis of the average speed of the user. This required continuous supervision during the measurements and manual determination of the beginning and the end of the measurement data recording process. During the presented research, RDM were performed by using the Symmetrical Double-Sided Two-Way Ranging (SDS-TWR) measurement method [26]. This allowed for the full automation of the system loss measurements. Time synchronization of both radio interfaces, using a dedicated signal line, allowed for obtaining the information on the system loss value and the corresponding distance between the receiver (Tx) and the receiver (Rx) at a given moment. Thus, all measurement data were recorded and analyzed regardless of the speed of the user in different phases of the scenario.
It is known that static scenarios are characterized by a smaller spread of propagation attenuation around the mean value and a smaller error in distance estimation using RDMs due to the absence of the small-scale fading component associated with the user's body motion, or lack of changes of its position. For this reason, it was decided to analyze only dynamic scenarios, which are much more complex in terms of radio wave propagation [25,27,28].
The measurements were carried out for two dynamic scenarios (SCNs): approaching to (APR), and departing from (DEP), the RN. These two scenarios correspond to a typical human behavior in indoor environments, and they are widely considered in WBAN research. They allow for the analysis of the measurement data independently of the side of the body on which the node is located, i.e., regardless of the antenna's placement on the left or right wrist. For the DEP scenario, the user started the motion at the distance of 2 m from the RN with the line-of-sight (LOS) conditions and walked towards the end point through the L-shaped corridor. After 6.5 m along the walking path there is a change between the zone with LOS conditions and the zone behind the corridor's corner with non-LOS (NLOS) conditions. The walk ended after travelling a distance of 14 m. In Figure 3, a floor plan of the investigated environment, which is a part of the passenger ferryboat's seventh deck, is shown.  The LOS and NLOS conditions were determined by the presence of the metallic walls between MN and RN. The APR scenario is reversed with respect to the DEP one, i.e., from the aspect of the walking direction and the starting and ending points of the route. For each scenario and MN placement, the measurements were repeated 20 times. During the measurements, there were no other people in the corridor nor inside the passenger cabins.

Analysis of the Results
During the measurements, 87,252 system loss values were collected with 40 ms interval between the consecutive measurements, for NB and UWB, for both motion scenarios and all three MN placements. It should be mentioned that the preliminary results, including the mean and the standard deviation system loss values are presented in [29,30], hence, they are not analyzed in detail in this The LOS and NLOS conditions were determined by the presence of the metallic walls between MN and RN. The APR scenario is reversed with respect to the DEP one, i.e., from the aspect of the walking direction and the starting and ending points of the route. For each scenario and MN placement, the measurements were repeated 20 times. During the measurements, there were no other people in the corridor nor inside the passenger cabins.

Analysis of the Results
During the measurements, 87,252 system loss values were collected with 40 ms interval between the consecutive measurements, for NB and UWB, for both motion scenarios and all three MN placements. It should be mentioned that the preliminary results, including the mean and the standard deviation system loss values are presented in [29,30], hence, they are not analyzed in detail in this article. The corresponding RDM were collected simultaneously. It should be pointed out that, in this article, the propagation attenuation is understood as the system loss [31], i.e., the difference in power delivered to the transmitting antenna and the power available at the terminals of the receiving one.
The preliminary analysis of the empirical system loss values for both frequency bands showed that the propagation model should include the estimation of the system loss separately for LOS and NLOS conditions, and separately for each MN placement. The results of the distance estimation by RDM indicated the significant differences in the radio wave propagation across the environment. In this environment the waves propagate along the corridor due to the metal structure and the waveguide effect, unlike in the typical indoor environments, where the radio waves can propagate along the direct line (even through the walls) between MN and RN [29,30].
The system loss (L SL(dB) ) is modeled by the general following expression: where: ∆L BS is the body shadowing component represented by a random variable, and • ∆L MF is the multipath component represented by a random variable.

Mean System Loss Component
The MSL component represents the decrease in the received signal power with distance due to absorption and dissipation of energy in the propagation medium, which is represented as follows: where: is MSL value at the reference distance d 0 (d 0 = 2 m for LOS and d 0 = 8 m for NLOS), • n is system loss exponent (SLE), and • d is the distance between Tx and Rx.
The MSL component was obtained by fitting the logarithmic distance function to the empirical data. Linear regression with the least squares method was used.
To analyze the accuracy of linear model, two well know measures were used, namely standard error of estimate (SEE), representing the standard deviation of empirical data around the regression line, and mean error (ME) [32].
The obtained values of the SLE (n) and the MSL (L S (d 0 )) at reference distance d 0 , both for 868 MHz and 6489 MHz bands, are presented in Table 1.
There are visible differences between the MSL parameters for various scenarios, placements of MN, and for different radio interfaces (NB and UWB). For each scenario, a significant differences in both parameters for the LOS and NLOS conditions were observed. The ME error was 0 dB, which means that the MSL model parameters are correctly estimated from the measurements.
In the 868 MHz NB channel the multipath propagation phenomenon and the waveguide effect in the corridor are clearly visible. The values of n for LOS conditions are in the range of [0.16; 0.34], being below the free-space value of 2 [33]. For NLOS, excluding the HE R scenario, the parameter n is in the range of [1.40; 3.80]. These differences are related to the high power level of the multipath components reaching the antenna compared to the power of the direct component [34]. The noticeably low value for the HE R scenario (0.16) can be explained by the unique environmental conditions, in which there is a strong reflected component from the opposite side of the corridor wall. Moreover, the polarizations of both antennas were matched and their mutual orientation changed to a small extent. As expected, the MSL at the reference distance is greater under NLOS conditions for most cases, but these differences are not significant. The SEE does not exceed the value of 3 dB for all the scenarios. Similar to the NB channel, there are also differences between the MSL model parameters for the UWB channel. The results show a significant difference in between the values of n for the LOS and NLOS conditions. Similar conclusions based on empirical data are presented in the available literature [14,35,36]. For the LOS conditions, the n values are in the range of [0.13; 0.85]. As expected, these values are below the free-space value, which is caused by the waveguide effect. It was also noticed that, similar to the NB interface, the highest value of n (0.85) was obtained for the DEP scenario and the TO F placement. However, the obtained values of n for the NLOS conditions are in the range of [1. 25; 3.46]. The two highest values of n were obtained for the MN mounted on the head (HE R ). Almost all scenarios are characterized by a higher MSL at the reference distance for the NLOS conditions relative to the LOS ones; the difference never exceeds 5.1 dB. The obtained SEE does not exceed the value of 3.1 dB.

Body Shadowing Component
According to Equation (1), one of the key elements of the system loss is represented by the component associated with the body shadowing ∆L BS . Due to the random nature of human motion, ∆L BS is also random, thus, it is represented by a random variable.
Obtaining the ∆L BS component is possible by a low-pass filtration of the empirical data. The method of this filtration depends on the analyzed measurement scenario. In static scenarios, the empirical data are averaged over time, yielding the composite component L SL (d), i.e.: where N t is the number of received signal samples obtained within a certain time window for the analyzed distance.
In dynamic scenarios, the filtration of the empirical data cannot be performed directly by using Equation (3) and requires spatial filtering due to the variable distance between the measurement nodes. Such filtration was realized with a linear moving average filter, which can be described by the following expression: where M is the length of the averaging window, which includes the number of samples corresponding to ten wavelengths (10λ). This is the selected value of the spatial filter resolution, determined from measurements [27]. It should be mentioned that the user's motion velocity was varying during the measurements, thus, the length of the filter was adjusted. These factors determined the development of the adaptive filter, whose length is selected on the basis of the estimated parameters of the user's motion during the measurements. The body shadowing component is obtained by subtracting the MSL (Equation (2)) from the composite component (4), i.e.: In order to model the large-scale fading component, the best fitting probability distribution was chosen. Due to the physics of the underlying phenomena, the losses affecting the received signal add up in the logarithmic domain. This yields a lognormal distribution of the large-scale fading component, as confirmed by the numerous studies available in the literature [37].
The lognormal probability distribution was fitted to the empirical data in linear domain (∆L BS ). This was done in MATLAB, using the fitdist function, which is the implementation of the maximum likelihood estimation.
In Figure 4, some examples of the lognormal probability density functions (PDFs) fitted to the empirical fading distribution for one of the measurement scenarios in NB channel were presented. The goodness-of-fit (GoF) metrics, i.e., the correlation coefficient (r) and the chi-square (χ 2 ) test were also evaluated. The correlation coefficient should be higher than the adopted acceptability threshold r TH = 0.90. The threshold value for the chi-square test, i.e., for the significance level α = 0.05, is χ 2 th = 28.87 for Rayleigh distribution, and χ 2 th = 27.58 for the others.  In order to visualize the obtained results in the best way, it was decided to present the distribution parameters, i.e., the average value µ and standard deviation σ in the logarithmic scale. In Table 2, the parameters of the fitted lognormal distribution are shown. The average value oscillates around 0 dB, which is in line with the results presented in the literature [28,29]. It also means that the adaptive low-pass filter (Equation (4)) was correctly implemented.  In order to visualize the obtained results in the best way, it was decided to present the distribution parameters, i.e., the average value µ and standard deviation σ in the logarithmic scale. In Table 2, the parameters of the fitted lognormal distribution are shown. The average value oscillates around 0 dB, which is in line with the results presented in the literature [28,29]. It also means that the adaptive low-pass filter (Equation (4)) was correctly implemented. For the NB channel the standard deviation is in the range of [1.2; 2.7] for almost all scenarios; the lower limit was obtained for LOS and the upper one for NLOS conditions. This range is wider than it is presented in the literature [28], i.e., [1.49; 1.52]. This is caused, among others factors, by the waveguide effect, which yields to a greater spread of the results for various MN placements.
All fitted distributions have passed the chi-square test. The value of the correlation coefficient in most cases is greater than the acceptance threshold. Only for the APR scenario and the MN placement on the chest (TO F ) can a significantly low (0.37) r parameter be observed. Therefore, this scenario was not taken into account during the analysis of the standard deviation.
For the UWB channel the observed standard deviation values are in the range of [1.3; 2.8], similarly as for the NB channel; the lower values were obtained for the LOS, while the higher ones were obtained for the NLOS conditions. In the majority of the analyzed scenarios, the standard deviation is higher for the UWB channel. Although both NB and UWB bands are significantly different in terms of frequency, occupied bandwidth, resistance to multipath effects, and frequency selectivity, it can be noticed that the values of the standard deviations are proportional to each other only if the differences between particular scenarios are considered.
The average value of σ equals 1.7 dB for NB and 2 dB for UWB channel. The obtained results show that the lognormal distribution, in general, properly represents the shadow fading component regardless of the MN placement and the propagation conditions.

Multipath Component
The system loss model for off-body communication (Equation (1)) also takes the small-scale fading component ∆L MF into account. It represents the impact of the user's motion and the multipath effect. Due to the random nature of this phenomenon, ∆L MF also has a random character and is represented by a random variable. Extraction of ∆L MF is possible by subtracting the value of the composite component L SL (d) (Equation (4)) from instantaneous values of the system loss L SL , i.e.: The methodology of determining the empirical distribution of the small-scale fading component and fitting the theoretical PDFs is very similar to that used for the body shadowing component, ∆L BS , and is also performed in the linear domain (∆L MF ).
Based on the literature dealing with the small-scale fading distributions, both for NB and UWB channels, the most frequently used ones were selected for fitting, i.e., lognormal, Rice, Rayleigh, Weibull, and Nakagami-m. For WBAN channels, the lognormal, Weibull, and Nakagami-m distributions were indicated as those that most accurately represent the small-scale fading amplitude distribution. Statistically, Rice and Rayleigh distributions are less often used for these channels [38].
In Figure 5, the exemplary results of the selected PDFs fitted to the empirical distribution of multipath fading is shown. Rice-and Rayleigh-fitted PDFs' distribution parameters were similar. This may be caused by the fact that, in a strong multipath environment like the investigated one, there is no dominant direct component, and the received signal is a sum of a large number of reflected components. Apparent asymmetry of the empirical distribution around the mean value makes the lognormal distribution the best fitting one. Apart from the subjective assessment, the GoF test has been performed using the correlation coefficient and the chi-square test. To present the results in a tabular form, a unified description of the P1 and P2 distributions' parameters was used, where:


for Rice distribution, P1 is the non-centrality parameter sR and P2 is the scale parameter σR,  for Nakagami-m distribution, P1 is the shape parameter m and P2 is the scale parameter Ω,  for Weibull distribution, P1 is the shape parameter αW and P2 is the scale parameter βW,  for Rayleigh distribution, P1 is a scale parameter σRAY, and  for lognormal distribution, P1 is the mean value µ and P2 is the standard deviation σ, both expressed in (dB). Tables 3-5 show the detailed list of fitted distributions' parameters for different MN placements and scenarios. The threshold value of the chi-square test is the same as for the body shadowing case.  To present the results in a tabular form, a unified description of the P 1 and P 2 distributions' parameters was used, where: • for Rice distribution, P 1 is the non-centrality parameter s R and P 2 is the scale parameter σ R , • for Nakagami-m distribution, P 1 is the shape parameter m and P 2 is the scale parameter Ω, • for Weibull distribution, P 1 is the shape parameter α W and P 2 is the scale parameter β W , • for Rayleigh distribution, P 1 is a scale parameter σ RAY , and • for lognormal distribution, P 1 is the mean value µ and P 2 is the standard deviation σ, both expressed in (dB). Tables 3-5 show the detailed list of fitted distributions' parameters for different MN placements and scenarios. The threshold value of the chi-square test is the same as for the body shadowing case.  For the majority of NB cases, the correlation coefficients satisfy the acceptance criteria. The best results were obtained for the lognormal distribution, although, in some cases, Nakagami-m also gives acceptable results. It may be observed that lognormal distribution fulfills the acceptability criteria for both GoF tests, regardless of the analyzed case. For the lognormal distribution, a significantly higher correlation coefficient was obtained compared to the other distributions. Thus, considering the data in the linear scale, the short-term fading is best described by the lognormal distributed random variable.
Similarly, all distributions fitted to the empirical UWB channel data are within the range of the acceptability of the chi-square test. It was also noticed that the value of the correlation coefficient in a few cases is greater than the selected acceptability threshold and typically oscillates close to this critical value. Based on the presented results, it can be assumed that the best results for the UWB channel were obtained for the lognormal and Nakagami-m distributions. However, it was decided to choose the lognormal distribution, which yields slightly better GoF test results. It is worth noting that, for the analyzed measurement scenarios, the Rice and Rayleigh distributions should not be used to describe the small-scale fading amplitude distribution, which confirms the conclusions from the research carried out so far [38].

Statistical Analysis of the System Loss Model
The developed propagation model, whose components are determined by Equations (2), (5), and (6), can be expressed as: where: is a random variable with lognormal distribution, which represents the power of the large-scale fading, due to the body shadowing effect, where µ B and σ B are the mean value and standard deviation, respectively; is a random variable with lognormal distribution, which represents the amplitude of the small-scale fading, caused by the multipath effect, where µ F and σ F are the mean value and standard deviation, respectively.
In order to analyze the accuracy of the developed propagation model, two statistical evaluation measures were used, i.e., the SEE and the ME [32]. In this case, the SEE is a measure of the model accuracy, and the ME determines the degree of dispersion of empirical data around the model. For the better fitting model, the SEE value is lower and the ME value is closer to 0 dB. It should be also noted that the model is considered to be sufficiently accurate if the SEE value does not exceed 8 dB [3,39].
Due to the two random factors included in Equation (7), it was decided to calculate both ME and SEE from 10,000 model simulations. In Table 6, the statistical analysis of the accuracy of the developed model for all scenarios is presented. Based on the obtained results, it was found that the developed model accurately describes the distribution of the system loss for the NB and UWB channels, operating at 868 MHz and 6489 MHz, respectively.
In Table 7, the general parameters of the channel model described by Equation (7) for LOS/NLOS conditions are presented. It is possible to use a generalized model in which the given ranges can be reduced to a single parameter value, by calculating the mean of its limit values. More detailed parameters of the model, which depend on the MN placements and scenarios, are given in Tables 1-5.

Conclusions
The growing demand for health monitoring devices, i.e., recording motion parameters or providing multimedia content, has contributed to the significant development of radio networks working within the human body or in its proximity. It applies directly to the evolution of WBANs, which is currently the major direction of modern 5G radio communication networks development, and (apart from medicine) are widely used in the entertainment, marketing, and security. To design such systems, it is very important to have an accurate radio channel model, especially when they are supposed to work in harsh environments.
In the article, the off-body NB and UWB system loss model for WBANs in a ferryboat environment is described. A mobile heterogeneous measurement equipment, which consists of three devices: miniaturized MN, stationary RN, and a DAS, was developed. The time-synchronized RDMs, implemented in the UWB interface, are a novel approach to system loss measurements and they enable one to automatize the measurement process and simplify the data analysis. Dynamic scenarios corresponding to a typical human behavior in an indoor environment were considered, where the lack of direct antenna visibility conditions, due to the presence of structural elements of the environment, were distinguished.
A detailed analysis of both radio channels' parameters was carried out. An analysis of the system loss distributions, SLE, MSL, large-scale fading corresponding to body shadowing and small-scale fading associated with the multipath phenomenon both for the NB and UWB channels was performed. The obtained SLE values are in the range of [0. 16; 3.8] for the NB channel, and [0.13; 3.46] for the UWB one. The lower values correspond to the LOS and the higher ones to the NLOS conditions. In addition, the MSL values at the reference distance for the NB channel, being in [64.7; 76.2], are much higher than in the UWB one, i.e., in [23.4; 30.7]. This is due to the transmission techniques, body influence, used antennas, and the impact of the environment.
The and [0.7; 1.7] for UWB channel, should be considered for radio link designing. A statistical analysis of the model parameters was also carried out and its good fit to the empirical data was demonstrated.
Future work will focus on determining the system loss model dependency on frequency. The behavior of channel impulse response, especially in the UWB channel, in accordance to the MN placement on the human body will be investigated.
Author Contributions: All authors contributed to developing the measurement equipment, designing and performing measurements, data analysis, scientific discussions, and writing the paper.