Normal forms and joint numerical ranges of doubly commuting matrices Original Research Article

A normal form is derived for finite sets of doubly commuting matrices, under simultaneous unitary similarity. The matrices need not be normal, but they commute with each other and with the adjoints of each other. The normal form is further used to study joint numerical ranges of doubly commuting matrices.