Logarithmic stability for a coefficient inverse problem of coupled Schrödinger equations