Super Multidimensional Scaling (SMDS)

• Super Multidimensional Scaling (SMDS) is a technique that extends classical MDS by integrating both distance and angle measurements to form a unified low-rank edge kernel.
• It features variants across complex, real, and quaternion domains, each tailored to optimize noise suppression and computational efficiency in localization tasks.
• SMDS is applied in wireless and rigid body localization, demonstrating robust performance even with high measurement noise and missing data.

Searching arXiv for papers on Super Multidimensional Scaling and closely related variants. I’ll look up the specific arXiv records and nearby papers to ground the article in the current literature. Super Multidimensional Scaling (SMDS) denotes an extension of classical multidimensional scaling in which pairwise distances are supplemented or reformulated to encode richer geometric structure. In the localization literature, SMDS augments distance-based reconstruction with angle measurements and organizes the resulting information through an edge-kernel construction, so that node coordinates can be recovered from a low-rank matrix representation. Recent work has pushed this line of development toward quaternion-domain formulations for three-dimensional localization and toward kernel-minor formulations for rigid body localization (Masuoka et al., 24 Apr 2025). The acronym is not fully uniform across the literature, however: in a distinct line of work on language-model interpretability, “SMDS” denotes “Supervised Multi-Dimensional Scaling,” a different method with different objectives (Tiblias et al., 1 Oct 2025).

1. Localization-oriented SMDS: core formulation

Within wireless localization and related inference problems, SMDS extends classical multidimensional scaling by incorporating angle measurements in addition to distances. In the formulation summarized by later rigid-body localization work, MDS reconstructs node positions from distance measurements, whereas SMDS fuses distance and angle information through complex-valued edge variables, with positions represented as complex numbers in $2$D (Führling et al., 1 Sep 2025).

For two nodes with location vectors in $\mathbb{R}^2$, the complex edge is written as

$$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$

where $d_p$ is the Euclidean distance and $\theta_p$ is the angle. Collecting all such edges into a vector $\mathbf{v}$ yields the kernel

$\mathcal{K} = \mathbf{v}^* \mathbf{v}^\intercal,$

whose entries encode the joint structure of distances and angles. In the ideal case this kernel is rank-$1$, so the edge vector can be recovered from the leading eigenpair through

$\hat{\mathbf{v}} = \sqrt{\lambda}\,\mathbf{u}.$

Coordinates are then obtained by solving a pseudo-inverse problem involving a coefficient matrix that specifies how edges depend on points (Führling et al., 1 Sep 2025).

This formulation makes the central methodological move of SMDS explicit: rather than treating angle information as an auxiliary correction to distance-only MDS, it embeds distances and directional information in a single algebraic object. This suggests that, in localization, SMDS is best understood as a low-rank edge-kernel framework rather than as a single fixed algorithm.

2. Quaternion-domain SMDS for three-dimensional localization

A major recent development is quaternion-domain super multidimensional scaling (QD-SMDS), introduced for low-complexity $3$D localization in wireless sensor networks. Conventional $\mathbb{R}^2$0D SMDS, when formulated in the real domain, represents node positions as vectors $\mathbb{R}^2$1, constructs a Gram edge kernel (GEK) from hybrid distance-angle measurements, and yields a GEK of rank $\mathbb{R}^2$2 in $\mathbb{R}^2$3D. Because low-rank truncation can suppress only the components outside that intrinsic rank, noise reduction is limited in the real-domain formulation (Masuoka et al., 24 Apr 2025).

QD-SMDS reformulates this construction in the quaternion domain. A $\mathbb{R}^2$4D position $\mathbb{R}^2$5 is mapped to the pure quaternion

$\mathbb{R}^2$6

Edge vectors are mapped similarly to quaternion edge variables $\mathbb{R}^2$7. Their products encode both inner-product and projected angular information: $\mathbb{R}^2$8 The additional terms capture angular or phase relations on the three orthogonal coordinate planes. The quaternion-domain GEK is then

$\mathbb{R}^2$9

which is rank-$$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$0 by construction (Masuoka et al., 24 Apr 2025).

