标签：right frac 平滑 曲线 cdots pi sum left
曲线平滑
1.曲线平滑步骤
第一步：创建频率域加权系数数组
\[{\rm{Gamma(i) = }}\frac{1}{{1 + par{{\left( {2 - 2\cos \frac{{i \cdot PI}}{N}} \right)}^2}}} \]其中\(par\)为平滑算法的平滑参数，\(par\)越大，平滑成都越狠，\(N\)为待平滑数组的总点数。
第二步：对待平滑数组\(X\)进行余弦正变换，转换到频率域,得到\(Xout\)
第三步：在频率域进行能量调整 \(xTemp(i)=Gamma(i)*Xout(i)\)
第四步：对调整后的频率域\(xTemp\)进行反余弦变换，得到平滑后的频率域一维数组。
2.一维余弦变换
2.1 傅里叶变换与反变换的公式为
\[F(u) = c\sum\limits_{n = 0}^{N - 1} f (n){e^{ - j\frac{{2\pi un}}{N}}},u = 0,1, \cdots ,N - 1 \]\[f(n) = d\sum\limits_{u = 0}^{N - 1} F (u){e^{j\frac{{2\pi un}}{N}}},n = 0,1, \cdots ,N - 1 \]2.2 DCT公式推导
基于公式
\[F(u) = c\sum\limits_{n = 0}^{N - 1} f (n){e^{ - j\frac{{2\pi un}}{N}}},u = 0,1, \cdots ,N - 1 \]令
\[f(n) = g(n){e^{ - j\frac{{p\pi u}}{{qN}}}},p,q \ne 0 \]其中\(g(n)\)通过个数为\(tN\left( {t \in {Z_{\ge 1}}} \right)\)的一个空间域（时域）实数采样序列：
\[g(0),g(1), \cdots ,g(N - 1),g(N), \cdots ,g(tN - 1),\forall n,g(n) \in R \]且当\(t>1\)时，\(\forall n \in \left\{ {N,N + 1, \cdots ,tN - 1} \right\}\),\(g(n)=0\),即在拥有前\(N\)个采样点的基础上新增\((t-1)N\)个零采样点，则有：
\[\begin{array}{l} F(u) = c\sum\limits_{n = 0}^{tN - 1} {\left[ {\left( {g(n){e^{ - j\frac{{p\pi u}}{{qN}}}}} \right) \cdot {e^{ - j\frac{{2\pi un}}{{tN}}}}} \right]} \$ \$ \\ = c\sum\limits_{n = 0}^{tN - 1} {\left[ {g(n) \cdot \left( {{e^{ - j\frac{{p\pi u}}{{qN}}}} \cdot {e^{ - j\frac{{2\pi un}}{{tN}}}}} \right)} \right]} \\ = c\sum\limits_{n = 0}^{tN - 1} g (n){e^{ - j\frac{{(pt + 2nq)u\pi }}{{qtN}}}}\$ \$ \\ = c\sum\limits_{n = 0}^{tN - 1} g (n)\left( {\cos \frac{{(pt + 2nq)u\pi }}{{qtN}} - i\sin \frac{{(pt + 2nq)u\pi }}{{qtN}}} \right)\$ \$ \\ = c\sum\limits_{n = 0}^{N - 1} g (n)\left( {\cos \frac{{(pt + 2nq)u\pi }}{{qtN}} - i\sin \frac{{(pt + 2nq)u\pi }}{{qtN}}} \right) \end{array} \]对\(F(u)\)取实部，得到计算\(C(u)\)的DCT公式
\[C(u)=\operatorname{Re}(F(u))=c \sum_{n=0}^{N-1} g(n) \cos \frac{(p t+2 n q) u \pi}{q t N}, u=0,1, \cdots, N-1 \]在DCT公式一般形式中，令
\[p=1, q=t=2 \]则有
\[C(u)=c \sum_{n=0}^{N-1} g(n) \cos \frac{(2+2 \times n \times 2) u \pi}{4 N}\\ =c \sum_{n=0}^{N-1} g(n) \cos \frac{(2 n+1) u \pi}{2 N}, u=0,1, \cdots, N-1 \]2.3IDCT公式推导
基于公式
\[f(n)=d \sum_{u=0}^{N-1} F(u) e^{j \frac{2 \pi n u}{N}}, n=0,1, \cdots, N-1 \]令
\[F(u)=C(u) e^{j \frac{p \pi u}{q N}} \]其中\(C(u)\)通过个数为\(tN\left( {t \in {Z_{\ge 1}}} \right)\)的一个频率域实数离散序列：
\[C(0), C(1), \cdots, C(N-1), C(N), \cdots, C(t N-1), \forall u, C(u) \in \mathbb{R} \]且当\(t>1\)时，\(\forall u \in\{N, N+1, \cdots, t N-1\}\)，\(C(u)=0\),相当于在拥有前N个离散点的基础上新增\((t-1)N\)个无效零函数点，则有：
\[f(n)=d \sum_{u=0}^{t N-1}\left[\left(C(u) e^{j \frac{p \pi u}{q N}}\right) \cdot e^{j \frac{2 \pi n u}{t N}}\right]\\ =d \sum_{u=0}^{t N-1}\left[C(u) \cdot\left(e^{j \frac{p \pi u}{q N}} \cdot e^{j \frac{2 \pi n u}{t N N}}\right)\right]\\ =d \sum_{u=0}^{t N-1} C(u) e^{j \frac{(p t+2 n g) u \pi}{q t N}}\\ =d \sum_{u=0}^{t N-1} C(u)\left(\cos \frac{(p t+2 n q) u \pi}{q t N}-i \sin \frac{(p t+2 n q) u \pi}{q t N}\right)\\ =d \sum_{u=0}^{N-1} C(u)\left(\cos \frac{(p t+2 n q) u \pi}{q t N}-i \sin \frac{(p t+2 n q) u \pi}{q t N}\right) \]对\(f(n)\)取实部，得到计算\(g(n)\)的IDCT公式
\[g(n)=\operatorname{Re}(f(n))=d \sum_{u=0}^{N-1} C(u) \cos \frac{(p t+2 n q) u \pi}{q t N}, n=0,1, \cdots, N-1 \]在IDCT公式一般形式中，令
\[p=1,q=t=2 \]则有
\[g(n)=d \sum_{u=0}^{N-1} C(u) \cos \frac{(1 \times 2+2 \times n \times 2) u \pi}{2 \times 2 \times N}\\ =d \sum_{u=0}^{N-1} C(u) \cos \frac{(2 n+1) u \pi}{2 N}, n=0,1, \cdots, N-1 \]标签：right,frac,平滑,曲线,cdots,pi,sum,left 来源： https://www.cnblogs.com/GeophysicsWorker/p/16470084.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。