离散正弦变换

             

预备知识 离散傅里叶变换

   由正弦级数($n = 1, 2, 3\dots$)

\begin{equation} f(x) = \sum_n C_n \sin\left(\frac{n\pi}{l}x\right) \end{equation}
\begin{equation} C_n = \frac2l \int_0^l f(x) \sin\left(\frac{n\pi}{l}x\right) \,\mathrm{d}{x} \end{equation}
不难推出正弦变换
\begin{equation} g(k) = \sqrt{\frac{2}{\pi}} \int_0^{\infty} f(x) \sin\left(kx\right) \,\mathrm{d}{x} \end{equation}
\begin{equation} f(x) = \sqrt{\frac{2}{\pi}} \int_0^{\infty} g(k) \sin\left(kx\right) \,\mathrm{d}{k} \end{equation}
注意这是一个正半轴的变换,且正反变换相同.

   正弦变换同样有采样定理,即若 $g(k)$ 的区间为 $[0, L_k]$,那么只需要取 $\Delta x = \pi/L_k$ 对 $f(x)$ 采样即可用以下插值公式精确还原 $f(x)$

\begin{equation} f(x) = \sum_{n=1}^\infty f(x_n)\frac{2x_n}{x+x_n} \operatorname{sinc} [\pi(x-x_n)/\Delta x] \end{equation}

1. 离散正弦变换

   把插值公式做正弦变换,得

\begin{equation} g(k) = \sqrt{\frac{2}{\pi}} \sum_{n=1}^\infty f(x_n) \sin\left(k x_n\right) \Delta x \end{equation}
现在假设 $f(x)$ 和 $g(k)$ 都只在 $[0, L_x]$ 和 $[0, L_k]$ 内,所以有
\begin{equation} \Delta x L_k = \Delta k L_x = N\Delta x\Delta k = \frac{L_xL_k}{N} = \pi \end{equation}
可得无损的离散正弦变换为
\begin{equation} g_q = \sum_{p = 1}^{N-1} f_p \sin\left(\pi pq/N\right) \end{equation}
\begin{equation} f_p = \sum_{q = 1}^{N-1} g_q \sin\left(\pi pq/N\right) \end{equation}
可以证明变换矩阵是对称的单位正交矩阵,所以逆矩阵就是矩阵本身.

致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择会员制,大量广告,内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 10 元,我们一个星期内就能脱离亏损, 并保证网站能在接下来的一整年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

         

© 小时科技 保留一切权利