import sympy as s

d=10
q=3

sigma = s.symbols("sigma")
W = s.MatrixSymbol("W",d,q)
Iq = s.Identity(q)
Id = s.Identity(d)

mat =Iq - W.T @ (sigma**2 * Id + W@W.T).I @ W