Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the management of heat in dielectric beams and membranes is an important part of the design in micro - electro - mechanical systems ( mems )  @xcite . while the physics of thermal transport in dielectric materials is well understood , the effect of fabrication processes on the surface physics is less clear and widely treated as a hidden variable in the evaluation of the device performance . this uncertainty can prolong the design - test cycle , where a parameter such as beam geometry is iterated until the target thermal conductance @xmath3 is obtained  @xcite .    in our application , the mems is a cryogenic detector known as the transition - edge sensor ( tes )  @xcite . figure  [ fig : q - band - detector ] shows a tes designed for measurements of the polarization state of the cosmic microwave background  @xcite . the tes is comprised of a superconducting moau bilayer deposited on a single - crystal silicon substrate . the silicon is etched around the bilayer and other electrical components to form a membrane that is supported by dielectric beams . the total thermal conductance @xmath4 of the beams is a critical parameter that determines the sensitivity ( @xmath5 ) , response time ( @xmath6 ) , and saturation power of the tes detector ( @xmath7 ) ( saturation energy in the case of a tes micro - calorimeter )  @xcite . the response time is also a function of the total heat capacity @xmath8 of the detector , which is the sum of the heat capacity of the silicon membrane and normal metal films . the thermal conductance from the membrane to the frame is effectively set by the ballistic - dominated thermal conductance of the short beam ( magnified ) . the long silicon beams are @xmath10 @xmath1m@xmath9 , and merely support tes bias and rf leads . the thermal conductance of the long silicon beams is diffusive - ballistic and sub - dominant to the conductance of the ballistic beam . the pd layer deposited on the membrane sets the detector heat capacity . , width=321 ]    at temperatures far below the debye temperature of the lattice ( @xmath11 k ) , anharmonic phonon collisions in the bulk of the dielectric are negligible , and the phonon - mean - free path @xmath0 is restricted by the cross sectional dimensions of the beam  @xcite . in state - of - the - art far - infrared and millimeter tes detectors operating at sub - kelvin temperatures  @xcite , the beam length @xmath12 is much larger than @xmath0 , and the thermal conductance is often described by the fourier law for phonon diffusion , @xmath13 , where @xmath14 and @xmath15 are the width and thickness of the beam , and @xmath16 is the bulk thermal conductivity of the dielectric . @xmath16 is determined empirically with the fabrication and testing of numerous beams of varying width and length , and the aspect ratio @xmath17 is tuned to achieve the desired conductance . however , this approach rarely provides a complete view of the thermal physics , and higher order effects such as surface roughness of the beams have been identified as the likely cause of the reported non - uniformity in @xmath3 across detector arrays  @xcite . a simple analysis can readily show the sensitivity of phonon propagation in long beams ( @xmath18 ) to the details of the surface physics . for a random gaussian surface , the roughness can be characterized by the ratio of the standard deviation @xmath19 of the surface features to the peak thermal wavelength @xmath20 of the incident phonon field . when @xmath21 , phonons are reflected specularly from the surface  @xcite . as the wavelength becomes much smaller than the scale size of the surface , @xmath22 , phonons are scattered over all angles , or diffused . assuming isotropic diffusion , it can be shown that for a rectangular beam , the phonon mean - free - path @xmath0 scales logarithmically with @xmath23 ( @xmath24 for @xmath25  @xcite ) . hence , in a beam where @xmath18 , phonons experience a large number of scattering events , @xmath26 , and the heat flux entering the beam at one end @xmath27 is reduced to @xmath28 at the exit port , where @xmath29 is the diffuse fraction . the derivative @xmath30 is a measure of the sensitivity of @xmath3 to changes in @xmath31 and proportional to the power - law function @xmath32 .    in practice , the statistics of the surface roughness can be heavily biased by the fabrication and is likely non - gaussian . dielectric beams can be roughened in various forms depending on the compatibility and properties of the etchants with the dielectric material . preliminary tests can shed light on the roughening steps , however , fabrication processes can be difficult to control in practice . details such as base pressure and process chamber conditioning during dry - etch steps can vary between batches of devices , especially in shared facilities . as a result , the characteristics of the roughness of a single surface can change significantly over spatial scales shorter than the beam length , and between wafers . furthermore , not all surfaces of a beam are rough to the same degree . these effects add to the unpredictable variability in @xmath31 and therefore @xmath3 in beams dominated by diffusive thermal transport . the work in this paper is focused on understanding the effects of fabrication on the thermal conductance in single - crystal silicon beams . the beam geometries explored span several orders of magnitude in @xmath33 . the results demonstrate that @xmath3 in short ( @xmath34 @xmath1 m , @xmath35 1 ) uniform beams with ballistic - dominated phonon transport is insensitive to the fabrication conditions and can be realized with precision in devices on different wafers . this approach is suitable for tes detector arrays in ground - based and air - borne telescopes , where the detector noise equivalent power requirements is greater than 10@xmath36 w/@xmath37 .    throughout this paper , the terms _ diffuse _ and _ ballistic _ are used to describe the scale size of the phonon transport in a silicon beam . in the presence of surface roughness , the thermal transport in a long beam with @xmath38 is dominated by diffuse reflections off the beam walls , and @xmath39 . for a short beam , ballistic transport dominates and @xmath40 . we emphasize that these limits bound the diffusive - ballistic conductance observed in practice .    in sec . [ sec : method ] , we briefly describe the device fabrication and readout hardware employed for the low - temperature measurements of @xmath3 and @xmath41 . in sec . [ sec : model ] , we describe the theory of radiative phonon transport in diffusive - ballistic beams . in sec . [ sec : diff - ball ] , we present results that illustrate the effect of fabrication on the diffusive - ballistic nature of the thermal conductance . in sec . [ sec : controlg ] , we describe the superior performance of ballistic beams for the control of @xmath3 of tes detectors . finally , we demonstrate that because of the low specific heat capacity of single - crystal silicon , the heat capacity of the detectors fabricated on this substrate can be effectively controlled using a normal metal film with high specific heat . the fabrication of the tes detectors reported in this paper is described in detail in ref . @xcite , but for completeness , we provide a brief description here . fabrication starts with a float - zone - refined p - type ( 100 ) boron doped silicon - on - insulator wafer . heat flow is primarily in the [ 110 ] direction . the single - crystal silicon layer is 5 @xmath1 m thick , with a room temperature resistivity of 5 k@xmath42-cm and a carrier density of 5@xmath4310@xmath44 @xmath45 . a moau bilayer is deposited by a combination of dc plasma sputter deposition and electron beam evaporation . the mo and au are photolithographically patterned to define the superconducting sensor . the au layer is etched by ion milling and the mo layer is etched by reactive ion etching ( rie ) . further metallic components are deposited on the device layer to enable the required operation of the pixel at rf frequencies , and to enable on - chip diagnostic capabilities such johnson noise thermometry as needed . au resistors are deposited by electron beam evaporation and patterned by lift - off . following an in - situ rf ar plasma cleaning of the wafer , a nb microstrip layer is sputter deposited , photolithographically patterned , and etched in fluorine - based rie chemistry . the silicon membrane is defined and reactive ion etched down to the buried oxide layer . the silicon wafer is then bonded to a pyrex handle wafer using a sacrificial wax . the backside of the silicon wafer is patterned by photolithography and the bulk silicon is etched by inductively - coupled plasma deep reactive - ion etching ( drie ) . the wax is then dissolved in solvent to release the devices . measurements were performed in a mechanically - cooled system that reaches 3 k at the coldest stage . further cooling of the detector stage down to 100 mk was achieved using an adiabatic demagnetization refrigerator . the devices were mounted and thermally anchored to a gold - plated copper package , and wirebonded to a superconducting quantum interference device ( squid ) cold stage amplifier for johnson noise current measurements . further amplification was achieved with a squid series array mounted at the 3 k stage . the temperature of the 100 mk cold stage was read with a calibrated ruo@xmath46 thermometer  @xcite with calibration traceable to the national institute of standards and technology . the precision of the temperature control was 32 @xmath1k rms at 150 mk . the room temperature electronics was comprised of a custom low - noise pre - amplifier , feedback electronics for flux - locked - loop operation , commercial low - noise high - gain amplifier and low - pass filter ( srs-650 ) , and commercial 16 bit digital - to - analog electronics ( ni pci-6120 ) for digitization of the amplified feedback signal . the method of johnson noise ( jn ) thermometry is frequently utilized for the characterization of the thermal properties of thin - films below 1 k  @xcite . at audio frequencies ( 0.5 - 20 khz ) and sub - kelvin temperatures , the current noise spectral density of a resistor @xmath47 in series with the squid input coil takes the form @xmath48^{-1}/r$ ] , where @xmath49 is the electron temperature , @xmath50 is the characteristic frequency of the circuit , and @xmath12 is the total inductance of the input coil and stray inductance in the input circuit loop . this equation directly relates measurements of @xmath51 to @xmath49 . for the test devices that are described in ref . @xcite , the jn thermometer and heater are 240@xmath43230@xmath430.27 @xmath1m@xmath9 au resistors deposited on top of interdigitated nb leads . the kapitza conductance at the au - silicon interface and the electron - phonon conductance in the au are approximately 10 nw / k each at 100 mk  @xcite . these conductances are a factor of 50 greater than the highest silicon beam conductance measured in our devices . we therefore assume the electron system in the thermometer is in thermal equilibrium with the phonons in the silicon membrane . the maximum expected temperature gradient across the membranes is less 0.5 mk across the 0.1 - 1 k temperature range . the random error in a temperature measurement from the jn current can be estimated from the radiometer equation @xmath52  @xcite , where @xmath53 khz is the bandwidth over which the current noise power spectral density was measured and @xmath54 s is the measurement time . hence , the fractional error throughout the temperature range of our measurements is less than 0.5% . measurements of the jn current as a function of bath temperature were used to fully calibrate the cryogenic and warm readout amplifiers . typical readout noise of 15 pa/@xmath37 referenced to the input circuit was achieved . the jn signal was 30 - 100 pa/@xmath37 in the range 0.1 - 1 k.    to measure @xmath3 of the silicon beams , a small signal dc heater power was applied to the membrane and the increase in temperature was recorded . the relationship between heater power and membrane temperature at a given bath temperature @xmath55 can be parameterized by the power law function @xmath56 , where @xmath57 is the heater power , @xmath58 is the membrane temperature , @xmath59 is a constant that is proportional to the number of phonon modes , and @xmath60 is the temperature exponent that is indicative of the temperature dependence of the phonon mean - free - path . the first order derivative yields the thermal conductance @xmath61 . in our measurements , the change in temperature was small , @xmath6210% , and the approximation @xmath63 was valid . we have also utilized tes current - voltage curves to determine the thermal conductance of detectors in which normal metal films for jn thermometry were omitted . when a tes is voltage - biased on the normal - to - superconducting phase transition , the silicon membrane temperature is fixed at the critical temperature of the moau bilayer , @xmath64 mk . the thermal conductance is then determined from a fit to the joule power dissipated in the tes as a function of bath temperature . in this case , @xmath65 and @xmath66 , where the parameters @xmath59 and @xmath60 are equivalent to those derived from jn thermometry in the same temperature range . in building a model of heat flow in a silicon beam , it is important to first identify the thermal baths that define the boundary conditions in the model . a thermal bath is ideally a blackbody that provides a matched termination for phonons at all wavelengths . in a detector such as that shown in fig . [ fig : q - band - detector ] , each silicon beam is connected to two thermal baths , the frame and the membrane . the electron system in the pd film is strongly coupled to the phonon system in the membrane . at 150 mk , the electron - phonon conductance in the pd is 30 nw / k  @xcite , and the kaptiza conductance between the pd and the silicon membrane is @xmath62 130 nw / k  @xcite . for comparison , the highest beam conductance measured is 240 pw / k at 150 mk . the silicon frame is effectively a blackbody since its lateral dimensions are much larger than the peak thermal wavelength . at 100 mk , the peak thermal wavelength is @xmath67 @xmath1 m in silicon , where @xmath68 and @xmath69 have the usual definitions and @xmath70 m s@xmath71 is the speed of sound for transverse waves . a phonon scattered from a beam into the silicon frame has a low probability of coherently backscattering into the beam . given the 5 @xmath1 m thickness of the silicon layer , and the thermal wavelengths above 100 mk ( @xmath72 @xmath1 m ) , the membrane and the frame can be modeled as three - dimensional multi - moded thermal baths , where the number of modes is given by the planck distribution and is much larger than unity . the beams behave as three - dimensional multi - moded acoustic waveguides since the widths ( @xmath73 @xmath1 m ) and 5 @xmath1 m thickness are much larger than @xmath74 at 0.1 - 1 k. phonons emitted by the baths couple radiatively to the silicon beams , propagate unimpeded in the bulk of the beams , and scatter at the beam surfaces . when the surfaces of the beams are perfectly smooth , phonons coupled to the beams are scattered elastically without any loss of coherence in the beams . the conductance in this case is determined by the beam areas and is referred to as ballistic . in the case of rough surfaces , phonons scatter diffusively with a large probability of backscatter . the scattering is elastic but the phonon wavevectors are randomized at the surface . if the surface features form cavities much larger than @xmath74 , the surface can be thought of as a blackbody . the conductance is thus boundary - limited and was first described theoretically by casimir  @xcite . an important assumption in casimir s formalism is the isotropic ( lambertian ) reemission of phonons by the beam surfaces .    in practice , the scattered pattern of phonons from a rough surface is strongly dependent on the surface characteristics . when the roughness can be described statistically , the reflection pattern is averaged over the ensemble of the surface realizations . for a surface with a gaussian correlation function , the roughness can be characterized by the correlation length @xmath75 and standard deviation @xmath19 of the surface features above a mean height . figure  [ fig : beampattern ] illustrates the scattered pattern as a function of @xmath76  @xcite . for a nearly smooth surface , @xmath77 , phonons are largely reflected in the specular direction . when @xmath78 , phonons incident on the surface are diffused in all directions , however , the flux is not lambertian . nonetheless , in numerous theoretical studies of boundary - limited scattering  @xcite , the reflected pattern is decomposed into a diffuse lambertian fraction @xmath29 and a specular fraction @xmath79 , as illustrated in fig . [ fig : beampattern ] . the diffuse reflection is modeled as inelastic scattering that redistributes the phonon wavevector , although strictly speaking all reflections are elastic at the surface facets that define the roughness . the surface can then be thought of as a blackbody with emittance @xmath29 . hence , a completely lossless but rough surface can be modeled as a lossy surface  @xcite . this approach is numerically less intensive than calculations of the reflection pattern based on the exact surface profile , or the statistics of the surface  @xcite .    ) and a very rough surface ( @xmath78 ) . the diffuse fraction of the realistic beam patterns are generally approximated as isotropic ( lambertian ) , which is relatively simple to integrate into models of thermal transport . , width=321 ]    our model of heat flow in a silicon beam has the following assumptions : ( i ) the ends of the beam constitute thermal baths with constant temperature and unity diffuse fraction . the baths continue to distances much larger than the beam length as shown in fig . [ fig : modeldiagram](a ) . this construct is numerically convenient from the perspective of the heat transport model presented , and is similar to that adopted in ref . ( ii ) the surfaces of the beam reflect phonons with a diffuse fraction @xmath29 that is isotropic , and a specular fraction given by @xmath79 . ( iii ) the gradient of the temperature profile across a beam is linearly dependent on @xmath29 and constant as a function of length . this assumption is valid in multi - moded beams , where the number of modes is much larger than unity and the diffuse phonon reflection can be treated as diffuse emission  @xcite . models that iteratively solve for the temperature profile in multi - moded beams show negligible deviations from a constant temperature gradient  @xcite . when the phonon mean - free - path is much smaller than the beam length , or @xmath29 is close to unity , phonon diffusion is the dominant term of the total heat flux , and the constant temperature gradient is an excellent approximation . the approximation also holds when @xmath29 approaches zero , or the mean - free - path is comparable to the beam length , since the temperature gradient is linearly dependent on @xmath29 ( the mathematical description of this detail is in the appendix ) . ( iv ) the phonon modes that propagate are two orthogonal transverse modes with a sound speed of @xmath80 m / s , and one dilatational mode with a sound speed of @xmath81 m / s . each mode has a thermal spectrum . ( v ) phonons propagate unimpeded in the crystalline bulk of the silicon beams , and only scatter at the beam surfaces . ( vi ) the diffuse flux is reemitted according to lambert s law of thermal radiation . ( vii ) all modes scatter diffusively with the same probability @xmath29 , and there is no spectral dependence for @xmath29 . it is possible to map @xmath82 by fitting measurements of @xmath4 and calculating the peak thermal wavelength as a function of @xmath58  @xcite . for a gaussian surface  @xcite , when @xmath77 , the specular power reflection coefficient is proportional to @xmath83 and the diffuse coefficient to @xmath84 ; when @xmath78 the reflection coefficient is independent of wavelength . hence , the thermal spectrum can be scattered very differently as temperature changes , and in recognizing the limitations of our model , we do not extrapolate the results to @xmath82 .    ) that extend to distances much larger than @xmath12 . ( b ) each surface element of the beam radiates phonons into all solid angles . the radiated flux from the infinitesimal area @xmath85 that reaches the shaded beam port is reduced depending on the diffuse fraction of the beam surfaces . ( c ) using the method of images , the path of specularly scattered phonons can be re - imaged to one vector . the components of the vector map the number of scattering events at each surface of the beam . , width=321 ]    figure  [ fig : modeldiagram](b ) illustrates the radiative model . the beam is divided into infinitesimal areas @xmath86 that radiate heat into all solid angles . the phonons are scattered several times before reaching the beam exit port and thermalizing in a bath . during each scattering event , the phonon flux is reduced by a factor @xmath87 that is the diffuse fraction of the surface @xmath88 . the nature of the problem is then to simultaneously calculate the flux that is viewed by @xmath89 from @xmath85 and count the number of scattering events from each surface of the beam . for the latter , the complex path of the phonons can be unfolded , or re - imaged , to a single vector that reaches an equivalent beam exit port as shown in fig . [ fig : modeldiagram](c ) . the components of the vector then equate to the number of scattering events from each surface . the details of the radiative model are described in the appendix  [ sec : app : a ] . we define an equivalent view factor for a rectangular waveguide @xmath90 , \label{eqn : fx - final}\end{aligned}\ ] ] where @xmath91 and @xmath92 are equivalent view factors between parallel and perpendicular surfaces of the beam respectively , and @xmath93 is the contribution of the direct ( unscattered ) flux from the surfaces of the beam to the cross sectional area at @xmath94 ( reference plane ) . @xmath95 is the diffuse fraction of the surface emitting the flux , and the index @xmath96 runs over all surfaces of the beam . the index @xmath97 runs over the surfaces perpendicular to @xmath96 . included in the terms @xmath91 and @xmath92 is a factor @xmath98 , where @xmath99 is the total number of scattering events from each surface as described in fig . [ fig : modeldiagram](c ) . for the baths , @xmath100 is obtained by setting @xmath101 in eq . [ eqn : fx - final ] . the total heat flux through the beam is @xmath102 , @xmath103dx_i ,   \label{eqn : px0}\end{aligned}\ ] ] where @xmath104 , and @xmath105 is an average group velocity . given the linear temperature profile @xmath106 , the approximation @xmath107 where @xmath58 is the average of the bath temperatures , and the relation @xmath108 , eq . [ eqn : px0 ] reduces to @xmath109 where @xmath110 is the net diffuse power from the beam surfaces , and @xmath111 is the net specular power that is exchanged between the baths . the conductance of the beam as a function of temperature when scaled by the maximum ( ballistic ) conductance @xmath112 is @xmath113 in the short beam limit , @xmath114 , and @xmath115 . for non - zero diffuse fractions and as the beam length becomes much larger than the cross - sectional dimensions , @xmath116 , and using @xmath117 where @xmath118 is the volume specific heat of the material , the mean free path takes the form @xmath119 .    when the diffuse fractions of the surfaces of the rectangular beam are different , all terms in the summation over the index @xmath96 in eq . [ eqn : fx - final ] must be evaluated . the numerical evaluation of the problem is reduced by a factor of two if there is symmetry in the diffuse fraction of the surfaces , @xmath120 and @xmath121 . for the silicon beams discussed in this paper , the top surface and sidewalls are diffuse , and the bottom surface is specular . it is important to note that the diffuse fraction of the top and sidewall surfaces are uncorrelated . the terms @xmath91 and @xmath92 in eq . [ eqn : fx - final ] are evaluated using monte - carlo integration with 10@xmath122 uniformly distributed random points averaged 10@xmath123 times . @xmath93 is evaluated with an adaptive quadrature routine . figure  [ fig : modelsimulation ] illustrates the accuracy of the model compared to an analytic expression for the mean - free - path @xmath0 in a rectangular beam  @xcite . at the diffuse limit ( @xmath39 ) , the conductance scales as @xmath124 , and at the ballistic limit , the conductance is independent of length . the expression @xmath125 is remarkably accurate over the range of beam length explored in fig . [ fig : modelsimulation ] . this fact highlights the importance of the phonon mean - free - path , that a single parameter can capture a significant fraction of the heat transport in the ballistic - diffusive limit . figure  [ fig : g - aol ] shows the thermal conductance of silicon beams as a function of beam area - to - length ratio . the cause of the variability in conductance ( a factor of @xmath625 for the same beam geometry ) is the roughening of the silicon beam surfaces during fabrication . the top surface of the beams are primarily roughened during a dry etch step that patterns a nb layer . the plasma chemistry , 75% cf4 , 15% o@xmath126 , etches silicon at a higher etch rate than nb . the exact processes that roughen the silicon are unclear , but the likely candidates are micro - masking due to non - volatile species in the plasma , and slower etch of native and in - situ nbo@xmath46 in the fluorine chemistry . the use of an etch stop such as 100 nm thick alo@xmath46 was successfully employed in early devices , however , these detectors exhibited very large uncontrolled heat capacity due to the amorphous nature of the etch stop  @xcite . test devices were also fabricated with a lift - off process that left the beam surfaces smooth . the conductance of these beams were higher , suggesting boundary - limited phonon propagation as the cause of the variability in conductance of the beams with rough surfaces . ( excluding the lift - off points ) , and is the simplest interpretation of the data that is intended to merely guide the eye . details of this data are summarized in table  [ tbl : data ] . , width=321 ]    we have explored the relationship between thermal conductance and surface roughness with an extensive dataset of atomic force microscopy ( afm ) images . non - destructive afm of a silicon beam surface a millimeter in length is not practical . instead , the silicon device layer was imaged on the handle wafer , in close proximity to each device . for a silicon surface with a gaussian roughness , @xmath127 nm and @xmath128 nm depending on the processing conditions . however , we have found that a large fraction of the roughened surfaces did not satisfy the mathematical assumptions required to apply the theory of phonon scattering by nearly - smooth gaussian surfaces , and the measured thermal conductance differed from expectations of such a theory , even when @xmath19 was much less than the typical phonon wavelength at the temperatures explored ( 0.1 - 1 k ) . we observed power spectral densities of surface height that were approximately power laws over the range of observable wavevector values . in this case , we found an empirical relationship between the power - spectral - density @xmath129 of the surface profile and the measured magnitude and temperature dependence of the phonon mean - free - path in our silicon beams , @xmath130 , where @xmath15 is the silicon beam thickness and @xmath131 is the phonon wavevector . as expected from beckmann s theoretical treatment of scattering from non - gaussian surfaces  @xcite , the observed relationship involves the correlation function of the local tilt angle of the surface facets , @xmath132 . it is challenging to precisely predict the phonon mean - free - path from information contained in the afm data due to practical sampling considerations . this occurs for several reasons : the statistics of the surface roughness are observed to vary significantly across the wafer and are sparsely sampled in @xmath133 @xmath1m@xmath123 patches set by the resolution and stability of the afm . in addition , for a multi - moded beam with thickness of order the width , details such as surface roughness of the beam sidewalls are significant , as is shown in sec . [ sec : controlg ] . non - destructive afm imaging of the sidewall geometry is extremely difficult to perform , and such a task becomes impractical when considering the number of silicon beams that constitutes a detector imaging array . however , we have found that afm images can provide qualitative insight to the thermal conductance in beams where the surface roughness relative to the thermal wavelength has a significant amplitude and is non - gaussian . suggests resonant phonon scattering from the rough surface . ( b ) afm images of 6.4@xmath436.4 @xmath1m@xmath123 patches representative of the top surface of the beams described in ( a ) . the scale size estimated from the peaks in @xmath134 is consistent with the roughness scales measured from these images . the rough surface is heavily pitted , with an average pit depth of 30 nm . the large distribution in pit size and shape is the likely cause of the observed broad resonance . since the scale size is much smaller for the smooth surfaces ( @xmath135 nm ) , the onset of the peaks occur at a much higher temperature ( shorter thermal wavelengths ) . the striping in the afm image of the smooth surface is due to residue and accumulation of particulates on the afm tip . beam dimension are @xmath136 mm , and @xmath137 @xmath1 m for devices labelled with squares , diamonds , and triangles respectively.,width=321 ]    figure  [ fig : cas - lim](a ) shows the conductance of smooth ( lift - off ) and roughened ( dry etch ) silicon beams when compared to the boundary - limited conductance in each case  @xcite . the phonon propagation transitions from diffuse to ballistic at 300 mk , where the peak thermal wavelength is @xmath138 nm . the minimum in conductance suggests the transition is resonant for the rough beam  @xcite . the inset in fig . [ fig : cas - lim](a ) shows the inverse phonon mean free path @xmath134 . the scale height of the scatterers @xmath139 can be estimated from the temperature at which @xmath134 peaks . using the rayleigh criterion , @xmath140  @xcite , where @xmath74 is the phonon wavelength at the peak temperature , and @xmath141 is the phonon angle of incidence to the surface . assuming an average over all angles , @xmath142 , the rayleigh criterion suggests @xmath143 nm for the rough beam and @xmath144 nm for the smooth beams . these estimates are consistent with the roughness measured from the afm images shown in figure  [ fig : cas - lim](b ) . the rough ( dry etch ) surface is heavily pitted , with an average pit depth of 30 nm and average size of 250 nm that is comparable to @xmath74 at 300 mk . the scale size is much smaller for the smooth ( lift - off ) surfaces , @xmath135 nm , which leads to the onset of resonance at a much higher temperature . in addition to phonon scattering by rough surfaces of a beam , organic residue such as photoresist on the surfaces can enhance diffuse scattering . thin layers of metals and non - metals hundreds of angstroms thick have been shown to induce diffuse transport in silicon  @xcite . in our devices , the reduction and variability in @xmath3 in a few cases were at least qualitatively attributed to residue accumulation on beams of order a millimeter in length . care must be taken when devices are cleaned since the cleaning process may introduce further residue , or roughness ( see example of silicon dipped in hf in ref . @xcite ) .    table  [ tbl : data ] summarizes the silicon beam geometries and conductances shown in fig . [ fig : g - aol ] . the exponent @xmath145 is inversely proportional to @xmath3 , indicating a transition from diffuse to ballistic transport . low @xmath3 devices have boundary - limited transport , and @xmath146 . high @xmath3 devices are approaching ballistic transport , and @xmath147 in the transition region . .conductance fit parameters @xmath148 and @xmath145 as a function of beam geometry , @xmath149 . the entries marked with an asterisk indicate beams with smooth top surfaces ( lift - off devices described in the text ) . [ cols="^,^,^,^,^ " , ]     [ tbl : data ] in the design of prototype tes detectors , the choice of the beam length and width was guided by the data in fig . [ fig : g - aol ] . several strategies were adopted in order to control the conductance and produce uniform detectors for the focal plane array described in ref . @xcite . to minimize the top surface area of the silicon beams roughened by the dry etch , the width of the nb leads were extended to cover a larger fraction of this surface , and the beam widths were defined with a combination of drie using the standard bosch process  @xcite followed by rie to remove the silicon down to a nb ground plane . the two processes are separated to reduce potential metallic contamination in the drie chamber while minimizing lateral silicon undercut from the isotropic rie . the bottom surface of the beams in contact with the nb ground plane are smooth , and scanning electron microscopy images of the sidewalls showed no significant features from the rie etching of the beams employed in test structures , fig . [ fig : sidewall - sem](a ) . the beam sidewalls were predicted to be nearly specular , with conductances comparable to the lift - off devices shown in fig . [ fig : g - aol ] . the top surface of the beams was expected to be slightly rough , @xmath150 , since the area exposed during the dry etch step was reduced by a factor of two . four beams 785 @xmath1 m in length and 13 @xmath1 m in width were integrated into the tes design . post fabrication , the beams had well - defined top surfaces , but the sidewalls were roughened by the drie steps , fig . [ fig : sidewall - sem](b ) . the total conductance was a factor of 5 below the predicted target range . in principle , the scale length and height of the sidewall features that include the drie scalloping as well as passivation residue present a very diffuse surface for phonons above 100 mk . figure  [ fig : longbeam - model ] compares the total conductance of the four beams to simulations based on the diffusive - ballistic model described in sec . [ sec : model ] . the diffuse fractions of the sidewall and top surface of the beams vary across a wafer and between wafers . the inferred diffuse fraction of the top surface of the beams is @xmath151 . a possible source for the moderate roughness is the 200 nm thick nb lead sputter - deposited on each beam . the lead and beam are generally considered as decoupled systems , and the nb conductance is estimated to be far below that of the silicon beam  @xcite . however , it is the diffusion process of the phonons in the silicon by the nb layer that is of importance and generally unknown . it has been shown that even in the presence of large acoustic mismatch , phonons in silicon can be diffused at surfaces with metallic layers tens to thousands of angstroms thick ( see fig .  8 in ref . @xcite for examples of phonon diffusion caused by al and ti films below 1 k ) . , devices marked with crosses , squares , diamonds , and circles have similar sidewall profiles . ( b ) to minimize the top surface area roughened by the dry etch of the nb lead layer , drie was used to define the beam width . the ridges at the top of the beam are scalloping from the drie cycles . possible sources for the bright features on the sidewall are leftover passivation layer , and charged silicon features . , width=321 ]    it is evident from the results in fig . [ fig : longbeam - model ] that a solution to the problem of conductance control is a design that is least sensitive to the effect of fabrication conditions such as surface roughening and residue . our approach was to reduce the length @xmath12 of a beam , which increases the phonon mean - free - path @xmath0 by reducing the interaction probability of the phonons with the rough sidewalls of the beam . the conductance is less sensitive to @xmath0 in the limit @xmath152 . figure  [ fig : shortbeam ] shows the simulated and measured total conductance when a ballistic beam 10 @xmath1 m long and 13 @xmath1 m wide is integrated into the detector as shown in fig .  [ fig : q - band - detector ] . given our knowledge of the surface roughness after the dry etch step ( fig . [ fig : longbeam - model ] ) , a diffuse fraction @xmath153 was expected and subsequently obtained in practice . m long beams shown in fig .  [ fig : q - band - detector ] . the measured conductance was a factor of @xmath625 below the design conductance ( hatched box ) . the cause of the variability is boundary - limited scattering . ( right ) simulations of the conductance as a function of beam wall diffuse fractions . the inferred diffuse fraction of the samples are moderately rough top surfaces , @xmath154 , and very rough sidewalls , @xmath155 . , width=321 ]    . the shaded and hatched regions are defined by the measured range of conductance of the long diffuse beams shown in fig . [ fig : longbeam - model ] . ( right ) the measured conductances are within the predicted range of the simulations ( light grey ) . , width=321 ]    the total conductance , which is largely set by the ballistic beam , was predicted prior to fabrication . the diffuse beams merely support microstrip and bias leads necessary for the operation of the detector . however , they contribute the largest fraction to the variance in the total conductance across a wafer and between wafers . the magnitude and variance in @xmath3 of these beams were estimated from the data in fig . [ fig : longbeam - model ] . the fractional standard deviation in total conductance is @xmath28% and within the requirement of the target instrument  @xcite . the range of conductances that can be accessed with a ballistic beam is suitable for low - noise bolometric detectors operating at ground - based or air - borne telescopes where the signal power is of order 1 - 3 pw / k . in the ideal design , a single ballistic beam sets the conductance of the detector and carries all signal and bias leads . this approach is very well suited to the fabrication of kilo - pixel arrays , especially when the signal is coupled to the detector with free - space absorbers . for space applications where the signal power is less than 1 pw , the ballistic beam can be a phononic structure with phonon bandgaps that can be engineered to reduce the conductance  @xcite . the heat capacity of the dielectric also plays an important role in the performance of thermal detectors . for a bolometer and calorimeter , @xmath41 affects the detector time constant , and in some cases anomalous heat capacity associated with the dielectric has been reported to increase thermodynamic noise of these detectors  @xcite .    the heat capacity of the silicon membranes were determined from @xmath156 , where the thermal time constant @xmath157 was measured from the decay in temperature of a silicon membrane to a pulse in heater power . compared to the debye specific heat  @xcite , @xmath158 at 200 mk and 0.7 at 1 k. the systematic difference between the measurements and the debye approximation is due to the large contribution to the total heat capacity of the devices from the normal metal au heater / thermometer and pd - au termination on the silicon membranes . nonetheless , the measured specific heat are extremely small . for the detectors shown in fig . [ fig : q - band - detector ] , the estimated silicon heat capacity is @xmath622 fj / k at 150 mk , which is advantageous since the overall detector heat capacity can be controlled with a normal metal that can be reliably deposited on the membrane . we have successfully used pd to achieve a heat capacity of 2 pj / k  @xcite . measurements indicate that the specific heat of pd below 1 k is 5.8 mj mol@xmath71 k@xmath159 , which is lower than the room temperature value of 9.4 mol@xmath71 k@xmath159  @xcite . we have successfully deposited 400 nm pd layers with a lift - off process . it is important to emphasize that there is no evidence for anomalous specific heat in the single - crystal silicon device layers . this is in contrast to the large excess specific heat of amorphous sin@xmath46 that is widely used in the fabrication of thermal detectors  @xcite . the thermal conductance of single - crystal silicon devices can be precisely controlled using a short beam with ballistic - dominated phonon transport . the conductance of the short beam is largely set by its area and is insensitive to the details of the fabrication conditions that can vary over time . this approach is superior to the conventional use of long beams for the control of thermal conductance . phonon transport in long beams is diffusive - ballistic and very sensitive to the detailed surface physics of the beams . the conductance is thus boundary - limited , and has been observed in beams with modest surface roughness ( 5 - 10 nm rms in surface height ) . resonant phonon scattering has also been observed at 300 mk in beams with highly pitted surfaces with pit depths of order 30 nm . boundary - limited scattering has contributed to a large variance , up to a factor of 5 , in the conductance of devices fabricated and tested during this work . in contrast , when a beam with ballistic - dominated thermal transport is integrated into the design , the uniformity in conductance is reduced to a fractional deviation of @xmath28% in devices fabricated across two wafers , and this variability is largely set by the conductance of diffuse beams that support electrical and microwave leads to the device . in addition , we have found no evidence of excess specific heat in single - crystal silicon membranes . hence , the total heat capacity of the detector can be effectively controlled with a normal metal film . for the transition - edge sensors ( tess ) described in this work , the heat capacity was determined by a pd layer 400 nm thick . hence , the strategies outlined for the precision control of thermal conductance and heat capacity are well suited to the fabrication of uniform large - format arrays of tess with sensitivities approaching 10@xmath36 w / hz@xmath160 . we gratefully acknowledge financial support from the nasa roses / apra program . k. rostem was supported by the nasa postdoctoral program at the goddard space flight center during the initial stages of this work . we thank samelys rodriguez for wirebonding and cryo - cable support . we thank our collaborators a. ali , j. appel , and t. a. marriage for their contribution to the data acquisition . we thank the national institute of standards and technology for providing the squid series array and time division multiplexer chips . we start from the usual definition of the infinitesimal view factor @xmath161  @xcite , @xmath162 where @xmath163 is the vector between the infinitesimal surfaces @xmath164 and @xmath165 , and @xmath166 are the angles between the normal vectors to @xmath167 and @xmath163 ( see fig .  [ fig : modeldiagram ] ) . we define an equivalent view factor for convenience , @xmath168 and note that the heat exchanged between the infinitesimal surfaces is @xmath169\,f_{da_i - da_j},\end{aligned}\ ] ] where @xmath170 are the temperatures of the surface elements . the total heat flux is the integral of @xmath171 over the beam length , @xmath172 where @xmath173 represents the surfaces of the waveguide , which in our case is a rectangular beam . the equivalent view factor between any two perpendicular surface elements of the beam is @xmath174 ^ 2 } \left [ \,f_m\prod_k ( 1-f_k)^{n_k}\right ] \right\ } \,da_i\,da_j , \\          & = & \mathcal{f}^\perp ( m ; x_i , y_i , y_j , z_j ) \,dx_i\,dy_i\,dy_j\,dz_j , \nonumber   \end{aligned}\ ] ] where the indices @xmath96 and @xmath88 sum over all surfaces of the beam . for parallel surface elements , @xmath175 ^ 2 } \left [ \,f_m\prod_k ( 1-f_k)^{n_k}\right ] \right\ } \,da_i\,da_j , \\          & = & \mathcal{f}^\parallel ( f_m ; x_i , y_i , y_j , z_j ) \,dx_i\,dy_i\,dx_j\,dy_j . \nonumber\end{aligned}\ ] ] @xmath99 is the number of specular reflections experienced by phonons reflecting off the element @xmath89 first , and then all other surfaces until a reference plane ( e.g. @xmath94 in fig . [ fig : modeldiagram](a ) ) is reached . @xmath99 can be determined from the geometry of the specular scattering as described in fig . [ fig : modeldiagram](c ) , and is the sum of scattering events from each surface . for example , if the radiating element is on the top surface , @xmath99 can be found for phonons reflecting off the top and bottom surfaces , and reflecting off the sidewalls only . in this way , a different diffuse fraction can be assigned to each surface . it is also useful to note that the linear dependence of the equivalent view factors on @xmath95 can be translated to a linear dependence of the temperature gradient on the diffuse fraction when eq . [ eqn : perp ] and  [ eqn : para ] are finally substituted into eq . [ eqn : px01 ] . the integral equations of the equivalent view factors over the surfaces of the beam are @xmath176 where the limits of the integrals depend on the surface @xmath96 from which the phonons were originally radiated . since a temperature gradient is only present in the principal direction of heat flow ( x - axis ) , all surface elements at the plane of constant @xmath177 are at the same temperature and radiate the same flux . hence , the equivalent view factor for the direct ( non - scattered ) heat flux @xmath93 from the surface to the reference plane can be calculated analytically from eq . [ eqn : vf ] by considering the view factor of a infinitesimal strip of width @xmath14 to the reference plane at @xmath94  @xcite .
Abstract:

