非線形の組み合わせ

　

ここでは、ある軸がスペクトル変数$X$ で線形にサンプリングされ、変数$S$ で表現される場合を考える。 スペクトル変数としては様々な組み合わせが考えられ、どの組み合わせも非線形でありうるが、どのスペクトル変数も $\nu, \lambda, \lambda_a, v$ のどれかと線形に結びつく。 これでもまだ多くの組み合わせがありうるので、$S$ と線形関係にある基本変数 $\nu, \lambda, \lambda_a, v$ のどれかを表す中間変数 $P$ を導入して計算方法を考えてみる。 $X = X(P)$ とその逆関係 $P = P(X)$ は基本変数 $\nu, \lambda, v$ の関係として次の左表のように表され、$S = S(P)$ とその逆関係の $P = P(S)$ は次の右表のような関係となる。

$\nu = \frac{c}{\lambda}$	$\nu = \nu_0\frac{c-v}{\sqrt{c^2-v^2}}$
$\lambda = \frac{c}{\nu}$	$\lambda = \lambda_0\frac{c+v}{\sqrt{c^2-v^2}}$
$v = c\frac{\nu_0^2 -\nu^2}{\nu_0^2 +\nu^2}$	$v = c\frac{\lambda^2 -\lambda_0^2}{\lambda^2 +\lambda_0^2}$

$\nu = \nu_0 (1-\frac{V}{c})$	$V = c\frac{\nu_0 - \nu}{\nu_0}$
$\nu = \frac{E}{h}$	$E = h\nu$
$\nu = c\kappa$	$\kappa = \frac{\nu}{c}$
$\lambda = \lambda_0 (1 + \frac{Z}{c})$	$Z = c\frac{\lambda - \lambda_0}{\lambda_0}$
$\lambda = \lambda_0 (1 + z)$	$z = \frac{\lambda - \lambda_0}{\lambda_0}$
$v = c\beta$	$\beta = \frac{v}{c}$

$S$ と $X$ の関係は中間変数である $P$ を介して $S(X) = S(P(X))$ と表される(逆関数は $X(S) = X(P(S))$)。 $S(P)$ は線形なので、$P$ と $X$ は常に異ならなければならない(そうでなければ $S(X)$ が線形になって非線形軸でなくなってしまう)。 このことから CTYPEka の最後の 3文字で使われる非線形アルゴリズムとしては次のようなものがあることになる。

コード	通常以下でサンプリングされる	以下で表現される
F2W	周波数	波長
F2V	周波数	見かけの速度
F2A	周波数	(空気中の)波長
W2F	波長	周波数
W2V	波長	見かけの速度
W2A	波長	(空気中の)波長
V2F	見かけの速度	周波数
V2W	見かけの速度	波長
V2A	見かけの速度	(空気中の)波長
A2F	(空気中の)波長	周波数
A2W	(空気中の)波長	波長
A2V	(空気中の)波長	見かけの速度
LOG	対数	どの座標タイプも可
GRI	検出器	どの座標タイプも可
GRA	検出器	どの座標タイプも可
TAB	通常ない	どの座標タイプも可

この 3文字コードの最初の 1文字は $X$ を表し、周波数(F)、波長(W)、(空気中の)波長(A)、見かけの速度(V) のどれかであり、3番目の文字は $P$ を表す。 例えば、'ZOPT-F2W' の場合 $X$ は周波数、$P$ は波長であり、この 2つの非線形変換は上左表により決まる。 最終的にスペクトル座標 $S$ として赤方偏移(ZOPT)を求めたいが、これは中間変数 $P$ (波長) と上右表の関係で結び付けられている。

線形にサンプリングされた $X$ からスペクトル座標 $S$ を求める一般的な方法は次のようになる。 $X$ は線形にサンプリングされているので、

$$X = X_r + w\frac{dX}{dw}$$ (31)

であり $dX/dw$ は定数である。この定数は、 $\frac{dS}{dw}\big\vert _r = 1$ という要求から決まる。これにより参照点では $S \approx S_r + w$ であることから、

$$\frac{dX}{dw} = \left. \frac{dP}{dS}\right\vert _r / \left. \frac{dP}{dX}\right\vert _r$$ (32)

と計算できる。 $S = S(P), P = P(X)$ が逆関数の $X = X(P), P = P(S)$ と同様にわかっているとすると、$S$ を $w$ の関数として表せて、

$$S(w) = S\left(P\left(X(P(S_r)) + w\left. \frac{dP}{dS}\right\vert _r /\left. \frac{dP}{dX}\right\vert _r\right)\right)$$ (33)

となる。 ここで $S_r$ は CRVALka で与えられる。 即ち、$S(w)$ を求めるには次の 3ステップを踏めばよい。

1. (31)を使って $w$ から $X$ を計算: $X_r = X(P(S_r))$ と $dX/dw$ は定数なので1度計算しておけばよい
2. $X$ から $P$ を計算: 適切な関係式を前ページの左表から適用する
3. $P$ から $S$ を計算: 適切な関係式を前ページの右表から適用する

原論文には具体的な計算例もあるので参照されたい。