A Subspace-based Robust Adaptive Capon Beamforming
|
|
- Hugh Welch
- 6 years ago
- Views:
Transcription
1 374 Progress In Electromgnetics Reserch Symposium 2006, Cmbridge, USA, Mrch A Subspce-bsed Robust Adptive Cpon Bemforg G. S. Lio, H. Q. Liu, nd J. Li Xidin University, Chin Abstrct Adptive bemforg suffers from performnce degrdtion in the presence of mismtch between the ctul nd presumed rry steering vector of the desired signl. This ide enlightens us, so we propose subspce pproch to dptive bemforg tht is robust to rry errors bsed on imizing MUSIC output power. The proposed method involves two steps, the first step is to estimte the ctul steering vector of the desired signl bsed on subspce technique, nd the second is to obtin optiml weight by utilizing the estimted steering vector. Our method belongs to the clss of digonl loding, but the optiml mount of digonl loding level cn be clculted precisely bsed on the uncertinty set of the steering vector. To obtin noise subspce needs eigen-decomposition tht hs hevy computtion lod nd knows the number of signls priori. In order to overcome this drwbck we utilized the POR (Power of R) technique tht cn obtin noise subspce without eigen-decomposition nd the number of signls priori. It is very interesting tht Li Jin s method is specil cse where m = 1, nd the proposed subspce pproch is the cse where m, so we obtined uniform frmework bsed on POR technique. This is lso n explntion why the performnce of the proposed subspce pproch excels tht of Li Jin s method. The excellent performnce of our lgorithm hs been demonstrted vi number of computer simultions. 1. Introduction Arry signl processing hs wide pplictions in rdr, communictions, sonr, coustics, seismology, nd medicine. One of the importnt tsks of rry processing is bemforg. The stndrd bemformers include the dely-nd-sum pproch, which is known to suffer from poor resolution nd high sidelobe problems. The Cpon bemformer dptively selects the weight vectors to imize the rry output power subject to the liner constrint tht the signl of interest (SOI) does not suffer from ny distortion [1]. The Cpon bemformer hs better resolution nd interference rejection cpbility thn the stndrd bemformer, provided tht the rry steering vector corresponding to the SOI is ccurtely known. In prctice, the knowledge of the SOI steering vector my be imprecise, the cse due to differences between the presumed signl steering vector nd the ctul signl steering vector. When this hppens, the Cpon bemforg my suppress the SOI s n interference, which result in rry performnce drsticlly reduced, especilly rry output signl-to-interference-plus-noise rtio (SINR) [4]. In the pst three decdes mny pproches hve been proposed to improve the robustness of the Cpon bemforg. Additionl liner constrints, including point nd derivtive constrints, hve been imposed to improve the robustness of the Cpon bemforg [2, 3]. However, for every dditionl liner constrints imposed, the bemformer loses one degree of freedom (DOF) for interference suppression. Moreover, these constrints re not explicitly relted to the uncertinty of the rry steering vector. Digonl loding (including its extended versions) hs been populr pproch to improve the robustness of the Cpon bemformer [4]. However, for most of the digonl loding methods, detering the digonl loding res n open problem. Recently there re some methods been proposed (for exmples, [5 7] nd reference therein) to this point. Mismtch between the presumed steering vector of the SOI nd the ctul one result in drsticlly reduced rry SINR, therefore if we cn estimte ctul steering vector of the SOI, robustness of the rry will be improved. In this pper, from the point of view of the subspce we propose novel robust Cpon bemformer, which involves two steps, the first step is to estimte ctul steering vector of SOI, nd the second is to clculte optiml weight by Cpon method. The rest of this pper is orgnized s follows. Section 2 contins bckground mteril. In section 3, the robust Cpon bemformer is developed. Computer simultion results illustrting the performnce of the robust Cpon bemformer re presented in Section 4. Finlly, Section 5 contins the conclusions. 2. Bckground 2.1. Signl Model We consider the stndrd nrrowbnd bemforg model in which set of M nrrowbnd plne wve signls, impinge on n rry of N sensors with hlf wvelength spcing, where M < N. The N 1 vector of received signls is given by
2 Progress In Electromgnetics Reserch Symposium 2006, Cmbridge, USA, Mrch x(t k ) = M 1 m=0 (θ m )s m (t k ) + n(t k ), k = 1,2,...,L (1) where s m (t k ),m = 0,...,M 1; k = 1,2,...,L re the source signls snpshots, (θ m ) = [1,e jπ sin θm,...,e jπ(n 1) sin θm ] T is the steering vector in the direction θ m, nd n(t k ), k = 1,2,...,L re the vectors contining dditive white noise smples, L is the number of the snpshots. Also, in this pper, the sources nd noise re ssumed to be sttisticlly uncorrelted. We ssume tht one of the signls is the desired signl, sy s 0 (t), nd tret the reing signls s interferences. Since s 0 (t) is uncorrelted with the noise nd interferences, the dt covrince mtrix hs the form, M 1 R = σ0(θ 2 0 ) H (θ 0 ) + σk(θ 2 k ) H (θ k ) + R n R s + R i+n (2) k=1 where R s = σ0(θ 2 0 ) H (θ 0 ), σi 2 = E{ s i (t k ) 2} is the power of ith signl, nd R i+n is the interference plus noise covrince mtrix. In prctice, the covrince mtrix R is estimted by ˆR = 1 L x n x H n (3) L where ll received signls hve zero mens nd L smples re independent. n= Cpon Bemforg The Cpon bemforg is s follows. Detere the N 1 vector w 0 tht is the solution to the following linerly constrined qudrtic imiztion problem, w wh Rw s.t.w H ā(θ 0 ) = 1 (4) where ā(θ 0 ) is presumed steering vector of the desired signl. Appling Lgrnge multiplier method results in the following solution, R 1 ā(θ 0 ) w 0 = ā H (θ 0 )R 1 (5) ā(θ 0 ) The rry men output power p 0 is 1 p 0 = ā H (θ 0 )R 1 ā(θ 0 ) The Cpon bemformer hs better resolution nd much better interference rejection cpbility thn the stndrd bemformer, provided tht the presumed rry steering vector of the SOI mtch ctul rry steering vector precisely. In prctice, the exct steering vector of the SOI is unvilble or its mesure/estimtion is imprecise, therefore, we only use the presumed ā(θ 0 ) insted of the ctul (θ 0 ) in the Cpon bemformer, nd the mismtch between the exct steering vector nd the presumed one my drsticlly degrde the performnce of the Cpon bemformer. The rry output SINR cn be written s, SINR = E[ wh 0 s 0 (t) 2 ] w H 0 R i+nw 0 = where σ0 2 = E( s 0 (t) ). Inserting (5) into (7) yields, ā SINR = σ0 2 H (θ 0 )R 1 i+n (θ 0) 2 ā H (θ 0 )R 1 i+nā(θ 0) where (θ 0 ) is the ctul steering vector, then (8) cn be rewritten s: SINR = σ 2 0 H (θ 0 )R 1 i+n (θ 0) w H 0 (6) σ0 w 2 H 0 (θ 0 ) 2 ( M 1 ) (7) σk 2(θ k) H (θ k ) + R n w 0 k=1 H (θ 0 )R 1 i+nā(θ 0) 2 ( H (θ 0 )R 1 i+n (θ 0))(ā H (θ 0 )R 1 i+nā(θ 0)) = SINR m cos 2 ((θ 0 ),ā(θ 0 );R 1 i+n ) (9) (8)
3 376 Progress In Electromgnetics Reserch Symposium 2006, Cmbridge, USA, Mrch where SINR m = σ0 2 H (θ 0 )R 1 i+n (θ 0) nd cos 2 ( ) is defined s, cos 2 H Zb 2 (,b;z) = ( H Z)(b H Zb) Clerly, 0 cos 2 (,b;z) 1. Therefore, rry output SINR is reduced due to mismtch between the presumed steering vector of the SOI nd its true vlue. In recent yers, digonl loding (DL) is populr pproch to improving the robustness of Cpon bemformer to the mismtch bove. In DL methods, the dt covrince ˆR is replced by ˆR+γI, where γ is positive constnt (see reference [4 6] for detils). The DL method proposed in [4] is used in Section 4 for comprisons. In the following section, novel robust bemforg is developed to llevite the effects of the steering vector mismtch on the SINR performnce of Cpon bemformer. 3. Robust Cpon Bemforg The robust bemforg problem we will del with in this pper cn be briefly stted s follows: Extend the Cpon bemformer so s to improve rry output SINR even only n imprecise knowledge of steering vector (θ 0 ) is vilble. To simplify the nottion, in wht follows, we sometimes omit the rgument θ of (θ) nd ā(θ). We ssume tht the only knowledge we hve bout (θ 0 ) is tht it belongs to the following uncertinty [5] [(θ 0 ) ā] H C 1 [(θ 0 ) ā] 1 (11) where C re given positive definite mtrix. As shown bove, rry performnce loses will occur in the presence of mismtch between the presumed nd ctul steering vectors of the SOI. If we estimte the ctul steering vector of the SOI s more precise s we cn, then performnce of the bemformer will be improved. The proposed robust Cpon bemforg is bsed on this ide. From subspce theory we know tht the ctul steering vector of desired signl is orthogonl to noise subspce, our pproch is bsed on the optimizing the projection of signl steering vector onto noise subspce. The steering vector is normed s 2 = H = N. To derive our robust Cpon bemformer, we use following constrined optimiztion H U n U H n (10) s.t.( ā) H C 1 ( ā) 1 (12) 2 = N where ā is known to us in dvnce, but hs error (mismtch to the ctul steering vector of the SOI). U n is the noise subspce, which is obtined by the eigen-decomposition of ˆR. To mke up the noise subspce, we ssume tht the number M, of plne wves impinging on the rry is known priori. We use this ssumption only for derivtions nd cncel it lter. Note tht we cn improve the estimtion ccurcy of the ctul steering vector of the SOI from (12), nd then obtin optiml weight w 0 by Cpon method. Without loss of generlity, we will consider solving (12) for the cse in which C = εi, (ε is user prmeter), then, (12) becomes H U n U H n s.t. ā 2 ε (13) 2 = N We use the Lgrnge multiplier methodology gin, which is bsed on the function L(,λ,µ) = H U n U H n + µ(2n ε ā H H ā) + λ( H N) (14) where µ 0, λ 0 re the Lgrnge multiplier. Hence, the unconstrined imiztion of (14) for fixed µ, λ, is given by Clerly, the optiml solution of is δl(,µ,λ) δ = 2U n U H n 2µā + 2λ = 0 (15) â = µ(u n U H n + λi) 1 ā (16)
4 Progress In Electromgnetics Reserch Symposium 2006, Cmbridge, USA, Mrch Inserting â into (14), imizing L(,λ,µ) with respect to µ gives δl(â,µ,λ) = 2N ε ā H â â H ā = 0 δµ (17) Then, we obtin 2N ε ˆµ = 2ā H (U n U H n + λi) 1 ā (18) Inserting ˆµ into (14), imizing Lgrnge function with respect to λ yields nd the following eqution cn be derived, δl(â, ˆµ,λ) δλ = â H â N = 0 (19) ā H (U n U H n + ˆλI) 2 ā [ā H (U n U H n + ˆλI) 1 ā] 2 = N (N ε 2 )2 (20) Then, the solution of ˆλ cn be obtined by some simple mnipultions. Substituting (18) into (16) yields â = (N ε 2 ) (U n U H n + ˆλI) 1 ā ā H (U n U H n + ˆλI) 1 ā To summrize, the proposed robust Cpon bemforg consists of following steps. The lgorithm: Step 1: Clculte dt covrince mtrix, i.e., ˆR = 1 L x n x H n L Step 2: Compute the eigen-decomposition of ˆR nd obtin the noise subspce U n. Step 3: Solve ˆλ in (20). Step 4: Use the ˆλ in Step 3 to clculte â = (N ε 2 ) (U n U H n + ˆλI) 1 ā ā H (U n U H n + ˆλI) 1 ā n=1 (21) (22) Step 5: Compute optiml weight by Cpon method, i.e., w 0 = α ˆR 1 1 â, α = (23) â H ˆR 1â The proposed robust bemforg belongs to the clss of digonl loding, but the optiml mount of digonl loding level cn be precisely clculted bsed on the uncertinty set of the steering vector. In the Section 4 computer simultion results demonstrte excellent performnce of the proposed lgorithm. In order to void eigen-decomposition nd knowing the number of signls priori, we use the POR pproch to obtin noise subspce. In [8], R is decomposed by EVD s R = [U s U n ] [ ] [ Λs + σvi 2 0 U H s 0 σvi 2 where Λ s = dig{δ1,...,δ 2 M 2 }, U s denotes the signl subspce. It pproximtes the noise subspce of R bsed on R m (m is positive integer). Accordingly {( σv 2m R m = U n U H σ 2 ) m } v n + U s dig δi 2 + U H σ2 s (25) v Clerly, (σ 2 v/(δ 2 i +σ2 v)) m is less thn 1 nd converge to zero for sufficiently lrge m. Theoreticlly, lim m σ 2m v R m = U n U H n. As result, we modify the criterion (12) nd consider the following POR cost function U H n ] (24) H ˆR m s.t. ā 2 ε (26) 2 = N
5 378 Progress In Electromgnetics Reserch Symposium 2006, Cmbridge, USA, Mrch By contrst, the (26) voids estimting tht dimension directly. Moreover, s m, the proposed the POR bemforg method in (26) converges to the subspce one in (12), nd it cn be shown tht the performnce of the POR method for finite m will converge to the subspce one through computer simultion. We compred our method with previous one in [6], where m = 1 in the section 4. () (b) Figure 1: Output SINR versus different SNR, pointing errors = 3, for () ε = 0.7, for (b) ε = 7. () (b) Figure 2: The Output SINR versus pointing errors for () ε = 0.7, for (b) ε = Computer Results Our motivtion of simultion is to demonstrte the performnce in the presence of some errors in the steering vector. In ll of the exmples considered below, we ssume uniform liner rry (ULA) with N = 20 sensors nd hlf-wvelength spcing is used. The sources emitted mutul independent nrrowbnd wveforms. All the results re chieved vi 50 Monte Crlo trils. In the first exmple, we consider the effect of the pointing error of the SOI on rry output SINR. The exct direction of rrivl of SOI is θ 0, of which ssumed vlue is θ 0 +, i.e., ā(θ 0 ) = (θ 0 + ). We ssume tht (θ 0 ) belongs to the uncertinty set (θ 0 ) ā(θ 0 ) 2 ε (27)
6 Progress In Electromgnetics Reserch Symposium 2006, Cmbridge, USA, Mrch where ε is user prmeter. Let ε 0 = (θ 0 ) ā(θ 0 ) 2. Then, choosing ε = ε 0. However, since is unknown in prctice, the ε we choose my be greter or less thn ε 0. To show tht the choice of ε is not criticl issue for our lgorithm, we will present simultion results with severl vlues of ε in eqution (21). In this exmple, the directions of the SOI nd n interference source re θ 0 = 30, θ 1 = 30, respectively. The ssumed direction of the SOI is θ 0 + = 33, which results exct ε 0 = The interference-to-noise rtio (INR) is 40 db. Figure 1 plots rry output SINR versus the SNR of the SOI when the number of snpshots is set to be L = 100. It is observed tht the proposed lgorithm (12) performs better thn other two lgorithms t ll input SNR. Also, since the error in steering vector of SOI is reltively lrge nd cnnot be negligible, the stndrd Cpon bemformer nd its digonl loding version suffer from severe performnce degrdtion when SNR increses. However, the proposed bemformer hs SINR loss of 5 db when SNR = 20 db. The proposed the POR method for different m over vrious input SNRs is lso illustrted in Figure 1. Obviously, the Output SINR for m = 2 nd m = 3 ll converge to subspce pproch (12), the counterprt for m = 1 [6] hs the lrge output SINR loss. Figure 2 shows the rry output SINR curve versus the pointing errors, in which SNR = 0 db, INR = 20 db, L = 100. In this figure, the excellent performnce chieved by the proposed lgorithm is observed, which shows the robustness to the pointing errors. It is noted tht, similr to other robust pproches, our method will worsen if there is/re strong interference sptilly closed to the SOI. The reson is tht for given uncertinty region (11), the solution of in optimiztion (12) is converge to the strong interference source. Also, it cn be seen tht the Output SINR of the proposed POR method increses s m increses, with m = 3 hs sme performnce with subspce one (12). 5. Conclusion In this pper, we discuss the performnce degrdtion due to the presence of steering vector uncertinty of the SOI, such s, direction of rrivl estimtion error, finite number of snpshots, nd rry response error, etc. A robust Cpon bemformer is developed by utilizing the orthogonlity between signl nd noise subspce. A more ccurte estimte of the ctul steering vector of the SOI is obtined vi constrined optimiztion, by which the optiml weight is computed ccording to Cpon bemforg. We hve shown tht the proposed lgorithm belongs to the clss of digonl loding pproches, nd the optiml mounts of digonl loding cn be precisely clculted. In order to void eigen-decomposition nd knowing the number of signls priori, we hve proposed POR-bsed robust bemforg scheme. It significntly outperforms the method proposed in [6] nd converge to the subspce one (12). The excellent performnce of our lgorithm hs been demonstrted vi number of computer simultions. Acknowledgment The work described in this pper ws supported by Ntionl Science Fund under grnt NFS REFERENCES 1. Cpon, J., High resolution frequency-wvenumber spectrum nlysis, Proc. IEEE, Vol. 57, , Aug Er, M. H. nd A. Cntoni, Derivtive constrints for brod-bnd element spce ntenn rry processor, IEEE Trns. Acoust., Speech, Signl Processing, Vol. ASSP-31, , Dec Buckley, K. M. nd L. J. Griffiths, An dptive generlized sidelobe cnceller with derivtive constrints, IEEE Trns. Antenns Propgt., Vol. AP-34, , Mr Crlson, B. D., Covrince mtrix estimtion errors nd digonl loding in dptive rrys, IEEE Trns. Aerospce nd Electronic System., Vol. 24, , Jul Stoic, P., Z.-S. Wng, nd J. Li, Robust cpon bemforg, IEEE Signl Processing Letters, Vol. 10, No. 6, , June Li, J., P. Stoic, nd Z.-S. Wng, Doubly constrined robust cpon bemforg, IEEE Trns. Signl Processing, Vol. 52, No. 9, , Sep Shhbzpnhi, S., A. B. Gershmn, Z.-Q. Luo, nd K. M. Wong, Robust dptive bemforg for generl-rnk signl models, IEEE Trns. Signl Processing, Vol. 51, No. 9, , Sep Xu, Z., P. Liu, nd X. Wng, Blind multiuser detection: from moe to subspce methods, IEEE Trns. Signl Processing, Vol. 52, No. 2, , Feb
ROBUST ADAPTIVE BEAMFORMER WITH FEASIBILITY CONSTRAINT ON THE STEERING VECTOR
ROBUST ADAPTIVE BEAMFORMER WITH FEASIBILITY CONSTRAINT ON THE STEERING VECTOR Wenyi Zhng nd Bhskr D. Ro Deprtment of Electricl nd Computer Engineering University of Cliforni, Sn Diego L Joll, CA 92093-0407,
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
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 informationDiscrete Least-squares Approximations
Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve
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 informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationMath Lecture 23
Mth 8 - Lecture 3 Dyln Zwick Fll 3 In our lst lecture we delt with solutions to the system: x = Ax where A is n n n mtrix with n distinct eigenvlues. As promised, tody we will del with the question of
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 informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
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 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 informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
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 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 informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationIntroduction to Hybrid Beamforming Techniques
Introduction to ybrid Bemforming Techniques Jmes Chen Advisor : Andy Wu Grdute Institute of Electronics Engineering Ntionl Tiwn University Tipei, Tiwn Mr 3, 205 2 Outline Introduction of Precoding Why
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 informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More information1 Linear Least Squares
Lest Squres Pge 1 1 Liner Lest Squres I will try to be consistent in nottion, with n being the number of dt points, nd m < n being the number of prmeters in model function. We re interested in solving
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 information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More information1 The Lagrange interpolation formula
Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x
More informationFredholm Integral Equations of the First Kind Solved by Using the Homotopy Perturbation Method
Int. Journl of Mth. Anlysis, Vol. 5, 211, no. 19, 935-94 Fredholm Integrl Equtions of the First Kind Solved by Using the Homotopy Perturbtion Method Seyyed Mhmood Mirzei Deprtment of Mthemtics, Fculty
More informationMatrices, Moments and Quadrature, cont d
Jim Lmbers MAT 285 Summer Session 2015-16 Lecture 2 Notes Mtrices, Moments nd Qudrture, cont d We hve described how Jcobi mtrices cn be used to compute nodes nd weights for Gussin qudrture rules for generl
More informationModule 6: LINEAR TRANSFORMATIONS
Module 6: LINEAR TRANSFORMATIONS. Trnsformtions nd mtrices Trnsformtions re generliztions of functions. A vector x in some set S n is mpped into m nother vector y T( x). A trnsformtion is liner if, for
More informationNotes on length and conformal metrics
Notes on length nd conforml metrics We recll how to mesure the Eucliden distnce of n rc in the plne. Let α : [, b] R 2 be smooth (C ) rc. Tht is α(t) (x(t), y(t)) where x(t) nd y(t) re smooth rel vlued
More informationJim Lambers MAT 169 Fall Semester Lecture 4 Notes
Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of
More informationLecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationElements of Matrix Algebra
Elements of Mtrix Algebr Klus Neusser Kurt Schmidheiny September 30, 2015 Contents 1 Definitions 2 2 Mtrix opertions 3 3 Rnk of Mtrix 5 4 Specil Functions of Qudrtic Mtrices 6 4.1 Trce of Mtrix.........................
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 informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More informationMath 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of
More informationS. S. Dragomir. 2, we have the inequality. b a
Bull Koren Mth Soc 005 No pp 3 30 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Abstrct Compnions of Ostrowski s integrl ineulity for bsolutely
More informationOrthogonal Polynomials and Least-Squares Approximations to Functions
Chpter Orthogonl Polynomils nd Lest-Squres Approximtions to Functions **4/5/3 ET. Discrete Lest-Squres Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationNumerical Linear Algebra Assignment 008
Numericl Liner Algebr Assignment 008 Nguyen Qun B Hong Students t Fculty of Mth nd Computer Science, Ho Chi Minh University of Science, Vietnm emil. nguyenqunbhong@gmil.com blog. http://hongnguyenqunb.wordpress.com
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 informationThe Velocity Factor of an Insulated Two-Wire Transmission Line
The Velocity Fctor of n Insulted Two-Wire Trnsmission Line Problem Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 Mrch 7, 008 Estimte the velocity fctor F = v/c nd the
More informationUnit #9 : Definite Integral Properties; Fundamental Theorem of Calculus
Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl
More informationdt. However, we might also be curious about dy
Section 0. The Clculus of Prmetric Curves Even though curve defined prmetricly my not be function, we cn still consider concepts such s rtes of chnge. However, the concepts will need specil tretment. For
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 informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationAdaptive beamforming for uniform linear arrays with unknown mutual coupling. IEEE Antennas and Wireless Propagation Letters.
Title Adaptive beamforming for uniform linear arrays with unknown mutual coupling Author(s) Liao, B; Chan, SC Citation IEEE Antennas And Wireless Propagation Letters, 2012, v. 11, p. 464-467 Issued Date
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 information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
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 INTRODUCTION TO UNIFORM CONVERGENCE. In the study of Fourier series, several questions arise naturally, such as: c n e int
A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE HANS RINGSTRÖM. Questions nd exmples In the study of Fourier series, severl questions rise nturlly, such s: () (2) re there conditions on c n, n Z, which ensure
More informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number 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 informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
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 information13: Diffusion in 2 Energy Groups
3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups
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 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 informationCalculus of Variations
Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function
More informationLecture Outline. Dispersion Relation Electromagnetic Wave Polarization 8/7/2018. EE 4347 Applied Electromagnetics. Topic 3c
Course Instructor Dr. Rymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mil: rcrumpf@utep.edu EE 4347 Applied Electromgnetics Topic 3c Wve Dispersion & Polriztion Wve Dispersion These notes & Polriztion
More informationLecture 21: Order statistics
Lecture : Order sttistics Suppose we hve N mesurements of sclr, x i =, N Tke ll mesurements nd sort them into scending order x x x 3 x N Define the mesured running integrl S N (x) = 0 for x < x = i/n for
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 informationWHEN finding the directions-of-arrival (DOA s) of narrowband
2166 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL 46, NO 8, AUGUST 1998 Modeling nd Estimtion of Ambiguities in Liner Arrys Athnssios Mniks, Member, IEEE, nd Christos Proukkis Abstrct In this pper, the
More informationNumerical Integration. 1 Introduction. 2 Midpoint Rule, Trapezoid Rule, Simpson Rule. AMSC/CMSC 460/466 T. von Petersdorff 1
AMSC/CMSC 46/466 T. von Petersdorff 1 umericl Integrtion 1 Introduction We wnt to pproximte the integrl I := f xdx where we re given, b nd the function f s subroutine. We evlute f t points x 1,...,x n
More informationEnergy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon
Phys463.nb 49 7 Energy Bnds Ref: textbook, Chpter 7 Q: Why re there insultors nd conductors? Q: Wht will hppen when n electron moves in crystl? In the previous chpter, we discussed free electron gses,
More informationW. We shall do so one by one, starting with I 1, and we shall do it greedily, trying
Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)
More informationRecitation 3: Applications of the Derivative. 1 Higher-Order Derivatives and their Applications
Mth 1c TA: Pdric Brtlett Recittion 3: Applictions of the Derivtive Week 3 Cltech 013 1 Higher-Order Derivtives nd their Applictions Another thing we could wnt to do with the derivtive, motivted by wht
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 information20 MATHEMATICS POLYNOMIALS
0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of
More informationMATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1
MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further
More information63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1
3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =
More informationEuler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )
Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationMain topics for the First Midterm
Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 2-3, Sections 4.1-4.8, nd Sections 5.1-5.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of the
More informationMultivariate problems and matrix algebra
University of Ferrr Stefno Bonnini Multivrite problems nd mtrix lgebr Multivrite problems Multivrite sttisticl nlysis dels with dt contining observtions on two or more chrcteristics (vribles) ech mesured
More information7.2 The Definite Integral
7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where
More informationINFORMATION THEORY AND ELECTROMAGNETISM: ARE THEY RELATED?
INFORMATION THEORY AND ELECTROMAGNETISM: ARE THEY RELATED? Sergey Loyk School of Informtion Technology nd Engineering (SITE), University of Ottw 161 Louis Psteur, Ottw, Ontrio, Cnd, K1N 6N5, emil: sergey.loyk@ieee.org
More informationUndergraduate Research
Undergrdute Reserch A Trigonometric Simpson s Rule By Ctherine Cusimno Kirby nd Sony Stnley Biogrphicl Sketch Ctherine Cusimno Kirby is the dughter of Donn nd Sm Cusimno. Originlly from Vestvi Hills, Albm,
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 informationCAAM 453 NUMERICAL ANALYSIS I Examination There are four questions, plus a bonus. Do not look at them until you begin the exam.
Exmintion 1 Posted 23 October 2002. Due no lter thn 5pm on Mondy, 28 October 2002. Instructions: 1. Time limit: 3 uninterrupted hours. 2. There re four questions, plus bonus. Do not look t them until you
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 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 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 informationChapter 1. Chapter 1 1
Chpter Chpter : Signls nd Systems... 2. Introduction... 2.2 Signls... 3.2. Smpling... 4.2.2 Periodic Signls... 0.2.3 Discrete-Time Sinusoidl Signls... 2.2.4 Rel Exponentil Signls... 5.2.5 Complex Exponentil
More informationQuadratic Forms. Quadratic Forms
Qudrtic Forms Recll the Simon & Blume excerpt from n erlier lecture which sid tht the min tsk of clculus is to pproximte nonliner functions with liner functions. It s ctully more ccurte to sy tht we pproximte
More informationDuality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.
Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we
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 informationIntegrals - Motivation
Integrls - Motivtion When we looked t function s rte of chnge If f(x) is liner, the nswer is esy slope If f(x) is non-liner, we hd to work hrd limits derivtive A relted question is the re under f(x) (but
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 informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk bout solving systems of liner equtions. These re problems tht give couple of equtions with couple of unknowns, like: 6 2 3 7 4
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 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 informationLecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.
Lecture 3 3 Solving liner equtions In this lecture we will discuss lgorithms for solving systems of liner equtions Multiplictive identity Let us restrict ourselves to considering squre mtrices since one
More informationLecture 20: Numerical Integration III
cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More informationA027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data
A07 Uncertinties in Locl Anisotropy Estimtion from Multi-offset VSP Dt M. Asghrzdeh* (Curtin University), A. Bon (Curtin University), R. Pevzner (Curtin University), M. Urosevic (Curtin University) & B.
More informationNumerical Integration
Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl
More informationModule 11.4: nag quad util Numerical Integration Utilities. Contents
Qudrture Module Contents Module 11.4: ng qud util Numericl Integrtion Utilities ng qud util provides utility procedures for computtion involving integrtion of functions, e.g., the computtion of the weights
More informationBases for Vector Spaces
Bses for Vector Spces 2-26-25 A set is independent if, roughly speking, there is no redundncy in the set: You cn t uild ny vector in the set s liner comintion of the others A set spns if you cn uild everything
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationFor the percentage of full time students at RCC the symbols would be:
Mth 17/171 Chpter 7- ypothesis Testing with One Smple This chpter is s simple s the previous one, except it is more interesting In this chpter we will test clims concerning the sme prmeters tht we worked
More informationDefinite integral. Mathematics FRDIS MENDELU
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the
More information