Let be nontrivial involutions, i.e., R=R−1≠±In and S=S−1≠±In. A matrix A∈Cn×n is referred to as (R,S)-symmetric if and only if RAS=A. The set of all (R,S)-symmetric matrices of order n is denoted by . Given a full column rank matrix X∈Cn×m, a matrix B∈Cm×m and a matrix A∗∈Cn×n. In structural dynamic model updating, we usually consider the sets and in the Frobenius norm sense, where the superscript H denotes conjugate transpose. Then we characterize the unique matrices and when R=RH and S=SH. By using the generalized singular value decomposition (GSVD), the necessary and sufficient conditions for the non-emptiness of S1 and the general representations of the elements in S1 and are derived, respectively. The analytical expressions of A∈S2 and are also obtained by using the GSVD, the canonical correlation decomposition (CCD) and the projection theorem. Finally, a corresponding numerical algorithm and some illustrated examples are presented.