... (-3 * j) que es simplemente -3 en un coseno, y (-1 * j) o simplemente -1 en un coseno
En base a esto, creo que tienes un poco de confusión sobre los números complejos.
La relación que quieres tener en cuenta es
\ $ exp (j \ theta) = \ cos (\ theta) + j \ sin (\ theta) \ $,
por lo que ninguno de tus ejemplos da un número real o imaginario puro.
Otra relación importante es cierta para cualquier exponencial:
\ $ x ^ y \ cdot x ^ z = x ^ {(y + z)} \ $.
Para responder a la pregunta específica, ya que k nunca aparece dentro del bucle for, es bastante fácil desenvolverlo todo:
zz = 10 * exp(-3*j) * exp(j*7833*pi*tt) * 2 * exp(-1*j) * exp(j*7833*pi*tt) * 2 * exp(-1*j) * exp(j*7833*pi*tt) * 2 * exp(-1*j) * exp(j*7833*pi*tt)
Ahora puedes reunir términos semejantes:
zz = 10 * 2 * 2 * 2 * exp(-3*j) * exp(-1*j) * exp(-1*j) * exp(-1*j) * exp(j*7833*pi*tt) * exp(j*7833*pi*tt) * exp(j*7833*pi*tt) * exp(j*7833*pi*tt)
Los factores reales no afectan la ganancia, así que ignora esos. Entonces tienes
zz = mag * exp(j * (-3 + -1 + -1 + -1 + 7833*pi*tt + 7833*pi*tt + 7833*pi*tt + 7833*pi*tt))
o
zz = mag * exp(j * (-6 + 4*7833*pi*tt) )
Ahora depende de lo que quieras decir por fase, ya que parece que estás representando una señal que varía en el tiempo. Si te refieres a la fase relativa en comparación con alguna otra señal con la misma frecuencia, lo que importa es el término -6 . Tu fase es -6 radianes. Si te refieres a la fase instantánea en cada muestra de tiempo, entonces es -6 + 31,332 \ $ \ pi \ $ radianes por unidad de tiempo.