comparación doble de matlab

Question

Estoy tratando de comparar una matriz de dobles a un doble escalar para la igualdad, pero la igualdad nunca se reconoce bajo ciertas circunstancias. Sospecho que esto tiene que ver con la forma en que se representa el doble (por ejemplo, 1.0 vs 1.00), pero no puedo resolverlo.

Por ejemplo, he generado una matriz compuesta de miles de valores dobles, el último pocos de los cuales en algún instante de tiempo se dan por

10.6000 
-11.0000 
10.2000 
22.6000 
3.4000 

Si comprobar la igualdad de 10.2 (o 10.2000) por el comando array==10.2 (o array=10.2000), devuelvo una matriz de 0. Si coloco los valores que se muestran en una matriz manualmente (por ejemplo, array=[10.6000 -11.0000 10.2000 22.6000 3.4000]), el comando es exitoso (es decir, array==10.2 devuelve 0 0 1 0 0). ¿Podría alguien explicar por qué la igualdad tiene éxito si ingreso los valores manualmente, pero falla si la matriz se genera en el contexto de un programa? Puedo rectificar el error de comparación utilizando una comparación aproximada en lugar de una exacta (por ejemplo, (array<10.20001) & (array>10.19999)), pero esto parece insatisfactorio.

Editar: Los valores en la matriz se generan mediante la suma o resta iterativa de un doble constante (por ejemplo, 0.2). El módulo de este conjunto por 0.2 debe estar por lo tanto en todas partes igual a 0. De hecho, el módulo de cada elemento es igual a cualquiera de 0 o 0.2, como se muestra a continuación para la secuencia anterior de los números en la matriz:

mod(array,0.2) 
... 
0.2000 
    0 
0.2000 
0.2000 
    0 

De nuevo, si los valores se colocan en una matriz de forma manual y el módulo se toma, se obtiene el valor esperado de todos 0 s.

Answer 1

La razón es que MATLAB truncó los números en la matriz para conservar solo 4 dígitos después del punto decimal cuando se visualiza. Es decir, el valor real de su matriz puede ser [10.600000000001, -10.99999999999, ...]. Tiene razón, esto se debe a la representación interna de los números de coma flotante en la computadora, lo que puede ocasionar pequeños errores en los cálculos.

Personalmente, creo que hay dos soluciones, una es aproximada como usted, mientras que la otra es redondear la matriz primero (por ejemplo, con this tool de FileExchange), y luego hacer la coincidencia exacta.

Answer 2

Algo probablemente esté en una sola precisión en alguna parte, y doble en otros lugares. La representación binaria de, p. 10.2 en cada uno es diferente ya que terminan después de un número diferente de bits. Por lo tanto son diferentes:

>> if (single(10.2) == 10.2) disp('honk'); end 
>> if (single(10.2) == single(10.2)) disp('honk'); end 
honk 

Tendrá que comprobar la igualdad dentro de alguna pequeña diferencia:

eps = 0.001; 
result = abs(array-10.2) < eps; 

Puede encontrar la precisión utilizado en una matriz mediante whos:

>> whos A 
    Name  Size   Bytes Class  Attributes 

    A   1x2     8 single  

Answer 3

Cree un archivo de función MATLAB que acepte los valores de un módulo (de 3 a 9; es decir, de Z3 a Z9) y generará el menor valor posible descrito por las condiciones del módulo.

Simulación de la muestra:

Z = [3 4 5]; % Modulo Z3, Z4, Z5 y

r = [2 1 4]; % El resto Valores

menos valor posible es 29.

Las entradas Z debe ser una matriz array ... donde en puede escribir cualquier número de 3 a 9 .... y usted puede escriba 3,4,5,6,7,8,9 en cualquier orden, en cualquier emparejamientos o agrupaciones ...

las entradas R debe ser igual al número de entradas z también ...

la salida debería producir el menor valor posible aunque las condiciones del módulo ...