A noise resilient distributed inference framework for body sensor networks
Common Phase Error due to Phase Noise in OFDM - Estimation and Suppression

COMMON PHASE ERROR DUE TO PHASE NOISE IN OFDM-ESTIMATION AND SUPPRESSIONDenis Petrovic,Wolfgang Rave and Gerhard FettweisV odafone Chair for Mobile Communications,Dresden University of Technology,Helmholtzstrasse18,Dresden,Germany{petrovic,rave,fettweis}@ifn.et.tu-dresden.deAbstract-Orthogonal frequency division multiplexing (OFDM)has already become a very attractive modulation scheme for many applications.Unfortunately OFDM is very sensitive to synchronization errors,one of them being phase noise,which is of great importance in modern WLAN systems which target high data rates and tend to use higher frequency bands because of the spectrum availability.In this paper we propose a linear Kalmanfilter as a means for tracking phase noise and its suppression.The algorithm is pilot based.The performance of the proposed method is investigated and compared with the performance of other known algorithms.Keywords-OFDM,Synchronization,Phase noise,WLANI.I NTRODUCTIONOFDM has been applied in a variety of digital commu-nications applications.It has been deployed in both wired systems(xDSL)and wireless LANs(IEEE802.11a).This is mainly due to the robustness to frequency selective fading. The basic principle of OFDM is to split a high data rate data stream into a number of lower rate streams which are transmitted simultaneously over a number of orthogonal subcarriers.However this most valuable feature,namely orthogonality between the carriers,is threatened by the presence of phase noise in oscillators.This is especially the case,if bandwidth efficient higher order modulations need to be employed or if the spacing between the carriers is to be reduced.To compensate for phase noise several methods have been proposed.These can be divided into time domain[1][2]and frequency domain approaches[3][4][5].In this paper we propose an algorithm for tracking the average phase noise offset also known as the common phase error(CPE)[6]in the frequency domain using a linear Kalmanfilter.Note that CPE estimation should be considered as afirst step within more sophisticated algorithms for phase noise suppression[5] which attempt to suppress also the intercarrier interference (ICI)due to phase noise.CPE compensation only,can however suffice for some system design scenarios to suppress phase noise to a satisfactory level.For these two reasons we consider CPE estimation as an important step for phase noise suppression.II.S YSTEM M ODELAn OFDM transmission system in the presence of phase noise is shown in Fig. 1.Since all phase noise sources can be mapped to the receiver side[7]we assume,without loss of generality that phase noise is present only at the front end of the receiver.Assuming perfect frequency and timing synchronization the received OFDM signal samples, sampled at frequency f s,in the presence of phase noise can be expressed as r(n)=(x(n) h(n))e jφ(n)+ξ(n).Each OFDM symbol is assumed to consist of a cyclic prefix of length N CP samples and N samples corresponding to the useful signal.The variables x(n),h(n)andφ(n)denote the samples of the transmitted signal,the channel impulse response and the phase noise process at the output of the mixer,respectively.The symbol stands for convolution. The termξ(n)represents AWGN noise with varianceσ2n. The phase noise processφ(t)is modelled as a Wiener process[8],the details of which are given below,with a certain3dB bandwidth∆f3dB.,0,1,2...m lX l=,0,1,2...m lR l=Fig.1Block diagram of an OFDM transmission chain.At the receiver after removing the N CP samples cor-responding to the cyclic prefix and taking the discrete Fourier transform(DFT)on the remaining N samples,the demodulated carrier amplitude R m,lkat subcarrier l k(l k= 0,1,...N−1)of the m th OFDM symbol is given as[4]:R m,lk=X m,lkH m,lkI m(0)+ζm,lk+ηm,lk(1)where X m,lk,H m,lkandηm,lkrepresent the transmitted symbol on subcarrier l k,the channel transfer function andlinearly transformed AWGN with unchanged variance σ2n at subcarrier l k ,respectively.The term ζm,l k represents intercarrier interference (ICI)due to phase noise and was shown to be a gaussian distributed,zero mean,randomvariable with variance σ2ICI =πN ∆f 3dB s[7].The term I m (0)also stems from phase noise.It does not depend on the subcarrier index and modifies all subcarriers of one OFDM symbol in the same manner.As its modulus is in addition very close to one [9],it can be seen as a symbol rotation in the complex plane.Thus it is referred to in the literature as the common phase error (CPE)[6].The constellation rotation due to CPE causes unaccept-able system performance [7].Acceptable performance can be achieved if one estimates I m (0)or its argument and compensates the effect of the CPE by derotating the received subcarrier symbols in the frequency domain (see Eq.(1)),which significantly reduces the error rate as compared to the case where no compensation is used.The problem of esti-mating the CPE was addressed by several authors [3][4][10].In [3]the authors concentrated on estimating the argument of I m (0)using a simple averaging over pilots.In [10]the argument of I m (0)was estimated using an extended Kalman filter,while in [4]the coefficient I m (0)itself was estimated using the LS algorithm.Here we introduce an alternative way for minimum mean square estimation (MMSE)[11]of I m (0)using a linear scalar Kalman filter.The algorithm is as [4]pilot based.III.P HASE N OISE M ODELFor our purposes we need to consider a discretized phase noise model φ(n )=φ(nT s )where n ∈N 0and T s =1/f s is the sampling period at the front end of the receiver.We adopt a Brownian motion model of the phase noise [8].The samples of the phase noise process are given as φ(n )=2πf c √cB (n )where f c is the carrier frequency,c =∆f 3dB /πf 2c [8]and B (n )represents the discretizied Brownian motion process,Using properties of the Brownian motion [12]the fol-lowing holds:B (0)=0and B (n +1)=B (n )+dB n ,n ∈N 0where each increment dB n is an independent random variable and dB n ∼√T s N (0,1).Noting that φ(n )=2πf c √cB (n )we can write the discrete time phase noise process equation asφ(n +1)=φ(n )+w (n )(2)where w (n )∼N (0,4π2f 2c cT s )is a gaussian randomvariable with zero mean and variance σ2w =4π2f 2c cT s .IV.CPE E STIMATION U SING A K ALMAN F ILTER Since all received subcarriers within one OFDM symbolare affected by the same factor,namely I m (0),the problem at hand can be seen as an example of estimating a constant from several noisy measurements given by Eq.(1)for which purpose a Kalman filter is well suited [11].For a Kalmanfilter to be used we need to define the state space model of the system.Define first the set L ={l 1,l 2,l 3,...l P }as a subset of the subcarrier set {0,1,...N −1}.Using Eq.(1)one can writeR m,l k =A m,l k I m,l k (0)+εm,l k(3)where A m,l k =X m,l k H m,l k and I m,l k (0)=I m (0)for all k =1,2...,P .Additional indexing of the CPE terms is done here only for convenience of notation.On the other hand one can writeI m,l k +1(0)=I m,l k (0).(4)Equations (3)and (4)are the measurement and processequation of the system state space model,where A m,l k represents the measurement matrix,while the process matrix is equal to 1and I m,l k (0)corresponds to the state of the system.The measuring noise is given by εm,l k which combines the ICI and AWGN terms in Eq.(1),the varianceof which for all l k equals σ2ε=(σ2ICI +σ2n ).The process noise equals zero.Note that the defined state space model is valid only for one OFDM symbol.For the state space model to be fully defined,knowledge of the A m,l k =X m,l k H m,l k is needed.Here we assume to have ideal knowledge of the channel.On the other hand we define the subset L to correspond to the pilot subcarrier locations within one OFDM symbol so that X m,q ,q ∈L are also known.We assume that at the beginning of each burst perfect timing and frequency synchronization is achieved,so that the phase error at the beginning of the burst equals zero.After the burst reception and demodulation,the demodulated symbols are one by one passed to the Kalman filter.For a Kalman filter initialization one needs for eachOFDM symbol an a priori value for ˆI m,l 1(0)and an a priori error variance K −m,1.At the beginning of the burst,when m =1,it is reasonable to adopt ˆI −1,l 1(0)=1.Within each OFDM symbol,say m th,the filter uses P received pilot subcarriers to recursively update the a priori value ˆI −1m,l 1(0).After all P pilot subcarriers are taken into account ˆI m,l P (0)is obtained,which is adopted as an estimate ofthe CPE within one OFDM symbol,denoted as ˆIm (0).The Kalman filter also provides an error variance of the estimateof I m,l P (0)as K m,P .ˆI m,l P(0)and K m,P are then used as a priori measures for the next OFDM symbol.The detailed structure of the algorithm is as follows.Step 1:InitializationˆI −m,l 1(0)=E {I −m,l 1(0)}=ˆI m −1(0)K −m,1=E {|I m (0)−ˆIm −1(0)|2}∼=E {|φm −ˆφm −1|2}=σ2CP E +K m −1,Pwhere σ2CP E =4π2N 2+13N +N CP ∆f 3dBf s(see [10]),K 0,P =0and φm =arg {I m (0)}.Repeat Step2and Step3for k=1,2,...,P Step2:a-posteriori estimation(update)G m,k=K−m,kH H m,lkH m,lkK−m,kH Hm,l k+(σ2ICI+σ2n)ˆIm,l k (0)=ˆI−m,l k(0)+G m,k[R m,lk−H m,l kˆI−m,l k(0)]K m,k=(1−G m,k H m,lk )K−m,kStep3:State and error variance propagationK−m,k+1=K m,k(5)ˆI−m,l k+1(0)=ˆI m,lk(0)Note that no matrix inversions are required,since the state space model is purely scalar.V.CPE C ORRECTIONThe easiest approach for CPE correction is to derotate all subcarriers l k of the received m th symbol R m,lkby φm=−arg{ˆI m(0)}.Unambiguity of the arg{·}function plays here no role since any unambiguity which is a multiple of2πrotates the constellation to its equivalent position in terms of its argument.The presented Kalmanfilter estimation algorithm is read-ily applicable for the decision feedback(DF)type of algo-rithm presented in[4].The idea there was to use the data symbols demodulated after thefirst CPE correction in a DFE manner to improve the quality of the estimate since that is increasing the number of observations of the quantity we want to estimate.In our case that would mean that after thefirst CPE correction the set L={l1,l2,l3,...l P}of the subcarriers used for CPE estimation,which previously corresponded to pilot subcarriers,is now extended to a larger set corresponding to all or some of the demodulated symbols. In this paper we have extended the set to all demodulated symbols.The Kalmanfilter estimation is then applied in an unchanged form for a larger set L.VI.N UMERICAL R ESULTSThe performance of the proposed algorithm is investigated and compared with the proposal of[4]which is shown to outperform other known approaches.The system model is according to the IEEE802.11a standard,where64-QAM modulation is used.We investigate the performance in AWGN channels and frequency selective channels using as an example the ETSI HiperLAN A-Channel(ETSI A). Transmission of10OFDM symbols per burst is assumed.A.Properties of an EstimatorThe quality of an estimation is investigated in terms of the mean square error(MSE)of the estimator for a range of phase noise bandwidths∆f3dB∈[10÷800]Hz.Table1 can be used to relate the phase noise bandwidth with other quantities.Figures2and3compare the MSE of the LS estimator from[4]and our approach for two channel types and both standard correction and using decision feedback. Note that SNRs are chosen such that the BER of a coded system after the Viterbi algorithm in case of phase noise free transmission is around1·10−4.Kalmanfilter shows better performance in all cases and seems to be more effective for small phase noise bandwidths. As expected when DF is used the MSE of an estimator is smaller because we are taking more measurements into account.Fig.2MSE of an estimator for AWGN channel.Fig.3MSE of an estimator for ETSI A channel.Table 1Useful relationsQuantitySymbolRelationTypical values for IEEE802.11aOscillator constant c [1radHz]8.2·10−19÷4.7·10−18Oscillator 3dB bandwidth ∆f 3dB [Hz]∆f 3dB =πf 2cc 70÷400Relative 3dB bandwidth ∆f 3dB ∆f car∆f 3dBfsN 2·10−4÷13·10−4Phase noise energy E PN [rad]E PN =4π∆f 3dB∆fcar0.0028÷0.016Subcarrier spacing∆f car∆f car =f s N312500HzB.Symbol Error Rate DegradationSymbol error rate (SER)degradation due to phase noise is investigated also for a range of phase noise bandwidths ∆f 3dB ∈[10÷800]Hz and compared for different correc-tion algorithms.Ideal CPE correction corresponds to the case when genie CPE values are available.In all cases simpleconstellation derotation with φ=−arg {ˆIm (0)}is used.Fig.4SER degradation for AWGN channel.In Figs.4and 5SER degradation for AWGN and ETSI A channels is plotted,respectively.It is interesting to note that as opposed to the ETSI A channel case in AWGN channel there is a gap between the ideal CPE and both correction approaches.This can be explained if we go back to Eq.(1)where we have seen that phase noise affects the constellation as additive noise.Estimation error of phase noise affects the constellation also in an additive manner.On the other hand the SER curve without phase noise in the AWGN case is much steeper than the corresponding one for the ETSI A channel.A small SNR degradation due to estimation errors will cause therefore large SER variations.This explains why the performance differs much less in the ETSI A channel case.Generally from this discussion a conclusion can be drawn that systems with large order of diversity are more sensitive to CPE estimation errors.Note that this ismeantFig.5SER degradation for ETSI A channel.not in terms of frequency diversity but the SER vs.SNR having closely exponential dependence.It can be seen that our approach shows slightly better performance than [4]especially for small phase noise bandwidths.What is also interesting to note is,that DF is not necessary in the case of ETSI A types of channels (small slope of SER vs.SNR)while in case of AWGN (large slope)it brings performance improvement.VII.C ONCLUSIONSWe investigated the application of a linear Kalman filter as a means for tracking phase noise and its suppression.The proposed algorithm is of low complexity and its performance was studied in terms of the mean square error (MSE)of an estimator and SER degradation.The performance of an algorithm is compared with other algorithms showing equivalent and in some cases better performance.R EFERENCES[1]R.A.Casas,S.Biracree,and A.Youtz,“Time DomainPhase Noise Correction for OFDM Signals,”IEEE Trans.on Broadcasting ,vol.48,no.3,2002.[2]M.S.El-Tanany,Y.Wu,and L.Hazy,“Analytical Mod-eling and Simulation of Phase Noise Interference in OFDM-based Digital Television Terrestial Broadcast-ing Systems,”IEEE Trans.on Broadcasting,vol.47, no.3,2001.[3]P.Robertson and S.Kaiser,“Analysis of the effects ofphase noise in OFDM systems,”in Proc.ICC,1995.[4]S.Wu and Y.Bar-Ness,“A Phase Noise SuppressionAlgorithm for OFDM-Based WLANs,”IEEE Commu-nications Letters,vol.44,May1998.[5]D.Petrovic,W.Rave,and G.Fettweis,“Phase NoiseSuppression in OFDM including Intercarrier Interfer-ence,”in Proc.Intl.OFDM Workshop(InOWo)03, pp.219–224,2003.[6]A.Armada,“Understanding the Effects of PhaseNoise in Orthogonal Frequency Division Multiplexing (OFDM),”IEEE Trans.on Broadcasting,vol.47,no.2, 2001.[7]E.Costa and S.Pupolin,“M-QAM-OFDM SystemPerformance in the Presence of a Nonlinear Amplifier and Phase Noise,”IEEE mun.,vol.50, no.3,2002.[8]A.Demir,A.Mehrotra,and J.Roychowdhury,“PhaseNoise in Oscillators:A Unifying Theory and Numerical Methods for Characterisation,”IEEE Trans.Circuits Syst.I,vol.47,May2000.[9]S.Wu and Y.Bar-ness,“Performance Analysis of theEffect of Phase Noise in OFDM Systems,”in IEEE 7th ISSSTA,2002.[10]D.Petrovic,W.Rave,and G.Fettweis,“Phase NoiseSuppression in OFDM using a Kalman Filter,”in Proc.WPMC,2003.[11]S.M.Kay,Fundamentals of Statistical Signal Process-ing vol.1.Prentice-Hall,1998.[12]D.J.Higham,“An Algorithmic Introduction to Numer-ical Simulation of Stochastic Differential Equations,”SIAM Review,vol.43,no.3,pp.525–546,2001.。
Interplay of friction and noise and enhancement of disoriented chiral condensate

