📄 dir_c.mht
字号:
From: <由 Microsoft Internet Explorer 5 保存>
Subject:
Date: Wed, 27 Sep 2006 22:40:20 +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/DIR.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>/* dir.c - IIR filtering in direct form */
double dir(M, a, L, b, w, v, x) /* usage: y =3D dir(M, a, L, =
b, w, v, x); */
double *a, *b, *w, *v, x; /* \(v,w\) are internal states =
*/
int M, L; /* denominator and numerator =
orders */
{
int i;
v[0] =3D x; /* current input sample */
w[0] =3D 0; /* current output to be =
computed */
for (i=3D0; i<=3DL; i++) /* numerator part */
w[0] +=3D b[i] * v[i];
for (i=3D1; i<=3DM; i++) /* denominator part */
w[0] -=3D a[i] * w[i];
for (i=3DL; i>=3D1; i--) /* reverse-order =
updating of \(v\) */
v[i] =3D v[i-1];
for (i=3DM; i>=3D1; i--) /* reverse-order =
updating of \(w\) */
w[i] =3D w[i-1];
return w[0]; /* current output sample */
}
</PRE></BODY></HTML>
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -