For a given random variable $X$ with mean $\mu_{\mathrm{old}}$ and standard deviation $\sigma_{\mathrm{old}}$ I would like to perform a transformation $g$ to obtain a new random variable $Y := g(X)$ which has a desired mean and standard deviation $\mu_{\mathrm{new}}$ and $\sigma_{\mathrm{new}}$ respectively without (if possible) changing other characteristics of the distribution (higher moments). A simple calculation shows that the transformation $g$ has the form
$$Y = g(X) := \frac{\sigma_{\mathrm{new}}}{\sigma_{\mathrm{old}}}X-\frac{\sigma_{\mathrm{new}}}{\sigma_{\mathrm{old}}}\mu_{\mathrm{old}} + \mu_{\mathrm{new}}$$
This fulfills the condition above since
$$E[Y] = \frac{\sigma_{\mathrm{new}}}{\sigma_{\mathrm{old}}}E[X]-\frac{\sigma_{\mathrm{new}}}{\sigma_{\mathrm{old}}}\mu_{\mathrm{old}} + \mu_{\mathrm{new}}=\frac{\sigma_{\mathrm{new}}}{\sigma_{\mathrm{old}}}\mu_{\mathrm{old}}-\frac{\mu_{\mathrm{new}}}{\sigma_{\mathrm{old}}}\mu_{\mathrm{old}} + \mu_{\mathrm{new}}=\mu_{\mathrm{new}}$$
and
$$\operatorname{Var}(Y) = \frac{\sigma^2_{\mathrm{new}}}{\sigma^2_{\mathrm{old}}}\operatorname{Var}(X)=\sigma^2_{\mathrm{new}}.$$
Does this normalization and scaling change other moments? I did some simulation in R and it seems no to be the case, at least for kurtosis and skewness. What about higher moments? Are there any reference for this?