The influence of the scalar unparticle and anomalous couplings at muon colliders in final states with multiple photons in the Randall–Sundrum model is evaluated in detail. The results indicate that with fixed collision energies, the total cross-sections for the production of multiple photons depend strongly on the polarization of the muon beams, the parameters of unparticle physics (the scaling dimension [Formula: see text], operator [Formula: see text], the energy scale [Formula: see text]) and also the strength of anomalous couplings. Numerical evaluation shows that the cross-sections for the production of four photons in finale states with the contribution of scalar anomalous couplings are much larger than that of the unparticle under the same conditions. In the Higgs–radion mixing, the cross-sections achieve the maximum value at the radion-dominated state, [Formula: see text][Formula: see text]GeV, in which the cross-section is much enhanced and can be measurable in current experiments.