The importance of this reformulation lies in the combination of compactness and denoising. By encoding distances together with projected angular relations in a single rank-$$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$1 quaternion Hermitian matrix, QD-SMDS maximizes the noise reduction effect of low-rank truncation. The same summary contrasts this with the real-domain $$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$2D GEK, whose rank is $$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$3, and with complex-domain SMDS, which achieves rank $$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$4 in $$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$5D but does not natively represent full $$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$6D geometry (Masuoka et al., 24 Apr 2025).

3. Recovery procedures, QSVD, and SVD-free variants

In QD-SMDS, localization proceeds by constructing the noisy quaternion GEK, applying quaternion singular value decomposition (QSVD), extracting the principal component, and recovering real coordinates from the estimated quaternion edge vector. The principal estimate is

$$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$7

where $$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$8 is the largest singular value and singular vector of the noisy GEK. The recovered edge vector is then mapped back to real-valued coordinates, using known anchor locations, a structure matrix, the Moore–Penrose pseudoinverse, and, when necessary, Procrustes alignment (Lukaj et al., 23 Jul 2025).

A later development introduces an SVD-free variant, QD-MRC-SMDS. Its derivation partitions the quaternion edge vector into anchor-anchor and anchor-target blocks, exploits the block structure of the quaternion GEK, and applies a linear filtering step described as maximal ratio combining (MRC). The resulting method avoids eigendecomposition and directly estimates node coordinates using only matrix multiplications within the quaternion domain. An iterative MRC refinement based on power iteration is also described, and is reported to converge typically in one step because of the underlying rank-$$v_p = (a_{x_j} - a_{x_i}) + j(b_{x_j} - b_{x_i}) = d_p(\cos\theta_p + j\sin\theta_p),$$9 structure; elsewhere in the same summary the iterative variant is described as needing typically $d_p$0–$d_p$1 iterations (Lukaj et al., 23 Jul 2025).

The localization literature therefore now contains several closely related SMDS variants.

Variant	Domain	Defining property
Conventional SMDS	Complex in $d_p$2D; real-domain formulation in $d_p$3D	GEK fuses distance and angle information
QD-SMDS	Quaternion domain	Rank-$d_p$4 GEK for $d_p$5D localization
QD-MRC-SMDS	Quaternion domain	SVD-free recovery via MRC and matrix multiplications
SMDS-based RBL variant	Complex-domain kernel minor	Uses anchor-anchor and target-target information to estimate anchor-target edges

A plausible implication is that the central design axis across these variants is not the choice of embedding space alone, but the extent to which the edge kernel can be made both information-rich and structurally low rank.

4. Performance regimes, robustness, and computational trade-offs

The empirical behavior of quaternion-domain SMDS is explicitly regime dependent. Two simulation settings recur in the summaries. In Scenario I, only mutual distances and angle differences of arrival are available. In Scenario II, azimuth and elevation angles are additionally measurable, for example via a planar antenna array (Masuoka et al., 24 Apr 2025).

The reported results distinguish sharply between low-noise and high-noise conditions. For small angular errors, such as $d_p$6 or $d_p$7, conventional SMDS can outperform QD-SMDS when distance error is large; the explanation given is that, when angular measurements are already accurate, the additional phase or angle information brought by QD-SMDS offers little added value (Masuoka et al., 24 Apr 2025). By contrast, for larger angular errors, reported as $d_p$8, QD-SMDS consistently outperforms SMDS, and the performance gap widens as angular error increases. When additional phase information is available in Scenario II, the improvement of QD-SMDS over SMDS becomes even more pronounced (Lukaj et al., 23 Jul 2025).

The later study also reports robustness under missing data. Because the quaternion GEK is rank-$d_p$9, matrix completion is more powerful than in the conventional rank-$\theta_p$0 SMDS setting; even with $\theta_p$1 missing GEK entries, QD-SMDS outperformed SMDS for moderate and large angle errors. The same study reports that QD-MRC-SMDS attains localization accuracy very close to full QD-SMDS, including under large measurement errors, but at much lower computational cost (Lukaj et al., 23 Jul 2025).

