Identify R2 with C and write A∈M2(R) as Ax=zx+wxˉ with z=2(a+d)+i(c−b), w=2(a−d)+i(b+c) for A=(acbd). Recall that A↦(z,w) is linear, that ∥A∥F2=2(∣z∣2+∣w∣2), and that the rotation Rα has coordinates (eiα,0). Show that R∈SO(2)min∥A−R∥F2=2((∣z∣−1)2+∣w∣2), and find the nearest rotation(s).
Hint
By linearity, A−Rα has coordinates (z−eiα,w), so ∥A−Rα∥F2=2(∣z−eiα∣2+∣w∣2). Minimise ∣z−eiα∣ over α.
Show solutionHide solution
Solution
Idea. In the coordinates (z,w) the Frobenius norm is Euclidean (up to the factor 2), so the nearest rotation to A is the nearest point of the circle {(eiα,0)} to the point (z,w).
Step 1: the squared distance.A−Rα has coordinates (z−eiα,w) since A↦(z,w) is linear. Hence ∥A−Rα∥F2=2(∣z−eiα∣2+∣w∣2).
Step 2: minimise over α. For any α, the reverse triangle inequality gives ∣z−eiα∣≥∣z∣−∣eiα∣=∣z∣−1. Equality holds iff eiα and z are positively proportional, i.e. eiα=z/∣z∣ when z=0. When z=0 we have ∣0−eiα∣=1=∣0−1∣ for all α, so the bound is attained for every α. Therefore αmin∥A−Rα∥F2=2((∣z∣−1)2+∣w∣2).
Nearest rotation. If z=0 the minimiser is unique: R=Rα with eiα=z/∣z∣, i.e. the rotation by argz (in matrices: R=∣z∣1(RezImz−ImzRez), the normalised "rotation-scaling part" 2(a+d)I+2(c−b)(01−10) of A). If z=0 every rotation is equally near (this happens iff A is "anti-conformal", of the form x↦wxˉ).
Check.A=2I: z=2, w=0; the formula gives 2(2−1)2=2, and indeed ∥2I−I∥F2=2. A=diag(1,−1): z=0, w=1; the formula gives 2(1+1)=4, and ∥diag(1,−1)−Rα∥F2=(1−cosα)2+sin2α+sin2α+(−1−cosα)2=4 for every α. A numerical minimisation over a fine grid of α on random matrices reproduces the formula.
Remark. For the other component of O(2), {z=0,∣w∣=1}, the same argument gives min∥A−F∥F2=2(∣z∣2+(∣w∣−1)2). Comparing the two minima, [(∣z∣−1)2+∣w∣2]−[∣z∣2+(∣w∣−1)2]=2(∣w∣−∣z∣), so a matrix is closer to SO(2) than to O(2)∖SO(2) iff ∣z∣>∣w∣, i.e. iff detA>0. This is the 2×2 case of the orthogonal Procrustes problem: the nearest orthogonal matrix is obtained from the SVD A=UΣVT as UVT.
Lesson. When a distance to a group is wanted, pick coordinates in which the group is a coordinate subset (here the circle in the z-plane) and the norm is Euclidean.