HW #4 Solutions: M697J Fall 2006
|
|
- Darren Dennis
- 6 years ago
- Views:
Transcription
1 HW #4 Solutios: M697J Fall (4.2 BHM) Eigevalue calculatios i two dimesios. Cosider the weighted Jacobi method applied to the model Poisso equatio i two dimesios o the uit square. Assume a uiform grid with h = i each directio. (a) Let v ij be the approximatio to the solutio at the grid poit (x i,y j ). Write the (i,j)th equatio of the correspodig discrete problem, where i,j. (b) Lettig A be the matrix of coefficiets for the discrete system, write the (i,j)th equatio for the eigevalue problem Av = λv. (c) Assume ad eigevector solutio of the form ( ) ikπ v ij = si si ( jlπ ), k,l. Usig sie additio rules, simplify this eigevalue equatio, cacel commo terms, ad show that the eigevalues are [ ( ) ( )] kπ lπ λ kl = 4 si 2 + si 2, k,l. 2 2 (d) As i the oe-dimesioal case, ote that the iteratio matrix of the weighted Jacobi method is give by P ω = I ωd A, where D correspods to the diagoal terms of A. Fid the eigevalues of P ω. (e) Usig a graphig utility, fid a suitable way to preset the two-dimesioal set of eigevalues (either a surface plot or multiple curves). Plot the eigevalues for ω = 2/3,4/5, ad = 6. (f) I each case, discuss the effect of the weighted Jacobi method o the low- ad high-frequecy modes. Be sure to ote that modes ca have high frequecies i oe directio ad low frequecies i the other directio. (g) What do you coclude about the optimal value of ω for the two-dimesioal problem? ANS: (a) We approximate the model problem u = f, u Γ = 0 by Dx 2v D2 yv = f, or v i,j + 2v i,j v i+,j h 2 + v i,j + 2v i,j v i,j+ h 2 = f i,j. (b) Note: give that there is o divisio by h 2 i the defitio of λ kl i part (c) it is ot icluded here. It may be that BHM meat the homogeeous problem i part (a). I either case, we have ( v i,j + 2v i,j v i+,j ) + ( v i,j + 2v i,j v i,j+ ) = λv i,j. (c) Due to liearity we ca treat each term o the left i (b) separately. The oe has ( ( ) v i,j + 2v i,j v i+,j = si (i )kπ + 2si ( ) ( )) ikπ si (i+)kπ si ( jlπ ).
2 Note that With this, ( ) (i ± )kπ si = si v i,j + 2v i,j v i+,j ( ) ikπ cos = ( 2cos ( kπ ( ) kπ ± cos ( ) ikπ si ( ) kπ. ) ( si ikπ ) ( + 2si ikπ )) ( ) si jlπ = 2 ( 2cos ( )) ( kπ si ikπ ) ( ) si jlπ = 4si 2 ( ) ( kπ 2 si ikπ ) ( ) si jlπ = λ x si ( ikπ ) si ( jlπ ), or, λ x = 4si 2 ( kπ 2). Repeatig these steps for the vi,j + 2v i,j v i,j+ term, ad addig, gives the result. (d) Agai, due to liearity, the kl th eigevalue of P ω is (P ω ) kl = (I ωd A) kl = ω [ ( ) kπ 4 λ kl = ω si 2 2 (e) (images o ext page!) ( )] lπ + si 2 2 for k,l.
3
4 2. (4.5 BHM) Two-dimesioal program. For the two-dimesioal problem, proceed agai i a modular way: (a) Write a subroutie that performs weighted Jacobi o a two-dimesioal grid. Withi the subroutie, it is easiest to refer to v ad f as two-dimesioal arrays. (b) Make the appropriate modificatios to the oe-dimesioal code to implemet biliear iterpolatio ad full weightig o a two-dimesioal grid. (c) Make the (mior) chages required i the mai program to create a two-dimesioal V-cycle program. Test this program o problems with kow exact solutios. For example, for fixed k ad l, take f(x,y) = C si (kπx) si (lπy) o the uit square (0 x,y ) where C is a costat. The C u(x,y) = π 2 k 2 + π 2 l 2 si(kπx) si (lπy). + σ (c) Preset results for k = 0, l = 20 takig M = N = 32. Use the optimal Jacobi weight ω = 4/5. (d) With your code ad the Jacobi weight ω = 4/5, reproduce the = 64 portio of the table o page 65 of the coursebook, makig sure to use the discrete L 2 orm h. Recall for d-dimesios, if h = h x = h x2 = = h xd the ( u h h = h d i ) {u h i } 2 /2 ( ) = h d/2 {u h i } 2 /2 = h d/2 u h 2. So we see, i this case the discrete L 2 orm is simply the stadard 2-orm, 2, scaled. Also ote that you ca use Matlab s orm to compute u h 2. ANS: The V-cycle code is available o the course website. Here are the results for part (d): cycle(s) res rat err rat e e e e e e e e e e e e e e e e e e e e-06 i
5 3. Cosider the 2-poit oe-dimesioal BVP { u + u = (π 2 si πx 2π cos πx)e x u(0) = u() = 0. (a) Write a MATLAB script to solve the problem by the FFT method, usig the Discrete Sie Trasform as implemeted by dst.m applied to the 2d order cetered FD scheme, assumig σ 0 is a costat, D 2 v i + σv i = f i, where D 2 = D + D. Assume a meshsize h = /2 p, where p is a positive iteger. For p = : 4, plot the exact solutio (u(x) vs. x) ad the umerical solutio (v i vs. x i ), icludig the boudary poits. The 4 plots should appear separately i oe figure, with axes labeled ad a title for each idicatig p. Ivestigate subplot i MATLAB for how to have multiple plots i a sigle figure widow. (b) For p = : 5 preset a table with the followig data - colum : h; colum 2: u h v h ; colum 3: u h v h /h 2 ; colum 4: cpu time; colum 5: (cpu time)/( log ), where h = /. Discuss the treds i each colum. Iclude a copy of your code. ANS: Here is the code for the D solver: fuctio [u,x] = poissd_fft(n,sigma,f_fuc) Poisso solve usig DST applied to: -u + \sigma u = f, u(0)=u()=0 h = /N; x = (0:N) *h; lam = 2*(-cos(x(2:N)*pi))/(h^2); f = eval(f_fuc); u = dst(f(2:n)); trasform u = u./(lam+sigma); Solve i Fourier space u = dst(u); trasform back u = [0; u; 0]; set BC Here is the code for the computatios ad plots, followed by the result: xx = 0:0.0:;xx=xx ; uu = si(pi*xx).*exp(xx); sigma = ; f = (pi^2*si(pi*x)-2*pi*cos(pi*x)).*exp(x) ; clf; for p=:4 [v,x] = poissd_fft(2^p,sigma,f); subplot(2,2,p),plot(xx,uu,x,v, * ),grid leged( u, v,2) axis( tight ),xlabel( x ),title([ p=,um2str(p)]) ed
6 p= p=2.5 u v.5 u v x p= x p=4.5 u v.5 u v x x (b) Here is the code. Note that the solver is called 50 times ad the cpu time averaged i order to get a accurate timig. sigma = ; f = (pi^2*si(pi*x)-2*pi*cos(pi*x)).*exp(x) ; h = zeros(5,); = zeros(5,); times = zeros(5,); err_if = zeros(5,); for p = :5 p tic; for j=:50 [v,x] = poissd_fft(2^p,sigma,f); ed stime = toc; h(p) = /(2^p); (p) = 2^p; times(p) = stime/50; u = si(pi*x).*exp(x); err_if(p) = max(abs(u-v)); ed format short e disp( ) disp( h if_err err/h^2 cputime cpu/(log) ) disp( ) disp( ) disp([h err_if err_if./(h.^2) times times./(.*log2())])
7 Here are the results. I executed the code twice so that FFTW was properly iitialized for the vector legths required. h if_err err/h^2 cputime cpu/(log) e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e-0.996e e e e e e e e e e e e e e e e e e e e e e-07 We ca see from the err/h 2 colum that the expected O(h 2 ) error is observed util h.2207e 04 sice util this poit err/h 2 approaches a costat. But the we lose accuracy sice roudoff error begis to domiate. Thus, while theoretically as h 0 we have covergece, floatig poit errors i the form of roudoff error evetually domiates! Note, that the cputime does ot exhibit the O( log ) scalig util is large. This is due to the fact that there are also O() computatios which are a sigificat portio of the total time util is large.
8 4. Cosider the 2D BVP o [0,] [0,] { u + σu = f(x,y) u Γ = g(x,y). (a) Write a MATLAB script to solve the problem by the FFT method, usig the Discrete Sie Trasform as implemeted by dst.m applied to the 2d order cetered FD scheme, assumig σ 0 is a costat, (D 2 x + D2 y )v i,j + σv i,j = f i,j, where D 2 = D + D. Your code should be geeral eough to hadle ihomogeeous BC (g 0) ad M N, i.e., /M = h x h y = /N. (b) Usig your code, for σ =, p = 4 : 0 ad M = N = 2 p, compute the approximate solutio give u(x,y) = (2π) 2 + (3π) 2 ex+siy si2πx si 3πy. I leave it to you to compute the appropriate f. Preset a table with the followig data - colum : h; colum 2: u h v h h ; colum 3: u h v h h /h 2 ; colum 4: cpu time; colum 5: (cpu time)/ 2 log, where h x = h y = /. Discuss the treds i each colum. Iclude a copy of your code. ANS: Here is the code for the 2D solver: fuctio [u,x,y] = poiss2d_fft(f,sigma) 2D Poisso eq -(u_xx+u_yy)+ \sigma u = f o [0,]x[0,] Dirichlet BC for u are stored i the first/last-row/colum of f [M N] = size(f); M = M-; N= N-; u = zeros(m+,n+); hx = /M; x = (0:hx:) ; lamx = 2*(-cos(x(2:M)*pi))/(hx^2); hy = /N; y = (0:hy:) ; lamy = 2*(-cos(y(2:N)*pi))/(hy^2); Fourier multiple matrix: lambda_x(i)+lambda_y(j) lamx_p_lamy = (repmat(lamx,,n-) + repmat(lamy,m-,)) ; set RHS ad adjust for possible o-zero Dirichlet BC v = f(2:m,2:n); v(:, ) = v(:, ) + f(2:m, )/(hy^2); v(:,n-) = v(:,n-) + f(2:m,n+)/(hy^2); v(,:) = v(,:) + f(,2:n)/(hx^2); v(m-,:) = v(m-,:) + f(m+,2:n)/(hx^2); solve, put v i u, ad set BC v = dst(dst(dst(dst(v) )./(lamx_p_lamy+sigma)) ); u(2:m,2:n) = v; u(,:) = f(,:); u(m+,:) = f(m+,:); u(:,) = f(:,); u(:,n+) = f(:,n+);
9 Here are the results followed by the discussio ad code: h err_h err_h/h^2 cpu cpu/((^2)log) code: e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e-07 sigma = ; p = 4:0; lep = legth(p); h = zeros(lep,); = zeros(lep,); T = zeros(lep,); err_h = zeros(lep,); trials = 5; fac = /((2*pi)^2 + (3*pi)^2); for i = :lep disp(um2str(p(i))) M = 2^p(i); hx = /M; x = 0:hx:; N = 2^p(i); hy = /N; y = 0:hy:; [yy xx]= meshgrid(y,x); u = fac*exp(xx+si(yy)).*si(2*pi*xx).*si(3*pi*yy); f = -fac*exp(xx+si(yy)).*si(3*pi*yy).*... -u_xx ((-4*pi^2)*si(2*pi*xx)+4*pi*cos(2*pi*xx))... -fac*exp(xx+si(yy)).*si(2*pi*xx).*... -u_yy (-9*pi^2*si(3*pi*yy)+6*pi*cos(3*pi*yy).*cos(yy)... -si(3*pi*yy).*si(yy)+si(3*pi*yy).*cos(yy).^2) + u ; f(,:) = u(,:); f(m+,:) = u(m+,:); f(:,) = u(:,); f(:,n+) = u(:,n+); tic; for j=:trials uu = poiss2d_fft(f,sigma); ed stime = toc; h(i) = hx; (i) = M; T(i) = stime/trials; err_h(i) = hx*orm(uu(:)-u(:)); ed format short e disp( ) disp( h err_h err_h/h^2 cpu cpu/((^2)log) ) disp( ) disp( ) disp([h err_h err_h./(h.^2) T T./((.^2).*log2())])
TMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationHomework Set 4 - SOLUTIONS
MATH 55 Scientific Computing II Spring 6 Homework Set 4 - SOLUTIONS. For N even define a discrete grid x h for [, π] by x k = k h where h = π/n and k N. Show e i( N +j)x k = e i( N j)x k for j =,,...,
More informationIterative Techniques for Solving Ax b -(3.8). Assume that the system has a unique solution. Let x be the solution. Then x A 1 b.
Iterative Techiques for Solvig Ax b -(8) Cosider solvig liear systems of them form: Ax b where A a ij, x x i, b b i Assume that the system has a uique solutio Let x be the solutio The x A b Jacobi ad Gauss-Seidel
More informationPROBLEMS AND SOLUTIONS 2
PROBEMS AND SOUTIONS Problem 5.:1 Statemet. Fid the solutio of { u tt = a u xx, x, t R, u(x, ) = f(x), u t (x, ) = g(x), i the followig cases: (b) f(x) = e x, g(x) = axe x, (d) f(x) = 1, g(x) =, (f) f(x)
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationFFTs in Graphics and Vision. The Fast Fourier Transform
FFTs i Graphics ad Visio The Fast Fourier Trasform 1 Outlie The FFT Algorithm Applicatios i 1D Multi-Dimesioal FFTs More Applicatios Real FFTs 2 Computatioal Complexity To compute the movig dot-product
More informationPUTNAM TRAINING, 2008 COMPLEX NUMBERS
PUTNAM TRAINING, 008 COMPLEX NUMBERS (Last updated: December 11, 017) Remark. This is a list of exercises o Complex Numbers Miguel A. Lerma Exercises 1. Let m ad two itegers such that each ca be expressed
More informationThe Jordan Normal Form: A General Approach to Solving Homogeneous Linear Systems. Mike Raugh. March 20, 2005
The Jorda Normal Form: A Geeral Approach to Solvig Homogeeous Liear Sstems Mike Raugh March 2, 25 What are we doig here? I this ote, we describe the Jorda ormal form of a matrix ad show how it ma be used
More informationChapter 10 Partial Differential Equations and Fourier Series
Math-33 Chapter Partial Differetial Equatios November 6, 7 Chapter Partial Differetial Equatios ad Fourier Series Math-33 Chapter Partial Differetial Equatios November 6, 7. Boudary Value Problems for
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationZ - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.
Z - Trasform The -trasform is a very importat tool i describig ad aalyig digital systems. It offers the techiques for digital filter desig ad frequecy aalysis of digital sigals. Defiitio of -trasform:
More informationChapter 9 - CD companion 1. A Generic Implementation; The Common-Merge Amplifier. 1 τ is. ω ch. τ io
Chapter 9 - CD compaio CHAPTER NINE CD-9.2 CD-9.2. Stages With Voltage ad Curret Gai A Geeric Implemetatio; The Commo-Merge Amplifier The advaced method preseted i the text for approximatig cutoff frequecies
More informationApplication 10.5B Rectangular Membrane Vibrations
Applicatio.5B Rectagular Membrae Vibratios Here we ivestigate the vibratios of a flexible membrae whose equilibrium positio is the rectagle x a, y b. Suppose it is released from rest with give iitial displacemet,
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationSNAP Centre Workshop. Basic Algebraic Manipulation
SNAP Cetre Workshop Basic Algebraic Maipulatio 8 Simplifyig Algebraic Expressios Whe a expressio is writte i the most compact maer possible, it is cosidered to be simplified. Not Simplified: x(x + 4x)
More informationAnalysis of a Numerical Scheme An Example
http://www.d.edu/~gtryggva/cfd-course/ Computatioal Fluid Dyamics Lecture 3 Jauary 5, 7 Aalysis of a Numerical Scheme A Example Grétar Tryggvaso Numerical Aalysis Example Use the leap-frog method (cetered
More informationBasic Iterative Methods. Basic Iterative Methods
Abel s heorem: he roots of a polyomial with degree greater tha or equal to 5 ad arbitrary coefficiets caot be foud with a fiite umber of operatios usig additio, subtractio, multiplicatio, divisio, ad extractio
More informationMATH 10550, EXAM 3 SOLUTIONS
MATH 155, EXAM 3 SOLUTIONS 1. I fidig a approximate solutio to the equatio x 3 +x 4 = usig Newto s method with iitial approximatio x 1 = 1, what is x? Solutio. Recall that x +1 = x f(x ) f (x ). Hece,
More informationEigenvalues and Eigenvectors
5 Eigevalues ad Eigevectors 5.3 DIAGONALIZATION DIAGONALIZATION Example 1: Let. Fid a formula for A k, give that P 1 1 = 1 2 ad, where Solutio: The stadard formula for the iverse of a 2 2 matrix yields
More informationComputational Fluid Dynamics. Lecture 3
Computatioal Fluid Dyamics Lecture 3 Discretizatio Cotiued. A fourth order approximatio to f x ca be foud usig Taylor Series. ( + ) + ( + ) + + ( ) + ( ) = a f x x b f x x c f x d f x x e f x x f x 0 0
More informationINTEGRATION BY PARTS (TABLE METHOD)
INTEGRATION BY PARTS (TABLE METHOD) Suppose you wat to evaluate cos d usig itegratio by parts. Usig the u dv otatio, we get So, u dv d cos du d v si cos d si si d or si si d We see that it is ecessary
More informationMath 142, Final Exam. 5/2/11.
Math 4, Fial Exam 5// No otes, calculator, or text There are poits total Partial credit may be give Write your full ame i the upper right corer of page Number the pages i the upper right corer Do problem
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationWhere do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway?
Where do eigevalues/eigevectors/eigeuctios come rom, ad why are they importat ayway? I. Bacgroud (rom Ordiary Dieretial Equatios} Cosider the simplest example o a harmoic oscillator (thi o a vibratig strig)
More informationDefinition of z-transform.
- Trasforms Frequecy domai represetatios of discretetime sigals ad LTI discrete-time systems are made possible with the use of DTFT. However ot all discrete-time sigals e.g. uit step sequece are guarateed
More informationSolution of EECS 315 Final Examination F09
Solutio of EECS 315 Fial Examiatio F9 1. Fid the umerical value of δ ( t + 4ramp( tdt. δ ( t + 4ramp( tdt. Fid the umerical sigal eergy of x E x = x[ ] = δ 3 = 11 = ( = ramp( ( 4 = ramp( 8 = 8 [ ] = (
More informationPolynomial Functions and Their Graphs
Polyomial Fuctios ad Their Graphs I this sectio we begi the study of fuctios defied by polyomial expressios. Polyomial ad ratioal fuctios are the most commo fuctios used to model data, ad are used extesively
More informationMAXIMALLY FLAT FIR FILTERS
MAXIMALLY FLAT FIR FILTERS This sectio describes a family of maximally flat symmetric FIR filters first itroduced by Herrma [2]. The desig of these filters is particularly simple due to the availability
More informationComplex Analysis Spring 2001 Homework I Solution
Complex Aalysis Sprig 2001 Homework I Solutio 1. Coway, Chapter 1, sectio 3, problem 3. Describe the set of poits satisfyig the equatio z a z + a = 2c, where c > 0 ad a R. To begi, we see from the triagle
More informationThe picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled
1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how
More informationTaylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH
Taylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. A ay poit i the eighbourhood of 0, the fuctio ƒ() ca be represeted by a power series of the followig form: X 0 f(a) f() f() ( ) f( ) ( )
More informationMath 220B Final Exam Solutions March 18, 2002
Math 0B Fial Exam Solutios March 18, 00 1. (1 poits) (a) (6 poits) Fid the Gree s fuctio for the tilted half-plae {(x 1, x ) R : x 1 + x > 0}. For x (x 1, x ), y (y 1, y ), express your Gree s fuctio G(x,
More informationEE Midterm Test 1 - Solutions
EE35 - Midterm Test - Solutios Total Poits: 5+ 6 Bous Poits Time: hour. ( poits) Cosider the parallel itercoectio of the two causal systems, System ad System 2, show below. System x[] + y[] System 2 The
More informationFinite Difference Derivations for Spreadsheet Modeling John C. Walton Modified: November 15, 2007 jcw
Fiite Differece Derivatios for Spreadsheet Modelig Joh C. Walto Modified: November 15, 2007 jcw Figure 1. Suset with 11 swas o Little Platte Lake, Michiga. Page 1 Modificatio Date: November 15, 2007 Review
More informationWeek 2: Probability review Bernoulli, binomial, Poisson, and normal distributions Solutions
Wee 2: Probability review Beroulli, biomial, Poisso, ad ormal distributios Solutios A Biomial distributio. To evaluate the mea ad variace of a biomial RV B with parameters, p), we will rely o the relatio
More informationx x x Using a second Taylor polynomial with remainder, find the best constant C so that for x 0,
Math Activity 9( Due with Fial Eam) Usig first ad secod Taylor polyomials with remaider, show that for, 8 Usig a secod Taylor polyomial with remaider, fid the best costat C so that for, C 9 The th Derivative
More informationAssignment Number 3 Solutions
Math 4354, Assigmet Number 3 Solutios 1. u t (x, t) = u xx (x, t), < x (1) u(, t) =, u(, t) = u(x, ) = x ( 1) +1 u(x, t) = e t si(x). () =1 Solutio: Look for simple solutios i the form u(x, t) =
More informationThe Z-Transform. (t-t 0 ) Figure 1: Simplified graph of an impulse function. For an impulse, it can be shown that (1)
The Z-Trasform Sampled Data The geeralied fuctio (t) (also kow as the impulse fuctio) is useful i the defiitio ad aalysis of sampled-data sigals. Figure below shows a simplified graph of a impulse. (t-t
More informationMath 210A Homework 1
Math 0A Homework Edward Burkard Exercise. a) State the defiitio of a aalytic fuctio. b) What are the relatioships betwee aalytic fuctios ad the Cauchy-Riema equatios? Solutio. a) A fuctio f : G C is called
More informationInverse Matrix. A meaning that matrix B is an inverse of matrix A.
Iverse Matrix Two square matrices A ad B of dimesios are called iverses to oe aother if the followig holds, AB BA I (11) The otio is dual but we ofte write 1 B A meaig that matrix B is a iverse of matrix
More information: Transforms and Partial Differential Equations
Trasforms ad Partial Differetial Equatios 018 SUBJECT NAME : Trasforms ad Partial Differetial Equatios SUBJECT CODE : MA 6351 MATERIAL NAME : Part A questios REGULATION : R013 WEBSITE : wwwharigaeshcom
More informationCalculus 2 - D. Yuen Final Exam Review (Version 11/22/2017. Please report any possible typos.)
Calculus - D Yue Fial Eam Review (Versio //7 Please report ay possible typos) NOTE: The review otes are oly o topics ot covered o previous eams See previous review sheets for summary of previous topics
More informationG-2 Applied Computational Aerodynamics. G-1 Program THINFOIL
G-2 Applied Computatioal Aerodyamics G-1 Program THINFOIL THINFOIL solves Laplace s Equatio by fiite differeces usig a variety of iteratio methods. It was writte by Valery Razgoyaev. The iteratio optios
More information6.3. Poisson's Equation
Iterative Methods for Liear Systems 267 6.3. Poisso's Equatio 6.3.1. Poisso's Equatio i Oe Dimesio We begi with a oe-dimesioal versio of Poisso's equatio, where f(x) is a give fuctio ad v(x) is the ukow
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationSolutions to home assignments (sketches)
Matematiska Istitutioe Peter Kumli 26th May 2004 TMA401 Fuctioal Aalysis MAN670 Applied Fuctioal Aalysis 4th quarter 2003/2004 All documet cocerig the course ca be foud o the course home page: http://www.math.chalmers.se/math/grudutb/cth/tma401/
More informationThe Advection-Diffusion equation!
ttp://www.d.edu/~gtryggva/cf-course/! Te Advectio-iffusio equatio! Grétar Tryggvaso! Sprig 3! Navier-Stokes equatios! Summary! u t + u u x + v u y = P ρ x + µ u + u ρ y Hyperbolic part! u x + v y = Elliptic
More informationA note on the modified Hermitian and skew-hermitian splitting methods for non-hermitian positive definite linear systems
A ote o the modified Hermitia ad skew-hermitia splittig methods for o-hermitia positive defiite liear systems Shi-Liag Wu Tig-Zhu Huag School of Applied Mathematics Uiversity of Electroic Sciece ad Techology
More informationLecture 8: October 20, Applications of SVD: least squares approximation
Mathematical Toolkit Autum 2016 Lecturer: Madhur Tulsiai Lecture 8: October 20, 2016 1 Applicatios of SVD: least squares approximatio We discuss aother applicatio of sigular value decompositio (SVD) of
More information4.3 Growth Rates of Solutions to Recurrences
4.3. GROWTH RATES OF SOLUTIONS TO RECURRENCES 81 4.3 Growth Rates of Solutios to Recurreces 4.3.1 Divide ad Coquer Algorithms Oe of the most basic ad powerful algorithmic techiques is divide ad coquer.
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More information7.1 Finding Rational Solutions of Polynomial Equations
Name Class Date 7.1 Fidig Ratioal Solutios of Polyomial Equatios Essetial Questio: How do you fid the ratioal roots of a polyomial equatio? Resource Locker Explore Relatig Zeros ad Coefficiets of Polyomial
More informationProblem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =
Mathematics Summer Wilso Fial Exam August 8, ANSWERS Problem 1 (a) Fid the solutio to y +x y = e x x that satisfies y() = 5 : This is already i the form we used for a first order liear differetial equatio,
More informationMath 257: Finite difference methods
Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...
More informationA New Recursion for Space-Filling Geometric Fractals
A New Recursio for Space-Fillig Geometric Fractals Joh Shier Abstract. A recursive two-dimesioal geometric fractal costructio based upo area ad perimeter is described. For circles the radius of the ext
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationKinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More informationHOMEWORK I: PREREQUISITES FROM MATH 727
HOMEWORK I: PREREQUISITES FROM MATH 727 Questio. Let X, X 2,... be idepedet expoetial radom variables with mea µ. (a) Show that for Z +, we have EX µ!. (b) Show that almost surely, X + + X (c) Fid the
More information( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.
Joh Riley 8 Jue 03 ANSWERS TO EXERCISES IN APPENDIX B Sectio B VECTORS AND SETS Exercise B-: Covex sets (a) Let 0 x, x X, X, hece 0 x, x X ad 0 x, x X Sice X ad X are covex, x X ad x X The x X X, which
More informationPC5215 Numerical Recipes with Applications - Review Problems
PC55 Numerical Recipes with Applicatios - Review Problems Give the IEEE 754 sigle precisio bit patter (biary or he format) of the followig umbers: 0 0 05 00 0 00 Note that it has 8 bits for the epoet,
More informationChapter 7: The z-transform. Chih-Wei Liu
Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability
More informationA Lattice Green Function Introduction. Abstract
August 5, 25 A Lattice Gree Fuctio Itroductio Stefa Hollos Exstrom Laboratories LLC, 662 Nelso Park Dr, Logmot, Colorado 853, USA Abstract We preset a itroductio to lattice Gree fuctios. Electroic address:
More informationRandom Models. Tusheng Zhang. February 14, 2013
Radom Models Tusheg Zhag February 14, 013 1 Radom Walks Let me describe the model. Radom walks are used to describe the motio of a movig particle (object). Suppose that a particle (object) moves alog the
More informationAssignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.
PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator  is defied by its actio
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationSolutions to Final Exam Review Problems
. Let f(x) 4+x. Solutios to Fial Exam Review Problems Math 5C, Witer 2007 (a) Fid the Maclauri series for f(x), ad compute its radius of covergece. Solutio. f(x) 4( ( x/4)) ( x/4) ( ) 4 4 + x. Sice the
More informationDiscrete-Time Systems, LTI Systems, and Discrete-Time Convolution
EEL5: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we begi our mathematical treatmet of discrete-time s. As show i Figure, a discrete-time operates or trasforms some iput sequece x [
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationBrief Review of Functions of Several Variables
Brief Review of Fuctios of Several Variables Differetiatio Differetiatio Recall, a fuctio f : R R is differetiable at x R if ( ) ( ) lim f x f x 0 exists df ( x) Whe this limit exists we call it or f(
More informationMath 113, Calculus II Winter 2007 Final Exam Solutions
Math, Calculus II Witer 7 Fial Exam Solutios (5 poits) Use the limit defiitio of the defiite itegral ad the sum formulas to compute x x + dx The check your aswer usig the Evaluatio Theorem Solutio: I this
More informationQuestion 1: The magnetic case
September 6, 018 Corell Uiversity, Departmet of Physics PHYS 337, Advace E&M, HW # 4, due: 9/19/018, 11:15 AM Questio 1: The magetic case I class, we skipped over some details, so here you are asked to
More informationSolutions. Number of Problems: 4. None. Use only the prepared sheets for your solutions. Additional paper is available from the supervisors.
Quiz November 4th, 23 Sigals & Systems (5-575-) P. Reist & Prof. R. D Adrea Solutios Exam Duratio: 4 miutes Number of Problems: 4 Permitted aids: Noe. Use oly the prepared sheets for your solutios. Additioal
More informationStatistical Pattern Recognition
Statistical Patter Recogitio Classificatio: No-Parametric Modelig Hamid R. Rabiee Jafar Muhammadi Sprig 2014 http://ce.sharif.edu/courses/92-93/2/ce725-2/ Ageda Parametric Modelig No-Parametric Modelig
More informationA) is empty. B) is a finite set. C) can be a countably infinite set. D) can be an uncountable set.
M.A./M.Sc. (Mathematics) Etrace Examiatio 016-17 Max Time: hours Max Marks: 150 Istructios: There are 50 questios. Every questio has four choices of which exactly oe is correct. For correct aswer, 3 marks
More informationR is a scalar defined as follows:
Math 8. Notes o Dot Product, Cross Product, Plaes, Area, ad Volumes This lecture focuses primarily o the dot product ad its may applicatios, especially i the measuremet of agles ad scalar projectio ad
More informationMATH 304: MIDTERM EXAM SOLUTIONS
MATH 304: MIDTERM EXAM SOLUTIONS [The problems are each worth five poits, except for problem 8, which is worth 8 poits. Thus there are 43 possible poits.] 1. Use the Euclidea algorithm to fid the greatest
More informationLesson 10: Limits and Continuity
www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals
More informationIN many scientific and engineering applications, one often
INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS, VOL 3, NO, FEBRUARY 07 5 Secod Degree Refiemet Jacobi Iteratio Method for Solvig System of Liear Equatio Tesfaye Kebede Abstract Several
More informationMatrix Theory, Math6304 Lecture Notes from November 27, 2012 taken by Charles Mills
Matrix Theory, Math6304 Lecture Notes from November 27, 202 take by Charles Mills Last Time (9/20/2) Gelfad s formula for spectral radius Gershgori s circle theorem Warm-up: Let s observe what Gershgori
More informationMachine Learning for Data Science (CS 4786)
Machie Learig for Data Sciece CS 4786) Lecture & 3: Pricipal Compoet Aalysis The text i black outlies high level ideas. The text i blue provides simple mathematical details to derive or get to the algorithm
More informationSection 1 of Unit 03 (Pure Mathematics 3) Algebra
Sectio 1 of Uit 0 (Pure Mathematics ) Algebra Recommeded Prior Kowledge Studets should have studied the algebraic techiques i Pure Mathematics 1. Cotet This Sectio should be studied early i the course
More informationPhysics 7440, Solutions to Problem Set # 8
Physics 7440, Solutios to Problem Set # 8. Ashcroft & Mermi. For both parts of this problem, the costat offset of the eergy, ad also the locatio of the miimum at k 0, have o effect. Therefore we work with
More informationHW #1 Solutions: M552 Spring 2006
HW #1 Solutions: M55 Spring 006 1. (1.4-Trefethen & Bau) Let f 1,...,f 8 be a set of functions defined on the interval [1,8] with the property that for any numbers d 1,...,d 8, there exists a set of coefficients
More informationMath 312 Lecture Notes One Dimensional Maps
Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,
More informationSalmon: Lectures on partial differential equations. 3. First-order linear equations as the limiting case of second-order equations
3. First-order liear equatios as the limitig case of secod-order equatios We cosider the advectio-diffusio equatio (1) v = 2 o a bouded domai, with boudary coditios of prescribed. The coefficiets ( ) (2)
More informationL 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!
aylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. At ay poit i the eighbourhood of =0, the fuctio ca be represeted as a power series of the followig form: X 0 f(a) f() ƒ() f()= ( ) f( ) (
More informationModified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations
Noliear Aalysis ad Differetial Equatios, Vol. 5, 27, o. 4, 57-7 HIKARI Ltd, www.m-hikari.com https://doi.org/.2988/ade.27.62 Modified Decompositio Method by Adomia ad Rach for Solvig Noliear Volterra Itegro-
More information10.6 ALTERNATING SERIES
0.6 Alteratig Series Cotemporary Calculus 0.6 ALTERNATING SERIES I the last two sectios we cosidered tests for the covergece of series whose terms were all positive. I this sectio we examie series whose
More informationThe Discrete Fourier Transform
The Discrete Fourier Trasform Complex Fourier Series Represetatio Recall that a Fourier series has the form a 0 + a k cos(kt) + k=1 b k si(kt) This represetatio seems a bit awkward, sice it ivolves two
More informationAnalysis of Experimental Measurements
Aalysis of Experimetal Measuremets Thik carefully about the process of makig a measuremet. A measuremet is a compariso betwee some ukow physical quatity ad a stadard of that physical quatity. As a example,
More informationMath 778S Spectral Graph Theory Handout #3: Eigenvalues of Adjacency Matrix
Math 778S Spectral Graph Theory Hadout #3: Eigevalues of Adjacecy Matrix The Cartesia product (deoted by G H) of two simple graphs G ad H has the vertex-set V (G) V (H). For ay u, v V (G) ad x, y V (H),
More informationMTH 142 Exam 3 Spr 2011 Practice Problem Solutions 1
MTH 42 Exam 3 Spr 20 Practice Problem Solutios No calculators will be permitted at the exam. 3. A pig-pog ball is lauched straight up, rises to a height of 5 feet, the falls back to the lauch poit ad bouces
More informationMark Lundstrom Spring SOLUTIONS: ECE 305 Homework: Week 5. Mark Lundstrom Purdue University
Mark udstrom Sprig 2015 SOUTIONS: ECE 305 Homework: Week 5 Mark udstrom Purdue Uiversity The followig problems cocer the Miority Carrier Diffusio Equatio (MCDE) for electros: Δ t = D Δ + G For all the
More informationMath 128A: Homework 1 Solutions
Math 8A: Homework Solutios Due: Jue. Determie the limits of the followig sequeces as. a) a = +. lim a + = lim =. b) a = + ). c) a = si4 +6) +. lim a = lim = lim + ) [ + ) ] = [ e ] = e 6. Observe that
More informationFrequency Domain Filtering
Frequecy Domai Filterig Raga Rodrigo October 19, 2010 Outlie Cotets 1 Itroductio 1 2 Fourier Represetatio of Fiite-Duratio Sequeces: The Discrete Fourier Trasform 1 3 The 2-D Discrete Fourier Trasform
More informationSession 5. (1) Principal component analysis and Karhunen-Loève transformation
200 Autum semester Patter Iformatio Processig Topic 2 Image compressio by orthogoal trasformatio Sessio 5 () Pricipal compoet aalysis ad Karhue-Loève trasformatio Topic 2 of this course explais the image
More informationChapter 1 Simple Linear Regression (part 6: matrix version)
Chapter Simple Liear Regressio (part 6: matrix versio) Overview Simple liear regressio model: respose variable Y, a sigle idepedet variable X Y β 0 + β X + ε Multiple liear regressio model: respose Y,
More informationSequences and Series of Functions
Chapter 6 Sequeces ad Series of Fuctios 6.1. Covergece of a Sequece of Fuctios Poitwise Covergece. Defiitio 6.1. Let, for each N, fuctio f : A R be defied. If, for each x A, the sequece (f (x)) coverges
More information