Approximate solution of one equation using a quadratic or cubic regression

 

Print!

Objectives To show how certain statistical knowledge can intervene in analysis
Public
Classes Final or BTS
Prérequis Knowledge of the method of least squares
Idea B.Egger

 

Unfolding:

That is to say to solve the equation x-cos(x)=0.

quadra1.gif (977 octets)

A fast study shows that this equation shows that this equation does not admit that a solution ranging between 0 and 1.

 

From some points of the curve representative of the function f(x)=x-cos(x), we will create polynomials of degree 2 or 3 passer by by the cloud formed by these points and will replace the preceding equation by polynomial equations.

Commencons by 4 points.

One re-enters the column c1 in the form of a continuation and the column c2 like image of c1 by F:

c1=seq......

c2=f(c1)

quadra2.gif (1334 octets)

 

Determination of a quadratic regression for the scatter plot considered.

One rage the result in y1.

quadra3.gif (1182 octets)    quadra4.gif (1168 octets)

quadra5.gif (1183 octets)

 

Even thing with a cubic regression. One arranges the result in y2.

quadra6.gif (2370 octets)

 

The " real " solution of the equation and approximate solutions (obtained much more quickly).

Only one disadvantage: other solutions parasitic, but easy to eliminate.

quadra7.gif (1415 octets)

Let us increase the number of points of the cloud, and remake calculations.

quadra8.gif (2115 octets)       quadra9.gif (2099 octets)

**time-out** one have improve the precision in increase the number of point, but it be not always the case... as it be possible also to find some situation for which from same cloud, the regression quadratic be good than the regression cubic.

We know that it is normal from the point of view of the statistics, but that can not be obvious for the pupils.

Let us examine the case of the function F defined by f(x) = ln(x) + 10x

 

quadra10.gif (2096 octets)

quadra11.gif (1503 octets)

quadra12.gif (1932 octets)