User:Egm6322.s12.team2/HW12
Contents
R 12.1 - Integration Points and Gauss-Legendre Quadrature[edit]
Given[edit]
Legendre series
Find[edit]
Compute the 1st non-zero coefficient in the Legendre series with increasing number of integration points until convergence to within 1% accuracy. Use the error formula (1) p. 45b-7 to estimate the necessary number of integration points to achieve the desired accuracy.
Solution[edit]
Withn the class context, non-zero coefficients have generally been solved by applying the Fourier-Legendre series and orthogonality to find non-zero coefficients. For team 2, the previous solutions of 9.1 and 10.7 will be expanded to solve this problem. The general pattern is described below.
Step 1: Define 1st Non-Zero Coefficient
The coefficient of a Fourier-Legendre series requires the definition below.
The coefficient is from the term for the Fourier-Legendre series of f(x).
Essentially, the series utilizes the Legendre polynomials and coefficients to define a function . The Legendre polynomials have the following orthogonality.
The problem does not explicitly give a function, f(x), for finding the 1st non-zero coefficient. Therefore, the team's approach is to utilize the function from 9.1 and 10.7 below.
The coefficients are then defined as follows.
From 9.1 and 10.7, the 1st non-zero coefficient of the chosen f(x) is
Step 2: Apply Qauss-Quadrature to Integral of Coefficient
The weighting function is
The error formula is shown below from 45-7b.
Step 3: Find Number of Integration Points
The following definition is used to find the number of integration points.
This can be rearranged to apply the error limit of 1% or 0.01.
Note that is the (2n)th order derivative of . Normally, this higher order derivative would require a rigorous analysis with numerical methods and sophisticated programming. But that can be minimized by noting that a derivative of and vice versa. So the following definition is utilized.
So the simplified error term is below.
The Matlab image and code are shown below.
Matlab code Matlab code http://en.wikiversity.org/wiki/File:PEA2_12_1_Matlab.PNG
http://en.wikiversity.org/wiki/File:PEA2_12_1_Error_Term.png
See images below. The data from the Excel shows that there are two integration points needed to get the error less than 0.01 or 1%. The error term derived was simply entered into MS exceo with and calculated for n = 1 to 20. At n =2 E < 0.01, so n = 2 will suffice.
http://en.wikiversity.org/wiki/File:Team2_12_1_Excel_Data.PNG
Egm6322.s12.team2.steele.m2 (talk) 20:32, 28 March 2012 (UTC)
R 12.2 - Alternative formula for weights of GL quad[edit]
Given[edit]
Using the following fomula
(12.2.1)
Find[edit]
Show that
(12.2.2)
and
(12.2.3)
Solution[edit]
Solved on my own |
From the 12.2.1, we have
(12.2.4)
Since ,(The proof is @ p.45-17) then
(12.2.5)
In terms of (5)p.45b-6
(12.2.6)
then substitute 12.2.5 into 12.2.6, we obtain
(12.2.3)
R 12.3 - Proof of Thean of Smynna[edit]
Given[edit]
(12.3.1)
(12.3.2)
(12.3.3)
Find[edit]
Plot the sequence versus the iteration number n
Solution[edit]
By using Eq(12.3.1) and Eq(12.3.2),
Finding
Plot the Graph
File:Ssskimmm1.png
Egm6322.s12.sungsik (talk) 21:08, 27 March 2012 (UTC)
R 12.4 - Computing the convergents of [edit]
Given[edit]
Continued fraction formula for :
(12.4.1)
First four convergents:
(12.4.2)
Find[edit]
Write a Matlab program to compute the convergents for i = 1,2,...10 and plot these convergents along with the line to visualize the convergence
Solution[edit]
Solved on my own |
This problem mainly deals with Matlab implementation. From Eq. 12.4.1 it is clear that the quotient of the terms in the denominator is always 2 due to the relation . It is seen in the Mtg 48e that the term on the right can be expressed as ; substituting in the first relation then produces the quotient of 2.
When implementing this code in Matlab, it is more convenient to 'build up' the residue term starting with the subtraction. As all contributions to the evaluation after are ignored, this value is simply 1/2. It also defines the final level of the denominator which must be computed. Thus, a for loop is constructed between [:-1:2] to iteratively develop the residue term. This is seen in the Matlab code in lines 11 - 14. The computation itself is quite simple; each successive iteration simply adds the quotient of 2 onto the previously evaluated derivative.
The convergents are computed in Line 19 by adding the first quotient of 1 onto each residue term. The final values are then plotted against the exact value (to machine precision) of and are seen to converge within 5% of the exact value as soon as the second convergent. It can also be seen that the first 4 terms match those given by Eq. 12.4.2.
Matt Shields 11:40,24 March 2012 (UTC)
R 12.5 Legendre polynomials and continued fraction[edit]
Given[edit]
Find[edit]
(i) Verify the last term for the nth approximant as shown above.
(ii) The Legendre polynomials and compare the results to those in R10.2 on p. 44b-6.
(iii) What can you observe about the factor in when generated using the contnued fraction compared to the factor of as obtained in R10.2?
(iv) Would the difference in these factors change the orthogonality of the polynomials generated by the continued fraction (1) p. 48-24b?
Solution[edit]
Part (i)
The solution requires an analysis of the continued fraction given for this problem.
There is a general pattern of continued fractions.
....
In this case
Part (ii)
By applying the formula from part(i), we get
Therefore, the 1st Legendre polynomial is
Next, we get the 2nd Legendre polynomial.
Therefore, the 2nd Legendre polynomial is
The 3rd Legendre polynomial is just found by the 3rd convergent.
The inverse to get the Legendre polynomial is . From this pattern, it is clear that for each Legendre polynomial the related factor is missing when compared to 10.6. To confirm this pattern, let us rely on these two sources.
http://en.wikipedia.org/wiki/Legendre_polynomials
http://digitalcommons.uconn.edu/cgi/viewcontent.cgi?article=1079&context=chem_educ
C.W. David of University of Connecticut above demonstrates that the Legendre polynomials are derived from continued fractions in such a way that the constant of the solution is not necessary since the derived polynomial fits the general Legendre differential equation as a solution.
Part (iii)
In part (ii) we found that the fractional term at the front of each Legendre component is absent in the solutions found by continued fractions. For example, the fraction (1/2) is missing from the 2nd and 3rd Legendre polynomials for continued fractions. The method of derivation based on continued fractions produce solutions based on the derivates of the Legendre differential equations and recursive insertion of solutions. The ratio of derivatives results in an integration to give the Legendre function.
Part (iv)
The orthogonal property is not dependent on the factor term. The orthogonality is based on this term below, so the solutions from continuing fractions with convergents are still orthonal. This solution is not dependent on the fractional term derived from Rodrigues formula. The absence of the fractional term does not affect orthongality, which is still valid for continued fractional solutions.
Egm6322.s12.team2.steele.m2 (talk) 20:30, 28 March 2012 (UTC)
Contributing Members[edit]
Problem Number | Assigned To | Solved By | Typed By | Proofread By | |
12.1 | Manuel Steele | Manuel Steele | Manuel Steele | Lang Xia | |
12.2 | Lang Xia | Lang Xia | Lang Xia | Manuel Steele | |
12.3 | sungsik kim | sungsik kim | sungsik kim | Matt Shields | |
12.4 | Matt Shields | Matt Shields | Matt Shields | Lang Xia | |
12.5 | Manuel Steele | Manuel Steele | Manuel Steele | sungsik kim |
References for Report 12[edit]
King, A.C., J. Billingham and S.R. Otto. "Differential Equations: Linear, Nonlinear, Ordinary, Partial." New York, NY: Cambridge University Press, 2003.
http://en.wikipedia.org/wiki/Legendre_polynomials
http://digitalcommons.uconn.edu/cgi/viewcontent.cgi?article=1079&context=chem_educ
Vu-Quoc, L. Class Lecture: Principles of Engineering Analysis. University of Florida, Gainesville, FL, (Meeting 45h) Mtg 45h Spring 2012.
Vu-Quoc, L. Class Lecture: Principles of Engineering Analysis. University of Florida, Gainesville, FL, (Meeting 48a) Mtg 48a Spring 2012.
Vu-Quoc, L. Class Lecture: Principles of Engineering Analysis. University of Florida, Gainesville, FL, (Meeting 48b) Mtg 48b Spring 2012.
Vu-Quoc, L. Class Lecture: Principles of Engineering Analysis. University of Florida, Gainesville, FL, (Meeting 48c) Mtg 48c Spring 2012.
Vu-Quoc, L. Class Lecture: Principles of Engineering Analysis. University of Florida, Gainesville, FL, (Meeting 48d) Mtg 48d Spring 2012.
Vu-Quoc, L. Class Lecture: Principles of Engineering Analysis. University of Florida, Gainesville, FL, (Meeting 48e) Mtg 48e Spring 2012.
http://mathworld.wolfram.com/Fourier-LegendreSeries.html