dtft_c.mht

来自「常用的DSP算法集锦」· MHT 代码 · 共 42 行

MHT
42
字号
From: <由 Microsoft Internet Explorer 5 保存>
Subject: 
Date: Wed, 27 Sep 2006 22:46:29 +0800
MIME-Version: 1.0
Content-Type: text/html;
	charset="gb2312"
Content-Transfer-Encoding: quoted-printable
Content-Location: http://www.oxbad.com/DSP/maindoc/arithmetic_source/c/DTFT.C
X-MimeOLE: Produced By Microsoft MimeOLE V6.00.2900.2962

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=3DContent-Type content=3D"text/html; charset=3Dgb2312">
<META content=3D"MSHTML 6.00.2900.2963" name=3DGENERATOR></HEAD>
<BODY><PRE>/* dtft.c - DTFT of length-L signal at a single frequency w =
*/

#include &lt;cmplx.h&gt;                                   /* complex =
arithmetic */

complex dtft(L, x, w)                                /* usage: =
X=3Ddtft(L, x, w); */
double *x, w;                                        /* \(x\) is =
\(L\)-dimensional */
int L;
{
       complex z, X;
       int n;

       z =3D cexp(cmplx(0, -w));                       /* set =
\(z=3De\sp{-j\om}\) */

       X =3D cmplx(0,0);                               /* initialize =
\(X=3D0\) */

       for (n=3DL-1; n&gt;=3D0; n--)
              X =3D cadd(cmplx(x[n], 0), cmul(z, X));

       return X;
}
</PRE></BODY></HTML>

⌨️ 快捷键说明

复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?