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