En trabajos recientes, se han propuesto esquemas de diferencias finitas exactas y no estándar para ecuaciones diferenciales lineales de retardo de primer orden escalares. El objetivo del presente trabajo es extender estos resultados anteriores a sistemas de ecuaciones diferenciales de retardo acopladas donde es un vector, y y son matrices reales conmutativas, en general no simultáneamente diagonalizables. Basándose en una expresión constructiva para la solución exacta de la ecuación vectorial, se obtiene un esquema exacto, y se proponen diferentes esquemas numéricos no estándar de orden creciente. Las propiedades de consistencia dinámica de los nuevos esquemas no estándar se ilustran con ejemplos numéricos, y se demuestran para una clase de métodos.