Accuracy Condition Numbers and Pivoting:
In this lecture we will discuss two separate issues of accuracy in solving linear systems. The firstpivoting is methods that ensure that Gaussian elimination proceeds as accurately as possible.
The second condition number is determining of how bad a matrix is. We will observe that if a matrixhas a bad condition number the solutions are unbalanced with respect to small changes in data.
The effect of rounding:
All computers store numbers as limited strings of binary floating point digits. This limits numbers toa preset number of significant digits and implies that after even the most basic calculations, roundingmust happen.
Consider the following exaggerated illustration- presume that our computer can only store 2 significantdigits and it is asked to do Gaussian elimination on:
Doing the exclusion exactly would produce:
However rounding to 2 digits our computer would store this as:
Backsolving this abridged system gives:
x1 = 0, x2 = 3.
This seems fine until you understand that backsolving the unrounded system gives
x1 = −1, x2 = 3.001.
Row Pivoting:
A way to fix the problem is to utilize pivoting which signifies to switch rows of the matrix. Asswitching rows of the augmented matrix now corresponds to switching the order of the equationsno harm is done:
Exact removal would produce:
Storing this result with only 2 important digits gives:
Now backsolving produces:
x1= −1, x2 = 3,
Which is precisely right.
The reason this worked is for the reason that 1 is bigger than 0.001. To pivot us switch rows therefore that the largestentry in a column is the one used to eliminate the others. In bigger matrices subsequent to each column iscompleted compare the diagonal element of the next column with all the entries below it. Switchit (as well as the entire row) with the one with greatest absolute value. For illustration in the followingmatrix the first column is finished as well as before doing the second column pivoting must occur since:
|− 2| >|1|:
Pivoting the 2nd as well as 3rd rows would produce:
Condition number:
In some systems problems take place even without rounding. Consider the subsequent augmented matrices:
Here we have the same A however two different input vectors:
b1 = (3/2, 1)′and b2 = (3/2, 5/6)′
which are appealing close to one another. You would expect afterwards that the solutions x1 and x2 wouldalso be close. Notice that this matrix doesn’t need pivoting. Eliminating precisely we get:
Now solving we discover:
x1 = (0, 3)′ and x2 = (1, 1)′
which aren’t close at all despite the fact that we did the calculations exactly. This pretence a newproblem some matrices are extremely sensitive to small changes in input data. The extent of thissensitivity is calculated by the condition number. The description of condition number is- considerall small changes δA as well as δb in A and b and the resulting change δx in the solution x. Then,
Put one more way changes in the input data get multiplied by the condition number to producechanges in the outputs. Therefore a high condition number is bad. It entails that small errors in theinput can cause large errors in the output.
In Mat lab enter:
> H = hilb(2)
which must result in the matrix above. Mat lab produces the condition number of a matrix withthe command:
>cond(H)
Therefore for this matrix small errors in the input can get magnified by 19 in the output. After that try thematrix:
> A = [ 1.2969 0.8648 ; .2161 .1441]>cond(A)
For this matrix little errors in the input can get magnified by 2.5×108 in the output! (We will observethis happen in the exercise.) This is obviously not extremely good for engineering where all measurements,constants as well as inputs are approximate.
Is there a resolution to the problem of bad condition numbers? Typically bad condition numbers in engineering contexts result from poor design. Therefore the engineering solution to bad conditioning is redesign.
Ultimately find the determinant of the matrix A above
>det(A)
Which will be little. If det(A) = 0 then the matrix is remarkable which is bad for the reason that it implies therewill not be a unique solution. The case here det(A) ≈ 0 is as well bad because it means the matrix isalmost singular. Even though det(A) ≈ 0 in general indicates that the condition number will be largethey are actually independent things. To observe this find the determinant and condition number ofthe matrix [1e-10,0;0,1e-10] and the matrix [1e+10,0;0,1e-10].
Latest technology based Matlab Programming Online Tutoring Assistance
Tutors, at the www.tutorsglobe.com, take pledge to provide full satisfaction and assurance in Matlab Programming help via online tutoring. Students are getting 100% satisfaction by online tutors across the globe. Here you can get homework help for Matlab Programming, project ideas and tutorials. We provide email based Matlab Programming help. You can join us to ask queries 24x7 with live, experienced and qualified online tutors specialized in Matlab Programming. Through Online Tutoring, you would be able to complete your homework or assignments at your home. Tutors at the TutorsGlobe are committed to provide the best quality online tutoring assistance for Matlab Programming Homework help and assignment help services. They use their experience, as they have solved thousands of the Matlab Programming assignments, which may help you to solve your complex issues of Matlab Programming. TutorsGlobe assure for the best quality compliance to your homework. Compromise with quality is not in our dictionary. If we feel that we are not able to provide the homework help as per the deadline or given instruction by the student, we refund the money of the student without any delay.
tutorsglobe.com mullikens scale assignment help-homework help by online electronegativity scales tutors
tutorsglobe.com linux assignment help-homework help by online computer programming tutors
Principles of segregation tutorial all along with the key concepts of Mendel's First Law, Punnett Squares, Locus, Homozygous or Heterozygous, Phenotypic Ratio, Genotypic Ratio and Monohybrid Cross
www.tutorsglobe.com offers psychology homework help, psychology assignment help, online psychology tutoring assistance, psychology solutions by online qualified and experienced psychology tutor's help.
Body organization of Astacus tutorial all along with the key concepts of Locomotion in astacus, digestive system and nutrition in astacus, Gaseous exchange in astacus, Excretory system in astacus
tutorsglobe.com portal of entry assignment help-homework help by online non-toxic determinants of virulence tutors
heat change tutorial all along with the key concepts of heat change under constant volume, enthalpy and enthalpy changes, relation between cp and cv of an ideal gas
General Characteristics of Microorganisms tutorial all along with the key concepts of functions of Microorganisms, Bacteria, Fungi, Viruses, Algae and Protozoa
Enzymes-Vitamins and Co-enzymes tutorial all along with the key concepts of Definition of Enzymes, Nomenclature of Enzymes, Categorization of Enzymes, General Properties of Enzymes, How Enzymes Work, Co-enzymes and Vitamins
The Energetic of Particle tutorial all along with the key concepts of a-Decay, Binding energy, disintegration energy of nuclei, Range of a -Particle, a-decay paradox, ß Decay, Electron capture or k-capture, ß-spectrum, Neutrino, ? - Decay
The advantages of the uniform costing high standards, cost comparision, healthy competition, automation of data, benefits to similar units, fair prices, decreased labour turnover.
One of the unique factors among the management accounting and financial accounting is that the management accounting does not comprise a unified structure.
tutorsglobe.com criticism of marginal theory assignment help-homework help by online marginal productivity theory of distribution tutors
One way in which long-term performance can be rewarded is through the granting of directors’ share options. This sort of reward though has provoked considerable controversy.
tutorsglobe.com pentose phosphate pathway assignment help-homework help by online plant physiology tutors
1958154
Questions Asked
3689
Tutors
1474818
Questions Answered
Start Excelling in your courses, Ask an Expert and get answers for your homework and assignments!!