On computational complexity, both SMDS and QD-SMDS are described as being dominated by an SVD of the GEK matrix, with similar scaling in the number of nodes. The quaternion version operates on quaternion-valued matrices, but is presented as low-complexity relative to Bayesian or convex optimization-based alternatives, and the SVD-free QD-MRC-SMDS further reduces complexity by avoiding the costly QSVD route (Masuoka et al., 24 Apr 2025).

5. Adaptation to rigid body localization

A separate development applies SMDS to rigid body localization (RBL). In this setting, a rigid body with landmark points is observed relative to anchors through measurements of distances and angles between the vehicle’s sensors and anchor landmark points. The proposed method uses a variation of SMDS in which only a minor part of the complex edge kernel is used, specifically the anchor-to-anchor and target-to-target information, while anchor-to-target edges are estimated (Führling et al., 1 Sep 2025).

The rigid body is modeled as

$\theta_p$2

where $\theta_p$3 is the conformation matrix, $\theta_p$4 is the rotation matrix, and $\theta_p$5 is the translation vector. The kernel is partitioned into blocks corresponding to anchor-anchor, anchor-target, and target-target interactions. The algorithm first uses the known anchor-anchor and target-target edges to initialize the unknown anchor-target edges, then refines those estimates iteratively through a Turbo-style update, reconstructs landmark coordinates, and finally solves for translation and rotation through a least-squares Procrustes problem (Führling et al., 1 Sep 2025).

The simulation setup reported for this method uses a $\theta_p$6 room, $\theta_p$7 anchors, $\theta_p$8 rigid-body landmark points, noisy gamma-distributed range measurements, and angular measurements with Tikhonov noise. Performance is evaluated in terms of mean squared error and the corresponding Cramér–Rao lower bound. In these experiments, the proposed full-SMDS approach yields the lowest MSE for both translation and rotation estimates and approaches the CRLB; SMDS using distance only outperforms plain MDS, especially when range errors are moderate; and in low-noise regimes SMDS is not always strictly better than MDS (Führling et al., 1 Sep 2025).

This adaptation is significant because it shows that SMDS can be deployed even when only a structured subset of the full kernel is available. That is a narrower measurement model than the fully connected edge setting often assumed in basic expositions.

6. Terminological divergence and relation to broader MDS research

The acronym “SMDS” is not unique to localization. In work on language-model interpretability, SMDS stands for “Supervised Multi-Dimensional Scaling,” a model-agnostic method for discovering feature manifolds in activations. There the objective is to learn a linear projection $\theta_p$9 such that Euclidean distances in the projected space match a label-derived hypothesis distance: $\mathbf{v}$0 A held-out stress metric is then used for model selection among candidate manifold geometries such as linear, circular, semicircular, log-linear, and cluster structures. The method is used to identify circles, lines, and clusters in temporal reasoning representations of LLMs (Tiblias et al., 1 Oct 2025).

This usage is methodologically distinct from localization-oriented Super Multidimensional Scaling. The shared acronym can create confusion, but the two literatures differ in input objects, optimization target, and intended interpretation: one addresses geometric localization from distance-angle measurements, while the other addresses supervised manifold discovery in latent representations.

A broader theoretical backdrop comes from work on multidimensional scaling on metric measure spaces. That literature generalizes MDS from finite point sets to compact metric spaces with Borel probability measures, replaces finite double-centering with an operator-theoretic kernel construction, and studies optimality through strain minimization in Hilbert–Schmidt norm. It also analyzes the MDS embedding of the geodesic circle $\mathbf{v}$1 into $\mathbf{v}$2 and into $\mathbf{v}$3, and discusses convergence of finite-sample embeddings to operator-theoretic limits (Adams et al., 2019). The summary associated with that work presents these extensions as relevant to “SMDS” in the sense of unifying and extending MDS beyond finite point clouds.

Finally, broader MDS research underscores the computational pressures that motivate low-complexity SMDS variants. Out-of-sample embedding methods for MDS have been developed to address the fact that common MDS approaches are at least quadratic in memory and computation, using landmarks together with either optimization-based or neural-network-based embedding for new points (Herath et al., 2021). A plausible implication is that the SVD-free and kernel-structured variants of SMDS belong to the same general computational trend: preserving the geometric strengths of MDS-derived methods while reducing the cost of large or noisy inference problems.