Improving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization
|
|
- Earl Leslie Nichols
- 6 years ago
- Views:
Transcription
1 1976 MONTHLY WEATHER REVIEW VOLUME 15 Improving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization PETER LYNCH Met Éireann, Dublin, Ireland DOMINIQUE GIARD CNRM/GMAP, Météo-France, Touloue, France VLADIMIR IVANOVICI Romanian Intitute of Meteorology and Hydrology, Bucharet, Romania 4 April 1996 and December 1996 ABSTRACT A method of diabatic initialization i decribed in which a nonrecurive digital filter i applied to both the backward (adiabatic) and forward (diabatic) tep. It ha clear advantage over previouly propoed cheme: the initialization require ignificantly le computation time and the reulting change in the analyzed field are conitently maller. Depite thee maller increment, the uppreion of puriou high-frequency noie i at leat a effective a for the earlier cheme. 1. Introduction Digital filtering ha been hown to be an effective mean of initializing data for numerical weather prediction. In Lynch and Huang (199, hereafter LH9) a imple filter wa applied to a equence of value centered on the initial time t 0; thee value were generated by two hort adiabatic model integration, one forward and one backward from t 0. It wa found that a total pan of 6 h wa required for effective elimination of high-frequency noie. Huang and Lynch (1993, hereafter HL93) howed, by mean of an optimal filter, that equally effective noie control could be achieved with a pan of only 3 h. They alo decribed a mean of incorporating diabatic effect: the backward adiabatic integration (of length 1.5 h) i followed by a diabatic forecat of twice the length, and the value generated by the diabatic integration are proceed by the filter. In Lynch and Huang (1994, hereafter LH94), recurive filter were applied to the initialization problem. Several cheme were invetigated; in one, deignated RAD recurive adiabatic-diabatic in LH94, the backward adiabatic integration wa proceed with a recurive filter, and the terminal output value wa ued to initiate a forward diabatic integration, to which re- Correponding author addre: Dr. Peter Lynch, Met Éireann, Dublin 9, Ireland. plynch@irmet.ie curive filtering wa again applied. The application of filtering to both the revere and forward integration i alo poible uing a nonrecurive filter. In thi note we examine the advantage of thi idea.. A imple conceptual model Let u conider a function of time x(t) governed by the ocillation equation dx i 0 x 0 (1) dt with the initial value x(0). The olution i x(0) exp(i 0 t). Suppoe thi olution i filtered by convolution with an impule repone function h(t): x (t) h x(t) h()x(t ) d. Subtituting the olution x(0) exp(i 0 t) into the integral, it follow immediately that x (t) H( 0 )x(t), where the frequency repone function H() i the Fourier tranform of h(t). Now let u ue a value of the filtered olution y(0) x (0) a an initial value for a olution of the inhomogeneou equation dy i y f(t). () 0 dt 1997 American Meteorological Society
2 AUGUST 1997 NOTES AND CORRESPONDENCE 1977 FIG. 1. Schematic illutration of original and modified filtering procedure. Thin line indicate backward adiabatic integration; thick line indicate forward diabatic integration. (a) Scheme of Huang and Lynch (1993): adiabatic tep i not filtered. (b) Modified cheme: both tep are filtered. The olution i eaily obtained by Laplace tranformation or otherwie. For inuoidal forcing, f(t) f 0 exp(i 1 t)( 0 1 ), it may be written ] 0[ exp(i0t) exp(i1t) y(t) y(0) exp(i0t) f. (3) i( ) If thi i filtered, the reult i ] 0[ i0t i1t H( 0)e H( 1)e ȳ(t) [H( 0)] x(t) f, (4) i( ) wherea tarting from the unfiltered value x(0), a in HL93, it would be ] 0[ i0t i1t H( 0)e H( 1)e ȳ(t) [H( 0)]x(t) f. (5) i( ) Let u conider that (1) correpond to the adiabatic model and () to the diabatic model, with phyical forcing repreented by f(t). Suppoe now that the adiabatic model i integrated and the olution x(t) i filtered to give x (t). Thi i ued to define initial data y(0) for an integration of the diabatic model (). The reulting olution y(t) i again filtered, producing ȳ(t). Then (4) how that the adiabatic part of the olution i effectively filtered twice, while the diabatic part i filtered once. If the amplitude f 0 of the forcing i ufficiently mall, a ingle filtering of the diabatic component hould be adequate to reduce it high frequency component to an acceptable level. The initialization cheme decribed in HL93 i illutrated chematically in Fig. 1a. The thin line repreent the revere adiabatic integration of duration T. Thi i not filtered; it terminal value, depicted by an open circle, i ued to initiate a forward diabatic integration (the thick line) that i then ubjected to filtering. Thi integration mut be of duration T to enure that the output of the ymmetric filter i valid at t 0. The modified cheme, to be analyzed here, i illutrated in Fig. 1b. The revere adiabatic integration (thin curve) i filtered, yielding an output valid at t T/. Thi value, depicted by the black pot, i ued to initiate a forward diabatic integration (thick curve) of duration T, centered on the initial time t 0. Conidering the adiabatic component of the olution, we ee that the HL93 cheme involve a ingle application of a filter of pan T; the modified cheme involve two application of a filter of pan T. Comparion between the cheme can be made by conidering the relationhip between the magnitude of the repone function, H T (), of the filter with the longer time pan and the quared repone function H T () of the filter with the horter time pan. The diabatic initialization cheme apply filtering to trajectorie that do not correpond, at time t 0, to the original analyi. Let u uppoe that both 0 and 1 are low frequencie, o that we may aume H( 0 ) H( 1 ) 1. Thu, the filtering doe not ubtantially affect thee component. The reult of a backward adiabatic integration of length T followed by a forward diabatic integration of the ame length may eaily be deduced: integrate (1) backward from x(0) to time t T to get x(t) x(0) exp(i 0 T); ue thi value to initiate a forward integration of () from t Tto t 0 to obtain 0[ ] exp(i0t) exp(i1t) y(0) x(0) f. (6) i( ) For hort time pan ( 0 T and 1 T mall), auming the filter repone i unity, we have ȳ(0) x(0) f 0 T. (7) Thi differ from the original value x(0) by the diabatic dicrepancy f 0 T, which i proportional to the amplitude f 0 of the diabatic forcing and alo to the time pan T. Since the pan of the new cheme i typically half that of the old cheme, the effect of thi dicrepancy i reduced. Of coure, the problem of diabatic dicrepancy diappear when uing a finalization or launch technique a decribed in Fillion et al. (1995) and LH94, repectively; but then the filtered field are not applicable at the initial time. 3. Some example of filter The filter ued in HL93 wa an optimal filter having pa-band and top-band edge with period p 15 h and 3 h (recall that the period i given by / t/). The total pan wa T 3 h and the time tep wa t 360. A hown in Lynch (1997), the Dolph Chebyhev filter with a top-band edge 3 h and the ame time pan (3 h) i a imple optimal filter with virtually the ame frequency repone. We will
3 1978 MONTHLY WEATHER REVIEW VOLUME 15 FIG.. Frequency repone (db) for Dolph filter with ripple ratio r 0.1. Solid line: filter order N M 1 37; 0 log H() plotted. Dahed line: filter order N M 1 19; 0 log H() plotted. FIG. 3. Amplitude repone for Dolph filter a in Fig. above. Solid line: filter order N M 1 37; H() plotted. Dahed line: filter order N M 1 19; H() plotted. therefore ue imple Dolph Chebyhev filter for the comparion below. We chooe a time tep t 300 o that M T/t 18 i an even number. Let u recall that the frequency repone of a Dolph Chebyhev filter with a top-band edge frequency i co H() rt M, (8) co where T M (x) i a Chebyhev polynomial of order M and the total time pan i T Mt. The ripple ratio r meaure the maximum repone in the top band ( ) and i determined by requiring H(0) 1: TM coh M coh (9) r co co (ee Lynch 1997). In Fig., two repone function are hown. The ordinate i the filter attenuation, 0 logh() (db). The olid curve i the repone of a Dolph Chebyhev filter (denote it H T ) of total pan T 3 h; the topband edge i 3 h and M 18. The dahed curve i the quared repone of a Dolph Chebyhev filter ( H T ) of total pan T 1.5 h; thu, M 9. The topband edge.5 h ha been choen by experiment o that the total attenuation in the top band wa about the ame in both cae: 1.3 db for the filter with the longer pan and 1. db for the quared repone of the filter with the horter time pan. A more ytematic mean of chooing will be decribed below. The amplitude of the frequency repone of the two filter are hown in Fig. 3 with the axi expanded for clarity: the repone are very imilar, but the pa band of the quared hort filter i omewhat broader. If a emi-lagrangian advection cheme i ued, relatively long time tep are poible. We conider a comparion between filter H T and H T having the ame top band and total pan a above, but with a time tep t 900. Thu, M 6 for H T and M 3 for H T. The frequency repone are hown in Fig. 4. Once again, FIG. 4. Frequency repone (db) for Dolph filter with ripple ratio r 0.1. Solid line: filter order N M 1 13; 0 log H() plotted. Dahed line: filter order N M 1 7; 0 log H() plotted.
4 AUGUST 1997 NOTES AND CORRESPONDENCE 1979 the attenuation in the top band i imilar for both filter (1.6 db in each cae) and the pa band of the quared hort filter i lightly broader. The repone in the pa band for filter with the longer and horter time tep were compared (by recaling the axi) and were found to be very imilar: thu, effective frequency dicrimination i achievable with the longer time tep. 4. Chooing the filter parameter Method of enuring that ingly and doubly applied filter have the ame level of damping for high frequencie will now be decribed. Let unprimed and primed quantitie refer to the filter with the longer and horter time pan, repectively. Firt, let u aume that the horter pan i exactly one-half of the longer one. We now determine the parameter for a filter with one-half of the pan, having a ripple ratio r r, o that it quared repone will have the ame high frequency attenuation a the longer filter. It top-band edge i determined by inverion of (9) with appropriate parameter value: [ ] co coh coh. (10) M r It wa found in the above example that wa omewhat larger than. It mut be acertained in each cae whether i mall enough to provide the required frequency electivity. Numerical tet uing (10) revealed that, in ome cae, the widening of the top band for the horter filter wa unacceptable. An alternative approach, which enure that both filter have the ame top-band edge, will now be decribed. Firt we note that, if and r are pecified, the required filter pan can be deduced from (9): we have T Mt with 1 coh (1/r) M. 1 coh 1 co Noting that, for mall, coh 1 (1/co ), amay be eaily verified by tandard Taylor expanion, and recalling that t/, we find that T coh, (11) r provided K or t K, which i alway true in practice. Now aume that the longer and horter pan are T Mt and T Mt, repectively, and that both filter have the ame top-band edge. The repone function are co H() rt M, co where r 1/T M (1/co /) for the longer filter, and co H() rt M, co where r 1/T M (1/co /) for the horter one. Since the filter with the horter pan i to be applied twice, we require r rto enure that the top-band attenuation i equivalent in both cae. Thu, from (11) we have T coh. (1) r Taking typical value 3 h and r 0.05 we find from (11) that the filter applied only once mut have a pan at leat equal to 3.5 h, wherea (1) implie that the filter applied twice achieve the ame high-frequency damping with a pan of only.08 h. Thi lead to a ubtantial reduction in computation time. 5. Application to a forecat model The new initialization cheme ha been evaluated by application to the limited-area pectral model ALADIN (partially decribed in Bubnová et al. 1995). Thi model i run quai-operationally at Météo-France and coupled with the global pectral NWP model ARPEGE. Horizontal repreentation of the variable i achieved by double Fourier erie. A nonhydrotatic verion of the model ha been developed, but the reult preented below are for the hydrotatic verion (application to the nonhydrotatic verion would imply involve filtering of the two additional prognotic variable). The gridize here i approximately 18.3 km and there are 7 vertical level. The time tep i t 450. It wa found to be neceary to ue a maller time tep, t 5, for the backward adiabatic integration to enure tability. The geographical area covered may be een later in Fig. 7. ALADIN doe not have it own analyi cheme. The initialized analyi produced by the global pectral model ARPEGE i tranformed to the reolution required for ALADIN. Thi tranformation introduce puriou noie that mut be removed by initialization on the limited domain. Four forecat were performed tarting from analye valid at 0000 UTC 1 January One wa from the tranformed analyi without any ubequent initialization. The other three forecat followed digital filtering initialization in three verion. 1) A Lanczo filter with pan 6 h and a cutoff c 6h.
5 1980 MONTHLY WEATHER REVIEW VOLUME 15 FIG. 5. The frequency repone of the filter ued in the three cheme decribed in the text. Solid line cheme 1, Lanczo filter applied once; dahed line cheme, Dolph Chebyhev filter applied once; dotted line cheme 3, Dolph Chebyhev filter applied twice. Note that for cheme 3 the quare of the repone function i plotted. FIG. 6. Evolution of the mean abolute urface preure tendency for the firt ix forecat hour, tarting from uninitialized data (olid) and from data filtered with cheme 1 (dotted), (dahed), and 3 (dot dahed). A backward integration of 3 h wa followed by a 6-h diabatic forecat. Filtering wa applied once only, to the forward integration. ) A Dolph Chebyhev filter with pan 4.5 h and topband edge 3 h. Filtering wa again applied once only, to the forward integration. 3) A Dolph Chebyhev filter with pan.5 h and topband edge 3 h. Filtering wa applied twice, to both the backward and forward integration. Thi i the new cheme. The frequency repone of the filter ued in the three cheme are hown in Fig. 5. Note that for the new cheme it i the quare of the repone function that i plotted ince the filter i applied twice. The Lanczo filter i preciely that ued in LH9 and i known to be effective. Note that the abcia i the period. Iti clear that the repone of the three filtering cheme are very imilar. The time pan ued here are longer than thoe conidered in ection 3 and allow a better damping of high frequencie, with a ripple ratio lower than 0.05 intead of 0.10 for Dolph Chebyhev filter. At the ame time, the attenuation of period longer than 1 h remain mall. Note that the contraint on the filter are not a trong in thi cae a might be required for a global model (a dicued in Fillion et al. 1995) or a limited-area model with it own aimilation cycle. The evolution of the abolute tendency of urface preure averaged over the forecat domain i hown in Fig. 6. The exceive noie preent in the uninitialized forecat (olid curve) i uccefully removed by all the filtering cheme, and they appear to be equally effective in thi regard. There i a mall amount of reidual noie but it i not of ufficient amplitude to caue concern. The average abolute urface preure tendency i an indicator of the level of noie in the external gravity wave component. It i alo neceary to tudy the impact of filtering on puriou internal gravity wave energy. For thi purpoe, the 500-hPa vertical velocity at the initial time wa examined for the four forecat (reult not hown). It wa found that large amplitude, mall-cale noie in the uninitialized run, particularly in the region of the Alp, wa effectively removed by all the filtering cheme. The change induced by the three filtering cheme were examined and compared. It wa found that the new cheme conitently led to maller change than either of the cheme involving ingle filtering. The impact on the 500-hPa temperature analyi i hown in Fig. 7. Figure 7a,b,c how, repectively, the increment due to filtering cheme 1,, and 3. It i clear that the change brought about by the new cheme are ignificantly maller than thoe of the alternative. Notwithtanding thi, the noie reduction i equally effective. It i reaonable to acribe the reduction in initialization increment for the new cheme to the reduced diabatic dicrepancy aociated with the horter pan ince, a we have een, the filter frequency repone i virtually identical for the three cheme. The aving in computation time for the new cheme may eaily be etimated. Let t B and c B be the time tep and cot per time tep for the backward adiabatic run, t and c the correponding value for the forward diabatic run, and T the total time pan. We may uppoe t B t/ and c B c/. Then the cot of the old cheme (ingle filtering) i
6 AUGUST 1997 NOTES AND CORRESPONDENCE 1981 FIG. 7. Impact on the 500-hPa temperature analyi: (a), (b), and (c) how the increment due to filtering cheme 1,, and 3. The contour interval i 0. K. T T 3Tc C cb c. tb t t If the filter pan of the new cheme i half a long, the cot i T T Tc C cb c, tb t t which reult in a aving of 33%. Thi wa the cae in the above example: the computation time (econd of
7 198 MONTHLY WEATHER REVIEW VOLUME 15 CPU uage) for cheme 1 wa 19 ; for cheme it wa 955 ; and for cheme 3 it wa 66. The reduction i of practical benefit in an operational context. The above reult demontrate that a ignificant gain in efficiency can be made by mean of the new cheme. Thi cheme reult in maller increment to the initial field but i equally effective in reducing high frequency noie in the forecat. The new cheme with Dolph Chebyhev filter ha been implemented in the operational ALADIN-LACE uite (focuing on central and eatern Europe) and in the preoperational ALADIN- FRANCE uite at Météo-France. It ha alo been implemented in the global ytem ARPEGE, in the framework of incremental initialization, o a to further reduce the diabatic dicrepancy and preerve tidal mode (a uggeted in LH94). In order to reaonably damp the large-cale gravity component introduced by analyi, a more elective Dolph Chebyhev filter ha been ued (with 5 h a the top-band edge). REFERENCES Bubnová, R., G. Hello, P. Bénard, and J.-F. Geleyn, 1995: Integration of the fully elatic equation cat in the hydrotatic preure terrain-following coordinate in the framework of the ARPE- GE/Aladin NWP ytem. Mon. Wea. Rev., 13, Fillion, L., H. L. Mitchell, H. Richie, and A. Staniforth, 1995: The impact of a digital filter finalization technique in a global data aimilation ytem. Tellu, 47A, Huang, X.-Y., and P. Lynch, 1993: Diabatic digital-filtering initialization: Application to the HIRLAM model. Mon. Wea. Rev., 11, Lynch, P., 1997: The Dolph Chebyhev window: A imple optimal filter. Mon. Wea. Rev., 15, , and X.-Y. Huang, 199: Initialization of the HIRLAM model uing a digital filter. Mon. Wea. Rev., 10, , and, 1994: Diabatic initialization uing recurive filter Tellu, 46A,
Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationLecture 10 Filtering: Applied Concepts
Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering
More informationDesign of Digital Filters
Deign of Digital Filter Paley-Wiener Theorem [ ] ( ) If h n i a caual energy ignal, then ln H e dω< B where B i a finite upper bound. One implication of the Paley-Wiener theorem i that a tranfer function
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationSampling and the Discrete Fourier Transform
Sampling and the Dicrete Fourier Tranform Sampling Method Sampling i mot commonly done with two device, the ample-and-hold (S/H) and the analog-to-digital-converter (ADC) The S/H acquire a CT ignal at
More informationLecture #9 Continuous time filter
Lecture #9 Continuou time filter Oliver Faut December 5, 2006 Content Review. Motivation......................................... 2 2 Filter pecification 2 2. Low pa..........................................
More informationA Simple Approach to Synthesizing Naïve Quantized Control for Reference Tracking
A Simple Approach to Syntheizing Naïve Quantized Control for Reference Tracking SHIANG-HUA YU Department of Electrical Engineering National Sun Yat-Sen Univerity 70 Lien-Hai Road, Kaohiung 804 TAIAN Abtract:
More information5.5 Application of Frequency Response: Signal Filters
44 Dynamic Sytem Second order lowpa filter having tranfer function H()=H ()H () u H () H () y Firt order lowpa filter Figure 5.5: Contruction of a econd order low-pa filter by combining two firt order
More informationDesign By Emulation (Indirect Method)
Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal
More informationMolecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions
Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,
More informationAnnex-A: RTTOV9 Cloud validation
RTTOV-91 Science and Validation Plan Annex-A: RTTOV9 Cloud validation Author O Embury C J Merchant The Univerity of Edinburgh Intitute for Atmo. & Environ. Science Crew Building King Building Edinburgh
More informationGain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays
Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,
More informationDigital Control System
Digital Control Sytem - A D D A Micro ADC DAC Proceor Correction Element Proce Clock Meaurement A: Analog D: Digital Continuou Controller and Digital Control Rt - c Plant yt Continuou Controller Digital
More informationChapter 4: Applications of Fourier Representations. Chih-Wei Liu
Chapter 4: Application of Fourier Repreentation Chih-Wei Liu Outline Introduction Fourier ranform of Periodic Signal Convolution/Multiplication with Non-Periodic Signal Fourier ranform of Dicrete-ime Signal
More informationComparing Means: t-tests for Two Independent Samples
Comparing ean: t-tet for Two Independent Sample Independent-eaure Deign t-tet for Two Independent Sample Allow reearcher to evaluate the mean difference between two population uing data from two eparate
More informationBogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationGiven the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is
EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -
More informationCodes Correcting Two Deletions
1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of
More informationLOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM. P.Dickinson, A.T.Shenton
LOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM P.Dickinon, A.T.Shenton Department of Engineering, The Univerity of Liverpool, Liverpool L69 3GH, UK Abtract: Thi paper compare
More informationRecent progress in fire-structure analysis
EJSE Special Iue: Selected Key Note paper from MDCMS 1 1t International Conference on Modern Deign, Contruction and Maintenance of Structure - Hanoi, Vietnam, December 2007 Recent progre in fire-tructure
More informationWhite Rose Research Online URL for this paper: Version: Accepted Version
Thi i a repoitory copy of Identification of nonlinear ytem with non-peritent excitation uing an iterative forward orthogonal leat quare regreion algorithm. White Roe Reearch Online URL for thi paper: http://eprint.whiteroe.ac.uk/107314/
More informationOnline supplementary information
Electronic Supplementary Material (ESI) for Soft Matter. Thi journal i The Royal Society of Chemitry 15 Online upplementary information Governing Equation For the vicou flow, we aume that the liquid thickne
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
More informationTHE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER
Proceeding of IMAC XXXI Conference & Expoition on Structural Dynamic February -4 Garden Grove CA USA THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER Yung-Sheng Hu Neil S Ferguon
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationChapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationControl Systems Analysis and Design by the Root-Locus Method
6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If
More informationChapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem
Chapter 5 Conitency, Zero Stability, and the Dahlquit Equivalence Theorem In Chapter 2 we dicued convergence of numerical method and gave an experimental method for finding the rate of convergence (aka,
More informationControl of Delayed Integrating Processes Using Two Feedback Controllers R MS Approach
Proceeding of the 7th WSEAS International Conference on SYSTEM SCIENCE and SIMULATION in ENGINEERING (ICOSSSE '8) Control of Delayed Integrating Procee Uing Two Feedback Controller R MS Approach LIBOR
More informationCompact finite-difference approximations for anisotropic image smoothing and painting
CWP-593 Compact finite-difference approximation for aniotropic image moothing and painting Dave Hale Center for Wave Phenomena, Colorado School of Mine, Golden CO 80401, USA ABSTRACT Finite-difference
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationSMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD
SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem
More information5. Fuzzy Optimization
5. Fuzzy Optimization 1. Fuzzine: An Introduction 135 1.1. Fuzzy Memberhip Function 135 1.2. Memberhip Function Operation 136 2. Optimization in Fuzzy Environment 136 3. Fuzzy Set for Water Allocation
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationDigital Control System
Digital Control Sytem Summary # he -tranform play an important role in digital control and dicrete ignal proceing. he -tranform i defined a F () f(k) k () A. Example Conider the following equence: f(k)
More informationAlternate Dispersion Measures in Replicated Factorial Experiments
Alternate Diperion Meaure in Replicated Factorial Experiment Neal A. Mackertich The Raytheon Company, Sudbury MA 02421 Jame C. Benneyan Northeatern Univerity, Boton MA 02115 Peter D. Krau The Raytheon
More informationDYNAMIC MODELS FOR CONTROLLER DESIGN
DYNAMIC MODELS FOR CONTROLLER DESIGN M.T. Tham (996,999) Dept. of Chemical and Proce Engineering Newcatle upon Tyne, NE 7RU, UK.. INTRODUCTION The problem of deigning a good control ytem i baically that
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationReal Sources (Secondary Sources) Phantom Source (Primary source) LS P. h rl. h rr. h ll. h lr. h pl. h pr
Ecient frequency domain ltered-x realization of phantom ource iet C.W. ommen, Ronald M. Aart, Alexander W.M. Mathijen, John Gara, Haiyan He Abtract A phantom ound ource i a virtual ound image which can
More informationSource slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis
Source lideplayer.com/fundamental of Analytical Chemitry, F.J. Holler, S.R.Crouch Chapter 6: Random Error in Chemical Analyi Random error are preent in every meaurement no matter how careful the experimenter.
More informationCorrection for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002
Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in
More informationSAMPLING. Sampling is the acquisition of a continuous signal at discrete time intervals and is a fundamental concept in real-time signal processing.
SAMPLING Sampling i the acquiition of a continuou ignal at dicrete time interval and i a fundamental concept in real-time ignal proceing. he actual ampling operation can alo be defined by the figure belo
More informationEE/ME/AE324: Dynamical Systems. Chapter 8: Transfer Function Analysis
EE/ME/AE34: Dynamical Sytem Chapter 8: Tranfer Function Analyi The Sytem Tranfer Function Conider the ytem decribed by the nth-order I/O eqn.: ( n) ( n 1) ( m) y + a y + + a y = b u + + bu n 1 0 m 0 Taking
More informationSIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. Solutions to Assignment 3 February 2005.
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuit II Solution to Aignment 3 February 2005. Initial Condition Source 0 V battery witch flip at t 0 find i 3 (t) Component value:
More informationRaneNote BESSEL FILTER CROSSOVER
RaneNote BESSEL FILTER CROSSOVER A Beel Filter Croover, and It Relation to Other Croover Beel Function Phae Shift Group Delay Beel, 3dB Down Introduction One of the way that a croover may be contructed
More informationSingular perturbation theory
Singular perturbation theory Marc R. Rouel June 21, 2004 1 Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly
More informationBy Xiaoquan Wen and Matthew Stephens University of Michigan and University of Chicago
Submitted to the Annal of Applied Statitic SUPPLEMENTARY APPENDIX TO BAYESIAN METHODS FOR GENETIC ASSOCIATION ANALYSIS WITH HETEROGENEOUS SUBGROUPS: FROM META-ANALYSES TO GENE-ENVIRONMENT INTERACTIONS
More informationStandard Guide for Conducting Ruggedness Tests 1
Deignation: E 69 89 (Reapproved 996) Standard Guide for Conducting Ruggedne Tet AMERICA SOCIETY FOR TESTIG AD MATERIALS 00 Barr Harbor Dr., Wet Conhohocken, PA 948 Reprinted from the Annual Book of ASTM
More informationCHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationClustering Methods without Given Number of Clusters
Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,
More informationProblem Set 8 Solutions
Deign and Analyi of Algorithm April 29, 2015 Maachuett Intitute of Technology 6.046J/18.410J Prof. Erik Demaine, Srini Devada, and Nancy Lynch Problem Set 8 Solution Problem Set 8 Solution Thi problem
More informationFinding the location of switched capacitor banks in distribution systems based on wavelet transform
UPEC00 3t Aug - 3rd Sept 00 Finding the location of witched capacitor bank in ditribution ytem baed on wavelet tranform Bahram nohad Shahid Chamran Univerity in Ahvaz bahramnohad@yahoo.com Mehrdad keramatzadeh
More informationPreemptive scheduling on a small number of hierarchical machines
Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,
More informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
Stochatic Optimization with Inequality Contraint Uing Simultaneou Perturbation and Penalty Function I-Jeng Wang* and Jame C. Spall** The John Hopkin Univerity Applied Phyic Laboratory 11100 John Hopkin
More informationProperties of Z-transform Transform 1 Linearity a
Midterm 3 (Fall 6 of EEG:. Thi midterm conit of eight ingle-ided page. The firt three page contain variou table followed by FOUR eam quetion and one etra workheet. You can tear out any page but make ure
More informationGNSS Solutions: What is the carrier phase measurement? How is it generated in GNSS receivers? Simply put, the carrier phase
GNSS Solution: Carrier phae and it meaurement for GNSS GNSS Solution i a regular column featuring quetion and anwer about technical apect of GNSS. Reader are invited to end their quetion to the columnit,
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationConvex Hulls of Curves Sam Burton
Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall
More informationBUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL
J. Autral. Math. Soc. Ser. B 4(999), 332 349 BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL J. LEE and J.-M. VANDEN-BROECK 2 (Received 22 April 995; revied 23 April 996)
More informationUSPAS Course on Recirculated and Energy Recovered Linear Accelerators
USPAS Coure on Recirculated and Energy Recovered Linear Accelerator G. A. Krafft and L. Merminga Jefferon Lab I. Bazarov Cornell Lecture 6 7 March 005 Lecture Outline. Invariant Ellipe Generated by a Unimodular
More informationSemiflexible chains under tension
Semiflexible chain under tenion B.-Y. Ha and D. Thirumalai Intitute for Phyical Science and Technology, Univerity of Maryland, College Park, Maryland 74 Received 16 September 1996; accepted 5 December
More informationS_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS
S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS by Michelle Gretzinger, Daniel Zyngier and Thoma Marlin INTRODUCTION One of the challenge to the engineer learning proce control i relating theoretical
More informationCHAPTER 6. Estimation
CHAPTER 6 Etimation Definition. Statitical inference i the procedure by which we reach a concluion about a population on the bai of information contained in a ample drawn from that population. Definition.
More informationSolutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam
BSc - Sample Examination Digital Control Sytem (5-588-) Prof. L. Guzzella Solution Exam Duration: Number of Quetion: Rating: Permitted aid: minute examination time + 5 minute reading time at the beginning
More informationLecture 7: Testing Distributions
CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting
More informationSpring 2014 EE 445S Real-Time Digital Signal Processing Laboratory. Homework #0 Solutions on Review of Signals and Systems Material
Spring 4 EE 445S Real-Time Digital Signal Proceing Laboratory Prof. Evan Homework # Solution on Review of Signal and Sytem Material Problem.. Continuou-Time Sinuoidal Generation. In practice, we cannot
More informationFOURIER-BASED METHODS FOR THE SPECTRAL ANALYSIS OF MUSICAL SOUNDS. Sylvain Marchand
FOUIE-BAS METHOS FO THE SPECTAL ANALYSIS OF MUSICAL SOUNS Sylvain Marchand Univerity of Bret, Lab-STICC CNS UM 628, 292 Bret, Brittany, France ABSTACT When dealing with muical ound, the hort-time Fourier
More informationThe continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.
1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker
More information[Saxena, 2(9): September, 2013] ISSN: Impact Factor: INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY
[Saena, (9): September, 0] ISSN: 77-9655 Impact Factor:.85 IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY Contant Stre Accelerated Life Teting Uing Rayleigh Geometric Proce
More informationAn Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem
An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...4
More informationChapter 13. Root Locus Introduction
Chapter 13 Root Locu 13.1 Introduction In the previou chapter we had a glimpe of controller deign iue through ome imple example. Obviouly when we have higher order ytem, uch imple deign technique will
More informationRoot Locus Diagram. Root loci: The portion of root locus when k assume positive values: that is 0
Objective Root Locu Diagram Upon completion of thi chapter you will be able to: Plot the Root Locu for a given Tranfer Function by varying gain of the ytem, Analye the tability of the ytem from the root
More informationMAHALAKSHMI ENGINEERING COLLEGE-TRICHY
DIGITAL SIGNAL PROCESSING DEPT./SEM.: CSE /VII DIGITAL FILTER DESIGN-IIR & FIR FILTER DESIGN PART-A. Lit the different type of tructure for realiation of IIR ytem? AUC APR 09 The different type of tructure
More informationModeling Local and Advective Diffusion of Fuel Vapors to Understand Aqueous Foams in Fighting Fires. Final Report
Modeling Local and Advective Diffuion of Fuel Vapor to Undertand Aqueou Foam in Fighting Fire Final Report Author: Andrew Brandon abrando@math.umd.edu Advior: Dr. Ramagopal Ananth Naval Reearch Laboratory,
More informationCalculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary
More informationPractice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions
For Quetion -6, rewrite the piecewie function uing tep function, ketch their graph, and find F () = Lf(t). 0 0 < t < 2. f(t) = (t 2 4) 2 < t In tep-function form, f(t) = u 2 (t 2 4) The graph i the olid
More informationμ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL =
Our online Tutor are available 4*7 to provide Help with Proce control ytem Homework/Aignment or a long term Graduate/Undergraduate Proce control ytem Project. Our Tutor being experienced and proficient
More informationAdvanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment
Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,
More informationonline learning Unit Workbook 4 RLC Transients
online learning Pearon BTC Higher National in lectrical and lectronic ngineering (QCF) Unit 5: lectrical & lectronic Principle Unit Workbook 4 in a erie of 4 for thi unit Learning Outcome: RLC Tranient
More informationIII.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES
III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SBSTANCES. Work purpoe The analyi of the behaviour of a ferroelectric ubtance placed in an eternal electric field; the dependence of the electrical polariation
More informationConvergence criteria and optimization techniques for beam moments
Pure Appl. Opt. 7 (1998) 1221 1230. Printed in the UK PII: S0963-9659(98)90684-5 Convergence criteria and optimization technique for beam moment G Gbur and P S Carney Department of Phyic and Atronomy and
More informationSuggestions - Problem Set (a) Show the discriminant condition (1) takes the form. ln ln, # # R R
Suggetion - Problem Set 3 4.2 (a) Show the dicriminant condition (1) take the form x D Ð.. Ñ. D.. D. ln ln, a deired. We then replace the quantitie. 3ß D3 by their etimate to get the proper form for thi
More informationA Bluffer s Guide to... Sphericity
A Bluffer Guide to Sphericity Andy Field Univerity of Suex The ue of repeated meaure, where the ame ubject are teted under a number of condition, ha numerou practical and tatitical benefit. For one thing
More informationChapter 2: Problem Solutions
Chapter 2: Solution Dicrete Time Proceing of Continuou Time Signal Sampling à 2.. : Conider a inuoidal ignal and let u ample it at a frequency F 2kHz. xt 3co000t 0. a) Determine and expreion for the ampled
More informationSimple Observer Based Synchronization of Lorenz System with Parametric Uncertainty
IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih
More informationSERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48)
Chapter 5 SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lecture 41-48) 5.1 Introduction Power ytem hould enure good quality of electric power upply, which mean voltage and current waveform hould
More informationMath Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK
ppendix 5 Scientific Notation It i difficult to work with very large or very mall number when they are written in common decimal notation. Uually it i poible to accommodate uch number by changing the SI
More informationOptimal Coordination of Samples in Business Surveys
Paper preented at the ICES-III, June 8-, 007, Montreal, Quebec, Canada Optimal Coordination of Sample in Buine Survey enka Mach, Ioana Şchiopu-Kratina, Philip T Rei, Jean-Marc Fillion Statitic Canada New
More informationCONTROL OF INTEGRATING PROCESS WITH DEAD TIME USING AUTO-TUNING APPROACH
Brazilian Journal of Chemical Engineering ISSN 004-6632 Printed in Brazil www.abeq.org.br/bjche Vol. 26, No. 0, pp. 89-98, January - March, 2009 CONROL OF INEGRAING PROCESS WIH DEAD IME USING AUO-UNING
More informationRobust Mould Level Control
5 American Control Conference June 8-1, 5. Portland, OR, USA ThA9.4 Robut Mould Level Control J. Schuurman, A. Kamperman, B. Middel, P.F.A van den Boch Abtract In the firt year of production ince, the
More informationThen C pid (s) S h -stabilizes G(s) if and only if Ĉpid(ŝ) S 0 - stabilizes Ĝ(ŝ). For any ρ R +, an RCF of Ĉ pid (ŝ) is given by
9 American Control Conference Hyatt Regency Riverfront, St. Loui, MO, USA June -, 9 WeC5.5 PID Controller Synthei with Shifted Axi Pole Aignment for a Cla of MIMO Sytem A. N. Gündeş and T. S. Chang Abtract
More informationAvoiding Forbidden Submatrices by Row Deletions
Avoiding Forbidden Submatrice by Row Deletion Sebatian Wernicke, Jochen Alber, Jen Gramm, Jiong Guo, and Rolf Niedermeier Wilhelm-Schickard-Intitut für Informatik, niverität Tübingen, Sand 13, D-72076
More informationLecture 8: Period Finding: Simon s Problem over Z N
Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing
More informationEfficient Methods of Doppler Processing for Coexisting Land and Weather Clutter
Efficient Method of Doppler Proceing for Coexiting Land and Weather Clutter Ça gatay Candan and A Özgür Yılmaz Middle Eat Technical Univerity METU) Ankara, Turkey ccandan@metuedutr, aoyilmaz@metuedutr
More informationUnified Design Method for Flexure and Debonding in FRP Retrofitted RC Beams
Unified Deign Method for Flexure and Debonding in FRP Retrofitted RC Beam G.X. Guan, Ph.D. 1 ; and C.J. Burgoyne 2 Abtract Flexural retrofitting of reinforced concrete (RC) beam uing fibre reinforced polymer
More informationTHE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY STABILIZING CONTROLLERS. Tatsuya Hoshikawa, Kou Yamada and Yuko Tatsumi
International Journal of Innovative Computing, Information Control ICIC International c 206 ISSN 349-498 Volume 2, Number 2, April 206 pp. 357 370 THE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY
More information