Research on Modeling and Compensating Method of Random Drift of MEMS Gyroscope
|
|
- Kathlyn Marsh
- 5 years ago
- Views:
Transcription
1 01 4th Interntionl Conference on Signl Processing Systems (ICSPS 01) IPCSIT vol. 58 (01) (01) IACSIT Press, Singpore DOI: /IPCSIT.01.V58.9 Reserch on Modeling nd Compensting Method of Rndom Drift of MEMS Gyroscope Yun Li +, Fengyng Dun nd Znping Li Deprtment of Control Engineering, Avition University, Chngchun, Chin Abstrct. Aiming t the performnce of Micro Electro Mech-nicl Systems (MEMS) gyroscope influenced by rndom drift lrgely, its error compensting method is nlyzed. Bsed on the time series nlysis method, the rw mesurement from the MEMS gyroscope is processed nd modeled by the Autoregressive (AR) model. With the AR model, time-vrying fctor dptive Klmn filter (AKF) is designed. The results show tht, fter filtering in the cse of constnt nd chnging ngulr rte, the stndrd devitions re 16% nd 70% of tht before filtering, nd the drift error cn be reduced by the filter effectively. Keywords: MEMS, rndom drift, AR model, AKF. 1. Introduction MEMS gyroscopes, with the dvntges of smll size, light weight nd low power consumption, re more nd more widely pplicted in low cost inertil system. As the precision level of MEMS gyroscopes re not high till now, even in low precision mesurement system, the rndom drift estimtion nd compenstion of gyroscopes re quite necessry. And the min wys of improving the system precision re the model identifiction nd filtering of the MEMS gyroscopes rndom drift [1]. The common methods of modeling of gyroscopes rndom drift lwys get high model order nd lrge mount of computtion, such s neurl networ [] nd wvelet nlysis [3], which pplicted in low cost system cn hrdly ensure rel-time property. Bsed on the time series nlysis theory, n AR model of gyroscope rndom drift is estblished. Then we process the gyroscope rndom drift with time-vrying fctor dptive Klmn filter, nd the method, which is helpful in improving rel-time property of the low cost system, is efficient with low computtion.. Error Model Modeling bsed on the time series nlysis theory requires collecting long-term mesurements, nd ensure tht the signl is zero mens vlue, sttionry nd norml [4]..1. Dt cquisition nd processing In this pper, we te STIM0 gyroscope [5] s reserch object. Under the condition of constnt temperture 5 C,we recorded 1 hours mesurements of gyroscope t uniform intervls of two seconds while it is sttic nd processed the output of Z-Axis. The rw dt contins constnt nd rndom components. So fter the men of it hs been removed, the rest is the signl of gyro rndom drift. During the long-term of smpling, the chnge of electricl prmeters nd outside environment could cuse the signl to generte trend, so it is necessry to te processing to ensure sttionry. + Corresponding uthor. E-mil ddress: @qq.com. 48
2 The liner trend in the signl cn be removed by the first or second order differentil tretment [6], [7]. Fig. 1 shows the signl fter it hs been processed. Fig.1. Processed dt of gyro rndom drift... Model Estblishment After processed, the sttic dt from the MEMS gyroscope turns into rndom error signl. The result of signl testing through reverse order testing method nd its third-order nd forth-order cumulnts show tht the signl hs stisfied modeling of time series nlysis theory requirements. Auto-Regressive nd Moving Averge (ARMA) Model [8] is bsic nd widely used time series model, nd it is defined s X ϕ X ϕ X ϕ X t 1 t 1 t m t m = θ θ θ t 1 t 1 t n t n X t is the time series { X t } t timet, m nd n re the orders of AR model nd Moving Averge (MA) model, ϕ i ( i= 1,,, m) θ j ( j = 1,,, n) is the AR prmeter, is the MA prmeter nd { } t is the stndrd devition of the sensor white noise. θ j = 0( j = 1,,, n) When, ARMA model reduces to AR model s Xt = t + ϕ1xt 1+ ϕxt + + ϕnxt n () ARMA model nd AR model hve different chrcteri-stics in uto-correltion function (ACF) nd prtil uto-correltion function (PACF). In AR model, the ACF is tiled nd the PACF is truncted, while in ARMA model both of ACF nd PACF re tiled. Therefore, the model of signl cn be identified by the chrcteristics of its ACF nd PACF. The ACF nd PACF of gyro rndom drift signl re presented s figure, which shows tht the PACF cuts off wheres the ACF hs jgged exponentil decy. Therefore the AR model is dopted to model the gyro rndom drift signl. After the prmeters of ech low order model hve been estimted by lest-squre method, the order determintion test of the AR model cn be conducted vi Aie Inform-tion Criterion (AIC) nd Byesin Informtion Criterion (BIC), which cn be normlly written s AIC( p) = N lnσ + p (3) BIC( p) = N lnσ + p ln N (4) p denotes the order of model, N denotes the number of Smpling points nd σ denotes is the vrince of residul t. And the test result is presented in Tble 1. As shown, by incresing the order of the AR model, AIC nd BIC vlues don t decrese remrbly. Although those vlues re the lower the better, the mount of computtion rises with the cubic of the order of AR model in the system. To further improve the system reltime property, the AR (1) model is used for the rndom drift error model. After the Model pln hs been finlized, it is necessry to perform the pplicbility test, testing whether { t } is white noise. By clculting, the one-step uto-correltion coeffic-ient of t is ρ,1 = nd the two-step cross-correltion coefficient between t nd X t is ρ X, = Both of the coefficients pproch zero, so { t } is considered to be white noise nd the model is pplicble. (1) 49
3 Fig.. Vlue of ACF nd PACF. Tble 1: Vlue of AIC nd BIC of Models ( 105) Model Criterio n AR(1) AR() AR(3) AR(4) AIC BIC The model is expressed in the form of X = ϕ X + t t 1 t NID σ (5) t ~ (0, ) 3. Filter Design After the drift error model hs been determined, we choose Klmn filer (KF) [8] to compenste the drift error. The liner stte eqution nd mesurement eqution re given s follows X = X +Γ W (6) φ, Z = HX + V (7) X is the system stte vector t time t φ,, 1 t is the stte trnsition mtrix between time point 1 nd t, Γ 1 H is the system noise distribution mtrix, is mesurement observtion mtrix, W 1 nd V re the rndom error of system nd mesurement nd they re per the following error sttistics Q EW [ ] = 0, CovW [, Wj] = Qδ j (8) E[ V] = 0, Cov[ V, Vj] = Rδ j (9) Cov[ W, V j] = 0 (10) nd R re vrinces of W nd V. The regulr KF estblished on this bsis cn effectively reduce the rndom drift error of gyro in the sttic or constnt ngulr rte stte, but in the stte of chnging ngulr rte the KF doesn t wor well. Through the nlysis we lerned tht the filter model chnged while the ngulr rte ws chnging. The more drmtic chnges in ngulr velocity, the pplicbility of the filter model is poorer, cuses the greter system error. But the direct reson for this is tht the old mesurement vlue influences the estimte of the current too much. Incresing the gin of current mesurement vlue in the stte estimtion eqution is n effective men to solve this problem. Bsed on the drift model before, AR(1) model tes the plce of KF stte eqution nd vrince of residul Q tes the plce of σ. H is set by 1. The AKF multiple updte steps my be performed One-step predicted stte estimte Updted stte estimte Filter gin Xˆ = ϕ Xˆ (11) / 1 1 X ˆ = X ˆ + K ( Z X ˆ ) (1) / 1 / 1 K P P R 1 = / 1 ( / 1+ ) (13) 50
4 One-step predicted estimte covrince Updted estimte covrince / 1 = sϕ P 1 + σ P (14) P = ( I K) P/ (15) 1 s s A fctor is dded into recursive eqution of regulr KF. The vlue selection of is not only to be ble to reduce the error under the condition of chnging ngulr rte, but lso to ensure the ccurcy in the constnt ngulr rte stte, nd expressed s follows ω0 s = (16) ω0 Δω Δω ω is the first order differentil of gyroscope mesurements t the current moment, 0 is the reference vlue obtined by experiments. Δω In the condition of constnt ngulr rte, pproches zero nd s pproches 1, the AKF doesn t chnge too much from KF. When the ngulr rte is chnging, s will djust the filter. 4. Simultion nd Experiment 4.1. Constnt ngulr rte experiment In order to verify the effectiveness of the model nd filter, use both regulr KF nd AKF methods to process the sttic dt t first. We set ϕ = nd X ˆ 0 = 0, P 0 [10] nd R [11] re set by 10 times of stndrd devition of rw dt nd one tenth of vrince. Filtering effect is shown in figure 3 nd figure 4, nd the error reduces obviously. The stndrd devition of error before nd fter filtering is shown in tble. The result shows tht the filter nd model estblished in the sttic condition re resonble. Compred with the regulr KF, the ccurcy of AKF doesn t reduce. Under the condition constnt ngulr rte, we get the sme result s it in the sttic condition. After processed by the two inds of filter, the stndrd devition of error is reduced to 16% of the rw dt s. Becuse of the limited spce, the detil description of the experiment is omitted. 4.. Chnging ngulr rte experiment In the sme environment condition, we te out prt of rw chnging ngulr rte dt (figure 5) from the MEMS gyroscope for experiment. Filtering effect is shown in figure 6 nd figure 7. The result shows tht filtering chnging ngulr rte dt with the regulr KF get poor effect nd the error vlue is chnging with the ngulr rte chnging, while the error of AKF isn t influenced by ngulr rte chnging lmost. The stndrd devition of error before nd fter filtering in the tble 3 shows tht AKF cn reduce the rndom drift error under the condition of chnging ngulr rte effectively. From bove, it shows tht bsed on the AR model before, regulr KF just cn reduce the rndom error in the cse of constnt ngulr rte, while AKF cn not only ensure the ccurcy in the constnt ngulr rte condition, but lso improve the filter effect when ngulr rte is chnging. Thus it cn be seen tht this lgorithm wored well. Fig.3. Klmn filter result of sttic experiment. Fig.4. Adptive Klmn filter result of sttic experiment. 51
5 Tble : Stndrd Devition of Error In Sttic Experiment Rw Dt Regulr KF AKF Stndrd devition (deg/s) Fig.5. Chnging ngulr rte dt of gyro. Fig.6. Klmn filter result of chnging ngulr rte experiment. 5. Conclusion Fig.7. Adptive Klmn filter result of chnging ngulr rte experiment. Tble 3: Stndrd Devition of Error In Chnging Angulr Rte Experiment Rw Dt Regulr KF AKF Stndrd devition (deg/s) In this pper, we nlysis the MEMS gyroscope output rw dte, extrct the rndom drift signl bsed on requirements of time series nlysis theory nd estblish n AR(1) model for filter. The results of experiments show tht in the constnt ngulr rte cse, both regulr KF nd AKF hve good performnces, but when the ngulr rte is chnging regulr KF loses efficcy, while AKF still performs well nd the stndrd devition of drift error reduces by 30%. The AKF lgorithm bsed on AR model, which is simple, effective nd high rel-time property, hs some prcticl vlue in MEMS gyroscope-bsed low cost integrted nvigtion system. 6. References [1] H. L. Chng, L. Xue, nd W. Qin, An integrted MEMS gyroscope rry with higher ccurcy output, Sensors, 008(8), pp
6 [] H. X. Lu, D. Z. Xi, nd B. L. Zhou, MEMS gyroscope's error modeling bsed on wvelet neurl networ of genetic lgorithms, Journl of Chinese Inertil Technology, 16(): 16-19, April 008. [3] J. H. Xing, L. Cong, nd H. L. Qin, Method for MEMS gyroscope dt processing bsed on wvelet denoising nd AR modeling, Computer Engineering nd Design, 010, 31(19): [4] A. Noureldin, T. B. Krmt, M. D. Eberts, nd A. El-Shfie. Performnce enhncement of MEMS-Bsed INS/GPS integrtion for low-cost nvigtion pplictions, IEEE Trnsctions On Vehiculr Technology, 58(3): , Mrch 009. [5] STIM0 Multi-Axis Gyro Module, NO: Sensonor Technologies AS, 010(3), pp. -4. [6] T. Meng, H. Wng, H. Li, G. H. He, nd Z. H. Jin, Error modeling nd filtering method for MEMS gyroscope, Systems Engineering nd Electronics. August 009, pp [7] J. Li, W. D. Zhng, nd J. Liu, Reserch on the ppliction of the time-seril nlysis bsed lmn filter in MEMS gyroscope rndom drift compenstion, Chinese Journl of Sensors nd Actutors, October 006, pp [8] S. Z. Yng, Y. Wu, nd J. P. Xun, Time series nlysis in engineering ppliction, Huzhong University of Science nd Technology Press, June. 007, pp [9] N. A. Crlson, Federted squre root filter for decentrlized prllel processes, IEEE Trnsctions on Aerospce nd Electronic System, 6(3):517-55, My [10] H. M. Qin, Q. X. Xi, X. T. Que, nd Q. Zhng, Algorithm for MEMS gyroscope bsed on Klmn filter, Journl of Hrbin Engineering University, 31(9):117-10, September 010. [11] X. S. Ji nd S. R. Wng, Reserch on the MEMS gyroscope rndom drift error, Journl of Astronutics, 7(4): , July
Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d
Interntionl Industril Informtics nd Computer Engineering Conference (IIICEC 15) Driving Cycle Construction of City Rod for Hybrid Bus Bsed on Mrkov Process Deng Pn1,, Fengchun Sun1,b*, Hongwen He1, c,
More informationPredict Global Earth Temperature using Linier Regression
Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationFig. 1. Open-Loop and Closed-Loop Systems with Plant Variations
ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses
More informationMonte Carlo method in solving numerical integration and differential equation
Monte Crlo method in solving numericl integrtion nd differentil eqution Ye Jin Chemistry Deprtment Duke University yj66@duke.edu Abstrct: Monte Crlo method is commonly used in rel physics problem. The
More informationChapter 9: Inferences based on Two samples: Confidence intervals and tests of hypotheses
Chpter 9: Inferences bsed on Two smples: Confidence intervls nd tests of hypotheses 9.1 The trget prmeter : difference between two popultion mens : difference between two popultion proportions : rtio of
More informationCBE 291b - Computation And Optimization For Engineers
The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn
More informationNon-Linear & Logistic Regression
Non-Liner & Logistic Regression If the sttistics re boring, then you've got the wrong numbers. Edwrd R. Tufte (Sttistics Professor, Yle University) Regression Anlyses When do we use these? PART 1: find
More informationCONTRIBUTION TO THE EXTENDED DYNAMIC PLANE SOURCE METHOD
CONTRIBUTION TO THE EXTENDED DYNAMIC PLANE SOURCE METHOD Svetozár Mlinrič Deprtment of Physics, Fculty of Nturl Sciences, Constntine the Philosopher University, Tr. A. Hlinku, SK-949 74 Nitr, Slovki Emil:
More informationMAC-solutions of the nonexistent solutions of mathematical physics
Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE
More informationSection 11.5 Estimation of difference of two proportions
ection.5 Estimtion of difference of two proportions As seen in estimtion of difference of two mens for nonnorml popultion bsed on lrge smple sizes, one cn use CLT in the pproximtion of the distribution
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationA Brief Review on Akkar, Sandikkaya and Bommer (ASB13) GMPE
Southwestern U.S. Ground Motion Chrcteriztion Senior Seismic Hzrd Anlysis Committee Level 3 Workshop #2 October 22-24, 2013 A Brief Review on Akkr, Sndikky nd Bommer (ASB13 GMPE Sinn Akkr Deprtment of
More informationA Signal-Level Fusion Model for Image-Based Change Detection in DARPA's Dynamic Database System
SPIE Aerosense 001 Conference on Signl Processing, Sensor Fusion, nd Trget Recognition X, April 16-0, Orlndo FL. (Minor errors in published version corrected.) A Signl-Level Fusion Model for Imge-Bsed
More informationConstruction and Selection of Single Sampling Quick Switching Variables System for given Control Limits Involving Minimum Sum of Risks
Construction nd Selection of Single Smpling Quick Switching Vribles System for given Control Limits Involving Minimum Sum of Risks Dr. D. SENHILKUMAR *1 R. GANESAN B. ESHA RAFFIE 1 Associte Professor,
More informationContinuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht
More informationA New Grey-rough Set Model Based on Interval-Valued Grey Sets
Proceedings of the 009 IEEE Interntionl Conference on Systems Mn nd Cybernetics Sn ntonio TX US - October 009 New Grey-rough Set Model sed on Intervl-Vlued Grey Sets Wu Shunxing Deprtment of utomtion Ximen
More information5 Probability densities
5 Probbility densities 5. Continuous rndom vribles 5. The norml distribution 5.3 The norml pproimtion to the binomil distribution 5.5 The uniorm distribution 5. Joint distribution discrete nd continuous
More informationChapter 5 : Continuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 216 Néhémy Lim Chpter 5 : Continuous Rndom Vribles Nottions. N {, 1, 2,...}, set of nturl numbers (i.e. ll nonnegtive integers); N {1, 2,...}, set of ll
More informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationA New Statistic Feature of the Short-Time Amplitude Spectrum Values for Human s Unvoiced Pronunciation
Xiodong Zhung A ew Sttistic Feture of the Short-Time Amplitude Spectrum Vlues for Humn s Unvoiced Pronuncition IAODOG ZHUAG 1 1. Qingdo University, Electronics & Informtion College, Qingdo, 6671 CHIA Abstrct:
More informationOnline Short Term Load Forecasting by Fuzzy ARTMAP Neural Network
Online Short Term Lod Forecsting by Fuzzy ARTMAP Neurl Network SHAHRAM JAVADI Electricl Engineering Deprtment AZAD University Tehrn Centrl Brnch Moshnir Power Electric Compny IRAN Abstrct: This pper presents
More informationAnalysis of the Effect of Time Delay on the Integrated GNSS/INS Navigation Systems
http://www.trnsnv.eu the Interntionl Journl on Mrine Nvigtion n Sfety of Se Trnsporttion Volume 7 Number June 13 DOI: 1.1716/11.7..6 Anlysis of the Effect of Time Dely on the Integrte / Nvigtion Systems
More informationTests for the Ratio of Two Poisson Rates
Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson
More informationDepartment of Physical Pharmacy and Pharmacokinetics Poznań University of Medical Sciences Pharmacokinetics laboratory
Deprtment of Physicl Phrmcy nd Phrmcoinetics Poznń University of Medicl Sciences Phrmcoinetics lbortory Experiment 1 Phrmcoinetics of ibuprofen s n exmple of the first-order inetics in n open one-comprtment
More informationData Assimilation. Alan O Neill Data Assimilation Research Centre University of Reading
Dt Assimiltion Aln O Neill Dt Assimiltion Reserch Centre University of Reding Contents Motivtion Univrite sclr dt ssimiltion Multivrite vector dt ssimiltion Optiml Interpoltion BLUE 3d-Vritionl Method
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationExpected Value of Function of Uncertain Variables
Journl of Uncertin Systems Vol.4, No.3, pp.8-86, 2 Online t: www.jus.org.uk Expected Vlue of Function of Uncertin Vribles Yuhn Liu, Minghu H College of Mthemtics nd Computer Sciences, Hebei University,
More informationLinear and Non-linear Feedback Control Strategies for a 4D Hyperchaotic System
Pure nd Applied Mthemtics Journl 017; 6(1): 5-13 http://www.sciencepublishinggroup.com/j/pmj doi: 10.11648/j.pmj.0170601.1 ISSN: 36-9790 (Print); ISSN: 36-981 (Online) Liner nd Non-liner Feedbck Control
More informationAcceptance Sampling by Attributes
Introduction Acceptnce Smpling by Attributes Acceptnce smpling is concerned with inspection nd decision mking regrding products. Three spects of smpling re importnt: o Involves rndom smpling of n entire
More informationMeasuring Electron Work Function in Metal
n experiment of the Electron topic Mesuring Electron Work Function in Metl Instructor: 梁生 Office: 7-318 Emil: shling@bjtu.edu.cn Purposes 1. To understnd the concept of electron work function in metl nd
More informationState space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies
Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response
More informationTravelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing
Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge
More informationThe steps of the hypothesis test
ttisticl Methods I (EXT 7005) Pge 78 Mosquito species Time of dy A B C Mid morning 0.0088 5.4900 5.5000 Mid Afternoon.3400 0.0300 0.8700 Dusk 0.600 5.400 3.000 The Chi squre test sttistic is the sum of
More informationEstimation of Binomial Distribution in the Light of Future Data
British Journl of Mthemtics & Computer Science 102: 1-7, 2015, Article no.bjmcs.19191 ISSN: 2231-0851 SCIENCEDOMAIN interntionl www.sciencedomin.org Estimtion of Binomil Distribution in the Light of Future
More information1B40 Practical Skills
B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need
More informationResearch Article Moment Inequalities and Complete Moment Convergence
Hindwi Publishing Corportion Journl of Inequlities nd Applictions Volume 2009, Article ID 271265, 14 pges doi:10.1155/2009/271265 Reserch Article Moment Inequlities nd Complete Moment Convergence Soo Hk
More informationOperations with Polynomials
38 Chpter P Prerequisites P.4 Opertions with Polynomils Wht you should lern: How to identify the leding coefficients nd degrees of polynomils How to dd nd subtrct polynomils How to multiply polynomils
More informationImproper Integrals. Type I Improper Integrals How do we evaluate an integral such as
Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph
More informationDefinition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim
Mth 9 Course Summry/Study Guide Fll, 2005 [1] Limits Definition of Limit: We sy tht L is the limit of f(x) s x pproches if f(x) gets closer nd closer to L s x gets closer nd closer to. We write lim f(x)
More informationIDENTIFYING ABBERANT SEGMENTS IN PERMANENT DOWNHOLE GAUGE DATA
IDENTIFYING ABBERANT SEGMENTS IN PERMANENT DOWNHOLE GAUGE DATA A REPORT SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR
More informationTemperature influence compensation in microbolometer detector for image quality enhancement
.26/qirt.26.68 Temperture influence compenstion in microolometer detector for imge qulity enhncement More info out this rticle: http://www.ndt.net/?id=2647 Astrct y M. Krupiński*, T. Sosnowski*, H. Mdur*
More information#6A&B Magnetic Field Mapping
#6A& Mgnetic Field Mpping Gol y performing this lb experiment, you will: 1. use mgnetic field mesurement technique bsed on Frdy s Lw (see the previous experiment),. study the mgnetic fields generted by
More informationf(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral
Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one
More informationConsequently, the temperature must be the same at each point in the cross section at x. Let:
HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
More informationHybrid Group Acceptance Sampling Plan Based on Size Biased Lomax Model
Mthemtics nd Sttistics 2(3): 137-141, 2014 DOI: 10.13189/ms.2014.020305 http://www.hrpub.org Hybrid Group Acceptnce Smpling Pln Bsed on Size Bised Lomx Model R. Subb Ro 1,*, A. Ng Durgmmb 2, R.R.L. Kntm
More informationQuantum Physics II (8.05) Fall 2013 Assignment 2
Quntum Physics II (8.05) Fll 2013 Assignment 2 Msschusetts Institute of Technology Physics Deprtment Due Fridy September 20, 2013 September 13, 2013 3:00 pm Suggested Reding Continued from lst week: 1.
More informationThe Atwood Machine OBJECTIVE INTRODUCTION APPARATUS THEORY
The Atwood Mchine OBJECTIVE To derive the ening of Newton's second lw of otion s it pplies to the Atwood chine. To explin how ss iblnce cn led to the ccelertion of the syste. To deterine the ccelertion
More informationQUADRATURE is an old-fashioned word that refers to
World Acdemy of Science Engineering nd Technology Interntionl Journl of Mthemticl nd Computtionl Sciences Vol:5 No:7 011 A New Qudrture Rule Derived from Spline Interpoltion with Error Anlysis Hdi Tghvfrd
More informationINVESTIGATION ON THE MODEL OF VORTEX-INDUCED
The Seventh Asi-Pcific Conference on Wind Engineering, November 8-1, 9, Tipei, Tiwn ABSTRACT INVESTIGATION ON THE MODEL OF VORTEX-INDUCED VIBRATIONS OF RECTANGULAR SUPER HIGH-RISE BUILDINGS Hi-Yng Wu 1
More informationMath 270A: Numerical Linear Algebra
Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner
More informationTHERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION
XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Deprtment 8.044 Sttisticl Physics I Spring Term 03 Problem : Doping Semiconductor Solutions to Problem Set # ) Mentlly integrte the function p(x) given in
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationApplication of Exp-Function Method to. a Huxley Equation with Variable Coefficient *
Interntionl Mthemticl Forum, 4, 9, no., 7-3 Appliction of Exp-Function Method to Huxley Eqution with Vrible Coefficient * Li Yo, Lin Wng nd Xin-Wei Zhou. Deprtment of Mthemtics, Kunming College Kunming,Yunnn,
More informationProbabilistic Stand Still Detection using Foot Mounted IMU
Probbilistic Stnd Still Detection using Foot Mounted IMU Jons Cllmer, Dvid Törnqvist nd Fredrik Gustfsson Division of Automtic Control Linköping University Linköping, Sweden {cllmer, tornqvist, fredrik}@isy.liu.se
More informationCredibility Hypothesis Testing of Fuzzy Triangular Distributions
666663 Journl of Uncertin Systems Vol.9, No., pp.6-74, 5 Online t: www.jus.org.uk Credibility Hypothesis Testing of Fuzzy Tringulr Distributions S. Smpth, B. Rmy Received April 3; Revised 4 April 4 Abstrct
More informationMinimum Energy State of Plasmas with an Internal Transport Barrier
Minimum Energy Stte of Plsms with n Internl Trnsport Brrier T. Tmno ), I. Ktnum ), Y. Skmoto ) ) Formerly, Plsm Reserch Center, University of Tsukub, Tsukub, Ibrki, Jpn ) Plsm Reserch Center, University
More information1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation
1 1.1. Liner Constnt Coefficient Equtions Section Objective(s): Overview of Differentil Equtions. Liner Differentil Equtions. Solving Liner Differentil Equtions. The Initil Vlue Problem. 1.1.1. Overview
More informationRecitation 3: More Applications of the Derivative
Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech
More informationReinforcement learning II
CS 1675 Introduction to Mchine Lerning Lecture 26 Reinforcement lerning II Milos Huskrecht milos@cs.pitt.edu 5329 Sennott Squre Reinforcement lerning Bsics: Input x Lerner Output Reinforcement r Critic
More informationSome remarks on GNSS integer ambiguity validation methods
Some remrks on GNSS integer mbiguity vlidtion methods T. Li nd J. Wng* Ambiguity resolution is n indispensble step in fst nd high precision Globl Nvigtion Stellite System (GNSS) bsed positioning. In generl,
More informationADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS
ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:
More informationLecture 23: Interpolatory Quadrature
Lecture 3: Interpoltory Qudrture. Qudrture. The computtion of continuous lest squres pproximtions to f C[, b] required evlutions of the inner product f, φ j = fxφ jx dx, where φ j is polynomil bsis function
More informationTime Truncated Two Stage Group Sampling Plan For Various Distributions
Time Truncted Two Stge Group Smpling Pln For Vrious Distributions Dr. A. R. Sudmni Rmswmy, S.Jysri Associte Professor, Deprtment of Mthemtics, Avinshilingm University, Coimbtore Assistnt professor, Deprtment
More informationTutorial 4. b a. h(f) = a b a ln 1. b a dx = ln(b a) nats = log(b a) bits. = ln λ + 1 nats. = log e λ bits. = ln 1 2 ln λ + 1. nats. = ln 2e. bits.
Tutoril 4 Exercises on Differentil Entropy. Evlute the differentil entropy h(x) f ln f for the following: () The uniform distribution, f(x) b. (b) The exponentil density, f(x) λe λx, x 0. (c) The Lplce
More informationChapter 0. What is the Lebesgue integral about?
Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous
More informationLECTURE NOTE #12 PROF. ALAN YUILLE
LECTURE NOTE #12 PROF. ALAN YUILLE 1. Clustering, K-mens, nd EM Tsk: set of unlbeled dt D = {x 1,..., x n } Decompose into clsses w 1,..., w M where M is unknown. Lern clss models p(x w)) Discovery of
More informationAn approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95
An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech
More informationINVESTIGATION OF MATHEMATICAL MODEL OF COMMUNICATION NETWORK WITH UNSTEADY FLOW OF REQUESTS
Trnsport nd Telecommuniction Vol No 4 9 Trnsport nd Telecommuniction 9 Volume No 4 8 34 Trnsport nd Telecommuniction Institute Lomonosov Rig LV-9 Ltvi INVESTIGATION OF MATHEMATICAL MODEL OF COMMUNICATION
More informationAn efficient fractional-pixel motion compensation based on Cubic convolution interpolation
Journl of Electricl nd Electronic Engineering 2014; 2(3): 47-54 Published online September 20, 2014 (http://www.sciencepublishinggroup.com/j/jeee) doi: 10.11648/j.jeee.20140203.11 ISSN: 2329-1613 (Print);
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationMath 360: A primitive integral and elementary functions
Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:
More informationMath 31S. Rumbos Fall Solutions to Assignment #16
Mth 31S. Rumbos Fll 2016 1 Solutions to Assignment #16 1. Logistic Growth 1. Suppose tht the growth of certin niml popultion is governed by the differentil eqution 1000 dn N dt = 100 N, (1) where N(t)
More informationBest Approximation in the 2-norm
Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion
More informationComposite Mendeleev s Quadratures for Solving a Linear Fredholm Integral Equation of The Second Kind
Globl Journl of Pure nd Applied Mthemtics. ISSN 0973-1768 Volume 12, Number (2016), pp. 393 398 Reserch Indi Publictions http://www.ripubliction.com/gjpm.htm Composite Mendeleev s Qudrtures for Solving
More informationStudent Activity 3: Single Factor ANOVA
MATH 40 Student Activity 3: Single Fctor ANOVA Some Bsic Concepts In designed experiment, two or more tretments, or combintions of tretments, is pplied to experimentl units The number of tretments, whether
More informationUNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction
Lesson : Logrithmic Functions s Inverses Prerequisite Skills This lesson requires the use of the following skills: determining the dependent nd independent vribles in n exponentil function bsed on dt from
More informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More informationThe EbayesThresh Package
The EbyesThresh Pckge Mrch 26, 2005 Title Empiricl Byes thresholding nd relted methods Version 1.3.0 Dte 2005-03-24 Author Bernrd W. Silvermn Mintiner Bernrd W. Silvermn
More informationNew data structures to reduce data size and search time
New dt structures to reduce dt size nd serch time Tsuneo Kuwbr Deprtment of Informtion Sciences, Fculty of Science, Kngw University, Hirtsuk-shi, Jpn FIT2018 1D-1, No2, pp1-4 Copyright (c)2018 by The Institute
More informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More informationSection 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40
Mth B Prof. Audrey Terrs HW # Solutions by Alex Eustis Due Tuesdy, Oct. 9 Section 5. #7,, 6,, 5; Section 5. #8, 9, 5,, 7, 3; Section 5.3 #4, 6, 9, 3, 6, 8, 3; Section 5.4 #7, 8,, 3, 5, 9, 4 5..7 Since
More information( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that
Arc Length of Curves in Three Dimensionl Spce If the vector function r(t) f(t) i + g(t) j + h(t) k trces out the curve C s t vries, we cn mesure distnces long C using formul nerly identicl to one tht we
More informationA SHORT NOTE ON THE MONOTONICITY OF THE ERLANG C FORMULA IN THE HALFIN-WHITT REGIME. Bernardo D Auria 1
Working Pper 11-42 (31) Sttistics nd Econometrics Series December, 2011 Deprtmento de Estdístic Universidd Crlos III de Mdrid Clle Mdrid, 126 28903 Getfe (Spin) Fx (34) 91 624-98-49 A SHORT NOTE ON THE
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More informationPartial Derivatives. Limits. For a single variable function f (x), the limit lim
Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles
More informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
More informationProbabilistic Investigation of Sensitivities of Advanced Test- Analysis Model Correlation Methods
Probbilistic Investigtion of Sensitivities of Advnced Test- Anlysis Model Correltion Methods Liz Bergmn, Mtthew S. Allen, nd Dniel C. Kmmer Dept. of Engineering Physics University of Wisconsin-Mdison Rndll
More informationu( t) + K 2 ( ) = 1 t > 0 Analyzing Damped Oscillations Problem (Meador, example 2-18, pp 44-48): Determine the equation of the following graph.
nlyzing Dmped Oscilltions Prolem (Medor, exmple 2-18, pp 44-48): Determine the eqution of the following grph. The eqution is ssumed to e of the following form f ( t) = K 1 u( t) + K 2 e!"t sin (#t + $
More informationResearch on the Quality Competence in Manufacturing Industry
Reserch on the Qulity Competence in Mnufcturing Industry Xioping M, Zhijun Hn Economics nd Mngement School Nnjing University of Science nd Technology Nnjing 210094, Chin Tel: 86-25-8431-5400 E-mil: hnzhij4531@sin.com
More informationReview of basic calculus
Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More information1. a) Describe the principle characteristics and uses of the following types of PV cell: Single Crystal Silicon Poly Crystal Silicon
2001 1. ) Describe the principle chrcteristics nd uses of the following types of PV cell: Single Crystl Silicon Poly Crystl Silicon Amorphous Silicon CIS/CIGS Gllium Arsenide Multijunction (12 mrks) b)
More informationAbstract inner product spaces
WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the
More informationResearch Article The Prediction for Shanghai Business Climate Index by Grey Model
Reserch Journl of Applied Sciences, Engineering nd echnology 7(4): 976-980, 04 DOI:0.906/rjset.7.69 ISSN: 040-7459; e-issn: 040-7467 04 Mwell Scientific Publiction Corp. Submitted: September 4, 03 Accepted:
More informationDeteriorating Inventory Model for Waiting. Time Partial Backlogging
Applied Mthemticl Sciences, Vol. 3, 2009, no. 9, 42-428 Deteriorting Inventory Model for Witing Time Prtil Bcklogging Nit H. Shh nd 2 Kunl T. Shukl Deprtment of Mthemtics, Gujrt university, Ahmedbd. 2
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationSection 14.3 Arc Length and Curvature
Section 4.3 Arc Length nd Curvture Clculus on Curves in Spce In this section, we ly the foundtions for describing the movement of n object in spce.. Vector Function Bsics In Clc, formul for rc length in
More information