Fourier Series Code. James K. Peterson. April 9, Department of Biological Sciences and Department of Mathematical Sciences Clemson University
|
|
- Darlene Nichols
- 5 years ago
- Views:
Transcription
1 Fourier Series Code James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April 9, 2018
2 Outline 1
3 We will need to approximate Fourier series expansions for arbitrary functions f when we solve the cable equation. All of these approximations require that we find inner products of the form < f, u n > for some function u n which is a sin or cos term. We must make sure our sin and cosine families are orthogonal so the Fourier coefficients can be accurate. If we have a familiy of functions which should be orthogonal, {g i } a standard test to see if our numerical inner product calculations are sufficiently accurate is to compute the N N matrix d = (< g i, g j >) which should be the essentially the N N identify matrix as the off diagonal entries should all be zero. However, if the inner product computations are inaccurate, off diagonal values need not be zero and our Fourier coefficients will be wrong.
4 1 f u n c t i o n c = i n n e r p r o d u c t ( f, g, a, b,n) h t ) f ( t ) g ( t ) ; 6 d e l x = ( b a ) /N; x = l i n s p a c e ( a, b,n+1) ; c = 0 ; f o r i =1:N c = c+ h ( x ( i ) ) d e l x ; 11
5 We can check how accurate we are with these inner product calculation by using them to do a Graham-Schmidt orthogonalization on the functions 1, t, t 2,..., t N on the interval [a, b]. The code to do this is in the function GrahamSchmidtTwo. This function returns the new orthogonal of length one functions g i and also prints out the matrix with terms < g i, g j > which should be an identity. The input NIP is the number of terms to use in the Riemann sum approximations to the inner product and N is size of the number of functions we perform GSO on.
6 f u n c t i o n g = GrahamSchmidtTwo ( a, b, N, NIP ) Perform Graham Schmidt O r t h o g o n a l i z a t i o n on a s e t o f f u n c t i o n s t ˆ 0,..., t ˆN 5 Setup f u n c t i o n h a n d l e s f = S e t U p F u n c t i o n s (N) ; g = c e l l (N+1,1) ; 10 n f = s q r t ( i n n e r p r o d u c t ( f {1}, f {1}, a, b, NIP ) ) ; g{1} x ) f {1}( x ) / n f ; d = z e r o s (N+1,N+1) ; f o r k =2:N+1 compute n e x t o r t h o g o n a l p i e c e 15 p h i x ) 0 ; f o r j = 1 : k 1 c = i n n e r p r o d u c t ( f {k }, g{ j }, a, b, NIP ) ; p h i x ) ( p h i ( x )+c g{ j }( x ) ) ; 20 p s i x ) ( f {k }( x ) p h i ( x ) ) ; n f = s q r t ( i n n e r p r o d u c t ( p s i, p s i, a, b, NIP ) ) ; g{k} x ) ( p s i ( x ) / n f ) ; 25 f o r i =1:N+1 f o r j =1:N+1 d ( i, j ) = i n n e r p r o d u c t ( g{ i }, g{ j }, a, b, NIP ) ; 30 d
7 Here is a typical run. g = GrahamSchmidtTwo ( 0, 2, 2, 5 0 ) ; This computes the GSO of the functions 1, t, t 2 on the interval [0, 2] using Riemann sum approximations with 50 points. The matrix < g i, g j > which is calculated is e e e e e e e e e+00 This is the 3 3 identify that we expect. Next, we do GSO on the functions 1, t,..., t 6. g = GrahamSchmidtTwo ( 0, 2, 6, 5 0 ) ;
8 This also generates the identity matrix we expect for < g i, g j > 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+00
9 To find the GSO of arbitary functions, use the GrahamScmidtThree function. f u n c t i o n g = GrahamSchmidtThree ( f, a, b, NIP ) Perform Graham Schmidt O r t h o g o n a l i z a t i o n on a s e t o f f u n c t i o n s f Setup f u n c t i o n h a n d l e s 5 N = l e n g t h ( f ) ; g = c e l l (N, 1 ) ; n f = s q r t ( i n n e r p r o d u c t ( f {1}, f {1}, a, b, NIP ) ) ; g{1} x ) f {1}( x ) / n f ; 10 d = z e r o s (N,N) ; f o r k =2:N compute n e x t o r t h o g o n a l p i e c e p h i x ) 0 ; f o r j = 1 : k 1 15 c = i n n e r p r o d u c t ( f {k }, g{ j }, a, b, NIP ) ; p h i x ) ( p h i ( x )+c g{ j }( x ) ) ; p s i x ) ( f {k }( x ) p h i ( x ) ) ; n f = s q r t ( i n n e r p r o d u c t ( p s i, p s i, a, b, NIP ) ) ; 20 g{k} x ) ( p s i ( x ) / n f ) ; check o r t h o g o n a l i t y f o r i =1:N f o r j =1:N 25 d ( i, j ) = i n n e r p r o d u c t ( g{ i }, g{ j }, a, b, NIP ) ; d
10 To use this, we have to setup a cell structure to hold our list of functions. For example, if f (t) = t 2, g(t) = sin(4t + 3) and h(t) = 1/(1 + t 2 ), we would do the GSO as follows: f t ) t. ˆ 2 ; g t ) s i n (4 t +3) ; h t ) 1./(1+ t. ˆ 2 ) ; F = c e l l ( 3, 1 ) ; F{1} = f ; F{2} = g ; F{3} = h ; G = GrahamSchmidtThree (F, 1,1,300) ; d =
11 We can then see the graphs of the original functions using the basic code T = l i n s p a c e ( 1,1,51) ; p l o t (T, F{1}(T),T, F{2}(T),T, F{3}(T) ) ; x l a b e l ( Time ) ; y l a b e l ( O r i g i n a l F u n c t i o n s ) ; t i t l e ( O r i g i n a l F u n c t i o n s on [ 1,1] ) ;
12 We can then plot the new orthonormal basis using the code T = l i n s p a c e ( 1,1,51) ; p l o t (T, G{1}(T),T, G{2}(T),T, G{3}(T) ) ; x l a b e l ( Time ) ; y l a b e l ( B a s i s F u n c t i o n s ) ; t i t l e ( B a s i s F u n c t i o n s on [ 1,1] ) ;
13 To compute the first N terms of the Fourier Cosine series of a function f on ( the interval [0, L], we first need a way to encode all the functions nπ cos L ). x We do this in the function SetUpCosines. f u n c t i o n f = S e t U p C o s i n e s ( L,N) Setup f u n c t i o n h a n d l e s 5 f = c e l l (N+1,1) ; f o r i =1:N+1 f { i } x ) c o s ( ( i 1) p i x /L ) ;
14 ( π This generates handles to the functions 1, cos L ), x cos ( forth ing with cos Nπ L x ). A similar function encodes the ( 2π L x ) and so corresponding sin functions we need for the first N terms of a Fourier Sine series. The function is called SetUpSines with code f u n c t i o n f = S e t U p S i n e s ( L,N) Setup f u n c t i o n h a n d l e s 5 f = c e l l (N, 1 ) ; f o r i =1:N f { i } x ) s i n ( i p i x /L ) ; ( π This generates handles to the functions sin L ), x sin ( Nπ forth ing with sin L ). x ( 2π L x ) and so
15 Let s check how accurate our innerproduct calculations are on the sin and cos terms. First, we modify the functions which set up the sin and cos functions to return functions of length one on the interval [0, L]. This is done in the functions SetUpOrthogCos and SetUpOrthogSin. f u n c t i o n f = SetUpOrthogCos ( L,N) Setup f u n c t i o n h a n d l e s 5 f = c e l l (N+1,1) ; f {1} x ) s q r t (1/ L ) ; f o r i =2:N+1 f { i } x ) s q r t (2/ L ) c o s ( ( i 1) p i x /L ) ; and 1 f u n c t i o n f = SetUpOrthogSin ( L,N) Setup f u n c t i o n h a n d l e s f = c e l l (N, 1 ) ; 6 f o r i =1:N f { i } x ) s q r t (2/ L ) s i n ( ( i ) p i x /L ) ;
16 Then, we can check the accuracy of the inner product calculations by computing the matrix d = (< f i, f j >). We do this with the new function CheckOrtho We input the function f and interval [a, b] points a and b and the number of terms to use in the Riemann sum approximation of the inner product, NIP. f u n c t i o n CheckOrtho ( f, a, b, NIP ) Perform Graham Schmidt O r t h o g o n a l i z a t i o n on a s e t o f f u n c t i o n s f 5 Setup f u n c t i o n h a n d l e s N = l e n g t h ( f ) ; f o r i =1:N 10 f o r j =1:N d ( i, j ) = i n n e r p r o d u c t ( f { i }, f { j }, a, b, NIP ) ; d 15 Let s try it with the cos functions. We use a small number of terms for the Riemann sums, NIP = 50 and compute < g i, g j > for the first 7 functions.
17 f = SetUpOrthogCos ( 5, 7 ) ; CheckOrtho ( f, 0, 5, 5 0 ) ; We do not get the identity matrix as we expect some of the off diagonal values are too large 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 e e e e e e e e e e e+00 Resetting NIP = 200, we find a better result.
18 1. 0 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 e e e e e e e e e e e+00 We can do even better by increasing NIP, of course. We encourage you to do these experiments yourself. We did not show the results for the sin terms, but they will be similar. We can then calculate the first N + 1 terms of a Fourier series or the first N terms of a Fourier Sine series using the functions FourierCosineApprox and FourierSineApprox, respectively.
19 Let s look at the Fourier Cosine approximation first. This function returns a handle to the function p which is the approximation p(x) = N n=0 < f, g i > g i (x) ( iπ where g i (x) = cos L ). x It also returns the Fourier cosine coefficients A 1 = 1 L < f, g 0 > through A N+1 = 2 L < f, g N >. We must choose how many points to use in our Riemann sum inner product estimates this is the input variable N. The other inputs are the function handle f and the length L. We also plot our approximation and the original function together so we can see how we did.
20 f u n c t i o n [ A, p ] = F o u r i e r C o s i n e A p p r o x ( f, L,M,N) p i s t h e Nth F o u r i e r A p p r o x i m a t i o n f i s t h e o r i g i n a l f u n c t i o n 5 M i s the number of Riemann sum terms i n the i n n e r product N i s the number of terms i n the approximation L i s t h e i n t e r v a l l e n g t h g e t t h e f i r s t N+1 F o u r i e r c o s a p p r o x i m a t i o n s 10 g = SetUpCosines ( L,N) ; g e t F o u r i e r C o s i n e C o e f f i c i e n t s A = z e r o s (N+1,1) ; A( 1 ) = i n n e r p r o d u c t ( f, g {1},0, L,M) /L ; f o r i =2:N+1 15 A( i ) = 2 i n n e r p r o d u c t ( f, g{ i },0, L,M) /L ; g e t Nth F o u r i e r C o s i n e A p p r o x i m a t i o n p x ) 0 ; f o r i =1:N+1 20 p x ) ( p ( x ) + A( i ) g{ i }( x ) ) ; x = l i n s p a c e ( 0, L, ) ; f o r i =1:101 y ( i ) = f ( x ( i ) ) ; 25 yp = p ( x ) ; f i g u r e s = [ Fourier Cosine Approximation with, i n t 2 s t r (N+1), term ( s ) ] ; p l o t ( x, y, x, yp ) ; 30 x l a b e l ( x axis ) ; y l a b e l ( y axis ) ; t i t l e ( s ) ;
21 We will test our approximations on two standard functions: the sawtooth curve and the square wave. To define these functions, we will use an auxiliary function splitfunc which defines a new function z on the interval [0, L] as follows { f (x) 0 x < L z(x) = 2, L g(x) 2 x L The arguments to splitfunc are the functions f and g, the value of x for we wish to find the output z and the value of L. f u n c t i o n z = s p l i t f u n c ( x, f, g, L ) i f x < L/2 z = f ( x ) ; e l s e z = g ( x ) ;
22 It is easy then to define a sawtooth and a square wave with the following code. The square wave, Sq has value H on [0, L 2 ) and value 0 on [ L 2, L]. In general, the sawtooth curve, Saw, is the straight line connecting the point (0, 0) to ( L 2, H) on the interval [0, L 2 ] and the line connecting ( L 2, H) to (L, 0) on the interval ( L 2, L). Thus, and Sq(x) = { H 0 x < L 0 L 2, 2 x L Saw(x) = { 2 L Hx 0 x < L 2, 2H 2 L Hx L 2 x L f u n c t i o n f = sawtooth ( L,H) f l e f t x ) 2 x H/L ; f r i g h t x ) 2 H 2 H x /L ; f x ) s p l i t f u n c ( x, f l e f t, f r i g h t, L ) ;
23 As an example, we build a square wave and sawtooth of height 10 on the interval [0, 10]. It is easy to plot these functions as well, although we won t show the plots yet. We will wait until we can compare the functions to their Fourier series approximations. However, the plotting code is listed below for your convenience. Note the plot must be set up in a for loop as the inequality checks in the function splitfunc do not handle a vector argument such as the x from a linspace command correctly. f x ) 1 0 ; g x ) 0 ; Saw = s awtooth ( 1 0, 1 0 ) ; Sq x ) s p l i t f u n c ( x, f, g, 1 0 ) ; x = l i n s p a c e ( 0, 1 0, ) ; f o r i =1:101 ysq ( i ) = Sq ( x ( i ) ) ; ysaw ( i ) = Saw ( x ( i ) ) ; p l o t ( x, ysq ) ; a x i s ( [. 1, ] ) ; p l o t ( x, YSaw) ;
24 We can then test the Fourier Cosine Approximation code on a saw function. f x ) 1 0 ; g x ) 0 ; Saw = s awtooth ( 1 0, 1 0 ) ; Sq x ) s p l i t f u n c ( x, f, g, 1 0 ) ; [ Acos, pcos ] = F o u r i e r C o s i n e A p p r o x ( Saw, 1 0, 1 0 0, 5 ) ;
25 We can then test the Fourier Cosine Approximation code on a square wave function. f x ) 1 0 ; g x ) 0 ; Saw = s awtooth ( 1 0, 1 0 ) ; Sq x ) s p l i t f u n c ( x, f, g, 1 0 ) ; [ Acos, pcos ] = F o u r i e r C o s i n e A p p r o x ( Sq, 1 0, 1 0 0, 5 ) ;
26 The code for FourierSineApprox is next. It is quite similar to the cosine approximation code and so we will say little about it. It returns the Fourier sine coefficients as A with ( A 1 = 2 L < f, g 1 > through A N = 2 L < f, g iπ N > where g i (x) = sin L ). x It also returns the handle to the Fourier sine approximation p(x) = N < f, g i ) > g i (x). n=1
27 f u n c t i o n [ A, p ] = F o u r i e r S i n e A p p r o x ( f, L,M,N) p i s t h e Nth F o u r i e r A p p r o x i m a t i o n f i s t h e o r i g i n a l f u n c t i o n M i s the number of Riemann sum terms i n the i n n e r product 5 N i s the number of terms i n the approximation L i s t h e i n t e r v a l l e n g t h g e t t h e f i r s t N F o u r i e r s i n e a p p r o x i m a t i o n s g = SetUpSines ( L,N) ; g e t F o u r i e r S i n e C o e f f i c i e n t s 10 A = z e r o s (N, 1 ) ; f o r i =1:N A( i ) = 2 i n n e r p r o d u c t ( f, g{ i },0, L,M) /L ; g e t Nth F o u r i e r S i n e A p p r o x i m a t i o n 15 p x ) 0 ; f o r i =1:N p x ) ( p ( x ) + A( i ) g{ i }( x ) ) ; x = l i n s p a c e ( 0, L, ) ; 20 f o r i =1:101 y ( i ) = f ( x ( i ) ) ; yp = p ( x ) ; f i g u r e 25 s = [ Fourier Sine Approximation with, i n t 2 s t r (N), term ( s ) ] ; p l o t ( x, y, x, yp ) ; x l a b e l ( x axis ) ; y l a b e l ( y axis ) ; t i t l e ( s ) ; 30
28 Let s test the approximation code on the square wave. using 22 terms. The sin approximations are not going to like the starting value at 10 as you can see in the figure below. The approximation is generated with the command f x ) 1 0 ; g x ) 0 ; Saw = s awtooth ( 1 0, 1 0 ) ; Sq x ) s p l i t f u n c ( x, f, g, 1 0 ) ; [ Asin, p s i n ] = F o u r i e r S i n e A p p r o x ( Sq, 1 0, 1 0 0, 2 2 ) ;
29 Here is another generic type of function to work with: a pulse. This is a function constant of an interval of the form [x 0 r, x 0 + r] and 0 off of that interval. We implement it in pulsefunc. f u n c t i o n z = p u l s e f u n c ( x, x0, r,h) i f x > x0 r && x < x0+r z = H; 5 e l s e z = 0 ; It is easy to use. Here is a pulse centered at 3 of radius 0.5 on the interval [0, 10]. P x ) p u l s e f u n c ( x, 3, 0. 5, 1 0 ) ; [ Asin, p s i n ] = F o u r i e r S i n e A p p r o x (P, 1 0, 1 0 0, 2 2 ) ;
30
31 We can combine the Fourier Sin and Fourier Cos Approximations by adding them and dividing by two. Here is the code with the plotting portion removed to make it shorter. f u n c t i o n [ A, B] = F o u r i e r A p p r o x ( f, L,M,N) 3 g = SetUpCosines ( L,N) ; B = z e r o s (N+1,1) ; B( 1 ) = i n n e r p r o d u c t ( f, g {1},0, L,M) /L ; f o r i =2:N+1 B( i ) = 2 i n n e r p r o d u c t ( f, g{ i },0, L,M) /L ; 8 p x ) 0 ; f o r i =1:N+1 p x ) ( p ( x ) + B( i ) g{ i }( x ) ) ; 13 h = SetUpSines ( L,N) ; A = z e r o s (N, 1 ) ; f o r i =1:N A( i ) = 2 i n n e r p r o d u c t ( f, h{ i },0, L,M) /L ; 18 q x ) 0 ; f o r i =1:N q x ) ( q ( x ) + A( i ) h{ i }( x ) ) ; 23 x = l i n s p a c e ( 0, L, ) ; f o r i =1:101 y ( i ) = f ( x ( i ) ) ; yp = 0. 5 ( p ( x )+q ( x ) ) ; 28
32 P x ) p u l s e f u n c ( x, 3, 0. 5, 1 0 ) ; [ A, B] = F o u r i e r A p p r o x (P, 1 0, 1 0 0, 2 2 ) ;
33 Homework Generate the Fourier Sine Approximation with varying number of terms and different values of NIP for the periodic square wave. Draw figures as needed Generate the Fourier Sine Approximation with varying number of terms and different values of NIP for pulses applied at various locations. Draw figures as needed Generate the Fourier Cosine Approximation with varying number of terms and different values of NIP for pulses applied at various locations. Draw figures as needed.
More Least Squares Convergence and ODEs
More east Squares Convergence and ODEs James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April 12, 219 Outline Fourier Sine and Cosine Series
More informationGeometric Series and the Ratio and Root Test
Geometric Series and the Ratio and Root Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 5, 2018 Outline 1 Geometric Series
More informationCable Convergence. James K. Peterson. May 7, Department of Biological Sciences and Department of Mathematical Sciences Clemson University
Cable Convergence James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University May 7, 2018 Outline 1 Fourier Series Convergence Redux 2 Fourier Series
More informationRiemann Integration. James K. Peterson. February 2, Department of Biological Sciences and Department of Mathematical Sciences Clemson University
Riemann Integration James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University February 2, 2017 Outline 1 Riemann Sums 2 Riemann Sums In MatLab 3 Graphing
More informationRiemann Integration. Outline. James K. Peterson. February 2, Riemann Sums. Riemann Sums In MatLab. Graphing Riemann Sums
Riemann Integration James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University February 2, 2017 Outline Riemann Sums Riemann Sums In MatLab Graphing
More informationFourier Sin and Cos Series and Least Squares Convergence
Fourier and east Squares Convergence James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University May 7, 28 Outline et s look at the original Fourier sin
More informationGeometric Series and the Ratio and Root Test
Geometric Series and the Ratio and Root Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 5, 2017 Outline Geometric Series The
More informationMore Series Convergence
More Series Convergence James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University December 4, 218 Outline Convergence Analysis for Fourier Series Revisited
More informationThe First Derivative and Second Derivative Test
The First Derivative and Second Derivative Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April 9, 2018 Outline 1 Extremal Values 2
More informationConvergence of Fourier Series
MATH 454: Analysis Two James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April, 8 MATH 454: Analysis Two Outline The Cos Family MATH 454: Analysis
More informationwhich arises when we compute the orthogonal projection of a vector y in a subspace with an orthogonal basis. Hence assume that P y = A ij = x j, x i
MODULE 6 Topics: Gram-Schmidt orthogonalization process We begin by observing that if the vectors {x j } N are mutually orthogonal in an inner product space V then they are necessarily linearly independent.
More informationThe First Derivative and Second Derivative Test
The First Derivative and Second Derivative Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 2017 Outline Extremal Values The
More informationLecture 5b: Starting Matlab
Lecture 5b: Starting Matlab James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University August 7, 2013 Outline 1 Resources 2 Starting Matlab 3 Homework
More informationRiemann Sums. Outline. James K. Peterson. September 15, Riemann Sums. Riemann Sums In MatLab
Riemann Sums James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 15, 2013 Outline Riemann Sums Riemann Sums In MatLab Abstract This
More informationPHYS 301 HOMEWORK #5
PHY 3 HOMEWORK #5 olutions On this homework assignment, you may use Mathematica to compute integrals, but you must submit your Mathematica output with your assignment.. For f x = -, - < x
More informationConvergence of Sequences
James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 5, 2018 Outline 1 2 Homework Definition Let (a n ) n k be a sequence of real numbers.
More informationSolving systems of ODEs with Matlab
Solving systems of ODEs with Matlab James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 20, 2013 Outline 1 Systems of ODEs 2 Setting Up
More informationAP Calculus Testbank (Chapter 9) (Mr. Surowski)
AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series
More informationFourier Sin and Cos Series and Least Squares Convergence
Fourier Sin and Cos Series and east Squares Convergence James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University December 4, 208 Outline Sin and Cos
More informationIntegration by Parts Logarithms and More Riemann Sums!
Integration by Parts Logarithms and More Riemann Sums! James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 16, 2013 Outline 1 IbyP with
More informationTaylor Polynomials. James K. Peterson. Department of Biological Sciences and Department of Mathematical Sciences Clemson University
James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 24, 2013 Outline 1 First Order Approximation s Second Order Approximations 2 Approximation
More informationVariation of Parameters
Variation of Parameters James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April 13, 218 Outline Variation of Parameters Example One We eventually
More informationConvergence of Sequences
Convergence of Sequences James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University February 12, 2018 Outline Convergence of Sequences Definition Let
More informationTaylor and Maclaurin Series
Taylor and Maclaurin Series MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Background We have seen that some power series converge. When they do, we can think of them as
More informationSolving Linear Systems of ODEs with Matlab
Solving Linear Systems of ODEs with Matlab James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 27, 2013 Outline Linear Systems Numerically
More informationChapter 17 : Fourier Series Page 1 of 12
Chapter 7 : Fourier Series Page of SECTION C Further Fourier Series By the end of this section you will be able to obtain the Fourier series for more complicated functions visualize graphs of Fourier series
More informationExtreme Values and Positive/ Negative Definite Matrix Conditions
Extreme Values and Positive/ Negative Definite Matrix Conditions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 016 Outline 1
More informationMath 24 Spring 2012 Sample Homework Solutions Week 8
Math 4 Spring Sample Homework Solutions Week 8 Section 5. (.) Test A M (R) for diagonalizability, and if possible find an invertible matrix Q and a diagonal matrix D such that Q AQ = D. ( ) 4 (c) A =.
More informationVector Functions & Space Curves MATH 2110Q
Vector Functions & Space Curves Vector Functions & Space Curves Vector Functions Definition A vector function or vector-valued function is a function that takes real numbers as inputs and gives vectors
More informationMathematical Induction Again
Mathematical Induction Again James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University January 12, 2017 Outline Mathematical Induction Simple POMI Examples
More informationMath 180 Written Homework Assignment #10 Due Tuesday, December 2nd at the beginning of your discussion class.
Math 18 Written Homework Assignment #1 Due Tuesday, December 2nd at the beginning of your discussion class. Directions. You are welcome to work on the following problems with other MATH 18 students, but
More informationAntiderivatives! Outline. James K. Peterson. January 28, Antiderivatives. Simple Fractional Power Antiderivatives
Antiderivatives! James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University January 28, 2014 Outline Antiderivatives Simple Fractional Power Antiderivatives
More informationMathematical Induction Again
Mathematical Induction Again James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University January 2, 207 Outline Mathematical Induction 2 Simple POMI Examples
More information1 angle.mcd Inner product and angle between two vectors. Note that a function is a special case of a vector. Instructor: Nam Sun Wang. x.
angle.mcd Inner product and angle between two vectors. Note that a function is a special case of a vector. Instructor: Nam Sun Wang Define angle between two vectors & y:. y. y. cos( ) (, y). y. y Projection
More informationDifferentiating Series of Functions
Differentiating Series of Functions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 30, 017 Outline 1 Differentiating Series Differentiating
More information3. Use absolute value notation to write an inequality that represents the statement: x is within 3 units of 2 on the real line.
PreCalculus Review Review Questions 1 The following transformations are applied in the given order) to the graph of y = x I Vertical Stretch by a factor of II Horizontal shift to the right by units III
More informationAntiderivatives! James K. Peterson. January 28, Department of Biological Sciences and Department of Mathematical Sciences Clemson University
! James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University January 28, 2014 Outline 1 2 Simple Fractional Power Abstract This lecture is going to talk
More informationGeneral Power Series
General Power Series James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University March 29, 2018 Outline Power Series Consequences With all these preliminaries
More informationThere are two things that are particularly nice about the first basis
Orthogonality and the Gram-Schmidt Process In Chapter 4, we spent a great deal of time studying the problem of finding a basis for a vector space We know that a basis for a vector space can potentially
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 203 Outline
More informationHölder s and Minkowski s Inequality
Hölder s and Minkowski s Inequality James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 1, 218 Outline Conjugate Exponents Hölder s
More informationINTRODUCTION TO REAL ANALYSIS II MATH 4332 BLECHER NOTES
INTRODUCTION TO REAL ANALYSIS II MATH 433 BLECHER NOTES. As in earlier classnotes. As in earlier classnotes (Fourier series) 3. Fourier series (continued) (NOTE: UNDERGRADS IN THE CLASS ARE NOT RESPONSIBLE
More informationMath Computer Lab 4 : Fourier Series
Math 227 - Computer Lab 4 : Fourier Series Dylan Zwick Fall 212 This lab should be a pretty quick lab. It s goal is to introduce you to one of the coolest ideas in mathematics, the Fourier series, and
More informationFIRST YEAR CALCULUS W W L CHEN
FIRST YER CLCULUS W W L CHEN c W W L Chen, 994, 28. This chapter is available free to all individuals, on the understanding that it is not to be used for financial gain, and may be downloaded and/or photocopied,
More informationUNIVERSITY OF SOUTHAMPTON. A foreign language dictionary (paper version) is permitted provided it contains no notes, additions or annotations.
UNIVERSITY OF SOUTHAMPTON MATH055W SEMESTER EXAMINATION 03/4 MATHEMATICS FOR ELECTRONIC & ELECTRICAL ENGINEERING Duration: 0 min Solutions Only University approved calculators may be used. A foreign language
More informationTaylor and Maclaurin Series. Approximating functions using Polynomials.
Taylor and Maclaurin Series Approximating functions using Polynomials. Approximating f x = e x near x = 0 In order to approximate the function f x = e x near x = 0, we can use the tangent line (The Linear
More informationProject One: C Bump functions
Project One: C Bump functions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 2, 2018 Outline 1 2 The Project Let s recall what the
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 2013 Outline
More informationMath 291-2: Final Exam Solutions Northwestern University, Winter 2016
Math 29-2: Final Exam Solutions Northwestern University, Winter 206 Determine whether each of the following statements is true or false f it is true, explain why; if it is false, give a counterexample
More informationLecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain.
Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain. For example f(x) = 1 1 x = 1 + x + x2 + x 3 + = ln(1 + x) = x x2 2
More informationConstrained Optimization in Two Variables
Constrained Optimization in Two Variables James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 17, 216 Outline Constrained Optimization
More informationDirchlet s Function and Limit and Continuity Arguments
Dirchlet s Function and Limit and Continuity Arguments James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University February 23, 2018 Outline 1 Dirichlet
More informationMore on Fourier Series
More on Fourier Series R. C. Trinity University Partial Differential Equations Lecture 6.1 New Fourier series from old Recall: Given a function f (x, we can dilate/translate its graph via multiplication/addition,
More informationEngg. Math. I. Unit-I. Differential Calculus
Dr. Satish Shukla 1 of 50 Engg. Math. I Unit-I Differential Calculus Syllabus: Limits of functions, continuous functions, uniform continuity, monotone and inverse functions. Differentiable functions, Rolle
More informationMath Real Analysis II
Math 4 - Real Analysis II Solutions to Homework due May Recall that a function f is called even if f( x) = f(x) and called odd if f( x) = f(x) for all x. We saw that these classes of functions had a particularly
More informationChapter 4 Sequences and Series
Chapter 4 Sequences and Series 4.1 Sequence Review Sequence: a set of elements (numbers or letters or a combination of both). The elements of the set all follow the same rule (logical progression). The
More information11.10a Taylor and Maclaurin Series
11.10a 1 11.10a Taylor and Maclaurin Series Let y = f(x) be a differentiable function at x = a. In first semester calculus we saw that (1) f(x) f(a)+f (a)(x a), for all x near a The right-hand side of
More information3 rd class Mech. Eng. Dept. hamdiahmed.weebly.com Fourier Series
Definition 1 Fourier Series A function f is said to be piecewise continuous on [a, b] if there exists finitely many points a = x 1 < x 2
More informationUniform Convergence Examples
Uniform Convergence Examples James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 13, 2017 Outline 1 Example Let (x n ) be the sequence
More information1 A complete Fourier series solution
Math 128 Notes 13 In this last set of notes I will try to tie up some loose ends. 1 A complete Fourier series solution First here is an example of the full solution of a pde by Fourier series. Consider
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationAlgebraic. techniques1
techniques Algebraic An electrician, a bank worker, a plumber and so on all have tools of their trade. Without these tools, and a good working knowledge of how to use them, it would be impossible for them
More informationMath 118B Solutions. Charles Martin. March 6, d i (x i, y i ) + d i (y i, z i ) = d(x, y) + d(y, z). i=1
Math 8B Solutions Charles Martin March 6, Homework Problems. Let (X i, d i ), i n, be finitely many metric spaces. Construct a metric on the product space X = X X n. Proof. Denote points in X as x = (x,
More information2016 Mathematics. Advanced Higher. Finalised Marking Instructions
National Qualifications 06 06 Mathematics Advanced Higher Finalised ing Instructions Scottish Qualifications Authority 06 The information in this publication may be reproduced to support SQA qualifications
More informationThursday 14 June 2012 Morning
Thursday 14 June 01 Morning A GCE MATHEMATICS (MEI) 4757 Further Applications of Advanced Mathematics (FP3) QUESTION PAPER *471575061* Candidates answer on the Printed Answer Book. OCR supplied materials:
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationMore On Exponential Functions, Inverse Functions and Derivative Consequences
More On Exponential Functions, Inverse Functions and Derivative Consequences James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University January 10, 2019
More informationConstrained Optimization in Two Variables
in Two Variables James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 17, 216 Outline 1 2 What Does the Lagrange Multiplier Mean? Let
More informationRemark 3.2. The cross product only makes sense in R 3.
3. Cross product Definition 3.1. Let v and w be two vectors in R 3. The cross product of v and w, denoted v w, is the vector defined as follows: the length of v w is the area of the parallelogram with
More informationMA5206 Homework 4. Group 4. April 26, ϕ 1 = 1, ϕ n (x) = 1 n 2 ϕ 1(n 2 x). = 1 and h n C 0. For any ξ ( 1 n, 2 n 2 ), n 3, h n (t) ξ t dt
MA526 Homework 4 Group 4 April 26, 26 Qn 6.2 Show that H is not bounded as a map: L L. Deduce from this that H is not bounded as a map L L. Let {ϕ n } be an approximation of the identity s.t. ϕ C, sptϕ
More informationHölder s and Minkowski s Inequality
Hölder s and Minkowski s Inequality James K. Peterson Deartment of Biological Sciences and Deartment of Mathematical Sciences Clemson University Setember 10, 2018 Outline 1 Conjugate Exonents 2 Hölder
More informationNote : This document might take a little longer time to print. more exam papers at : more exam papers at : more exam papers at : more exam papers at : more exam papers at : more exam papers at : more
More informationFourier and Partial Differential Equations
Chapter 5 Fourier and Partial Differential Equations 5.1 Fourier MATH 294 SPRING 1982 FINAL # 5 5.1.1 Consider the function 2x, 0 x 1. a) Sketch the odd extension of this function on 1 x 1. b) Expand the
More informationUnit 10 Parametric and Polar Equations - Classwork
Unit 10 Parametric and Polar Equations - Classwork Until now, we have been representing graphs by single equations involving variables x and y. We will now study problems with which 3 variables are used
More informationConsequences of Continuity
Consequences of Continuity James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 4, 2017 Outline Domains of Continuous Functions The Intermediate
More informationFourier Series for Functions Defined on Arbitrary Limited Intervals with Polynomial Expansion
American Review of Mathematics and Statistics December 2016, Vol. 4, No.2, pp. 10-17 ISSN: 2374-2348 (Print), 2374-2356 (Online) Copyright The Author(s).All Rights Reserved. Published by American Research
More informationBC Exam 1 - Part I 28 questions No Calculator Allowed - Solutions C = 2. Which of the following must be true?
BC Exam 1 - Part I 8 questions No Calculator Allowed - Solutions 6x 5 8x 3 1. Find lim x 0 9x 3 6x 5 A. 3 B. 8 9 C. 4 3 D. 8 3 E. nonexistent ( ) f ( 4) f x. Let f be a function such that lim x 4 x 4 I.
More informationMath 201 Assignment #11
Math 21 Assignment #11 Problem 1 (1.5 2) Find a formal solution to the given initial-boundary value problem. = 2 u x, < x < π, t > 2 u(, t) = u(π, t) =, t > u(x, ) = x 2, < x < π Problem 2 (1.5 5) Find
More informationWeighted SS(E) = w 2 1( y 1 y1) y n. w n. Wy W y 2 = Wy WAx 2. WAx = Wy. (WA) T WAx = (WA) T b
6.8 - Applications of Inner Product Spaces Weighted Least-Squares Sometimes you want to get a least-squares solution to a problem where some of the data points are less reliable than others. In this case,
More informationElementary maths for GMT
Elementary maths for GMT Linear Algebra Part 1: Vectors, Representations Algebra and Linear Algebra Algebra: numbers and operations on numbers 2 + 3 = 5 3 7 = 21 Linear Algebra: tuples, triples... of numbers
More informationnatural frequency of the spring/mass system is ω = k/m, and dividing the equation through by m gives
77 6. More on Fourier series 6.. Harmonic response. One of the main uses of Fourier series is to express periodic system responses to general periodic signals. For example, if we drive an undamped spring
More informationCHAPTER ONE FUNCTIONS AND GRAPHS. In everyday life, many quantities depend on one or more changing variables eg:
CHAPTER ONE FUNCTIONS AND GRAPHS 1.0 Introduction to Functions In everyday life, many quantities depend on one or more changing variables eg: (a) plant growth depends on sunlight and rainfall (b) speed
More informationFourier Analysis Fourier Series C H A P T E R 1 1
C H A P T E R Fourier Analysis 474 This chapter on Fourier analysis covers three broad areas: Fourier series in Secs...4, more general orthonormal series called Sturm iouville epansions in Secs..5 and.6
More informationTic, Toc: Pendulum Motion
Tic, Toc: Pendulum Motion Activity 25 Pendulum motion has long fascinated people. Galileo studied pendulum motion by watching a swinging chandelier and timing it with his pulse. In 1851 Jean Foucault demonstrated
More informationThis Week. Professor Christopher Hoffman Math 124
This Week Sections 2.1-2.3,2.5,2.6 First homework due Tuesday night at 11:30 p.m. Average and instantaneous velocity worksheet Tuesday available at http://www.math.washington.edu/ m124/ (under week 2)
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 3 2, 5 2 C) - 5 2
Test Review (chap 0) Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Solve the problem. ) Find the point on the curve x = sin t, y = cos t, -
More informationOverview of Fourier Series (Sect. 6.2). Origins of the Fourier Series.
Overview of Fourier Series (Sect. 6.2. Origins of the Fourier Series. Periodic functions. Orthogonality of Sines and Cosines. Main result on Fourier Series. Origins of the Fourier Series. Summary: Daniel
More informationMidterm 1 Review. Distance = (x 1 x 0 ) 2 + (y 1 y 0 ) 2.
Midterm 1 Review Comments about the midterm The midterm will consist of five questions and will test on material from the first seven lectures the material given below. No calculus either single variable
More informationIn this chapter you will learn how to use MATLAB to work with lengths, angles and projections in subspaces of R and later in certain linear spaces.
Chapter 5. Introduction In this chapter you will learn how to use MATLAB to work with lengths angles and n projections in subspaces of R and later in certain linear spaces.. Lengths and Angles Recall from
More informationMath 1272 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Exam ) This fraction appears in Problem 5 of the undated-? exam; a solution can be found in that solution set. (E) ) This integral appears in Problem 6 of the Fall 4 exam;
More informationMath 116 Final Exam April 21, 2016
Math 6 Final Exam April 2, 206 UMID: Instructor: Initials: Section:. Do not open this exam until you are told to do so. 2. Do not write your name anywhere on this exam. 3. This exam has 4 pages including
More informationMatrices and Vectors
Matrices and Vectors James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 11, 2013 Outline 1 Matrices and Vectors 2 Vector Details 3 Matrix
More informationVector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms.
Vector Spaces Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. For each two vectors a, b ν there exists a summation procedure: a +
More informationDO NOT OPEN THIS BOOKLET UNTIL YOU ARE TOLD TO DO SO.
AP Calculus AB Exam SECTION I: Multiple Choice 016 DO NOT OPEN THIS BOOKLET UNTIL YOU ARE TOLD TO DO SO. At a Glance Total Time 1 hour, 45 minutes Number of Questions 45 Percent of Total Score 50% Writing
More informationUniform Convergence Examples
Uniform Convergence Examples James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 13, 2017 Outline More Uniform Convergence Examples Example
More information2009 GCE A Level Solution Paper 1
2009 GCE A Level Solution Paper i) Let u n = an 2 + bn + c. u = a + b + c = 0 u 2 = 4a + 2b + c = 6 u 3 = 9a + 3b + c = 5 Using GC, a =.5, b = 8.5, c = 7. u n =.5n 2 8.5n + 7. (ii) Let y =.5n 2 8.5n +
More informationHW3 - Due 02/06. Each answer must be mathematically justified. Don t forget your name. 1 2, A = 2 2
HW3 - Due 02/06 Each answer must be mathematically justified Don t forget your name Problem 1 Find a 2 2 matrix B such that B 3 = A, where A = 2 2 If A was diagonal, it would be easy: we would just take
More informationConsequences of Continuity
Consequences of Continuity James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 4, 2017 Outline 1 Domains of Continuous Functions 2 The
More informationUniversity of Leeds, School of Mathematics MATH 3181 Inner product and metric spaces, Solutions 1
University of Leeds, School of Mathematics MATH 38 Inner product and metric spaces, Solutions. (i) No. If x = (, ), then x,x =, but x is not the zero vector. So the positive definiteness property fails
More information10.2-3: Fourier Series.
10.2-3: Fourier Series. 10.2-3: Fourier Series. O. Costin: Fourier Series, 10.2-3 1 Fourier series are very useful in representing periodic functions. Examples of periodic functions. A function is periodic
More information