a r X i v :h e p -p h /0011003v 1 1 N o v 2000Interplay of friction and noise and enhancement of disoriented chiral condensateA.K.Chaudhuri ∗Variable Energy Cyclotron Centre 1-AF,Bidhan Nagar,Calcutta-700064Using the Langevin equation for the linear σmodel,we have investigated the effect of friction and noise on the possible disoriented chiral condensate formation.Friction and noise are supposed to suppress longwavelength oscillations and growth of disoriented chiral condensate domains.Details simulation shows that for heavy ion collisions,interplay of friction and noise occur in such a manner that formation of disoriented chiral condensate domains are enhanced.25.75.+r,12.38.Mh,11.30.RdIn recent years there is much excitement about the possibility of formation of disoriented chiral condensate (DCC).The idea was suggested by Rajagopal and Wilczek [1]and also by Bjorken and others [2].In hadron-hadron or in heavy ion collisions,a macroscopic region of space-time may be created within which the chiral order parameter is not oriented in the same direction in the internal O (4)=SU (2)×SU (2)space as the ordinary vacuum.The misaligned condensate has the same quark content and quantum numbers as do pions and essentially constitute a classical pion field.The system will finally relaxes to the true vacuum and in the process can emit coherent pions.Possibility of producing classical pion fields in heavy ion collisions had been discussed earlier by Anslem [3].DCC formation in hadronic or in heavy ion collisions can lead to the spectacular events that some portion of the detector will be dominated by charged pions or by neutral pions only.In contrast,in a general event all the three pions (π+,π−and π0)will be equally well produced.This may then be the natural explanation of the so called Centauro events [4].Microscopic physics governing DCC phenomena is not well known.It is in the regime of non-perturbative QCD as well as nonlinear phenomena,theoretical understanding of both of which are limited.One thus uses some effective field theory like linear σmodel with various approximations to simulate the chiral phase transition [5–7].In the linear sigma model chiral degrees of freedom are described by the the real O(4)field Φ=(σ,π→).Because of the isomorphism between the groups O (4)and SU (2)×SU (2),the later being the appropriate group for two flavour QCD,linear sigma model can effectively model the low energy dynamic of QCD [8].Explicit simulation with linear sigma model,indicate that DCC depends critically on the initial field configurations.With quench like initial condition DCC domains of 4-5fm in size can form [7].Initial conditions other than quench lead to much smaller domain size.Quench scenario assume that the effective potential governing the evolution of long wave length modes immediately after the phase transition at T c turns to classical one at zero temperature.It can happen only in case of very rapid cooling and expansion of the fireball.In heavy ion collisions quench like initial conditions are unlikely.Very recently,effect of external media on possible DCC is being investigated [9–16].Indeed,in heavy ion collisions,even if some region is created where chiral symmetry is restored,that region will be continually interacting with surrounding medium (mostly pions).The surrounding medium can be conveniently modelled by a white noise source and one can use Langevin equation for linear σmodel to simulate the DCC formation under the influence of external media.Recently it has been shown that in φ4model,hard modes can be integrated out to obtain a Langevin type of equations for the soft modes [9].Biro and Greiner [10]using a Langevin equation for the linear σmodel,investigated the interplay of friction and white noise on the evolution and stability of zero mode pion fields.In general friction and noise reduces the amplification of zero modes.But in some trajectories,large amplification may occur [10].We also obtain similar results [13].It is popular expectation that friction and noise will reduce the DCC formation probability.The expectation was found to be true in one dimensional calculations with zero modes only [10,13].To see how far this expectation is valid when higher modes are included,in the present letter we have investigated DCC formation in 3+1dimension.Two scenarios were considered.In scenario I,we use the equation of motion of linear σmodel with quenched initial condition to simulate DCC formation without friction and noise.In scenario II,we solve the Langevin equation for linear σmodel,with the same quenched initial condition to simulate DCC formation in presence of friction and parative study of these two scenarios reveal that certainly the noise and friction affect DCC formation probability.However contrary to one dimensional calculations,large DCC formation is seen in scenario II rather than in scenario I.Indeed,it will be shown that large DCC domains are formed in scenario II only,not in scenario I.This effect is particular to heavy ion collisions as explained below.The Langevin equation for linear σmodel,at temperature T can be written as,[∂2τ+ηi )∂∂x 2−∂2τ2∂2ion collisions.ησ,πandζσ,πare the friction coefficients and the white noise for theσandπfields.We note that if the friction and the noise terms are dropped from eq.1the resulting equation is for the scenario I,i.e.just the equation of motion of linearσmodel.The noise sourceζand the frictionηare not independent.They are related by thefluctuations-dissipation relation. We use white noise source with zero average and correlation as demanded from thefluctuation-dissipation relation,<ζa(τ,x,y,Y)>=0(2a) <ζa(τ,x,y,Y)ζb(τ′,x′,y′,Y′)>=2Tη1N )2m2σ1−4m2π1−exp(−m2σ/2mπT)1N )2N−11−4m2π4T(3b)Solution of eq.1require initialfields configurations.They were distributed according to a random Gaussian with,<σ>=(1−f(r))fπ(4a)<πi>=0(4b)<σ2>−<σ>2=<π2i>−<πi>2=f2π/6f(r)(4c)<˙σ>=<˙πi>=0(4d)<˙σ2>=<˙π2>=4f2π/6f(r)(4e)The interpolation functionf(r)=[1+exp(r−r0)/Γ)]−1(5) separates the central region from the rest of the system.We have taken r0=6.4fm andΓ=.7fm.The initialfield configurations corresponds to quench like condition[1]and as told in the beginning are unlikely to be obtained in a heavy ion collisions.We still use it as they are the most favourable initial conditions to produce DCC like phenomena. In the simulation results presented below,the initial time and temperature were assumed to beτi=1fm/c and T i=200 MeV.Effect of expansion was included through the cooling law,T(τ)=T iτii,j|π(i)||π(j)|(7)where the sum is taken over those grid points i and j such that the distance between i and j is r.Infig.1,we have compared the correlation function in scenario I and II.Initially there is no correlation length beyond the lattice spacing of1fm.Correlations starts to develop at later time.It increases for about7fm/c,then decreases again. Interestingly,larger correlation length is obtained in the scenario II,than in scenario I.Thus at7fm/c,correlation length in scenario I is only∼2fm,while that in scenario II is∼4fm.Increased correlation in scenario II is contrary to popular expectation that friction and noise will reduce correlation.Enhancement of correlation length in scenario II,with friction and noise is corroborated in the pionfield distribution also.Infig.2,we have shown the xy contour plot of theπ2component,at rapidity Y=0.Field distribution atτi=1 fm/c and after7fm/c of evolution are shown.Initially there is no correlation.Domain like structure is seen both in scenario I and II,after7fm/c evolution.The positive and negative components of theπ2separates out.Here again, much larger domains are formed in scenario II than in scenario I.It may be noted that domain like structure seen in one of the component ofπfield do not necessarily convert into physical domains.They are indication of larger correlation length only.Physical domain should contain either charged or neutral pion only.Thus in physical domain neutral to total pion ratio should differ considerably from the isospin symmetric value of1/3.Assuming that the pion density is proportional to the amplitude’s square,infig.3,we have shown the contour plot of the neutral to total pion ratio,at rapidity Y=0.π23R3(x,y,Y,τ)== π23dxdyNπtotalin one component of theπfield,there is no(large)physical domain formation in scenario I.On the contrary,large physical domains are seen to be formed in scenario II after7fm/c of evolution.We have also studied the rapidity distribution of the neutral to pion ratio in both the scenarios.Throughout the rapidity range,the ratiofluctuate about the isospin symmetric value of1/3,fluctuations being more in scenario II than in scenario I.However,it was also noted thatfluctuations will be considerably less when integrated over time.Rapidity distribution of neutral to total pion ratio may not be able to single out the DCC events.for7fm/c.shown.。
Teleportation and spin squeezing utilizing multimode entanglement of light with atoms

a r X i v :q u a n t -p h /0511174v 1 17 N o v 2005Teleportation and spin squeezing utilizing multimode entanglement of light withatomsK.Hammerer 1,E.S.Polzik 2,3,J.I.Cirac 11Max-Planck–Institut f¨u r Quantenoptik,Hans-Kopfermann-Strasse,D-85748Garching,Germany2QUANTOP,Danish Research Foundation Center for Quantum Optics,DK 2100Copenhagen,Denmark3Niels Bohr Institute,DK 2100Copenhagen,Denmark We present a protocol for the teleportation of the quantum state of a pulse of light onto the collective spin state of an atomic ensemble.The entangled state of light and atoms employed as a resource in this protocol is created by probing the collective atomic spin,Larmor precessing in an external magnetic field,offresonantly with a coherent pulse of light.We take here for the first time full account of the effects of Larmor precession and show that it gives rise to a qualitatively new type of multimode entangled state of light and atoms.The protocol is shown to be robust against the dominating sources of noise and can be implemented with an atomic ensemble at room temperature interacting with free space light.We also provide a scheme to perform the readout of the Larmor precessing spin state enabling the verification of successful teleportation as well as the creation of spin squeezing.PACS numbers:03.67.Mn,32.80.QkI.INTRODUCTIONQuantum teleportation -the disembodied transport of quantum states -has been demonstrated so far in sev-eral seminal experiments dealing with purely photonic [1]or atomic [2]systems.Here we propose a protocol for the teleportation of a coherent state carried initially by a pulse of light onto the collective spin state of ∼1011atoms.This protocol -just as the recently demonstrated direct transfer of a quantum state of light onto atoms [3]-is particularly relevant for long distance entanglement distribution,a key resource in quantum communication networks [4].Our scheme can be implemented with just coherent light and room-temperature atoms in a single vapor cell placed in a homogeneous magnetic field.Existing proto-cols in Quantum Information (QI)with continuous vari-ables of atomic ensembles and light [4]are commonly de-signed for setups where no external magnetic field is ap-plied such that the interaction of light with atoms meets the Quantum non-demolition (QND)criteria [5,6].In contrast,in all experiments dealing with vapor cells at room-temperature [3,7]it is,for technical reasons,ab-solutely essential to employ magnetic fields.In experi-ments [3,7]two cells with counter-rotating atomic spins were used to comply with both,the need for an exter-nal magnetic field and the one for an interaction of QND character.So far it was believed to be impossible to use a single cell in a magnetic field to implement QI protocols,since in this case -due to the Larmor precession -scat-tered light simultaneously reads out two non-commuting spin components such that the interaction is not of QND type.In this paper we do not only show that it is well possible to make use of the quantum state of light and atoms created in this setup but we demonstrate that -for the purpose of teleportation [8,9]-it is in factbetter to do so.As compared to the state resulting from the common QND interaction the application of an external magnetic field enhances the creation of correlations between atoms and light,generating more and qualitatively new,multimode type of entanglement.The results of the paper can be summarized as follows:(i )Larmor precession in an external magnetic field enhances the creation of entanglement when a collective atomic spin is probed with off-resonant light.The resulting entanglement involves multiple modes and is stronger as compared to what can be achieved in a comparable QND interaction.(ii )This type of entangled state can be used as a resource in a teleportation protocol,which is a simple generalization of the standard protocol [8,9]based on Einstein-Podolsky-Rosen (EPR)type of entanglement.For the experimentally accessible parameter regime the teleportation fidelity is close to optimal.The protocol is robust against imperfections and can be implemented with state of the art technique.(iii )Homodyne detection of appropriate scattering modes of light leaves the atomic state in a spin squeezed state.The squeezing can be the same as attained from a comparable QND measurement of the atomic spin [10,11].The same scheme can be used for atomic state read-out of the Larmor precessing spin,necessary to verify successful teleportation.We would like to note that it was shown recently in [12]that the effect of a magnetic field can enhance the capacity of a quantum memory in the setup of two cells.Teleportation in the setup of a single cell without mag-netic field was addressed in [13].The paper is organized as follows:The three points above are presented in sections II,III and IV,in this order.Some of the details in the calculations of sections III and IV are moved to appendices B and C.2II.INTERACTIONWe consider an ensemble of N at Alkali atoms with total ground state angular momentum F ,placed in a constant magnetic field causing a Zeeman splitting of Ωand ini-tially prepared in a fully polarized state along x .The collective spin of the ensemble is then probed by an offresonant pulse which propagates along z and is linearly polarized along x .Thorough descriptions of this inter-action and the final state of light and atoms after the scattering can befoundin [14,15,16,17,18]and espe-cially in [19,20,21,22]for the specific system we have in mind.We derive the final state here with a special focus on the effects of Larmor precession and light prop-agation in order to identify the light modes which are actually populated in the scattering process.In appendix A we show that the interaction is ade-quately described by a HamiltonianH =H at +H li +V,H at =Ω√N ph N at F a 1σΓ/2A ∆where N ph is the over-all number of photons in the pulse,a 1is a constant characterizing the ground state’s vector polarizability,σis the scattering cross section,Γthe decay rate,∆the detuning and A the effective beam cross section.Changing to a rotating frame with respect to H at by defining X I (t )=exp(−iH at t )X exp(iH at t )and evaluat-ing the Heisenberg equations for these operators yields the following Maxwell-Bloch equations∂t X I (t )=κT cos(Ωt )p (0,t ),(2a)∂t P I (t )=κT sin(Ωt )p (0,t ),(2b)(∂t +c∂z )x (z,t )=κcT[cos(Ωt )P I (t )−sin(Ωt )X I (t )]δ(z ),(∂t +c∂z )p (z,t )=0,where ∂t (z )denotes the partial derivative with respect to t (z ).These equations have a clear interpretation.Light noise coming from the field in quadrature with the classi-cal probe piles up in both,the X and P spin quadrature,but it alternately affects only one or the other,changing with a period of 1/Ω.Conversely atomic noise adds to the in phase field quadrature only and the signal comes alternately from the X and P spin quadrature.The out of phase field quadrature is conserved in the interaction.To solve this set of coupled equations it is convenient to introduce a new position variable,ξ=ct −z ,to elim-inate the z dependence.New light quadratures defined by ¯x (ξ,t )=x (ct −ξ,t ),¯p (ξ,t )=p (ct −ξ,t )also have a simple interpretation:ξlabels the slices of the pulse moving in and out of the ensemble one after the other,starting with ξ=0and terminating at ξ=cT .The Maxwell equations now read∂t ¯p (ξ,t )=0,(2c)∂t ¯x (ξ,t )=κcT[cos(Ωt )P I (t )−sin(Ωt )X I (t )]δ(ct −ξ).(2d)The solutions to equations (2a,2b,2c)areX I (t )=X I (0)+κT t 0d τcos(Ωτ)¯p (cτ,0),(3a)P I (t )=P I (0)+κT t 0d τsin(Ωτ)¯p (cτ,0),(3b)¯p (ξ,t )=¯p (ξ,0)(3c)and the formal solution to (2d)is¯x (ξ,t )=¯x (ξ,0)+(3d)+κT [cos(Ωξ/c )P I (ξ/c )−sin(Ωξ/c )X I (ξ/c )].As mentioned before,both atomic spin quadratures are affected by light but,as is evident from the solutions for X (t ),P (t ),they receive contributions from different and,in fact,orthogonal projections of the out-of-phase field.As we will show in the following,the corresponding projections of the in-phase field carry in turn the signal of atomic quadratures after the interaction.It is there-fore convenient to explicitly introduce operators for these modes [21].We define a cosine component before the in-teractionp in c = T Td τcos(Ωτ)¯p (cτ,0),(4a)x in c =TTd τcos(Ωτ)¯x (cτ,0)(4b)and a sine component p in s ,x ins with cos(Ωτ)replaced by sin(Ωτ).In frequency space these modes consist of spec-tral components at sidebands ωc ±Ωand are closely re-lated to the sideband modulation modes introduced in3[26]for the description of two photon processes.It is eas-ily checked that these modes are asymptotically canon-ical,[x in c ,p in c ]=[x in s ,p ins ]=i [1+O (n −10)]≃i ,and inde-pendent,[x in c ,p in s ]=O (n −10)≃0,if we assume n 0≫1for n 0=ΩT ,the pulse length measured in periods of Larmor precession.In terms of these modes the atomic state after the in-teraction X out =X I (T ),P out =P I (T )is given by X out =X in +κ2p inc ,P out =P in +κ2p ins .(5a)The final state of cosine (sine)modes is de-scribed by x out c(s),p outc(s),defined by equations (4)with ¯x (cτ,0),¯p (cτ,0)replaced by ¯x (cτ,T ),¯p (cτ,T )respec-tively.Since the out-of-phase field is conserved we have triviallyp out c =p in c ,p out s =p in s .(5b)Deriving the corresponding expressions for the cosine andsine components of the field in phase,x out c ,x outs ,raises some difficulties connected to the back action of light onto itself.This effect can be understood by noting that a slice ξof the pulse receives a signal of atoms at a time ξ/c [see equation (3d)]which,regarding equations (3a,3b),in turn carry already the integrated signal of all slices up to ξ.Thus,mediated by the atoms,light acts back on itself.The technicalities in the treatment of this effect are given in appendix B where we identify relevant ”back action modes”,x c ,1,p c ,1,x s ,1,p s ,1,in terms of which one can express x out c =x in c +κ2P in +κ√2 2p in s ,1,(5c)x out s=x in s −κ2X in−κ√2 2p in c ,1.(5d)The last two terms in both lines represent the effect of back action,part of which involves the already defined cosine and sine components of the field in quadrature.The remaining part is subsumed in the back action modes which are again canonical and independent from all other modes.Equations (5)describe the final state of atoms and the relevant part of scattered light after the pulse has passed the atomic ensemble and are the central result of this section.Treating the last terms in equations (5c,5d)as noise terms,it is readily checked by means of the separa-bility criteria in [27]that this state is fully inseparable,i.e.it is inseparable with respect to all splittings be-tween the three modes.For the following teleportation protocol the relevant entanglement is the one between atoms and the two light modes.Figure 1shows the von Neumann entropy E vN of the reduced state of atoms in its dependence on the coupling strength κand in com-parison with the entanglement created without magnetic field in a pure QND interaction of atoms and light.The amount of entanglement is significantly enhanced.FIG.1:Von Neumann Entropy of the reduced state of atoms versus coupling strength kappa for the state of equation 5(full line)and for the state generated without magnetic field in a pure QND interaction (dashed line)with the same cou-pling strength.Application of a magnetic field significantly enhances the amount of light-atom entanglement.III.TELEPORTATION OF LIGHT ONTOATOMSIn this section we will show how the multimode entan-glement between light and atoms generated in the scat-tering process can be employed in a teleportation proto-col which is a simple generalization of the standard proto-col for continuous variable teleportation using EPR-type entangled states [8,9].We first present the protocol and evaluate its fidelity and then analyze its performance un-der realistic experimental conditions.A.Basic protocolFigure 2depicts the basic scheme which,as usually,consists of a Bell measurement and a feedback operation.Input The coherent state to be teleported is encoded in a pulse which is linearly polarized orthogonal to the classical driving pulse and whose carrier frequency lies at the upper sideband,i.e.at ωc +Ω.The pulse envelope has to match the one of the classical pulse.As is shown in appendix B,canonical operators y,q with [y,q ]=i de-scribing this mode can conveniently be expressed in terms of cosine and sine modulation modes,analogous to equa-tions (4),defined with respect to the carrier frequency.One findsy =12(y s +q c ),q =−12(y c −q s ).(6)A coherent input amounts to having initially ∆y 2=∆q 2=1/2and an amplitude y , q with mean photon number n ph =( y 2+ q 2)/2.Bell measurement This input is combined at a beam splitter with the classical pulse and the scattered light.At the ports of the beam splitter Stokes vector compo-nents S y and S z are measured by means of standard po-larization measurements.Given the classical pulse in xFIG.2:Scheme for teleportation of light onto atoms:As de-scribed in section II,a classical pulse(linearly polarized along x)propagating along the positive z direction is scattered offan atomic ensemble contained in a glass cell and placed in a constant magneticfield B along x.Classical pulse and scat-tered light(linearly polarized along y)are overlapped with a with a coherent pulse(linearly polarized along z)at beam splitter B S.By means of standard polarization measurements Stokes vector components S y and S z are measured at one and the other port respectively,realizing the Bell measurement. The Fourier components at Larmor frequencyΩof the corre-sponding photocurrents determine the amount of conditional displacement of the atomic spin which can be achieved by ap-plying a properly timed transverse magneticfield b(t).See section III A for details.polarization this amounts to a homodyne detection of in-and out-of-phasefields of the orthogonal polarization component.The resulting photocurrents are numerically demodulated to extract the relevant sine and cosine com-ponents at the Larmor frequency[20].Thus one effec-tively measures the commuting observables˜x c=12 x out c+y c,˜x s=12 x out s+y s ,(7)˜q c=12 p out c−q c ,˜q s=12 p out s−q s .Let the respective measurement results be given by ˜Xc,˜X s,˜Q c and˜Q s.Feedback Conditioned on these results the atomic state is then displaced by an amount˜X s−˜Q c in X and −˜X c−˜Q s in P.This can be achieved by means of two fast radio-frequency magnetic pulses separated by a quar-ter of a Larmor period.In the ensemble average thefinal state of atoms is simply given byXfin=X out+˜x s−˜q c,Pfin=P out−˜x c−˜q s.(8)This description of feedback is justified rigorously in ap-pendix C.Relating these expressions to input operators,wefind by means of equations(5),(6)and(7)Xfin= 1−κ√2 2p in c+12x in s−16 κ2 P in−121−κ√√2 2p in s,1+q.(9b)This is the main result of this section.Teleportationfidelity Taking the mean of equations (9)with respect to the initial state all contributions due to input operators and back action modes vanish such that Xfin = y and Pfin = q .Thus,the am-plitude of the coherent input light pulse is mapped on atomic spin quadratures as desired.In order to prove faithful teleportation also the variances have to be con-served.It is evident from(9)that thefinal atomic spin variances will be increased as compared to the coherent input.These additional terms describe unwanted excess noise and have to be minimized by a proper choice of the couplingκ.As afigure of merit for the telepor-tation protocol we use thefidelity,i.e.squared over-lap,of input andfinal state.Given that the means are transmitted correctly thefidelity is found to be F= 2 (1+2(∆Xfin)2)(1+2(∆Pfin)2) −1/2.The variances of thefinal spin quadratures are readily calculated tak-ing into account that all modes involved are independent and have initially a normalized variance of1/2.In this way a theoretical limit on the achievablefidelity can be derived depending solely on the coupling strengthκ.In figure3we take advantage of the fact that the amount of entanglement between light and atoms is a monotonously increasing function ofκsuch that we can plot thefi-delity versus the entanglement.This has the advantage that we can compare the performance of our teleporta-tion protocol with the canonical one[8,9]which uses a two-mode squeezed state of the same entanglement as a resource and therefore maximizes the teleportationfi-delity for the given amount of entanglement.No physical state can achieve a higherfidelity with the same entan-glement.This follows from the results of[30]where it was shown that two-mode squeezed states minimize the EPR variance(and therefore maximize the teleportation fidelity)for given entanglement.The theoreticalfidelity achievable in our protocol is maximized forκ≃1.64cor-responding to F≃.77.But also for experimentally more feasible values ofκ≃1can thefidelity well exceed the classical limit[28,29]of1/2and,moreover,compari-son with the values achievable with a two-mode squeezed state shows that our protocol is close to optimal.FIG.3:(a)Theoretical limit on the achievablefidelity F ver-sus entanglement between atoms and light measured by the von Neumann entropy E vN of the reduced state of atoms.The grey area is unphysical.For moderate amounts of entangle-ment our protocol is close to optimal.(b)Coupling strengthκversus entanglement.The dashed lines indicate the maximal fidelity of F=.77which is achieved forκ=1.64.B.Noise effects and Gaussian distributed input Under realistic conditions the teleportationfidelity will be degraded by noise effects like decoherence of the atomic spin state,light absorption and reflection losses and also because the coupling constantκis experimen-tally limited to valuesκ≃1.On the other hand the classicalfidelity bound to be beaten will be somewhat higher than1/2since the coherent input states will nec-essarily be drawn according to a distribution with afinite width in the mean photon number¯n.In this section we analyze the efficiency of the teleportation protocol un-der these conditions and show that it is still possible to surpass any classical strategy for the transmission and storage of coherent states of light[28,29].During the interaction atomic polarization decays due to spontaneous emission and collisional relaxation.In-cluding a transverse decay thefinal state of atoms is given byX out= √βf X,(10a) P out= √βf P.(10b)as follows from the discussion in appendix A.βis the atomic decay parameter and f X,f P are Langevin noise operators with zero mean.Their variance is experimen-tally found to be close to the value corresponding to a coherent state such that f2X = f2P =1/2.Light absorption and reflection losses can be taken into account in the same way asfinite detection efficiency.For example the statistics of measurement outcome˜X s will not stem from the signal mode˜x s alone but rather from the noisy mode√ǫfx,swhereǫis the photon loss parameter and f x,s is a Langevin noise operator of zero mean and variance f2x,s =1/2.Analogous expres-sions have to be used for the measurements of˜x c,˜q s and ˜q c which will be adulterated by Langevin terms f x,c,f q,s and f q,c respectively.In principle each of the measure-ment outcomes can be fed back with an independently chosen gain but for symmetry reasons it is enough to distinguish gain coefficients g x,g q for the measurement outcomes of sine and cosine components of x and q re-spectively.Including photon loss,finite gain and atomic decay,as given in(10),equations(8),describing thefinal state of atoms after the feed back operation,generalize toXfin= βf X+g x √ǫf x,s (11a)−g q √ǫf q,c ,Pfin= βf P−g x √ǫf x,c (11b)−g q √ǫf q,s .For non unit gains a given coherent amplitude( y , q ) will not be perfectly teleported onto atoms and the cor-respondingfidelity will be degraded by this mismatch according toF( y , q )=2[1+2(∆Xfin)2][1+2(∆Pfin)2]·exp−( y − Xfin )21+2(∆Pfin)2 .If the input amplitudes are drawn ac-cording to a Gaussian distribution p( y , q )=exp[−( y 2+ q 2)/2¯n]/2π¯n with mean photon number¯n the averagefidelity[with respect to ( y , q )]is readily calculated.The exact expression in terms of initial operators can then be derived by means of equations(5),(6),(7)and(11)but is not particularly enlightening.Infigure4we plot the averagefidelity, optimized with respect to gains g x,g q,in its dependence on the atomic decayβfor various values of photon loss ǫ.We assume a realistic valueκ=0.96for the coupling constant and a mean number of photons¯n=4for the distribution of the coherent input.For feasible values ofβ,ǫ 0.2the averagefidelity is still well above the classical bound on thefidelity[28,29].This proves that the proposed protocol is robust against the dominating noise effects in this system.The experimental feasibility of the proposal is illus-trated with the following example.Consider a sam-ple of N at=1012Cesium atoms in a glass cell placed in a constant magneticfield along the x-direction caus-ing a Zeeman splitting ofΩ=350kHz in the F=4 ground state multiplet.The atoms are pumped into m F=4and probed on the D2(F=4→F′=3,4,5)FIG.4:(a)Averagefidelity achievable in the presenceof atomic decayβ,reflection and light absorption lossesǫ=8%,12%,16%,couplingκ=0.96and Gaussian dis-tributed input states with mean photon number¯n=4.Thefi-delity benchmark is in this case5/9(dashed line).(b)Respec-tive optimal values for gains g x(solid lines)and g q(dashedlines).transition.The classical pulse contains an overall num-ber of N ph=2.51013photons,is detuned to the blueby∆=1GHz,has a duration T=1ms and can have an effective cross section of A≃6cm2due to thermalmotion of atoms.Under these conditions the tensor polarizability can be neglected(∆/ωhfs≃10−1).Also n0=ΩT=350justifies the use of independent scatter-ing modes.The couplingκ≃1and the depumping ofground state populationη≃10−1as desired.IV.SPIN SQUEEZING AND STATE READ-OUTIn this section we present a scheme for reading out either of the atomic spin components X,P by means of a probe pulse interacting with the atoms in the one way as described in section II.The proposed scheme allows one,on the one hand,to verify successful receipt of the coherent input subsequent to the teleportation protocol of section III and,on the other hand,enables to generate spin squeezing if it is performed on a coherent spin state. It is well known[15,31]and was demonstrated exper-imentally[10,11]that the pure interaction V,as given in equation(1),can be used to perform a QND measure-ment of either of the transverse spin components.At first sight this seems not to be an option in the scenario under consideration since the local term H at,account-ing for Larmor precession,commutes with neither of the spin quadratures such that the total Hamiltonian does not satisfy the QND criteria[5,6].As we have shown in section II Larmor precession has two effects:Scattered light is correlated with both transverse components andsuffers from back action mediated by the atoms.Thus, in order to read out a single spin component one has to overcome both disturbing effects.Our claim is that this can be achieved by a simultane-ous measurement of x outc,p outs,p outs,1or x outs,p outc,p outc,1 if,respectively,X or P is to be measured.In the follow-ing we consider in particular the former case but every-thing will hold with appropriate replacements also for a measurement of P.As shown infigure5the set of observablesx outc,p outs,p outs,1can be measured simultaneously by a measurement of Stokes component S y after aπ/2rota-tion is performed selectively on the sine component of the scattered light.The cosine component of the corre-sponding photocurrent will give an estimate of x outcandthe sine component of p outs.Multiplying the photocur-rent’s sine component by the linear function defining the back action mode,equation(B1),will give in additionan estimate of p outs,1.Note that thefield out of phase is conserved in the interaction such thatp outs,1=p in s,1,p outc,1=p in c,1,(12) i.e.the results will have shot noise limited variance.It is then evident from equation(5c)that the respective photocurrents together with an a priori knowledge ofκare sufficient to estimate the mean X .The conditional variances after the indicated measure-ments are∆X2|{x out c,p out s,p out s,1}=(∆X in)222,(13b)corresponding to a pure state.Obviously the variance in X is squeezed by a factor(1+κ2/2)−1.Note that the squeezing achieved in a QND measurement without magneticfield but otherwise identical parameters is given by(1+κ2)−1.From this we conclude that the quality of the estimate for X ,as measured f.e.by input-output coefficients known from the theory of QND measurements [5,6],can be the same as in the case without Larmor precession albeit only for a higher couplingκ. Equations(13)are conveniently derived by means of the formalism of correlation matrices [32].For the operator valued vector R=(X,P, x c,p c,x s,p s,x c,1,p c,1,x s,1,p s,1)equations(5),(12) and(B2)define via R out=S(κ) R in a symplectic linear transformation S(κ).The contributions of p in c,2andp in s,2to x outs,1and x outc,1as given in(B2)are treated as noise and do not contribute to the symplectic trans-formation S but enter the input-output relation for the correlation matrix as an additional noise term as follows.The correlation matrix is as usually defined byγi,j=tr{ρ(R i R j+R j R i)}.The initial state is then an10×10identity matrix and thefinal state is γout=S(κ)S(κ)T+γnoise where the diagonal matrixFIG.5:Scheme for spin measurement:After the scattering a π/2rotation is performed on the scattered light modulated at the Larmor frequency such as to affect only the sine(cosine) component.Standard polarization measurement of S y and appropriate postprocessing allows to read out the mean of X(P),leaving the atoms eventually in a spin squeezed state.γnoise=diag[0,0,0,0,0,0,1,0,1,0](κ/2)4/15accountsfor noise contributions due correlations to second order back action modes c.f.equations(B2).In order to evaluate the atomic variances after a measurement ofx out c ,p outs,p outs,1the correlation matrixγout is split upinto blocks,γout= A C C T Bwhere A is the2×2subblock describing atomic variances. Now,the state A′after the measurement can be found by evaluating[32]A′=A−limx,n→∞C1(∆−iΓ/2) a01+ia1 F×(A3)8 with real dimensionless coefficients a j of order unity andΓthe excited states’decay rate.The non-hermitian partof the resulting Hamilton operator describes the effectof light absorption and loss of ground state populationdue to depumping in the course of interaction.In thefollowing we will focus on the coherent interaction and,for the time being,take into account only the hermitiancomponent.The effects of light absorption and atomicdepumping are treated below.Coherent interaction Since scattering of light occurspredominantly in the forward direction[22]it is legiti-mate to adopt a one dimensional model such that the(negative frequency component of the)electricfield prop-agating along z is given byE(−)(z,t)=E(−)(z) ey+E(−)(z,t) e xE(−)(z)=ρ(ωc) b dωa†(ω)e−ikzE(−)(z,t)=ρ(ωc)ωc/4πǫ0Ac and A denotes the pulse’scross sectional area,N ph the overall number of photonsin the pulse and T its duration.We restrict thefieldin x polarization to the classical probe pulse,since onlythe coupling of atoms to the y polarization is enhancedby the coherent probe.Furthermore we implicitly assumefor the classical pulse a slowly varying envelope such thatit arrives at z=0at t=0and is then constant fora time bining this expression for thefield withexpressions(A2)and(A3)for the atomic polarizability inequation(A1)yieldsV=−i κ4πJTbdω d zj(z) a(ω)e−i[(k c−k)z−ωc t]−h.c.where we defined a dimensionless coupling constant κ=N ph Ja1σΓ/2A∆withσthe scattering cross sec-tion on resonance.We now definefield quadratures for spatially localized modes[24,25]asx(z)=14πbdω a(ω)e−i(k c−k)z+h.c. ,(A4a)p(z)=−i4πbdω a(ω)e−i(k c−k)z−h.c. (A4b)with commutation relations[x(z),p(z′)]=icδ(z−z′) where the delta function has to be understood to have a width on the order of c/b.Since we assumed thatΩ≪b, the time it takes for such a fraction of the pulse to cross the ensemble is much smaller than the Larmor period 1/Ω.During the interaction with one of these spatiallylocalized modes the atomic state does not change appre-ciable and we can simplify the interaction operator to V= κ(JT)−1/2J z p(0)where J z= i F(i)z and we as-sumed that the ensemble is located at z=0and changedto a frame rotating at the carrier frequencyωc.A last approximation concerns the description of theatomic spin state.Initially the sample is prepared in acoherent spin state with maximal polarization along x, i.e.in the eigenstate of J x with maximal eigenvalue J. We can thus make use of the Holstein-Primakoffapprox-imation[23]which allows to describe the spin state as a Gaussian state of a single harmonic oscillator.The first step is to express collective step up/down operators (along x),J±=J y±iJ z,in terms of bosonic creation and annihilation operators,[b,b†]=11,asJ+=√11−b†b/2J b,J−=√11−b†b/2J. It is easily checked that these operators satisfy the correct commutation relations[J+,J−]=2J x if one identifies J x=J−b†b.The fully polarized initial state thus corresponds to the ground state of the harmonic oscillator.Note that this map-ping is exact.Under the condition that b†b ≪J one can approximate J+≃√2Jb†and therefore J z≃−i 2and P=−i(b−b†)/√。
Bosch 蒸汽干衣机商品说明书

1234Display panel / ButtonsCongratulations – you have chosen to buy a modern, high-quality Bosch domestic appliance.A distinctive feature of the condensation dryer with heat pump and automatic cleaning of the heat exchanger is its low energy consumption.Every dryer which leaves our factory is carefully checked to ensure that it functions correctly and is in perfect condition.Should you have any questions, our after-sales service will be pleased to help.Environmentally-responsible disposalDispose of packaging in an environmentally-responsible manner. This appliance is designated according to European directive 2002/96/EC which governs waste electrical and electronic equipment - WEEE).The directive provides an EU-wide framework for the return and recycling of used appliances.For further information about our products, accessories, spare parts and services, please visit: DryingÏ÷Spin speed Ï÷Spin speed exclusively to display the anticipated drying time (Ð÷Adjust drying level ContentsPageʋPreparation . . . . . . . . . . . . . . . . . . . . . . . . . . .2ʋControl panel . . . . . . . . . . . . . . . . . . . . . . . . . .2ʋDrying . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3/4ʋDisplay panel and buttons . . . . . . . . . . . . . . . .5ʋNotes on laundry . . . . . . . . . . . . . . . . . . . . . . .6ʋOverview of programmes . . . . . . . . . . . . . . . .7ʋImportant Information . . . . . . . . . . . . . . . . . . .8ʋConsumption rates. . . . . . . . . . . . . . . . . . . . . .9ʋNormal noises . . . . . . . . . . . . . . . . . . . . . . . .10ʋWhat to do if ... / After-sales service. . . . .11/12ʋSafety instructions . . . . . . . . . . . . . . . . . . . . .13Sort the laundryThe fluff filter consists of two parts.Clean the inner and outer fluff filter every 1.Open the door. Remove fluff from thedoor/door area.back together and reinsert into the appliance.6Labelling of fabricsFollow the manufacturer's care information.hc Drying tips––Remove the laundry and switch off Start/Pause buttonTwo-part fluff filterAir inletClean fluff filters reduce energy consumption.Depending on the selected programme, it may not be possible to select individual options button in the centre of the programme selector.Select the programme by turning the outer ring on the programme selector (the selector can be turned in both directions).Change the functions in the display using the buttons C/board dry extra,Dryeren Instruction manualDisplay Congratulations – you have chosen to buy a modern, high-qualityBosch domestic appliance.A distinctive feature of the condensation dryer with heat pump andautomatic cleaning of the heat exchanger is its low energyconsumption.Every dryer which leaves our factory is carefully checked to ensurethat it functions correctly and is in perfect condition.Should you have any questions, our after-sales service will bepleased to help.Environmentally-responsible disposalDispose of packaging in an environmentally-responsible manner.This appliance is designated according to European directive2002/96/EC which governs waste electrical and electronicequipment - WEEE).The directive provides an EU-wide framework for the return andrecycling of used appliances.For further information about our products, accessories, spare partsand services, please visit: Contents PageʋPreparation . . . . . . . . . . . . . . . . . . . . . . . . . . .2ʋControl panel . . . . . . . . . . . . . . . . . . . . . . . . . .2ʋDrying . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3/4ʋDisplay panel and buttons . . . . . . . . . . . . . . . .5ʋNotes on laundry . . . . . . . . . . . . . . . . . . . . . . .6ʋOverview of programmes . . . . . . . . . . . . . . . .7ʋImportant Information . . . . . . . . . . . . . . . . . . .8ʋConsumption rates. . . . . . . . . . . . . . . . . . . . . .9ʋNormal noises . . . . . . . . . . . . . . . . . . . . . . . .10ʋWhat to do if ... / After-sales service. . . . .11/12ʋSafety instructions . . . . . . . . . . . . . . . . . . . . .13button in the centre of the programme selector.Select the programme by turning the outer ring on the programmselector (the selector can be turned in both directions).1234DryingSort the laundryThe fluff filter consists of two parts.Clean the inner and outer fluff filter every 1.Open the door. Remove fluff from thedoor/door area.back together and reinsert into the appliance.6Remove the laundry and switch off Start/Pause buttonTwo-part fluff filterAir inletClean fluff filters reduce energy consumption.Change the functions in the display using the buttons C/board dry extra,Ï÷Spin speedÏ÷Spin speedexclusively to display the anticipated drying time (Ð÷Adjust drying level Labelling of fabricsFollow the manufacturer's care information. hc–Drying tips––Depending on the selected programme, it may not be possible to select individual optionsDryeren Instruction manual。
iProx无接近感触感应传感器安装说明说明书

iProx inductive proximity sensors—mounting close togetherDescriptionThe iProxீ is a powerful family of inductive proximity sensors featuring high sensing performance right out of the box. What makes the iProx unique from other inductive sensors is the ability to extensively customize the operating characteristics to suit a particular application.Mounting sensors close togetherWhen mounting iProx sensors close together, it is necessary to take into consideration problems that can be causedby two or more sensors communicating with each other (also known as “cross-talk”). This problem can arise when two or more sensors are mounted side by side (as shown in Figure 1) or facing each other (as shown in Figure 2).Figure 1. Proximity Sensors Mounted Side by SideFigure 2. Proximity Sensors Mounted Facing Each Other Standard inductive proximity sensors have a similar frequency and will interfere with each other when operated close together. Until iProx, your best solution may have been to buy special sensors designed to operate on different frequencies. The disadvantage of this solution is that your choice of sensing range and body style is usually very limited.The DC versions of iProx have three different noise immunity settings, while the AC versions have two settings. See Figure 3, Figure 4, Figure 5, and Figure 6 for optimum noise immunity settings based upon center-to-center spacing between mounted sensors. Note that in some cases, more than one noise immunity setting is available. In this case,the operator can choose the combinations of noise immunity settings most desirable for the application. The iProx breaks many of the traditional rules of inductive proximity sensors. For instance, it is possible to mount the sensors so that the sensing fields overlap, so long as the proper noise immunity settings are chosen.Figure 3. DC iProx Side-by-Side Configuration (S1)T able 1. DC iProx Side-by-Side Configuration (S1)Diameter (D)Sideby SideHigh NoiseImmunityFactoryDefault12 mm Shielded0–19 mmat 50 Hz0 mm–Infinityat 10 Hz30 mm–Infinityat 580 HzUnshielded0–40 mmat 50 Hz0 mm–Infinityat 10 Hz100 mm–Infinityat 300 Hz18 mm Shielded0–35 mmat 50 Hz0 mm–Infinityat 10 Hz80 mm–Infinityat 390 HzUnshielded0–75 mmat 50 Hz0 mm–Infinityat 10 Hz130 mm–Infinityat 150 Hz30 mm Shielded0–65 mmat 50 Hz0 mm–Infinityat 10 Hz110 mm–Infinityat 240 HzUnshielded0–75 mmat 50 Hz0 mm–Infinityat 10 Hz130 mm–Infinityat 90 Hz2Technical Data TD05301002EEffective June 2013iProx inductive proximity sensors—mounting close togetherEATON Figure 4. DC iProx Facing Configuration (S 2)T able 2. DC iProx Facing Configuration (S 2)Diameter (D)Side by SideHigh Noise ImmunityFactory Default12 mm Shielded 0–25 mm at 50 Hz 0 mm–Infinity at 10 Hz50 mm–Infinity at 580 HzUnshielded0–55 mm at 50 Hz 0 mm–Infinity at 10 Hz120 mm–Infinity at 300 Hz18 mm Shielded 0–45 mm at 50 Hz 0 mm–Infinity at 10 Hz100 mm–Infinity at 390 HzUnshielded0–90 mm at 50 Hz 0 mm–Infinity at 10 Hz160 mm–Infinity at 150 Hz30 mm Shielded 0–80 mm at 50 Hz 0 mm–Infinity at 10 Hz130 mm–Infinity at 240 HzUnshielded0–90 mm at 50 Hz 0 mm–Infinity at 10 Hz 160 mm–Infinity at 90 HzFigure 5. AC iProx Side-by-Side Configuration (S 1)T able 3. AC iProx Side-by-Side Configuration (S 1) ቢDiameter (D)Factory DefaultHigh Noise Immunity12 mm Shielded 0–19 mm/30 mm–Infinity at 30 Hz 0 mm–Infinity at 10 HzUnshielded0–40 mm/50 mm–Infinity at 30 Hz 0 mm–Infinity at 10 Hz18 mm Shielded 0–35 mm/80 mm–Infinity at 30 Hz 0 mm–Infinity at 10 HzUnshielded0–75 mm/130 mm–Infinity at 30 Hz 0 mm–Infinity at 10 Hz30 mm Shielded 0–65 mm/110 mm–Infinity at 30 Hz 0 mm–Infinity at 10 HzUnshielded0–75 mm/130 mm–Infinity at 30 Hz0 mm–Infinity at 10 HzቢThese specifications may not meet final product specifications.Figure 6. AC iProx Facing Configuration (S 2)T able 4. AC iProx Facing Configuration (S 2) ቢDiameter (D)Factory DefaultHigh Noise Immunity12 mm Shielded 0–25 mm/50 mm–Infinity at 30 Hz 0 mm–Infinity at 10 HzUnshielded0–55 mm/120 mm–Infinity at 30 Hz 0 mm–Infinity at 10 Hz18 mm Shielded 0–45 mm/100 mm–Infinity at 30 Hz 0 mm–Infinity at 10 HzUnshielded0–90 mm/160 mm–Infinity at 30 Hz 0 mm–Infinity at 10 Hz30 mm Shielded 0–80 mm/130 mm–Infinity at 30 Hz 0 mm–Infinity at 10 HzUnshielded0–90 mm/160 mm–Infinity at 30 Hz0 mm–Infinity at 10 Hzቢ These specifications may not meet final product specifications.ote:N There is a correlation between noise immunity and operatingfrequency. When setting the sensor for side-by-side sensing or high noise immunity, the operating frequency of the sensor will be reduced. Refer to the iProx Programming Software User Guide (P50228) for details.3Technical Data TD05301002EEffective June 2013iProx inductive proximity sensors—mounting close togetherEATON Setting iProx for side-by-side or facing sensingTo set up iProx sensors for side-by-side or face-to-face operation, just follow the simple procedure below. This process requires that the iProx programming software be installed on your computer. Consult the iProx Programming Software User Guide (P50228) for detailed installation and operating instructions.1. Connect the programming device to your computer. Connectthe remote programmer (E59RP1) to your computer’s serial or USB port. 2. Launch that the iProx programming software from the Start menu.3. Affix the magnetic puck to the face of the iProx sensor. To ensureproper mounting, it may be necessary to remove the mounting nuts from the sensor.4. Ensure that the sensor(s) you intend to program are powered on.5. The iProx programming software should automatically detect thesensor(s). See the “Connection Status” information displayed at the bottom of the “Getting Started” window to confirm that the software has detected the sensor(s). If the software does not detect the connected sensor, you may need to select a wider COM port range in the software settings. For more troubleshooting information, consult the iProx Programming Software User Guide (P50228) by clicking on the “Help” window and selecting “Contents/FAQs.”6. From the “Getting Started” window, click “Configure iProx” to modify the parameters of your iProx sensor.7. Under “Step 1: Select Sensor to Modify,” click the drop-downmenu and select the sensor you want to modify. It may take a few seconds to communicate with the sensor. Once this is complete, the rest of the screen should become enabled, allowing you to modify the parameters of the selected sensor.8. Under the “Response Time/Noise Immunity” section at the topright of the screen, drag the slider until you see the side-by-side value enabled.9. Click the “Program” button at the bottom of the screen. Repeatthis process for each sensor.Side-by-side vs. array sensingAlthough this procedure allows sensors to be operated close to other sensors in a side-by-side configuration (see Figure 7), Eaton does not recommend that iProx sensors be used in an array (Figure 8). It is possible that cross-ommunication can occur between sensors in an array. If the requirements of your application call for this arrangement, please contact Eaton’s Sensor ApplicationEngineering Department at 1-800-426-9184. They will work with youto find a solution.Figure 7. Sensors Oriented Side by SideFigure 8. Sensors Oriented in an ArrayTechnical Data TD05301002E Effective June 2013iProx inductive proximity sensors—mounting close togetherEaton is a registered trademark.All other trademarks are propertyof their respective owners.Eaton1000 Eaton BoulevardCleveland, OH 44122United States© 2013 EatonAll Rights ReservedPrinted in USAPublication No. TD05301002E / Z13430 June 2013。
振动噪声专业词汇

1 振动信号的时域、频域描述振动过程(Vibration Process)简谐振动(Harmonic Vibration)周期振动(Periodic Vibration)准周期振动(Ouasi-periodic Vibration)瞬态过程(Transient Process)随机振动过程(Random Vibration Process)各态历经过程(Ergodic Process)确定*过程(Deterministic Process)振幅(Amplitude)相位(Phase)初相位(Initial Phase)频率(Frequency)角频率(Angular Frequency)周期(Period)复数振动(Complex Vibration)复数振幅(Complex Amplitude)峰值(Peak-value)平均绝对值(Average Absolute Value)有效值(Effective Value,RMS Value)均值(Mean Value,Average Value)傅里叶级数(FS,Fourier Series)傅里叶变换(FT,Fourier Transform)傅里叶逆变换(IFT,Inverse Fourier Transform) 离散谱(Discrete Spectrum)连续谱(Continuous Spectrum)傅里叶谱(Fourier Spectrum)线*谱(Linear Spectrum)幅值谱(Amplitude Spectrum)相位谱(Phase Spectrum)均方值(Mean Square Value)方差(Variance)协方差(Covariance)自协方差函数(Auto-covariance Function)互协方差函数(Cross-covariance Function)自相关函数(Auto-correlation Function)互相关函数(Cross-correlation Function)标准偏差(Standard Deviation)相对标准偏差(Relative Standard Deviation)概率(Probability)概率分布(Probability Distribution)高斯概率分布(Gaussian Probability Distribution) 概率密度(Probability Density)集合平均(Ensemble Average)时间平均(Time Average)功率谱密度(PSD,Power Spectrum Density)自功率谱密度(Auto-spectral Density)互功率谱密度(Cross-spectral Density)均方根谱密度(RMS Spectral Density)能量谱密度(ESD,Energy Spectrum Density)相干函数(Coherence Function)帕斯瓦尔定理(Parseval''s Theorem)维纳,辛钦公式(Wiener-Khinchin Formula)2 振动系统的固有特*、激励与响应振动系统(Vibration System)激励(Excitation)响应(Response)单自由度系统(Single Degree-Of-Freedom System) 多自由度系统(Multi-Degree-Of- Freedom System) 离散化系统(Discrete System)连续体系统(Continuous System)刚度系数(Stiffness Coefficient)自由振动(Free Vibration)自由响应(Free Response)强迫振动(Forced Vibration)强迫响应(Forced Response)初始条件(Initial Condition)固有频率(Natural Frequency)阻尼比(Damping Ratio)衰减指数(Damping Exponent)阻尼固有频率(Damped Natural Frequency)对数减幅系数(Logarithmic Decrement)主频率(Principal Frequency)无阻尼模态频率(Undamped Modal Frequency)模态(Mode)主振动(Principal Vibration)振型(Mode Shape)振型矢量(Vector Of Mode Shape)模态矢量(Modal Vector)正交* (Orthogonality)展开定理(Expansion Theorem)主质量(Principal Mass)模态质量(Modal Mass)主刚度(Principal Stiffness)模态刚度(Modal Stiffness)正则化(Normalization)振型矩阵(Matrix Of Modal Shape)模态矩阵(Modal Matrix)主坐标(Principal Coordinates)模态坐标(Modal Coordinates)模态分析(Modal Analysis)模态阻尼比(Modal Damping Ratio)频响函数(Frequency Response Function)幅频特* (Amplitude-frequency Characteristics)相频特* (Phase frequency Characteristics)共振(Resonance)半功率点(Half power Points)波德图(BodéPlot)动力放大系数(Dynamical Magnification Factor)单位脉冲(Unit Impulse)冲激响应函数(Impulse Response Function)杜哈美积分(Duhamel’s Integral)卷积积分(Convolution Integral)卷积定理(Convolution Theorem)特征矩阵(Characteristic Matrix)阻抗矩阵(Impedance Matrix)频响函数矩阵(Matrix Of Frequency Response Function) 导纳矩阵(Mobility Matrix)冲击响应谱(Shock Response Spectrum)冲击激励(Shock Excitation)冲击响应(Shock Response)冲击初始响应谱(Initial Shock Response Spectrum)冲击剩余响应谱(Residual Shock Response Spectrum) 冲击最大响应谱(Maximum Shock Response Spectrum) 冲击响应谱分析(Shock Response Spectrum Analysis)3 模态试验分析模态试验(Modal Testing)机械阻抗(Mechanical Impedance)位移阻抗(Displacement Impedance)速度阻抗(Velocity Impedance)加速度阻抗(Acceleration Impedance)机械导纳(Mechanical Mobility)位移导纳(Displacement Mobility)速度导纳(Velocity Mobility)加速度导纳(Acceleration Mobility)驱动点导纳(Driving Point Mobility)跨点导纳(Cross Mobility)传递函数(Transfer Function)拉普拉斯变换(Laplace Transform)传递函数矩阵(Matrix Of Transfer Function)频响函数(FRF,Frequency Response Function)频响函数矩阵(Matrix Of FRF)实模态(Normal Mode)复模态(Complex Mode)模态参数(Modal Parameter)模态频率(Modal Frequency)模态阻尼比(Modal Damping Ratio)模态振型(Modal Shape)模态质量(Modal Mass)模态刚度(Modal Stiffness)模态阻力系数(Modal Damping Coefficient)模态阻抗(Modal Impedance)模态导纳(Modal Mobility)模态损耗因子(Modal Loss Factor)比例粘*阻尼(Proportional Viscous Damping)非比例粘*阻尼(Non-proportional Viscous Damping)结构阻尼(Structural Damping,Hysteretic Damping) 复频率(Complex Frequency)复振型(Complex Modal Shape)留数(Residue)极点(Pole)零点(Zero)复留数(Complex Residue)随机激励(Random Excitation)伪随机激励(Pseudo Random Excitation)猝发随机激励(Burst Random Excitation)稳态正弦激励(Steady State Sine Excitation)正弦扫描激励(Sweeping Sine Excitation)锤击激励(Impact Excitation)频响函数的H1 估计(FRF Estimate by H1)频响函数的H2 估计(FRF Estimate by H2)频响函数的H3 估计(FRF Estimate by H3)单模态曲线拟合法(Single-mode Curve Fitting Method)多模态曲线拟合法(Multi-mode Curve Fitting Method)模态圆(Mode Circle)剩余模态(Residual Mode)幅频峰值法(Peak Value Method)实频-虚频峰值法(Peak Real/Imaginary Method)圆拟合法(Circle Fitting Method)加权最小二乘拟合法(Weighting Least Squares Fitting method) 复指数拟合法(Complex Exponential Fitting method)1.2 振动测试的名词术语1 传感器测量系统传感器测量系统(Transducer Measuring System)传感器(Transducer)振动传感器(Vibration Transducer)机械接收(Mechanical Reception)机电变换(Electro-mechanical Conversion)测量电路(Measuring Circuit)惯*式传感器(Inertial Transducer,Seismic Transducer)相对式传感器(Relative Transducer)电感式传感器(Inductive Transducer)应变式传感器(Strain Gauge Transducer)电动力传感器(Electro-dynamic Transducer)压电式传感器(Piezoelectric Transducer)压阻式传感器(Piezoresistive Transducer)电涡流式传感器(Eddy Current Transducer)伺服式传感器(Servo Transducer)灵敏度(Sensitivity)复数灵敏度(Complex Sensitivity)分辨率(Resolution)频率范围(Frequency Range)线*范围(Linear Range)频率上限(Upper Limit Frequency)频率下限(Lower Limit Frequency)静态响应(Static Response)零频率响应(Zero Frequency Response)动态范围(Dynamic Range)幅值上限Upper Limit Amplitude)幅值下限(Lower Limit Amplitude)最大可测振级(Max.Detectable Vibration Level)最小可测振级(Min.Detectable Vibration Level)信噪比(S/N Ratio)振动诺模图(Vibration Nomogram)相移(Phase Shift)波形畸变(Wave-shape Distortion)比例相移(Proportional Phase Shift)惯*传感器的稳态响应(Steady Response Of Inertial Transducer)惯*传感器的稳击响应(Shock Response Of Inertial Transducer)位移计型的频响特* (Frequency Response Characteristics Vibrometer)加速度计型的频响特* (Frequency Response Characteristics Accelerometer) 幅频特*曲线(Amplitude-frequency Curve)相频特*曲线(Phase-frequency Curve)固定安装共振频率(Mounted Resonance Frequency)安装刚度(Mounted Stiffness)有限高频效应(Effect Of Limited High Frequency)有限低频效应(Effect Of Limited Low Frequency)电动式变换(Electro-dynamic Conversion)磁感应强度(Magnetic Induction,Magnetic Flux Density)磁通(Magnetic Flux)磁隙(Magnetic Gap)电磁力(Electro-magnetic Force)相对式速度传(Relative Velocity Transducer)惯*式速度传感器(Inertial Velocity Transducer)速度灵敏度(Velocity Sensitivity)电涡流阻尼(Eddy-current Damping)无源微(积)分电路(Passive Differential (Integrate) Circuit)有源微(积)分电路(Active Differential (Integrate) Circuit)运算放大器(Operational Amplifier)时间常数(Time Constant)比例运算(Scaling)积分运算(Integration)微分运算(Differentiation)高通滤波电路(High-pass Filter Circuit)低通滤波电路(Low-pass Filter Circuit)截止频率(Cut-off Frequency)压电效应(Piezoelectric Effect)压电陶瓷(Piezoelectric Ceramic)压电常数(Piezoelectric Constant)极化(Polarization)压电式加速度传感器(Piezoelectric Acceleration Transducer)中心压缩式(Center Compression Accelerometer)三角剪切式(Delta Shear Accelerometer)压电方程(Piezoelectric Equation)压电石英(Piezoelectric Quartz)电荷等效电路(Charge Equivalent Circuit)电压等效电路(Voltage Equivalent Circuit)电荷灵敏度(Charge Sensitivity)电压灵敏度(Voltage Sensitivity)电荷放大器(Charge Amplifier)适调放大环节(Conditional Amplifier Section)归一化(Uniformization)电荷放大器增益(Gain Of Charge Amplifier)测量系统灵敏度(Sensitivity Of Measuring System)底部应变灵敏度(Base Strain Sensitivity)横向灵敏度(Transverse Sensitivity)地回路(Ground Loop)力传感器(Force Transducer)力传感器灵敏度(Sensitivity Of Force Transducer)电涡流(Eddy Current)前置器(Proximitor)间隙-电压曲线(Voltage vs Gap Curve)间隙-电压灵敏度(V oltage vs Gap Sensitivity)压阻效应(Piezoresistive Effect)轴向压阻系数(Axial Piezoresistive Coefficient)横向压阻系数(Transverse Piezoresistive Coefficient)压阻常数(Piezoresistive Constant)单晶硅(Monocrystalline Silicon)应变灵敏度(Strain Sensitivity)固态压阻式加速度传感器(Solid State Piezoresistive Accelerometer) 体型压阻式加速度传感器(Bulk Type Piezoresistive Accelerometer) 力平衡式传感器(Force Balance Transducer)电动力常数(Electro-dynamic Constant)机电耦合系统(Electro-mechanical Coupling System)2 检测仪表、激励设备及校准装置时间基准信号(Time Base Signal)李萨茹图(Lissojous Curve)数字频率计(Digital Frequency Meter)便携式测振表(Portable Vibrometer)有效值电压表(RMS V alue V oltmeter)峰值电压表(Peak-value V oltmeter)平均绝对值检波电路(Average Absolute Value Detector)峰值检波电路(Peak-value Detector)准有效值检波电路(Quasi RMS V alue Detector)真有效值检波电路(True RMS Value Detector)直流数字电压表(DVM,DC Digital V oltmeter)数字式测振表(Digital Vibrometer)A/D 转换器(A/D Converter)D/A 转换器(D/A Converter)相位计(Phase Meter)电子记录仪(Lever Recorder)光线示波器(Oscillograph)振子(Galvonometer)磁带记录仪(Magnetic Tape Recorder)DR 方式(直接记录式) (Direct Recorder)FM 方式(频率调制式) (Frequency Modulation)失真度(Distortion)机械式激振器(Mechanical Exciter)机械式振动台(Mechanical Shaker)离心式激振器(Centrifugal Exciter)电动力式振动台(Electro-dynamic Shaker)电动力式激振器(Electro-dynamic Exciter)液压式振动台(Hydraulic Shaker)液压式激振器(Hydraulic Exciter)电液放大器(Electro-hydraulic Amplifier)磁吸式激振器(Magnetic Pulling Exciter)涡流式激振器(Eddy Current Exciter)压电激振片(Piezoelectric Exciting Elements)冲击力锤(Impact Hammer)冲击试验台(Shock Testing Machine)激振控制技术(Excitation Control Technique)波形再现(Wave Reproduction)压缩技术(Compression Technique)均衡技术(Equalization Technique)交越频率(Crossover Frequency)综合技术(Synthesis Technique)校准(Calibration)分部校准(Calibration for Components in system)系统校准(Calibration for Over-all System)模拟传感器(Simulated Transducer)静态校准(Static Calibration)简谐激励校准(Harmonic Excitation Calibration)绝对校准(Absolute Calibration)相对校准(Relative Calibration)比较校准(Comparison Calibration)标准振动台(Standard Vibration Exciter)读数显微镜法(Microscope-streak Method)光栅板法(Ronchi Ruling Method)光学干涉条纹计数法(Optical Interferometer Fringe Counting Method)光学干涉条纹消失法(Optical Interferometer Fringe Disappearance Method)背靠背安装(Back-to-back Mounting)互易校准法(Reciprocity Calibration)共振梁(Resonant Bar)冲击校准(Impact Exciting Calibration)摆锤冲击校准(Ballistic Pendulum Calibration)落锤冲击校准(Drop Test Calibration)振动和冲击标准(Vibration and Shock Standard)迈克尔逊干涉仪(Michelson Interferometer)摩尔干涉图象(Moire Fringe)参考传感器(Reference Transducer)3 频率分析及数字信号处理带通滤波器(Band-pass Filter)半功率带宽(Half-power Bandwidth)3 dB 带宽(3 dB Bandwidth)等效噪声带宽(Effective Noise Bandwidth)恒带宽(Constant Bandwidth)恒百分比带宽(Constant Percentage Bandwidth)1/N 倍频程滤波器(1/N Octave Filter)形状因子(Shape Factor)截止频率(Cut-off Frequency)中心频率(Centre Frequency)模拟滤波器(Analog Filter)数字滤波器(Digital Filter)跟踪滤波器(Tracking Filter)外差式频率分析仪(Heterodyne Frequency Analyzer) 逐级式频率分析仪(Stepped Frequency Analyzer)扫描式频率分析仪(Sweeping Filter Analyzer)混频器(Mixer)RC 平均(RC Averaging)平均时间(Averaging Time)扫描速度(Sweeping Speed)滤波器响应时间(Filter Response Time)离散傅里叶变换(DFT,Discrete Fourier Transform) 快速傅里叶变换(FFT,Fast Fourier Transform)抽样频率(Sampling Frequency)抽样间隔(Sampling Interval)抽样定理(Sampling Theorem)抗混滤波(Anti-aliasing Filter)泄漏(Leakage)加窗(Windowing)窗函数(Window Function)截断(Truncation)频率混淆(Frequency Aliasing)乃奎斯特频率(Nyquist Frequency)矩形窗(Rectangular Window)汉宁窗(Hanning Window)凯塞-贝塞尔窗(Kaiser-Bessel Window)平顶窗(Flat-top Window)平均(Averaging)线*平均(Linear Averaging)指数平均(Exponential Averaging)峰值保持平均(Peak-hold Averaging)时域平均(Time-domain Averaging)谱平均(Spectrum Averaging)重叠平均(Overlap Averaging)栅栏效应(Picket Fence Effect)吉卜斯效应(Gibbs Effect)基带频谱分析(Base-band Spectral Analysis)选带频谱分析(Band Selectable Sp4ctralAnalysis)细化(Zoom)数字移频(Digital Frequency Shift)抽样率缩减(Sampling Rate Reduction)功率谱估计(Power Spectrum Estimate)相关函数估计(Correlation Estimate)频响函数估计(Frequency Response Function Estimate) 相干函数估计(Coherence Function Estimate)冲激响应函数估计(Impulse Response Function Estimate) 倒频谱(Cepstrum)功率倒频谱(Power Cepstrum)幅值倒频谱(Amplitude Cepstrum)倒频率(Quefrency)4 旋转机械的振动测试及状态监测状态监测(Condition Monitoring)故障诊断(Fault Diagnosis)转子(Rotor)转手支承系统(Rotor-Support System)振动故障(Vibration Fault)轴振动(Shaft Vibration)径向振动(Radial Vibration)基频振动(Fundamental Frequency Vibration)基频检测(Fundamental Frequency Component Detecting) 键相信号(Key-phase Signal)正峰相位(+Peak Phase)高点(High Spot)光电传感器(Optical Transducer)同相分量(In-phase Component)正交分量(Quadrature Component)跟踪滤波(Tracking Filter)波德图(Bode Plot)极坐标图(Polar Plot)临界转速(Critical Speed)不平衡响应(Unbalance Response)残余振幅(Residual Amplitude)方位角(Attitude Angle)轴心轨迹(Shaft Centerline Orbit)正进动(Forward Precession)同步正进动(Synchronous Forward Precession)反进动(Backward Precession)正向涡动(Forward Whirl)反向涡动(Backward Whirl)油膜涡动(Oil Whirl)油膜振荡(Oil Whip)轴心平均位置(Average Shaft Centerline Position)复合探头(Dual Probe)振摆信号(Runout Signal)电学振摆(Electrical Runout)机械振摆(Mechanical Runout)慢滚动向量(Slow Roll Vector)振摆补偿(Runout Compensation)故障频率特征(Frequency Characteristics Of Fault)重力临界(Gravity Critical)对中(Alignment)双刚度转子(Dual Stiffness Rotor)啮合频率(Gear-mesh Frequency)间入简谐分量(Interharmonic Component)边带振动(Side-band Vibration)三维频谱图(Three Dimensional Spectral Plot)瀑布图(Waterfall Plot)级联图(Cascade Plot)阶次跟踪(Order Tracking)阶次跟踪倍乘器(Order Tracking Multiplier)监测系统(Monitoring System)适调放大器(Conditional Amplifier) 趋势分析(Trend Analysis)倒频谱分析(Cepstrum Analysis)直方图(Histogram)确认矩阵(Confirmation Matrix)通频幅值(Over-all Amplitude)幅值谱(Amplitude Spectrum)相位谱(Phase Spectrum)报警限(Alarm Level)。
001 (ISSCC tutorial)Noise Analysis in Switched-Capacitor Circuits

© 2011 IEEE
IEEE International Solid-State Circuits Conference
© 2011 IEEE
Thermal Noise Power
• Nyquist showed that
PSD ( f ) = 4kT
• The total average noise power of a resistor in a certain frequency band is therefore
– Examples: Audio systems, wireless transceivers, sensor interfaces
• Electronic noise directly trades with power dissipation and speed • Electronic noise is a major concern in modern technologies with reduced VDD
• The noise of a MOSFET operating in the triode region is approximately equal to that of a resistor • In the saturation region, the thermal noise can be modeled using a drain current source with power spectral density
• We can model the noise using an equivalent voltage or current generator
2 vn
= Pn ⋅ R = 4kT ⋅ R ⋅ Δf
Introduction to HF Radio Propagation

Introduction to HF Radio Propagation1. The Ionosphere1.1 The Regions of the IonosphereIn a region extending from a height of about 50 km to over 500 km, some of the molecules of the atmosphere are ionised by radiation from the Sun to produce an ionised gas. This region is called the ionosphere, figure 1.1.Ionisation is the process in which electrons, which are negatively charged, are removed from (or attached to) neutral atoms or molecules to form positively (or negatively) charged ions and free electrons. It is the ions that give their name to the ionosphere, but it is the much lighter and more freely moving electrons which are important in terms of high frequency (HF: 3 to 30 MHz) radio propagation. Generally, the greater the number of electrons, the higher the frequencies that can be used.During the day there may be four regions present called the D, E, F1 and F2 regions. Their approximate height ranges are:• D region 50 to 90 km;• E region 90 to 140 km;•F1 region 140 to 210 km;•F2 region over 210 km.During the daytime, sporadic E (section 1.6) is sometimes observed in the E region, and at certain times during the solar cycle the F1 region may not be distinct from the F2 region but merge to form an F region. At night the D, E and F1 regions become very much depleted of free electrons, leaving only the F2 region available for communications; however it is not uncommon for sporadic E to occur at night.Only the E, F1, sporadic E when present, and F2 regions refract HF waves. The D region is important though, because while it does not refract HF radio waves, it does absorb or attenuate them (section 1.5).The F2 region is the most important region for high frequency radio propagation as:•it is present 24 hours of the day;•its high altitude allows the longest communication paths;•it usually refracts the highest frequencies in the HF range.The lifetime of electrons is greatest in the F2 region which is one reason why it is present at night. Typical lifetimes of electrons in the E, F1 and F2 regions are 20 seconds, 1 minute and 20 minutes, respectively.Because the F1 region is not always present and often merges with the F2 region, it is not normally considered when examining possible modes of propagation. Throughout this report, discussion of the F region refers to the F2 region.Figure 1.1 Day and night structure of the ionosphere.1.2 Production and Loss of ElectronsRadiation from the Sun causes ionisation in the ionosphere. Electrons are produced when this radiation collides with uncharged atoms and molecules, figure 1.2. Since this process requires solar radiation, production of electrons only occurs in the daylight hemisphere of the ionosphere.When a free electron combines with a charged ion a neutral particle is usually formed,figure 1.3. Essentially, loss is the opposite process to production. Loss of electrons occurs continually, both day and night.500400300200100Figure 1.2 Production.Figure 1.3 Loss.1.3 Observing the IonosphereThe most important feature of the ionosphere in terms of radio communications is its ability to refract radio waves. However, only those waves within a certain frequency range will be refracted. The range of frequencies refracted depends on a number of factors (section 1.4). Various methods have been used to investigate the ionosphere,and the most widely used instrument for this purpose is the ionosonde, figure 1.4. Note that many references to ionospheric communications speak of reflection of the wave. It is, however, a refraction process.An ionosonde is a high frequency radar which sends very short pulses of radio energy vertically into the ionosphere. If the radio frequency is not too high, the pulses are refracted back towards the ground. The ionosonde records the time delay between transmission and reception of the pulses. By varying the oscillation frequency of the pulses, a record is obtained of the time delay at different frequencies.Frequencies less than about 1.6 MHz are interfered with by AM broadcast stations. As the frequency is increased, echoes appear first from the lower E region and subsequently, with greater time delay, from the F1 and F2 regions. Of course, at night echoes are returned only from the F2 region and possibly sporadic E since the other regions have lost most of their free electrons.Positively charged ionFree electronPositively charged ionFigure 1.4 Ionosonde operation.Today, the ionosphere is "sounded" not only by signals sent up at vertical incidence. Oblique sounders send pulses of radio energy obliquely into the ionosphere (the transmitter and receiver are separated by some distance). This type of sounder can monitor propagation on a particular circuit and observations of the various modes being supported by the ionosphere can be made. Backscatter ionosondes rely on echoes reflected from the ground and returned to the receiver, which may or may not be at the same site as the transmitter. This type of sounder is used for over-the-horizon radar. 1.4 Ionospheric VariationsThe ionosphere is not a stable medium that allows the use of one frequency over the year, or even over 24 hours. The ionosphere varies with the solar cycle, the seasons, the circuit and during any given day. So, a frequency which may provide successful propagation now, may not do so an hour later.1.4.1 Variations due to the Solar CycleThe Sun goes through a periodic rise and fall in activity which affects HF communications; solar cycles vary in length from 9 to 14 years. At solar minimum, only the lower frequencies of the HF band will be supported by the ionosphere, while at solar maximum the higher frequencies will successfully propagate, figure 1.5. This is because there is more radiation being emitted from the Sun at solar maximum, producing more electrons in the ionosphere which allows the use of higher frequencies.Figure 1.5 The relationship between solar cycles and E and F region frequencies at Townsville. Dotted vertical lines indicate start of each year. Note also theseasonal variationsThere are other consequences of the solar cycle. Around solar maximum there is a greater likelihood of large solar flares occurring. Flares are huge explosions on the Sun which emit radiation that ionises the D region causing increased absorption of HF waves. Since the D region is present only during the day, only those communication paths which pass through daylight will be affected. The absorption of HF waves travelling via the ionosphere after a flare has occurred is called a short wave fade-out (section 3.1). Fade-outs occur instantaneously and affect lower frequencies the most. Lower frequencies are also the last to recover. If it is suspected or confirmed that a fade-out has occurred, it is advisable to try using a higher frequency, if possible. The duration of fade-outs can vary between about 10 minutes to over an hour depending on the duration and intensity of the flare.1.4.2 Seasonal VariationsE region frequencies are greater in summer than winter, figure 1.5. However, the variation inF region frequencies is more complicated. In both hemispheres, F region noon frequencies generally peak around the equinoxes (March and September). Around solar minimum the summer noon frequencies are, as expected, generally greater than those in winter, but around solar maximum winter frequencies tend to be higher than those in summer. In addition, frequencies around the equinoxes (March and September) are higher than those in summer or winter for both solar maximum and1.4.3 ° latitude at H i g h e s t u s a b l e f r e q u e n c y (M H z )1.4.4 Daily VariationsFrequencies are normally higher during the day and lower at night, Figure 1.7. With dawn, solar radiation causes electrons to be produced in the ionosphere and frequencies increase reaching their maximum around noon. During the afternoon,frequencies begin falling due to electron loss and with darkness the D, E and F1regions disappear. So, communication during the night is by the F2 region and absorption of radio waves is lower. Through the night, frequencies gradually decrease,reaching their minimum just before dawn.Figure 1.7 E and F layer frequencies for a Singapore to Ho Chi Minh circuit sometime in a solar cycle.1.5 Variations in AbsorptionAbsorption was discussed in section 1.4.1 when describing how solar flares can cause disruptions or degradations to communication paths which pass through daylight.Absorption in the D region also varies with the solar cycle, being greatest around solar maximum. Signal absorption is also greater in summer and during the middle of the day, figure 1.8. There is a variation in absorption with latitude, with more absorption occurring near the equator and decreasing towards the poles. Lower frequencies are absorbed to a larger extent, so it is advisable to use as high a frequency as possible,particularly during the day when absorption is greatest.1012024********1618202224Time (UT)Frequency(MHz)Local time (hours)V e r t i c a l a b s o r p t i o n (d B )Figure 1.8 Daily and seasonal variations in absorption at Sydney, 2.2 MHz.Around the polar regions absorption can affect communications quite dramatically at times. Sometimes high energy protons ejected from the Sun during large solar flares will move down the Earth's magnetic field lines and into the polar regions. These protons can cause increased absorption of HF radio waves as they pass through the D region. This increased absorption may last for as long as 10 days and is called a Polar Cap Absorption event (PCA), section 3.2.1.6 Sporadic ESporadic E may form at any time during the day or night. It occurs at altitudes between 90 to 140 km (the E region), and may be spread over a large area or be confined to a small region. It is difficult to know where and when it will occur and how long it will persist. Sporadic E can have a comparable electron density to the F region. This implies that it can refract comparable frequencies to the F region. Sometimes a sporadic E layer is transparent and allows most of the radio wave to pass through it to the F region, however, at other times the sporadic E layer obscures the F region totally and the signal does not reach the receiver (sporadic E blanketing). If the sporadic E layer is partially transparent, the radio wave is likely to be refracted at times from the F region and at other times from the sporadic E. This may lead to partial transmission of the signal or fading, figure 1.9.Sporadic E in the low and mid-latitudes occurs mostly during the daytime and early evening, and is more prevalent during the summer months. At high latitudes, sporadic E tends to form at night.Signal partially transmitted by sporadic E and partially by the Flayer. Any received signal may be weak or fade.Wave passes through sporadic E and is refracted by the F layer. Figure 1.9 Some possible paths when sporadic E is present. The ground reflection will depend on the strength of the signal and other factors such as groundtype.1.7 Spread FSpread F occurs when the F region becomes diffuse due to irregularities in that region, which scatter the radio wave. The received signal is the superposition of a number of waves refracted from different heights and locations in the ionosphere at slightly different times. At low latitudes, spread F occurs mostly during the night hours and around the equinoxes. At mid-latitudes, spread F is less likely to occur than at low and high latitudes. Here it is more likely to occur at night and in winter. At latitudes greater than about 40°, spread F tends to be a night time phenomenon, appearing mostly around the equinoxes, while around the magnetic poles, spread F is often observed both day and night. At all latitudes there is a tendency for spread F to occur when there is a decrease in F region frequencies. That is, spread F is often associated with ionospheric storms (section 3.3).2. HF Communications2.1 Types of HF PropagationHigh Frequency (3 to 30 MHz) radio signals can propagate to a distant receiver, Figure 2.1, via the:•ground wave: near the ground for short distances, up to 100 km over land and 300 km over sea. Attenuation of the wave depends on antenna height, polarisation, frequency, ground types, terrain and/or sea state;•direct or line-of-sight wave: this wave may interact with the earth-reflected wave depending on terminal separation, frequency and polarisation;•sky wave: refracted by the ionosphere, all distances.2.2 Frequency Limits of Sky WavesNot all HF waves are refracted by the ionosphere, there are upper and lower frequency bounds for communications between two terminals. If the frequency is too high, the wave will penetrate the ionosphere, if it is too low, the strength of the signal will be lowered due to absorption in the D region. The range of usable frequencies will vary:•throughout the day;•with the seasons;•with the solar cycle;•from place to place;depending on the ionospheric region used for communications. While the upper limit of frequencies varies mostly with these factors, the lower limit is also dependent on receiver site noise, antenna efficiency, transmitter power, E layer screening (section 2.6) and absorption by the ionosphere.Figure 2.1 Types of HF propagation.2.3 The Usable Frequency RangeFor any circuit there is a Maximum Usable Frequency (MUF) which is determined by the state of the ionosphere in the vicinity of the refraction area(s) and the length of the circuit. The MUF is refracted from the area of maximum electron density of a region. Therefore, frequencies higher than the MUF for a particular region will penetrate thatregion. During the day it is possible to communicate via both the E and F layers usingdifferent frequencies. The highest frequency supported by the E layer is the EMUF, while that supported by the F layer is the FMUF.The F region MUF in particular varies during the day, seasonally and with the solar cycle. The data collected over the years displays a range of frequencies observed and the IPS predictions mirror this. A range of F region MUFs is provided in the predictions and this range extends from the lower decile MUF (called the Optimum Working Frequency, OWF), through the median MUF to the upper decile MUF. These MUFs have a 90%, 50% and 10% chance of being supported by the ionosphere, respectively. IPS predictions usually cover a period of one month, so the OWF should provide successful propagation 90% of the time or 27 days of the month. The median MUF should provide communications 50% or 15 days of the month and the upper decile MUF 10% or 3 days of the month. The upper decile MUF is the highest frequency of the range of MUFs and is most likely to penetrate the ionosphere, figure 2.2.Figure 2.2 Range of usable frequencies. If the frequency, f, is close to the ALF then the wave may suffer absorption in the D region. If the frequency is abovethe EMUF then propagation is via the F region. Above the FMUF the wave is likely to penetrate the ionosphere.The chances of successful propagation discussed above rely on the monthly prediction of solar activity being correct. Sometimes unforeseen events occur on the Sun resulting in the monthly predictions being inaccurate. The role of the Australian Space Forecast Centre (ASFC) at IPS is to provide corrections to the monthly predictions, warning customers of changes in communication conditions.The D region does not allow all frequencies to be used since the lower the frequency the more likely it is to be absorbed. The Absorption Limiting Frequency (ALF) is provided as a guide to the lower limit of the usable frequency band. The ALF is significant only for circuits with refraction points in the sunlit hemisphere. At night, the ALF falls to zero, allowing frequencies which are not usable during the day to successfully propagate.2.4 Hop LengthsThe hop length is the ground distance covered by a radio signal after it has been refracted once from the ionosphere and returned to Earth, figure 2.3. The upper limit of the hop length is set by the height of the ionosphere and the curvature of the Earth. For E and F region heights of 100 km and 300 km, the maximum hop lengths with an elevation angle of 4°, are 1800 km and 3200 km, respectively. Distances greater than these will require more than one hop. For example, a distance of 6100 km will require a minimum of 4 hops by the E region and 2 hops via the F region. More hops would be required with larger antenna elevation angles.Figure 2.3 Hop lengths based upon an antenna elevation angle of 4° and heights forthe E and F layers of 100 km and 300 km, respectively.If the elevation angle of the antenna is able to be altered easily or there is a choice of antennas, it may be possible to use a certain number of hops such that an opponent's position does not coincide with a reflection point on the Earth. This will however, be affected also by the time of day and what frequencies are available. However, in most cases the directional properties of the antenna are such that this is not possible.2.5 Propagation ModesThere are many paths by which a sky wave may travel from a transmitter to a receiver.The mode by a particular layer which requires the least number of hops between the transmitter and receiver is called the first order mode. The mode that requires one extra hop is called the second order mode. For a circuit with a path length of 5000 km, the first order F mode has two hops (2F), while the second order F mode has three hops (3F). The first order E mode has the same number of hops as the first order F mode. If this results in a hop length of greater than 2050 km, which corresponds to an elevation angle of 0°, then the E mode is not possible. This also applies to the second order E mode. Of course, the E region modes will only be available on daylight circuits.Simple modes are those propagated by one region, say the F region. IPS predictions are made only for these simple modes, Figure 2.4. More complicated modes consisting of combinations of refractions from the E and F regions, ducting and chordal modes are also possible, figure 2.5.E region(100 km)F region(300 km)Figure 2.4 Examples of simple propagation modes.Chordal modes and ducting involve a number of refractions from the ionosphere without intermediate reflections from the Earth. There is a tendency to think of the regions of the ionosphere as being smooth, however, the ionosphere undulates and moves, with waves passing through it which may affect the refraction of the signal. The ionospheric regions may tilt and when this happens chordal and ducted modes may occur. Ionospheric tilting is more likely near the equatorial anomaly, the mid-latitude trough and in the sunrise and sunset sectors. When these types of modes do occur,signals can be strong since the wave spends less time traversing the D region and being attenuated during ground reflections.Figure 2.5 Other propagation modes.Because of the high electron density of the daytime ionosphere in the vicinity of 15° of the magnetic equator (near the equatorial anomaly), transequatorial paths can use these enhancements to propagate on higher frequencies. Any tilting of the ionosphere may result in chordal modes, producing good signal strength over long distances.Ducting may result if tilting occurs and the wave becomes trapped between refracting regions of the ionosphere. This is most likely to occur in the equatorial ionosphere, near the auroral zone and mid-latitude trough. Disturbances to the ionosphere, such as travelling ionospheric disturbances (section 2.9), may also account for ducting and chordal mode propagation.E regionF region2.6 E Layer ScreeningFor daytime communications via the F region, the lowest usable frequency via the one hop F mode is dependent upon the presence of the E region. If the operating frequency for the 1F mode is below the two hop EMUF, then the signal is unlikely to propagate via the F region due to screening by the E region, figure 2.6. This is also because the antenna elevation angles of the 1F and 2E modes are similar.Figure 2.6 E layer screening occurs if communications are required by the 1F modeand the operating frequency is close to or below the EMUF for the 2Emode. Note the paths through the D region for each wave.A sporadic E layer may also screen a wave from the F region. Sometimes sporadic E can be quite transparent, allowing most of the wave to pass through it. At other times it will partially screen the F region leading to a weak or fading signal, while at other times sporadic E can totally obscure the F region with the possible result that the signal does not arrive at the receiver, figure 1.9 (section 1.6).2.7 Frequency, Range and Elevation AngleFor oblique propagation, there are three dependent variables:• frequency;• range or path length;• antenna elevation angle.The diagrams below illustrate the changes to the ray paths when each of these is fixed in turn.E regionF regionD regionFigure 2.7: Elevation angle fixed:• As the frequency is increased toward the MUF, the wave is refracted higher in the ionosphere and the range increases, paths 1 and 2;• At the MUF for that elevation angle, the maximum range is reached, path 3;• Above the MUF, the wave penetrates the ionosphere, path 4.Figure 2.7 Elevation angle fixed.Figure 2.8: Path length fixed (point-to-point circuit):• As the frequency is increased towards the MUF, the wave is refracted from higher in the ionosphere. To maintain a circuit of fixed length, the elevation angle must therefore be increased, paths 1 and 2;• At the MUF, the critical elevation angle is reached, path 3. The critical elevation angle is the elevation angle for a particular frequency, which if increased, would cause penetration of the ionosphere;• Above the MUF, the ray penetrates the ionosphere, path 4.Figure 2.8 Path length fixed.Figure 2.9: Frequency fixed:•At low elevation angles the path length is greatest, path 1;•As the elevation angle is increased, the path length decreases and the ray is refracted from higher in the ionosphere, paths 2 and 3;•If the frequency will return when sent vertically up into the ionosphere, then the skip distance is zero. However, if this is not the case, then as the elevation angle is increased, the range decreases. If the elevation angle is increased beyond the critical elevation angle for that frequency then the wave penetrates the ionosphere and there is an area around the transmitter within which no sky wave communications can be received, path 4. To communicate within the skip zone, the frequency must be lowered.Figure 2.9 Frequency fixed.2.8 Skip ZonesA propagation path will consist of high and low angle rays corresponding to the wave propagating from the transmit antenna at a range of angles. The high ray is transmitted at a high angle from the antenna. These rays may travel by different paths through the ionosphere. The edge of the skip zone corresponds to the high and low rays taking the same path through the ionosphere, with the resulting signal often being stronger. Within the skip zone the signal fades due to the waves penetrating the ionosphere.In figure 2.9, path 3 corresponds to the high and low rays taking the same path through the ionosphere, this corresponds to the MUF. As the frequency is increased toward the MUF, the height of refraction of the low angle ray increases and the height of refraction of the high angle ray decreases until they are both refracted at the same point in the ionosphere.Skip zones can often be used to advantage if it is desired that communications are not heard by a particular receiver. Selecting a different frequency will alter the size of the skip zone and if the receiver is within the skip zone and out of reach of the ground wave, then it is unlikely that it will receive the communications. However, factors such as sidescatter, where reflection from the Earth outside the skip zone results in the wavetransmitting into the skip zone may affect the reliability of this.Skip zones vary in size during the day, with the seasons, and with solar activity. During the day, solar maximum and around the equinoxes, skip zones generally are smaller in area. The ionosphere during these times has increased electron density and so is able to support higher frequencies.2.9 FadingMultipath fading results from dispersion of the signal by the transmitting antenna. A number of modes propagate which have variations in phase and amplitude. These waves may interfere with each other if they reach the receiver, figure 2.10.Figure 2.10 Multipath fading. The signal may travel by a number of paths which, if they arrive at the receiver and are of similar amplitude, may interfere and cause fading.Disturbances known as Travelling Ionospheric Disturbances, TIDs, may cause a region to be tilted, resulting in the signal being focussed or defocussed, figure 2.11. Fading periods of the order of 10 minutes or more can be associated with these structures. TIDs travel horizontally at 5 to 10 km/minute with a well defined direction of travel affecting higher frequencies first. Some originate in auroral zones following an event on the Sun and these may travel large distances. Others originate in weather disturbances. TIDs may cause variations in phase, amplitude, polarisation and angle of arrival of a wave.Polarisation fading results from changes to the polarisation of the wave along the propagation path. The receiving antenna is unable to receive parts of the signal; this type of fading can last for a fraction of a second to a few seconds.Skip fading can be observed around sunrise and sunset particularly, when the operating frequency is close to the MUF, or when the receiving antenna is positioned close to the boundary of the skip zone. At these times of the day, the ionosphere is unstable and the frequency may oscillate above and below the MUF causing the signal to fade in and out. If the receiver site is close to the skip zone boundary, as the ionosphere fluctuates, the skip zone boundary also fluctuates.Figure 2.11 Focussing and defocussing effects caused by tilting and travelling ionospheric disturbances (TIDs).2.10 NoiseRadio noise arises from internal and external origins. Internal or thermal noise is generated in the receiving system and is usually negligent when compared to external sources of noise. External radio noise originates from natural (atmospheric and galactic) and man-made (environmental) sources.Atmospheric noise, which is caused by thunderstorms, is normally the major contributor to radio noise in the HF band and will especially degrade circuits passing through the day-night terminator. Atmospheric noise is greatest in the equatorial regions of the world and decreases with increasing latitude. Its effect is also greater on lower frequencies, hence it is usually more of a problem around solar minimum and at night when lower frequencies are needed.Galactic noise arises from our galaxy. Receive antennas with higher frequencies more likely to be affected by this type of noise.Man-made noise includes ignition noise, neon signs, electrical cables, power transmission lines and welding machines. This type of noise depends on the technology used by the society and its population. Interference may be intentional, such as jamming, due to propagation conditions or the result of others working on the same frequency.Man-made noise tends to be vertically polarised, so selecting a horizontally polarised antenna may help in reducing noise. Using a narrower bandwidth, or a directional receiving antenna (with a lobe in the direction of the transmitting source and a null in the direction of the unwanted noise source), will also aid in reducing the effects of noise. Selecting a site with a low noise level and determining the major noise sources are important factors in establishing a successful communications system.2.11 VHF and 27 MHz PropagationVHF and 27 MHz are used for line-of-sight or direct wave communication, for example, ship-to-ship or ship-to-shore. The frequency bands are divided into channels and one channel is usually as good as the next. This is in contrast to medium frequency (MF: 300 kHz to 3 MHz) and HF where the choice of a frequency channel may be crucial for good communications.Because VHF and 27 MHz operate mainly by line-of-sight, it is important to mount the antenna as high as possible and free from obstructions. Shore stations are usually on the tops of hills to provide maximum range, but even the highest hills do not provide。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
With recent advances in low power wireless sensing technologies, the concept of BSN has shown significant strength in continuous monitoring of patients under their natural physiological status [1]. Due to the diversity of the environment and physiological conditions they may experience, understanding of the context within which the signals are collected plays an important role for the accurate prediction of adverse events. Reliable detection of patient activity, however, requires the use of a large number of context sensors around the body. This can potentially introduce a significant burden to power consumption and bandwidth utilisation. Directly sending all the sensory data to the centralised processing unit requires extensive transmission power and a high bandwidth at the central processing unit. Decreasing the transmission range and required bandwidth will greatly reduce the power consumption and prolong the life span of the sensors. Clustering data transmission among neighbourhood can also alleviate the problem of data collision.
5) Arc list
Current network
The arc A-C is not
added since B-C are A B C D E
barely correlated
6) Adding hidden nodes
R
H
H
A BCDE
Fig. 1. A diagrammatic overview of the construction of a BN with hidden nodes with complete-linkage clustering, where only the top arc is considered in each step.
D
D – decision (root) node
H – Hidden node
S1-S8 – Sensor nodes
H
H
S1
S2
S3
S4
S5
S6
S7
S8
Hidden nodes indicate redundant sensors Fig. 2. The structure of the BN with hidden nodes learned from the training dataset.
CurC,D)
(A,C)
Current network A B CDE
4) Arc list (A,C)
Current network
The arc C-D is not
added since C-E are A B C D E
not correlated
1) Arc list (A,B)
(D,E) (C,D) (A,C) (A,E) (B,D) (B,E) (A,D) (B,C) (C,E)
Current network A B CDE
Filter out arcs with low correlations
2) Arc list (D,E)
(C,D) (A,C)
To assess the overall performance of the proposed method, an ETH reference dataset was used. This consisted of sensor data obtained from accelerometers in performing eight different activities [2]. After extracting the temporal features from the dataset, reference data is constructed by selecting six representative features and two highly correlated features. From the reference data set, a BN is obtained by learning the structure from the training dataset, and hidden nodes were inserted to represent the dependency among correlated children [4], as shown in Fig.2.
13
2. Model Description
The main idea of introducing noise resilience to the sensing architecture is the introduction of hidden nodes to a traditional BN to resolve the violation of the independence assumption when there is a considerable amount of dependency among child nodes. This allows the transformation of a multiply connected network into a singly connected network. Complete-link clustering based on a pairwise-dependency measure, such as Pearson’s correlation coefficient or L1 metric, is used to form clusters of correlated child nodes as shown in Fig. 1. In practice, a close-form representation for deriving link matrices is difficult for a network with hidden nodes. This is because the hidden nodes are unobserved variables and the node statistics cannot be directly calculated. We used instead a backward propagation method, proposed by Kwoh et. al, is used for parameter learning [3].
* Department of Computing, Imperial College London, South Kensington Campus 180 Queen’s Gate, London, SW7 2AZ, United Kingdom. {st01, rck, benlo, dfg, gzy}@ /
A NOISE RESILIENT DISTRIBUTED INFERENCE FRAMEWORK FOR BODY SENSOR NETWORKS
Surapa Thiemjarus, Rachel King, Benny Lo, Duncan Gilles, Guang-Zhong Yang*
Abstract
Existing research has shown that inferencing with message passing is potentially useful for distributed sensing systems. However, resilience to communication error and node failures is major obstacle to overcome. Paskin et. al proposed a robust message passing algorithm for reasoning in a junction tree model [5]. By converting a standard multiply connected Bayesian Network (BN) [6] into a cluster tree and combining nodes into a clique, the problem of non-convergence and incorrect update of the posterior probabilities due to the loopy feedback in a multiply connected model can be avoided. Furthermore, existing research in BSN has also highlighted need for built-in redundancies in the sensor network for dealing with motion artefact and node failures. The purpose of this paper is to propose a Bayesian framework that permits distributed inferencing with a high level of noise tolerance.