1. Introduction
Due to the complexity of the internal structure of the human body and the influence of the external background environment, when a microwave radiometer is applied to measure the temperature of human tissue, the temperature measurement accuracy is often difficult to meet the requirements [
1,
2]. Therefore, it is necessary to provide calibration for the human body temperature microwave radiometer through a precise blackbody calibration target to improve the temperature measurement accuracy. In addition, according to the relationship between the frequency of a microwave and its penetration depth in the human body [
3,
4], only by determining the appropriate operating frequency band can the microwave radiometer correctly receive the radiation brightness temperature of the human tissue to be measured. This paper specifically measures the temperature of human multi-layer skin tissue and selects the Ku-Band as the working frequency band of the blackbody calibration target [
5]. The design principle and design method of the precise blackbody calibration target in the Ku-Band are introduced below.
The blackbody calibration target is to convert the energy with a physical temperature of
into microwave energy with a brightness temperature of
at an emissivity close to 1 and radiate it out. To achieve the best performance of the blackbody calibration target, it is necessary to: (1) optimize the electromagnetic characteristics of the blackbody calibration target to make the emissivity the highest [
6]; and (2) analyze the physical temperature characteristics of the blackbody calibration target to ensure it has good temperature uniformity, stability and a small temperature gradient [
7,
8].
According to Kirchhoff’s radiation law [
9,
10,
11], the sum of the reflectivity and emissivity of a microwave blackbody calibration target is 1. Designing a microwave blackbody calibration target with high emissivity means designing a structure with a high electromagnetic absorption rate and a low electromagnetic reflectivity [
12,
13]. Therefore, the analysis of electromagnetic scattering characteristics can provide important information for the design and measurement of a high-emissivity blackbody calibration target. There have been some studies which focused on the analysis and calculation of blackbody scattering characteristics. In recent years, Finite Difference Time Domain (FDTD) is often used for establishing an electromagnetic wave scattering model of blackbody [
14,
15,
16]. Then, for an infinite cone array, the reflectivity can be solved based on the Radar Cross Section (RCS), while for a finite cone array, the reflectivity can be solved by integrating the full-space scattering power. Other studies have pointed out that the cone array is usually a periodic structure. The reflectivity of infinite periodic array structures can be evaluated more accurately using Floquet mode [
17,
18,
19,
20].
The temperature characteristics of the blackbody calibration target are independent of electromagnetic characteristics, and the most common analysis method is to use finite element software to conduct a thermal simulation analysis of the temperature field of the blackbody calibration target [
21,
22,
23]. Based on the analysis, the factors affecting the uniformity of the temperature field can be discovered. The shape, material and heating conditions of the cavity structure can then be designed to give the blackbody calibration target a good temperature uniformity.
Although the internal structure of the blackbody calibration target has a great influence on the temperature uniformity, the control of its temperature is also a key point. The cavity of a blackbody radiation target needs to obtain a uniform temperature distribution, which requires a higher standard for temperature control. At present, the temperature control system in the black body calibration source mostly uses the traditional PID control algorithm, its parameters often need to be adjusted by experienced experts for many tests and debugging, the process is troublesome and time-consuming and laborious. Some researchers have combined the traditional PID control with fuzzy control [
24,
25], which can give rise to the full advantages of fuzzy control in dynamic control with good rapidity and make the stability of the system meet the requirements through the fine control of PID. With the in-depth study of intelligent algorithms and the good self-organization ability of neural networks, the research of PID algorithms based on neural networks in the field of temperature control has gradually increased [
26,
27]. Compared with traditional PID control, this method has significant advantages related to its adaptability, robustness and control quality.
The main work of this paper includes: (1) Due to the periodicity of the blackbody structure, this paper first uses the Floquet mode to study the scattering characteristics of an infinite array [
28]. Then, according to the spatial symmetry, this paper proposes to use a two-dimensional scattering model to explore the two-dimensional pattern on the XOZ plane to further analyze the scattering characteristics of the finite array. The above methods greatly improve the computing efficiency and reduce the requirements for hardware devices. (2) The scattering properties and temperature uniformity of calibration targets are discussed separately in the existing literature. This paper notes that the shape and structure of the coated cone array will affect both its scattering characteristics and the temperature uniformity. Under the premise of satisfying both the emissivity and the temperature uniformity, a new design method is determined to find a set of design parameters to optimize the performance of the calibration target. (3) Several advanced PID algorithms are analyzed and compared, and an optimal algorithm is selected to be applied to the blackbody calibration target to obtain faster and more accurate temperature control. (4) We complete the processing and installation of the blackbody calibration target, test the emissivity and temperature uniformity of the equipment, and finally, introduce an uncertainty analysis.
The overall design process of the blackbody calibration target is shown in
Figure 1.
4. Research of Temperature-Control Algorithms
In this section, the blackbody calibration target is taken as the control object, and the temperature-control algorithms include PID, PSO-PID, fuzzy-PID and BP-PID.
The temperature-control system of the blackbody calibration target is approximately equivalent to the first-order inertial element plus the delay element. The general mathematical form of the transfer function is:
This paper divides the temperature change situations into two types: heating and cooling. The transfer functions in the two cases are set as Equation (18) and Equation (19):
For the BP-PID algorithm, to optimize the parameters of the PID algorithm using the BP neural network, the number of hidden layers of the BP network, the learning rate of the weights, and the momentum factor can be adjusted. The learning rate has the greatest impact on the temperature-control accuracy, and the momentum factor will affect whether the temperature-control process produces oscillations. To prevent the output result from failing to converge, it is necessary to control the size of PID control increment Uk obtained after each training of the BP neural network and make the temperature control result converge to the set temperature value. In addition, for the transfer function of cooling, it is necessary to set the control rate of the PID parameter as a negative number to achieve the cooling effect. The experimental results show that when the hidden layer is 8, the learning rate is 0.0005, and the momentum factor is 0.3. Thus, an ideal heating effect can be obtained. However, it cannot be guaranteed that the cooling process can achieve an ideal effect every time. Since the weight matrix of the BP neural network is randomly assigned each time, it is easy for the BP network to fall into a local minimum.
The PSO-PID algorithm is based on the traditional PID algorithm in addition to the objective function of the particle swarm optimization algorithm and the objective loss function with time weight, and sets the three parameters P, I, D as the three position variables Kp, Ki and Kd. The number of iterations of particle swarm optimization are set to 10, the size of particle swarm optimization to 100 and the inertia factor to 0.6. Then, the PID optimization range is set. With the increase in the number of iterations, the particle swarm optimization algorithm dynamically changes the PID parameters so that the fitness of the objective function gradually decreases, reaching the effect of convergence to the set temperature value.
For the Fuzzy-PID algorithm, the number of output discrete samples in the Fuzzy controller is adjusted. The larger the value is, the faster the PID converges, although there will be a marginal decline effect. Compared with traditional PID regulation, Fuzzy-PID regulation is more intelligent, gradually refined and has a faster response speed and less overshoot.
Through the above analysis of the parameter tuning methods of three PID optimization algorithms, the effects of four temperature-control algorithm models are compared. The overall algorithm structure block diagram is shown in
Figure 12.
To simulate the heating condition, this paper sets the target temperature as 30 °C and the initial temperature as 0°C. The heating curve is shown in
Figure 13. For another cooling condition, this paper sets the target temperature as −30 °C and the initial temperature as 0 °C. The cooling curve is shown in
Figure 14.
As shown in
Figure 13, BP-PID has the best performance due to its small overshoot and the fastest response speed. Compared with traditional PID, PSO-PID and Fuzzy-PID algorithm have some improvements, but there is still a certain gap compared to the BP-PID algorithm.
As shown in
Figure 14, the comparison of the above four algorithms that Fuzzy-PID has the largest overshoot, PSO-PID algorithm has small overshoot but the convergence speed is very slow, while BP-PID has the fastest convergence and small overshoot.
Considering the overshoot, convergence time, convergence accuracy and other indicators of the heating and cooling process, the final blackbody calibration source uses the BP-PID temperature-control algorithm. However, in the process of multiple tests, the temperature-control effect of BP-PID may be unstable. Therefore, when using this algorithm, the parameters of the model should be set reasonably to obtain the best effect.
5. Processing and Testing
5.1. Processing of Blackbody Calibration Target
To ensure that the above performance indicators of the calibration source simulation can be achieved, the engineering prototype of the calibration source needs to be precision machined. The main components of the prototype are described according to the processing process of the blackbody calibration source. As shown in
Figure 15a, the brass matrix forms the cone array.
Figure 15b shows the coated cone, whose surface is made of domestic P-0001 microwave-absorbing material with a density of 4.2 g/cm
3 with an attenuation of −35 dB/cm at 10 GHz frequency. During coating, the coating and curing agent is mixed at a ratio of 100:3.5, and then the cone surface is painted and cured several times to reach the required coating thickness. Next, the coated cone is placed in an oven at 120 °C for 20 min to achieve full curing.
Figure 15c coated cone array. The coated cone is fixed on the bottom plate of the thermal insulation chamber with M4 screws, and the array layout is designed according to the simulation.
Figure 15d shows the temperature controller AI-516PD2G. The internal temperature-control algorithm is implemented using the BP-PID algorithm. The relevant parameters of the algorithm are set through the serial port command protocol. The measurement range is 0–1300 °C, the resolution is 0.1 °C, the maximum allowable error is
3.9 °C and the accuracy class is 0.3 °C.
Figure 15e shows the calibration source cavity structure. The outer side of the cavity is surrounded by thermal insulation materials. The shell is made of 6061 aluminum alloy. The entire cavity is fixed on the temperature-equalizing plate of the thermostatic heating table with screws, so as to realize the temperature control of the cone array coated in the cavity.
Figure 15f shows the complete blackbody calibration source.
5.2. Electromagnetic Emissivity Test Scenario and Results
The practical testing scene includes the arched test system, vector network analyzer, 1–18 GHz double-ridged horn antenna and test cable. The diagram of the arched test system is shown in
Figure 16, and the practical testing scene is shown in
Figure 17.
The basic principle of the arch method is the dual antenna power method. A pair of antennas with the same frequency and polarization direction are used as the transmitting antenna and the receiving antenna, respectively. The transmitting antenna and the receiving antenna are placed on the same side of the tested sample. The maximum gain direction always points to the tested sample. The transmitting antenna generates excitation to the tested sample, and its reflected signal is received by the receiving antenna. The reflectivity of the tested sample is:
where
Pt,
Pr are the transmitting power of the transmitting antenna and receiving power of the receiving antenna,
Gt,
Gr are the gain of the transmitting antenna and the gain of the receiving antenna,
F represents the test frequency and
L represents distance from the phase center of the transmitting antenna to the phase center of the receiving antenna through the reflection path of the material.
The arch test system is used to locate the transmitting and receiving antennas, so that the maximum gain direction of the transmitting and receiving antennas always points to the tested sample. In addition, the absorbing material is laid on the arch frame ground. The steps to follow when using the arch test system for testing include: (1) Select the antenna with the same frequency (one transmitting antenna and one receiving antenna), and connect the antenna with the port of the vector network analyzer using high-performance test cables; (2) Install the antenna on the arch frame to ensure that the two antennas are aligned, and raise them to the top of the arch frame; (3) Place the pyramidal absorbing material on the test bench; (4) Set the frequency band of the network analyzer. Then select the S11 parameters for testing and reset the test parameters to zero; (5) Place the tested blackbody calibration on the metal plate and record the parameters.
Reflectivity and emissivity are tested and the results are shown in
Table 2. It can be seen that, in the frequency band, the emissivity
of the blackbody calibration target is greater than 0.998.
5.3. Temperature Distribution Test Scenarios and Results
At the ambient temperature of 24 °C, an electronic thermometer, K-type thermocouple, thermal conductive silicone grease and other materials were used to test the thermostatic heating table. The device was set at 100 °C, and the target temperature was reached in 3 min and 25 s. The test restarted after 5 min of stabilization. The minimum measured data was 97 °C and the maximum was 102 °C at the initial time. The device was then placed at a constant temperature for 30 min and tested again, with a minimum temperature of 99 °C and a maximum of 101 °C. The temperature test points of the thermostatic heating table are shown in
Figure 18, and the specific data are shown in
Table 3.
5.4. Uncertainty Analysis
The main sources of the uncertainty of the brightness temperature of the blackbody calibration target are: the uncertainty introduced by the temperature measurement of the blackbody calibration target ; the uncertainty introduced by the calculation of the effective emissivity of the blackbody cavity at variable temperature ; and the uncertainty introduced by the ambient temperature .
The uncertainty of the temperature measurement of the blackbody calibration target includes: the uncertainty introduced by the calibration of temperature sensor ; the uncertainty introduced by the stability of temperature sensor ; the uncertainty introduced by the electric measuring instrument ; the uncertainty of heat conduction in the cavity bottom ; the uncertainty introduced by temperature control stability and the uncertainty introduced by the uniformity of the cavity bottom temperature . Among them, is mainly determined by the BP-PID temperature-control algorithm used in this paper.
The uncertainty
introduced by the temperature measurement of the blackbody calibration target can be calculated by the following equation:
The combined standard uncertainty is:
The temperature sensor used in this paper is a K-type thermocouple, and the uncertainty introduced by its calibration is 0.1 K. According to the test, the long-term stability of the K-type thermocouple is 1.3 K, and the introduced uncertainty K.
The electrical measuring instrument used in this paper is AI-516PD2G, and the uncertainty introduced by it is 0.3 K. The uncertainty obtained through the simulation is 0.3 K, because the proportion of the area with a temperature difference between the cone surface and the bottom that is less than 0.3K exceeds 90%. The uncertainty is very small and can be ignored. The uncertainty can be obtained from the multi-point distribution experiment and has a value of 0.44 K.
The effective emissivity of the blackbody calibration target is , and therefore the calculation result of is 0.0008. Similarly, the influence of ambient temperature on brightness temperature is small enough that it can be ignored.
Finally, through calculation, we achieve the combined standard uncertainty = 0.9729.
6. Conclusions and Outlooks
The paper describes the basic principle of the calibration of microwave temperature measurement system, discusses the theoretical basis of the blackbody calibration target design, and summarizes the method of designing a calibration target based on electromagnetics simulation. Then, the scattering characteristics and temperature-distribution characteristics of the blackbody calibration target are analyzed theoretically. Based on the finite element method in the simulation software Comsol, the corresponding scattering model and heat transfer model are established. Combined with the simulation results of the two models, the parameter optimization results of the blackbody calibration source are given. Through comparative experiments, it is proven that the method can effectively improve the blackbody calibration source in terms of its emissivity and temperature uniformity. To achieve faster and more accurate temperature control, this paper compares three advanced PID algorithms—BP-PID, PSO-PID, and Fuzzy PID—by simulation. BP-PID has the advantages of a small overshoot and fast response in both heating and cooling conditions. Therefore, this algorithm is selected as the temperature-control algorithm of the blackbody calibration source. After completing the above design, this paper realized the processing of blackbody calibration source, and completed the emissivity test based on the arch test system. At the same time, the multi-point distribution method is used to simulate the temperature uniformity. The test results showed that the emissivity of the designed blackbody calibration target is higher than 0.998 in Ku-band, and the temperature uniformity is good.
In conclusions, the design of the Ku-Band blackbody calibration target has been completed. However, there are still some deficiencies that need to be addressed.
The blackbody calibration target prototype should be further tested, and the influence of emissivity and temperature-distribution uniformity on the overall calibration accuracy of the microwave radiometer should be given.
The heat transfer model is established on the premise that the temperature of the bottom surface of the cone array is consistent and uniform, and there will be temperature differences on the plane in practice. Moreover, the effect of air convection is not considered. Later, the definition of air convection will be added to the modeling, and the temperature-distribution gradient of the bottom surface of the calibration target will be set.