Answer: we report on the diffusive - ballistic thermal conductance of multi - moded single - crystal silicon beams measured below 1 k. it is shown that the phonon mean - free - path @xmath0 is a strong function of the surface roughness characteristics of the beams . this effect is enhanced in diffuse beams with lengths much larger than @xmath0 , even when the surface is fairly smooth , 5 - 10 nm rms , and the peak thermal wavelength is 0.6 @xmath1 m . resonant phonon scattering has been observed in beams with a pitted surface morphology and characteristic pit depth of 30 nm . hence , if the surface roughness is not adequately controlled , the thermal conductance can vary significantly for diffuse beams fabricated across a wafer . in contrast , when the beam length is of order @xmath0 , the conductance is dominated by ballistic transport and is effectively set by the beam area . we have demonstrated a uniformity of @xmath28% in fractional deviation for ballistic beams , and this deviation is largely set by the thermal conductance of diffuse beams that support the micro - electro - mechanical device and electrical leads . in addition , we have found no evidence for excess specific heat in single - crystal silicon membranes . this allows for the precise control of the device heat capacity with normal metal films . we discuss the results in the context of the design and fabrication of large - format arrays of far - infrared and millimeter wavelength cryogenic detectors .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: cognitive radio ( cr ) is a promising approach to efficiently utilize the scarce rf spectrum resource . in this paradigm , knowledge about spectrum occupancy in time , frequency , and space that is both accurate and timely is crucial in allowing cr networks to opportunistically use the spectrum , and avoid interference to a primary user ( pu ) . in particular , information about pu location will enable several key capabilities in cr networks including improved spatio - temporal sensing , intelligent location - aware routing , as well as aiding spectrum policy enforcement . the localization problem in cr networks is in general different from localization in other applications such as wireless sensor networks ( wsn ) and global positioning system ( gps ) , in which the target to be localized cooperates with the localization devices . in contrast , a pu does not communicate directly with the crs during the localization process . this scenario is referred to as non - cooperative localization and does not allow the use of conventional techniques such as time - of - arrival ( toa ) and time - delay - of - arrival ( tdoa ) ranging . several techniques have been proposed in the literature for estimating the position of active pus in a cr network . these schemes include both range - based @xcite and range - free @xcite algorithms . range - free positioning schemes , such as centroid localization @xcite , have attracted a lot of interest because of their simplicity and robustness to changes in wireless propagation properties such as path loss . this characteristic makes them suitable candidates for systems requiring coarse - grained , but reliable and cost - effective techniques . on the other hand , range - based techniques such as multi - lateration @xcite , offer better estimates , but require accurate information about the target s distance that can only be achieved either through cooperation with the target or a very precise knowledge of the path loss model .    in this paper , we focus on analyzing the performance of a particular , low complexity , coarse - grained localization algorithm , referred to as weighted centroid localization ( wcl ) @xcite . in this technique , pu location is approximated as the weighted average of all secondary user positions within its transmission range . strictly speaking , wcl is not an entirely range - free technique because it requires additional information aside from simple connectivity , namely received signal strength ( rss ) measurements . wcl and wcl - like methods are suitable for pu localization because of their low computational complexity . they also do not require cooperation from the pu and only use the readily available rss information . in addition , wcl always converges to a solution , as opposed to some range - based techniques , and do not require estimation of the path loss exponent beforehand . recent papers on wcl have proposed improvements on the basic algorithm to enhance its performance for specific scenarios @xcite . all previous wcl studies reported in the literature have been analyzed either via simulation or field experiments @xcite . to the best of our knowledge , there has been no analysis of the probability distribution of the localization error for even the simplest version of wcl in the presence of shadowing in both wsn and cr literature . we present the first theoretical analysis of the error distribution of wcl in this paper . the presented analytical framework models varying levels of shadowing , including both independent and correlated shadowing environments , and also takes inaccuracy of sensor node positioning into consideration . using this analysis in conjunction with numerical simulations we observe and quantify the robustness of wcl to variable channel conditions and node placements . we also investigate variable node participation as a simple modification that can improve wcl performance .    in the latter part of this paper we propose and investigate a distributed cluster - based implementation of wcl in the context of cr networks . previous work primarily assumes a centralized fusion node aggregating data from other nodes . in the context of wsn where wcl was first utilized @xcite , the primary motivation was node self - localization where a node applies wcl using measurements from beacon nodes in close proximity to it . this approach naturally dictates a centralized system where the node to be localized acts as the fusion center and only requires minimal interaction with the beacons . in the case of pu localization in cr networks , however , where the primary motivation is estimating the position of a non - cooperative node , there no longer exists any particular node that has all the requisite information to perform wcl . in this case a fusion node is arbitrarily chosen . this approach does not scale very well in terms of communication overhead and transmit power , which limits its application to low node densities or small network size . there is therefore a need to create distributed implementations of wcl for this technique to be practically applied in cr localization . a standard approach for distributing control in sensor networks is clustering . several algorithms exist in the literature for creating either disjoint @xcite or overlapping @xcite clusters . in @xcite , a distributed technique based on the concept of mean - shift @xcite . however , the technique requires repeatedly forming a virtual cluster until it converges to the target location . this introduces additional communication overhead and also requires several iterations to converge . we present an implementation which saves total transmit power needed for localization process communication , and investigate its accuracy as well as its effect on per node computational complexity and total transmit power consumption . the rest of the paper is organized as follows . the system model is introduced in section [ sec : systemmodel ] . the theoretical analysis of centralized wcl accuracy is presented in section [ sec : analysis ] , while a practical and distributed implementation of wcl is proposed and analyzed in section [ sec : distributed ] . numerical results for various scenarios used to evaluate performance of wcl algorithms are presented in section [ sec : evaluation ] . finally , the paper is concluded in section [ sec : conclusion ] . this section presents the general framework for our analysis . we first review the assumptions about the system setting and propagation model , and then present the basic wcl algorithm . we adopt the pu - centric model @xcite in which the pu transmitter is assumed to be at the center area of a circle . the rationale behind this choice is that in a practical scenario , sensor nodes that obtain enough received power ( i.e. above sensitivity threshold ) approximately form a circle around the primary transmitter . thus , sensors within this circle are used to jointly estimate the position of the pu . an implicit assumption in this model is that the pu is surrounded by enough nodes so that edge effects can be disregarded . assume @xmath0 sensor nodes in a circle of fixed radius @xmath1 . the number of sensors within the circle is jointly affected by several factors , such as transmit power of the pu ( which effectively determines @xmath1 ) , node density , node selection methodology , etc . @xmath0 is a design parameter for wcl deployment , the choice of which will have significant influence on the localization accuracy , as will be shown in our analysis . the pu located at the center has a coordinate defined as @xmath2^t=[0 , 0]^t$ ] . for the @xmath3th sensor node , we define its 2-dimensional location as @xmath4^t$ ] where @xmath5 . sensors are placed independently in both dimensions . we assume sensors obtain information about their own locations , which is modeled as imperfect . the measured location of the @xmath3th sensor is distributed as @xmath6 , where @xmath7 defines standard deviation of the error in sensor node self - localization . we further assume the self - localization errors are independent among sensors and are independent from received powers . we adopt the following channel model for received signal strength used in wcl . the received power of the @xmath3th node from the primary transmitter , @xmath8 , is given by @xmath9 in which the first two terms characterize path loss and the last term describes shadowing effect . the path loss is characterized by three constants : reference power @xmath10 , reference distance @xmath11 and path loss exponent @xmath12 . denote collection of shadowing variables as @xmath13 $ ] , then its distribution is characterized by @xmath14 . we consider the following shadowing cases    _ = _ s^2 * i_n , & i.i.d case , + \{_}_ij = _ i-_j /x_c , & correlated case , [ shadowing ] *    where @xmath15 is the correlation distance @xcite within which the shadowing effects of nodes are correlated , and @xmath16 is the identity matrix with dimension @xmath0 .    the theoretical analysis presented in section [ sec : analysis ] is based on the centralized architecture , which will also give us an upper bound on the performance . in each localization period , all sensors involved in the wcl will send their rss measurements to the fusion center , where an estimate will be made about the coordinate of the primary transmitter . the considered algorithm is a modification of relative - span wcl proposed in @xcite . in this scheme , the weights are guaranteed to be non - negative in db scale . the 2-dimensional estimated location of pu is formed as @xmath17 the weight factor of node @xmath3 is defined as @xmath18 , @xmath19 where @xmath20 is a constant calculated from average received power of the node on the border with some margin , to guarantee the probability that receive power of all nodes in the area below @xmath20 is sufficiently small , say 1% . @xmath21 is the maximum received power among nodes , and @xmath22 is the span of received power . the localization error is then given by @xmath23^t$ ] , where @xmath24 and @xmath25 are the 1-dimensional errors along the x and y - axis respectively . finally , the performance of the wcl scheme is evaluated by distance error , which is given by @xmath26 this section presents the analysis of wcl error using probability theory , ultimately leading to analytical expressions for the error distribution , first of wcl applied to 1-dimensional coordinates , and finally the 2-dimensional wcl error given by ( [ error - dist ] ) .        first we derive the 1-dimensional localization error for the i.i.d . shadowing case . we define a constant @xmath27 for each node @xmath3 as @xmath28 one can verify that @xmath29 . the location estimation in one dimension , say , x - axis , is naturally given by @xmath30 based on our i.i.d assumption among all @xmath8 s , one can verify that @xmath31 for the statistical property of @xmath32 , since @xmath33 s are also i.i.d and are independent of @xmath8 s , the mean of @xmath32 is given by @xmath34 . the variance of the @xmath3th term of @xmath32 is given by    rcl [ variance a ] [ ( p_i - p _ ) x_i ] & = & [ ( p_i - p_)^2 ] [ x_i^2 ] - [ p_i - p_]^2 [ x_i]^2 + & = & _ l^2 _ i^2 + _ s^2 _ i^2 + _ l^2 _ s^2 .    thus the variance of @xmath32 is derived as @xmath35 . from central limit theorem , we approximate @xmath32 as a single gaussian variable , i.e. @xmath36 , which is verified by numerical simulations that are not presented here due to space limit . since @xmath32 and @xmath37 are correlated due to summation over the same set of variables , i.e. @xmath8 s , @xmath24 is a ratio of two correlated gaussian variables . in order to obtain statistical properties of @xmath24 , we first calculate correlation coefficient @xmath38 of @xmath32 and @xmath37 , which is given by @xmath39/(\sigma_a \sigma_b )        = \operatorname{\mathbb{e}}\left [ \sum_{i=1}^{n } \sum_{j=1}^{n } x_i s_i s_j \right]/(\sigma_a \sigma_b ) . \label{correlation - coeff1}\ ] ] note that    = _ i _ s^2 , & @xmath40 , + 0 , & @xmath41 , [ correlation ]    thus the correlation coefficient is obtained as @xmath42 where @xmath43^t$ ] , and @xmath44 is the @xmath45th norm of the vector .    from the previous analysis we know that @xmath24 is the ratio of two dependent gaussian random variables with non - zero means for which the closed - form pdf is available in ( * ? ? ? * eqn 7.14 ) . based on the pdf of @xmath24 , we calculate its mean @xmath46 and variance @xmath47 by numerical integration . alternatively , we also calculate @xmath46 and @xmath47 through approximation results of @xcite , which assume @xmath24 follows normal distribution . the approximation results are given by    rcl m__p & & ( m_a / m_b ) + _ b^2 m_a / m_b^3 - _ ab _ a _ b /m_b^2 + _ _ p^2 & & _ b^2 m_a^2/m_b^4 + _ a^2/m_b^2 - 2 _ ab _ a _ b m_a /m_b^3 . [ approx ]    the gaussian approximation ( [ approx ] ) is verified via numerical simulations . specifically , the mean and variance calculated from integration of the closed - form pdf are compared with the approximated mean and variance , for different node numbers . our results show that for @xmath48 , the approximation error for the means is essentially zero and less than 3% for the variances .    up to this stage , we are able to calculate the mean and variance of the estimation error in the x - axis as @xmath49 , @xmath50 . the same analysis applies for calculating the mean and variance of the estimation error along the y - axis , namely , @xmath51 and @xmath52 . note that the estimation errors in x - axis and y - axis are correlated , therefore the statistical properties of the 2-dimensional error , @xmath53 , can not be obtained directly .      in this case , since @xmath8 s in ( [ est x ] ) are correlated , thus the distribution of @xmath24 needs to be evaluated again . since @xmath37 in ( [ est x ] ) is the linear combination of correlated gaussian variables , it is still gaussian distributed @xcite , i.e. @xmath54 . thus , @xmath55 , and @xmath56 is given by @xmath57 , where @xmath58 is the @xmath59th element of the covariance matrix of @xmath60^{t}$ ] , defined as @xmath61                   = \operatorname{\mathbb{e}}[s_i s_j ] = \sigma_s^2 \lambda_{ij } , \label{covariance_b}\ ] ] where @xmath62 . for the distribution of @xmath32 , it follows directly that @xmath63 . its variance is given by @xmath64 , where @xmath65 is the @xmath59th element of the covariance matrix of @xmath66^{t}$ ] , which can be represented by @xmath67[(p_j - p_{min})x_j- x_j \mu_j]\ } \ieeenonumber \\ = \operatorname{\mathbb{e}}[x_i s_i x_j s_j ] , \label{covariance_a}\ ] ] where @xmath68 $ ] is given by     = _ s^2 ( _ i^2 + _ l^2 ) , & @xmath40 , + _ i _ j _ s^2 _ ij , & @xmath41 . [ correlation ]    based on central limit theorem , we approximate @xmath32 as a gaussian random variable , i.e. @xmath69 . the approximation is validated via simulations , and it applies for @xmath70 .    the correlation coefficient @xmath38 is given by @xmath71/(\sigma_a \sigma_b )                                 = ( \sigma_s^2 \sum_{i=1}^{n }                                     \sum_{j=1}^{n } \overline{x}_i \lambda_{ij } ) / ( \sigma_a \sigma_b ) . \label{correlation coeff}\ ] ] substituting the expressions for @xmath72 and @xmath73 into ( [ correlation coeff ] ) , after some simple manipulation , @xmath74 where @xmath75 and @xmath76 is a @xmath77 vector of all @xmath78 s . since the location estimate @xmath24 is formed by @xmath79 , we use the closed - form pdf given by ( * ? ? ? * eqn 7.14 ) to calculate @xmath46 and @xmath47 by numerical integration . gaussian approximation can also be used similarly to i.i.d shadowing case in sec . [ subsubsec : iidshadowing ] .      in our previous analysis , we obtained the statistical properties of 1-dimensional location estimation error for both i.i.d . and correlated shadowing case . from gaussian approximation of 1-dimensional localization error ( [ approx ] ) , we know that @xmath80 , where @xmath81^t$ ] , @xmath82,\ ] ] and @xmath83 is the correlation coefficient between 1-dimensional errors , the calculation of which is presented in the appendix . the calculation of localization error @xmath53 involves computing 2-norm of @xmath84 . since the two variables in @xmath84 are dependent gaussian r.v.s with different variances , we perform the following transformations . first we form two independent gaussian variables from @xmath84 , by the standard technique of decorrelation @xcite , @xmath85^t$ ] , where @xmath86 is the orthogonal eigenvalue decomposition matrix of @xmath87 , which satisfies @xmath88 , where @xmath89 is a diagonal matrix . therefore , one can verify that ( * ? ? ? * eqn 2.1 - 159 ) @xmath90 . then by the orthogonal - invariant property of 2-norm @xcite , we know @xmath91 .    in order to obtain the pdf of @xmath92 , let us first consider pdf of @xmath93 . @xmath94 and @xmath95 are two chi - square variables with one degree of freedom . then the pdf of @xmath96 is the convolution of these two chi - square pdfs , denoted as @xmath97 . the pdf of the localization error @xmath98 is given by @xmath99 . as a result , the closed - form pdf of @xmath53 , for given sensor node positions , is given by    rcl p_e_l(x ) & = & ( - ) + & & _ i=0^ _ l=0^ , x 0 [ pdf_2d ]    where @xmath100 , @xmath101 and @xmath102 , @xmath103 are means and variances of @xmath104 and @xmath105 , respectively , and @xmath106 and @xmath107 are gamma function and @xmath108th - order modified bessel function of the first kind , respectively . the mean and variance of @xmath53 , denoted as @xmath109 and @xmath110 , are calculated using their definitions , by numerical integration over the above pdf . note that the above calculations are based on a particular fixed node placement . if the nodes are randomly placed within the area according to some distribution , say uniform distribution , integrating @xmath111 times over all node coordinates gives the theoretical average performance . based on the theoretical analysis , wcl requires large number of nodes in order to achieve a high level of accuracy . this causes large communication overhead in terms of number of messages and transmit power consumption when implemented in a centralized manner . we now present a practical distributed implementation of wcl to address this problem in the context of cr networks . clustering is commonly used in implementation of distributed algorithms . the challenge in designing these algorithms is to reduce the amount of information exchanged among clusters . the proposed algorithm accomplishes this by splitting the task into two phases . first , the cluster with the highest average rss is selected . this cluster then acts as the coordinator for the second phase which performs the actual wcl algorithm . nodes , @xmath112 , in the cr network are assumed to be organized into clusters , @xmath113 , using the clustering algorithm presented in @xcite . this particular algorithm is based on geographic boundaries and generates clusters which have a predefined shape ( hexagonal in this case ) . the regularity in the shape of the cluster is beneficial for the first step of the distributed algorithm since wcl works best for a radially symmetric distribution of nodes . however , the proposed algorithm would work for any arbitrary clustering algorithm as long as nodes are clustered based on proximity . the pu is assumed to be randomly placed in the entire area and can thus be in any of these clusters . each cluster selects one node , @xmath114 , to be the cluster head . the cluster head maintains a table of all adjacent clusters , @xmath115 , and the location of its member nodes , @xmath116 . once a pu is detected by a node , using an existing pu detection method such as energy detection @xcite , it transmits its rss measurement to the corresponding cluster head which then participates in the distributed wcl algorithm outlined in the next section . using a predetermined fusion rule , such as a majority vote , the cluster head makes a decision on the presence of the pu . all clusters that decide that the pu is present , participate in the distributed wcl algorithm detailed in the next subsection . these clusters form a set defined as @xmath117 . if in addition , some prior knowledge of the pu transmit power is available , a situation typical of a cr network coexisting with a licensed incumbent network , then a predetermined threshold can be used to further reduce the number of clusters in @xmath117 . the threshold can be chosen based on the average node spacing to get a certain desired average number of nodes . we investigate the effect of this number on wcl accuracy in sec . [ subsec : variableparticipation ] .      the first phase , which selects the head cluster , is shown in algorithm [ alg : selection ] . the objective is to select the cluster with the highest average rss , @xmath118 . this is done in a distributed manner by comparing the average rss between neighboring clusters . however , directly implementing this procedure would require a large amount of inter - cluster communication .    to solve this issue , a cluster only communicates its average rss to the neighboring cluster , @xmath119 , in the direction of the pu relative to the cluster s centroid .    although techniques for estimating the angle of arrival of a transmitter could be beneficial , these methods require cluster heads to have multiple antennas which increases complexity of sensors . in order to get an approximation of the pu direction , we propose a metric , @xmath120 , shown in line 5 of algorithm [ alg : selection ] . we get @xmath120 by subtracting the wcl result of all nodes within the same cluster from the geometric centroid of those nodes . this vector approximates the rss gradient within the cluster . it corresponds to an approximate direction of the pu since on average the rss in ( [ received power ] ) decreases as distance to pu increases . therefore , the negative of the gradient is directed towards the point with maximum rss . once the cluster with the highest average rss is selected , algorithm [ alg : mwcl ] is performed . in this step , the final wcl calculation is done using nodes within a given radius , @xmath121 , of the selected cluster s strongest node ( sn ) , @xmath122 . @xmath121 is adaptively chosen to be the minimum of the distance to the closest map edge , @xmath123 , and the cluster radius @xmath124 . this reduces the border effect , which is the tendency of the location estimate to be biased towards the center when the pu is close to the map border , by ensuring a radially symmetric node distribution even when the pu is close to the network edge . in order to improve the localization accuracy , additional nodes from adjacent clusters are also included in the wcl calculation . using the location of @xmath122 and @xmath121 , these adjacent clusters can reduce the amount of information exchanged by only transmitting the location and rss measurements of nodes that fall within this circle . in the case of large clusters , each cluster could also choose to calculate a weighted centroid using its own nodes and transmit the result of this calculation only . since wcl is essentially an averaging algorithm , this technique allows even algorithm [ alg : mwcl ] to be distributed . @xmath125 @xmath126 @xmath127 \mathbf{l}(n_j ) } { \sum_{n_j\in c_i}(\operatorname{p}(n_j ) - \min(\operatorname{p}(n_j))}$ ] @xmath128 @xmath129 select @xmath130 proceed to algorithm [ alg : mwcl ]    @xmath131 @xmath132 @xmath133 @xmath134 \mathbf{l}(n_j ) } { \sum_{n_j \in c^*_{wcl}}\left[\operatorname{p}(n_j ) - \min(\operatorname{p}(n_j)\right]}$ ]    * function *  @xmath135 * return * @xmath136 ,   @xmath137 * end function *     the power efficiency and computational complexity of this distributed approach could further be improved by taking into account a practical mobility model of the pu . there has been a lot of work such as @xcite , which discussed how to incorporate past observations of the target position and a mobility model to improve localization . however , applying these methods introduces complexity to wcl and defeats the purpose of such a low complexity scheme .    intuitively , assuming the pu continuously transmits during successive wcl estimations , a pu would most likely be in the same cluster as it was in the previous wcl calculation , or it could move to an adjacent cluster . by considering this fact , algorithm [ alg : selection ] can further limit the clusters involved in choosing the head cluster , which is the most power consuming phase , by only taking the previous cluster head and its adjacent clusters as @xmath117 . we analyze the total transmit power and computational complexity of cwcl and the proposed dwcl algorithm . in order to do so , we first need to investigate the transmit power on a single link of length @xmath138 . the minimum transmit power is obtained from ( [ received power ] ) as @xmath139 , where @xmath140 is the minimum received power to correctly receive the message . for the computation complexity estimates , we assume the number of operations ( ops ) for addition , substraction and comparison is 1 ; for multiplication and division is 10 .      in cwcl , the total number of messages is the sum of reporting messages from nodes to the fusion center , thus , @xmath0 . for the transmit power consumption , we assume the fusion center is roughly at the center of a circle area of radius @xmath1 , and nodes are placed according to a uniform distribution . in this case , one can verify that , @xmath141 , the distance from the @xmath3th node to the fusion center satisfies @xmath142 , where @xmath143 is the uniform distribution from @xmath32 to @xmath37 . thus the average total power consumption is expressed as @xmath144 = \operatorname{\mathbb{e } } [ p_{r,\min } \sum_{i=1}^{n } ( d_i^{\gamma}10^{-s_i/10 } ) ] \approx n p_{r,\min } \operatorname{\mathbb{e}}[d_i^{\gamma } ] \operatorname{\mathbb{e}}[10^{-s_i/10}],\ ] ] where we use statistical average to replace sum of realizations .    note that @xmath145 = \operatorname{\mathbb{e}}[{(d_i^2)}^{\gamma/2}]$ ] , which is the @xmath146th moment of an uniform variable . in order to calculate @xmath147 $ ] , we need the distribution of @xmath148 , which is a function of central gaussian variable @xmath149 . thus we can derive the desired pdf as @xcite @xmath150    finally , the average total power consumption is given by @xmath151 = n p_{r,\min } \operatorname{\mathbb{e}}[d_i^{\gamma } ] \int_0^{\infty } x p(x ) d x.\ ] ]    from the definition of wcl estimate , number of ops of cwcl is @xmath152 . assume that : @xmath153 is the average number of nodes in one cluster ; @xmath154 is the average number of clusters ; @xmath155 is the number of neighboring clusters for each cluster ; @xmath156 is the average percentage of clusters who decide they are the head cluster in step 7 of algorithm 1 ( clusters that have to perform step 8) . the total number of messages for algorithm 1 is @xmath157 . to calculate transmit power , in step 1 ) we use the same procedure as cwcl ( with a smaller radius ) ; step 2 ) and 3 ) are communication among cluster heads , we use the average distance between centers of two clusters as the estimate of the average distance among cluster heads . in algorithm 2 , the total number of messages is @xmath158 and the calculation of power consumption for this process is the same as step 2 ) of algorithm 1 .    in terms of computational complexity , the number of ops for algorithm 1 is @xmath159 while the number of ops for algorithm 2 is @xmath160 . in this section , we evaluate the performance of cwcl algorithm and its proposed distributed implementation . for cwcl algorithm , using the theoretical framework presented in section [ sec : analysis ] , we consider impact of sensor node positioning error , different levels of shadowing @xmath161 and different correlation distances @xmath162 . the effect of node placement is investigated under random grid and uniform distribution . prior work on wcl uses mean localization error as the performance metric of interest . in @xcite , a slightly modified performance metric where the mean error is normalized by the transmission range @xmath1 of pu is adopted . simulations in these prior work have shown that the accuracy of wcl improves as the density of nodes is increased . this result is intuitive since the node spacing decreases as the density of nodes is increased . in this work , we introduce a novel performance metric where the mean localization error @xmath163 is normalized with respect to the average node spacing @xmath164 . the motivation for choosing this metric is to eliminate the effect of geometry scaling which intuitively improves accuracy .        in order to isolate the effect of shadowing on wcl performance we consider a fixed grid placement of nodes with perfect positioning ( i.e. @xmath165 ) and fixed pu position . the normalized mean localization error for this scenario is shown in fig .  [ fig : fixgridsquare_shadow_mean ] . for given shadowing variance , as the number of nodes increases , there is a slight improvement in the accuracy , which is a consequence of averaged shadowing over many nodes . note that this gain saturates at approximately 200 nodes or more . we also observe that wcl is quite robust to shadowing , since for an additional @xmath166 of shadowing from @xmath167 to @xmath168 , normalized error increases only 5% . therefore wcl is well suited for pu localization in severely shadowed environments . the results for mean localization error under correlated shadowing are shown in fig . [ fig : fixgridsquare_corr_mean ] . assuming fixed value of @xmath169 , and different levels of @xmath162 normalized to @xmath164 are presented . in highly correlated shadowing such as @xmath170 , increasing the number of nodes always worsens the error . for the high node density case , e.g. 400 nodes , the average error is as high as 47% of @xmath164 . we conclude that when all nodes are included in the wcl calculation , correlation has a degrading effect to accuracy . further investigation shows that the correlation distance resulting in the worst performance is proportional to the number of nodes involved .      the impact of positioning accuracy of sensing nodes is presented in fig . [ fig : imperfect location knowledge ] with a fixed grid placement of nodes in an i.i.d . shadowing environment with @xmath171 . analysis shows that errors in sensor node positioning decrease the wcl accuracy but the impact is marginal . for example , @xmath7 of @xmath172 m increases the error for about @xmath78% .    ) as a function of variance in the participating nodes positioning error , @xmath7 . for this scenario , nodes are placed on a grid with equal spacing ( fixed grid placement ) in an area with @xmath173 m . pu is at the center of the area . ] the previous subsections considered fixed node placement and pu location . realistically , node placement can be random and the pu also appears at random locations relative to the sensor nodes . next we consider two cases of random node placement : random grid and uniformly distributed placement .    in the random grid scenario , the sensor nodes are still on a grid , but the pu can appear uniformly within the center box . a grid distribution of nodes can be practically implemented if nodes are explicitly deployed for the purpose of acting as a dedicated localization sensor network . evenly placed beacons have been shown to provide better accuracy for wcl in wsns @xcite . results for random grid placement in fig . [ fig : randomgridanduniformrandomplacement ] show approximately 10% performance loss compared to fixed pu location scenarios . the uniformly distributed placement has been used in @xcite to more accurately model scenarios where the sensor node positions are not predefined , such as in mobile ad - hoc networks . the results for uniformly distributed placement are also shown in fig . [ fig : randomgridanduniformrandomplacement ] . it is observed that randomness in the node placement can increase the error as much as three times compared to the fixed node placement .    ) with varying correlation distance normalized to the node spacing . results of two placements , random grid ( rg ) and uniform random ( ur ) , are shown in an area with @xmath173 m . for both placements , @xmath174 . ]      the results given so far assumed an ideal circular coverage area for the pu . works such as @xcite indicated that in practical scenarios , the circular model can not be guarranteed due to several non - idealities such as fading , shadowing , and interference in the system . intuitively , the degree - of - irregularity ( _ doi _ ) parameter determines how closely the coverage area can be approximated by a circle . _ = 0 means the transmission range is perfectly circular . the model used in our simulation is the _ model briefly described in @xcite where the transmission range is a correlated gaussian r.v . centered at @xmath175 with variance @xmath176 .    in fig . [ fig : varydoi ] the effect of _ doi _ is investigated in the presence of uncorrelated medium level shadowing ( @xmath177 ) . the result shows that _ doi _ between 0% and 40% increases the error by 4% . _ doi _ is not a significant impairment for wcl because nodes closest to the pu , where the effect of irregularity is negligible , are weighted more than nodes closer to the edge of the coverage area . however , a slight increase in the normalized mean error can be observed as more nodes are used at high doi ( 30% or 40% ) . the denser the node distribution , the more nodes are included in wcl , but the more nodes there are , the more deviation from a perfect circular coverage coverage area , leading to lower accuracy .     with varying degree - of - irregularity ( doi ) . for this scenario , nodes are uniformly and independently distributed ( uniform random placement ) in an area with average @xmath175 . ]      from fig . [ fig : randomgridanduniformrandomplacement ] we observe that in correlated shadowing increasing the number of nodes also increases the error . one solution to this problem is to reduce the number of participating nodes used in wcl for this scenario . this could be achieved by ignoring nodes farther away from the pu to reduce the number of nodes resulting in less degradation .    in fig . [ fig : varyingparticipation ] we investigate varying levels of node participation for a fixed average node spacing , @xmath164 . in this case the total number of nodes is 100 , and the density remains the same while the number of involving nodes is varied . to select which nodes are included in the wcl calculation , nodes are sorted according to their rss values . the nodes with high received power are included in wcl . as a result , increasing the participation for the uncorrelated case always improves the accuracy . thus , 100% participation is optimal . however , even with a little correlation , this number is greatly reduced and about 10% of the total number of nodes is optimal . another selection scheme was previously proposed in @xcite where only nodes whose rss is within 15% of the stronges node s rss are used . the localization accuracy of wcl is also compared to other techniques suitable for cr scenarios . this includes lateration , centroid and strongest node ( sn ) . lateration is a range - based technique that uses least - squares estimation using distance estimates to the target . the centroid algorithm , on the other hand , is a range - free technique which calculates the average of all node positions in range of the pu . finally , the sn technique simply selects the node with the highest rss and uses its location as an estimate of the pu position . results in fig . [ fig : comaprison with other techiques ] show that wcl performs better under higher @xmath161 compared to lateration . although this seems counter - intuitive at first , the result can be explained by taking note that range - based techniques require an accurate estimate of the distance to the target to provide accurate results . in the case of non - interactive localization , where the only available information is the rss from the pu , shadowing induces a very large error in the range estimates . under low shadowing ( @xmath178 db ) , the lateration method is more accurate than wcl . at 0 db shadowing ( not shown in the figure ) , the range - based technique will achieve perfect localization . this result supports the claim that wcl is a viable method in the context of non - interactive localization because of its robustness . the centroid algorithm performs slightly worse than wcl because it does not utilize rss information . however , the performance difference becomes smaller at higher level of shadowing .    . for this scenario , nodes are uniformly and independently distributed ( uniform random placement ) in an area with @xmath175 ( uncorrelated shadowing ) . note that the centroid technique is unaffected by @xmath179 . ] [ fig : dwcl compare ] presents simulation results that compare normalized mean error of cwcl , sn and dwcl under different shadowing environments . sn is included in the comparison because it forms the basis of the head cluster selection for dwcl so it acts as an upperbound to the error . we observe that the proposed dwcl performance is between cwcl and sn algorithms . it is also interesting to note that the accuracy of dwcl degrades relative to cwcl and approaches the performance of sn when shadowing is high . the reason behind this is that for cluster head selection , dwcl is dependent on the sn algorithm . the accuracy of the sn degrades at large shadowing shown in fig . [ fig : comaprison with other techiques ] . [ fig : fixgrid_nodepower_theo ] compares the transmit power consumption for cwcl and dwcl with different number of clusters . in the simulation we set @xmath180 which is the mean received power at the border of the area of @xmath175 , assuming @xmath181 dbm , @xmath182 m and @xmath183 in ( [ received power ] ) . we observe that dwcl with 16 clusters consumes about @xmath184dbm less power for each node compared with cwcl . the analysis also shows that transmit power of each node is reduced as the number of node increases , since the average communication distance is getting smaller in dwcl . [ fig : fixgrid_cpu_theo ] compares computational complexity estimated in terms of number of ops . it shows that distributed algorithms require about three times more operations than cwcl , however total transmit power is reduced . in this paper we present the first analytical framework for the performance of wcl for localization of pu transmitters in cr networks . we derived the closed - form expression for the error distribution of wcl parameterized by node density , node placement , shadowing variance , correlation distance and errors in sensor nodes positions . using the theoretical results in conjunction with numerical simulations , the robustness of wcl against various physical conditions is investigated and quantified . based on the above analysis , we provide guidelines in terms of number of nodes and node placement required to achieve certain localization accuracy using wcl and propose variable node participation for improved wcl performance . to facilitate the deployment of wcl in an energy efficient context , a distributed cluster - based implementation of wcl is also proposed , which does not assume a fusion center that collects all rss information from nodes . the distributed method achieves comparable accuracy with its centralized counterpart , and greatly reduces total power consumption .        for notational convenience we define @xmath186 and form the vector @xmath187^{t}$ ] , so that @xmath188 with @xmath189 . note that this definition includes both i.i.d and correlated shadowing case . the first term in ( [ dimensional correlation ] ) then becomes @xmath190 = \operatorname{\mathbb{e}}\left[\frac{\mathbf{x}^{t } \mathbf{q } \mathbf{q}^{t } \mathbf{y}}{\mathbf{1}^{t } \mathbf{q } \mathbf{q}^{t }     \mathbf{1}}\right].\ ] ]    define three new variables @xmath191 , @xmath192 and @xmath193 , and their vector form @xmath194^{t}$ ] . now the first term in ( [ dimensional correlation ] ) becomes @xmath195 $ ] . as shown in sec . ii , @xmath196 , and @xmath197 . @xmath198 can be approximated as a single gaussian variable as shown in sec . [ subsubsec : iidshadowing ] since it is the same as @xmath32 in ( [ est x ] ) . we denote its distribution as @xmath199 , similarly @xmath200 . the distribution of @xmath201 is given as @xmath202 , where @xmath203 and @xmath204 . as a result @xmath205 , where @xmath206^t$ ] is the mean vector and @xmath207 is the covariance matrix . to derive @xmath207 we find @xmath208 as @xmath209 - m_1 m_2                           = \operatorname{\mathbb{e}}[(\sum_{i=1}^{n } x_i q_i)(\sum_{j=1}^{n } y_j q_j ) ] - m_1 m_2 = \overline{\textbf{x}}^{t } \mathbf{\omega } \overline{\textbf{y } } - m_1 m_2\ ] ] where in the last the equation we used the fact @xmath210 = \overline{x}_i \overline{y}_j \mathbf{\omega}_{ij}$ ] . similarly we calculate all entries of @xmath207 .    in order to simplify the calculation of ( [ reform ] ) , we perform the following transformation to @xmath211 , @xmath212 , thus @xmath213 . clearly @xmath214 , where @xmath215 is the @xmath3th column of @xmath216 . then we have @xmath217 = \operatorname{\mathbb{e}}[\frac{(m_1 + \boldsymbol{\omega}_{1}^{t } \widetilde{\textbf{t}})(m_2 + \boldsymbol{\omega}_{2}^{t } \widetilde{\textbf{t}})}{(m_3 + \boldsymbol{\omega}_{3}^{t } \widetilde{\textbf{t}})^2}].\ ] ]    since the expectation in ( [ results ] ) involves multiplication and division of three correlated random variables , the calculation can not be done directly . in order to reduce the number of variables involved , first we form @xmath218 $ ] , then perform qr factorization to reduce dimensionality , obtaining @xmath219 , in which @xmath86 is orthogonal and @xmath220 is upper triangular @xmath221.\ ] ]      note that since @xmath225 , denominator of ( [ results ] ) can further be expressed as @xmath226 \widetilde{\textbf{t}}\right)^{2 }      \doteq \left(m_3 + r_{11 } \widetilde{t}_{1}\right)^{2 } , \label{den2_new}\ ] ] where we use the fact that @xmath227 is the first column of @xmath220 ; and @xmath228^{t}$ ] . similarly , numerator of ( [ results ] ) can also be further expressed as      therefore , the expectation in ( [ reform ] ) becomes @xmath229     = \frac { r_{12 } r_{13 } } { r_{11}^2 } \operatorname{\mathbb{e}}\left[\frac{\widetilde{t}_{1}^2 + \left(\frac{m_1}{r_{12 } } + \frac{m_2}{r_{13}}\right ) \widetilde{t}_{1 } + \left(\frac{m_1 m_2 + r_{22 } r_{23}}{r_{12 } r_{13}}\right)}{\left(\widetilde{t}_{1 } + \frac{m_3}{r_{11}}\right)^{2}}\right ] . \label{exp1_new}\ ] ] after factorization ( [ exp1_new ] ) becomes @xmath230 = \frac { r_{12 } r_{13 } } { r_{11}^2 } \!\!\left(\!\!1 + \frac{\left(\frac{m_1}{r_{12 } } + \frac{m_2}{r_{13}}\right)-2 \frac{m_3}{r_{11}}}{\frac{m_3}{r_{11 } } } + \frac{\left(\frac{m_1 m_2 + r_{22 } r_{23}}{r_{12 } r_{13}}\right)+\left({\frac{m_3}{r_{11}}}\right)^2-\left(\frac{m_1}{r_{12 } } + \frac{m_2}{r_{13}}\right ) \frac{m_3}{r_{11}}}{1+\left({\frac{m_3}{r_{11}}}\right)^2 } \right).\ ] ] substituting ( [ corr - coeff_new ] ) into ( [ dimensional correlation ] ) gives the result . j.  wang , p.  urriza , y.  han , and d.  abri , `` performance analysis of weighted centroid algorithm for primary user localization in cognitive radio networks , '' in _ proc . asilomar conference on signals , systems , and computers _ , pacific grove , ca , usa , nov . 710 , 2010 . h.  chen , q.  shi , r.  tan , h.  poor , and k.  sezaki , `` mobile element assisted cooperative localization for wireless sensor networks with obstacles , '' _ ieee trans . wireless commun . _ , vol .  9 , no .  3 , pp . 956 963 , mar . 2010 s.  liu , y.  chen , w.  trappe , and l.  j. greenstein , `` non - interactive localization of cognitive radios based on dynamic signal strength mapping , '' in _ proc . conf . on wireless on - demand network systems and services _ , snowbird , ut , usa , feb . 24 , 2009 . r.  behnke , j.  salzmann , r.  grossmann , d.  lieckfeldt , d.  timmermann , and k.  thurow , `` strategies to overcome border area effects of coarse grained localization , '' in _ proc . 6th workshop on positioning , navigation and communication _ , hannover , germany , eu , mar . 19 , 2009 . q.  zhou , x.  li , and y.  xu , `` mean shift based collaborative localization with dynamically clustering for wireless sensor networks , '' in _ proc . conf . on communications and mobile computing _ , vol .  2 , yunnan , china , jan . 68 , 2009 . m.  k. simon , _ probability distributions involving gaussian random variables : a handbook for engineers , scientists and mathematicians_. 1em plus 0.5em minus 0.4emnew york , ny , usa : springer - verlag , 2006 . r.  chandra , j.  padhye , a.  wolman , and b.  zill , `` a location - based management system for enterprise wireless lans , '' in _ proc . symposium on networked systems design and implementation _ , cambridge , ma , usa , apr . 1113 , 2007 .
Abstract:

Answer: information about primary transmitter location is crucial in enabling several key capabilities in cognitive radio networks , including improved spatio - temporal sensing , intelligent location - aware routing , as well as aiding spectrum policy enforcement . compared to other proposed non - interactive localization algorithms , the weighted centroid localization ( wcl ) scheme uses only the received signal strength information , which makes it simple to implement and robust to variations in the propagation environment . in this paper we present the first theoretical framework for wcl performance analysis in terms of its localization error distribution parameterized by node density , node placement , shadowing variance , correlation distance and inaccuracy of sensor node positioning . using this analysis , we quantify the robustness of wcl to various physical conditions and provide design guidelines , such as node placement and spacing , for the practical deployment of wcl . we also propose a power - efficient method for implementing wcl through a distributed cluster - based algorithm , that achieves comparable accuracy with its centralized counterpart .