clean
nettoie une matrice (arrondi à zéro des termes très petits)
Séquence d'appel
B = clean(A) B = clean(A, epsa) B = clean(A, epsa, epsr)
Arguments
- A
scalaire, vecteur, matrice ou hypermatrice de nombres décimaux ou de nombres complexes (en encodage dense ou creux), ou de polynômes ou fractions rationnelles à coefficients réels ou complexes.
- epsa,epsr
nombres réels positifs : tolérances. Par défaut = 10-10 pour les deux.
Description
Cette fonction met à 0 les termes de A
de valeur absolue < epsa, ou de valeur relative
< epsr (par rapport à max(abs(A)), %inf et %nan exclus)
Pour des polynômes ou fractions rationnelles, clean() est appliquée de la même façon
à l'ensemble des coefficients de tous les éléments de A.
Exemples
--> m = grand(3,5,"uin",1,9) .* 10.^grand(3,5,"uin",-16,1) m = 5.000D-08 5. 0.8 7.000D-08 4.000D-11 0.00007 7. 4.000D-09 6.000D-13 2.000D-16 2.000D-08 1.000D-14 1. 8. 0.7 --> clean(m) ans = 5.000D-08 5. 0.8 7.000D-08 0. 0.00007 7. 4.000D-09 0. 0. 2.000D-08 0. 1. 8. 0.7 --> clean(m, 1e-7) ans = 0. 5. 0.8 0. 0. 0.00007 7. 0. 0. 0. 0. 0. 1. 8. 0.7
Pour des fractions rationnelles :
--> w = [3+2*x, 2+x; 3+x, 2-x]/5
w =
0.6 +0.4x 0.4 +0.2x
0.6 +0.2x 0.4 -0.2x
--> r = w * inv(w)
r =
2
-3.140D-16 + x -1.110D-15 + 3.331D-16x
--------------- ------------------------
2
x 1.3333333x + x
2 2
1.110D-15 - 2.220D-16x - 5.551D-17x -1.110D-15 + 1.3333333x + x
------------------------------------ ----------------------------
2 2
1.3333333x + x 1.3333333x + x
--> clean(r)
ans =
1 0
-- --
1 1
0 1
-- --
1 1
Voir aussi
- %eps — epsilon (floating-point relative accuracy)
- cond — conditionnement
- nearfloat — get previous or next floating-point number
- number_properties — determine floating-point parameters
| Report an issue | ||
| << ceil | Virgule flottante | double >> |