开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Linux/Unix编程 > 查看源码
voicin.c
资源名称：一个类似QQ的软件源码.rar　[点击查看]
上传用户：zslianheng
上传日期：2013-04-03
资源大小：946k
文件大小：28k
源码类别：

Linux/Unix编程

开发平台：

Visual C++

voicin.c：源码内容
  1. /*
  3. $Log: voicin.c,v $
  4. Revision 1.1.1.1  2002/07/28 05:23:51  freeman_yong
  5. lpc10 codec
  7.  * Revision 1.2  1996/08/20  20:45:00  jaf
  8.  * Removed all static local variables that were SAVE'd in the Fortran
  9.  * code, and put them in struct lpc10_encoder_state that is passed as an
  10.  * argument.
  11.  *
  12.  * Removed init function, since all initialization is now done in
  13.  * init_lpc10_encoder_state().
  14.  *
  15.  * Revision 1.1  1996/08/19  22:30:14  jaf
  16.  * Initial revision
  17.  *
  19. */
  21. #ifdef P_R_O_T_O_T_Y_P_E_S
  22. extern int voicin_(integer *vwin, real *inbuf, real *lpbuf, integer *buflim, integer *half, real *minamd, real *maxamd, integer *mintau, real *ivrc, integer *obound, integer *voibuf, integer *af, struct lpc10_encoder_state *st);
  23. /* comlen contrl_ 12 */
  24. /*:ref: vparms_ 14 14 4 6 6 4 4 6 4 4 4 4 6 6 6 6 */
  25. #endif
  27. /*  -- translated by f2c (version 19951025).
  28.    You must link the resulting object file with the libraries:
  29. -lf2c -lm   (in that order)
  30. */
  32. #include "f2c.h"
  34. /* Common Block Declarations */
  36. extern struct {
  37.     integer order, lframe;
  38.     logical corrp;
  39. } contrl_;
  41. #define contrl_1 contrl_
  43. /****************************************************************************/
  45. /*  VOICIN Version 52 */
  47. /* $Log: voicin.c,v $
  48. /* Revision 1.1.1.1  2002/07/28 05:23:51  freeman_yong
  49. /* lpc10 codec
  50. /*
  51.  * Revision 1.2  1996/08/20  20:45:00  jaf
  52.  * Removed all static local variables that were SAVE'd in the Fortran
  53.  * code, and put them in struct lpc10_encoder_state that is passed as an
  54.  * argument.
  55.  *
  56.  * Removed init function, since all initialization is now done in
  57.  * init_lpc10_encoder_state().
  58.  *
  59.  * Revision 1.1  1996/08/19  22:30:14  jaf
  60.  * Initial revision
  61.  * */
  62. /* Revision 1.10  1996/03/29  17:59:14  jaf */
  63. /* Avoided using VALUE(9), although it shouldn't affect the function of */
  64. /* the code at all, because it was always multiplied by VDC(9,SNRL), */
  65. /* which is 0 for all values of SNRL.  Still, if VALUE(9) had an initial */
  66. /* value of IEEE NaN, it might cause trouble (I don't know how IEEE */
  67. /* defines Nan * 0.  It should either be NaN or 0.) */
  69. /* Revision 1.9  1996/03/29  17:54:46  jaf */
  70. /* Added a few comments about the accesses made to argument array VOIBUF */
  71. /* and the local saved array VOICE. */
  73. /* Revision 1.8  1996/03/27  18:19:54  jaf */
  74. /* Added an assignment to VSTATE that does not affect the function of the */
  75. /* program at all.  The only reason I put it in was so that the tracing */
  76. /* statements at the end, when enabled, will print a consistent value for */
  77. /* VSTATE when HALF .EQ. 1, rather than a garbage value that could change */
  78. /* from one call to the next. */
  80. /* Revision 1.7  1996/03/26  20:00:06  jaf */
  81. /* Removed the inclusion of the file "vcomm.fh", and put its contents */
  82. /* into this file.  It was included nowhere else but here. */
  84. /* Revision 1.6  1996/03/26  19:38:09  jaf */
  85. /* Commented out trace statements. */
  87. /* Revision 1.5  1996/03/19  20:43:45  jaf */
  88. /* Added comments about which indices of OBOUND and VOIBUF can be */
  89. /* accessed, and whether they are read or written.  VOIBUF is fairly */
  90. /* messy. */
  92. /* Revision 1.4  1996/03/19  15:00:58  jaf */
  93. /* Moved the DATA statements for the *VDC* variables later, as it is */
  94. /* apparently illegal to have DATA statements before local variable */
  95. /* declarations. */
  97. /* Revision 1.3  1996/03/19  00:10:49  jaf */
  98. /* Heavily commented the local variables that are saved from one */
  99. /* invocation to the next, and how the local variable FIRST is used to */
  100. /* avoid the need to assign most of them initial values with DATA */
  101. /* statements. */
  103. /* A few should be initialized, but aren't.  I've guessed initial values */
  104. /* for two of these, SFBUE and SLBUE, and I've convinced myself that for */
  105. /* VOICE, the effects of uninitialized values will die out after 2 or 3 */
  106. /* frame times.  It would still be good to choose initial values for */
  107. /* these, but I don't know what reasonable values would be (0 comes to */
  108. /* mind). */
  110. /* Revision 1.2  1996/03/13  16:09:28  jaf */
  111. /* Comments added explaining which of the local variables of this */
  112. /* subroutine need to be saved from one invocation to the next, and which */
  113. /* do not. */
  115. /* WARNING!  Some of them that should are never given initial values in */
  116. /* this code.  Hopefully, Fortran 77 defines initial values for them, but */
  117. /* even so, giving them explicit initial values is preferable. */
  119. /* WARNING!  VALUE(9) is used, but never assigned a value.  It should */
  120. /* probably be eliminated from the code. */
  122. /* Revision 1.1  1996/02/07 14:50:28  jaf */
  123. /* Initial revision */
  126. /****************************************************************************/
  128. /*        Voicing Detection (VOICIN) makes voicing decisions for each half */
  129. /*  frame of input speech.  Tentative voicing decisions are made two frames*/
  130. /*   in the future (2F) for each half frame.  These decisions are carried */
  131. /*   through one frame in the future (1F) to the present (P) frame where */
  132. /*   they are examined and smoothed, resulting in the final voicing */
  133. /*   decisions for each half frame. */
  134. /*        The voicing parameter (signal measurement) column vector (VALUE) */
  135. /*   is based on a rectangular window of speech samples determined by the */
  136. /*  window placement algorithm.  The voicing parameter vector contains the*/
  137. /*  AMDF windowed maximum-to-minimum ratio, the zero crossing rate, energy*/
  138. /*   measures, reflection coefficients, and prediction gains.  The voicing */
  139. /*  window is placed to avoid contamination of the voicing parameter vector*/
  140. /*   with speech onsets. */
  141. /*        The input signal is then classified as unvoiced (including */
  142. /*   silence) or voiced.  This decision is made by a linear discriminant */
  143. /*   function consisting of a dot product of the voicing decision */
  144. /*   coefficient (VDC) row vector with the measurement column vector */
  145. /*  (VALUE).  The VDC vector is 2-dimensional, each row vector is optimized*/
  146. /*   for a particular signal-to-noise ratio (SNR).  So, before the dot */
  147. /*   product is performed, the SNR is estimated to select the appropriate */
  148. /*   VDC vector. */
  149. /*        The smoothing algorithm is a modified median smoother.  The */
  150. /*  voicing discriminant function is used by the smoother to determine how*/
  151. /*   strongly voiced or unvoiced a signal is.  The smoothing is further */
  152. /*   modified if a speech onset and a voicing decision transition occur */
  153. /*   within one half frame.  In this case, the voicing decision transition */
  154. /*  is extended to the speech onset.  For transmission purposes, there are*/
  155. /*   constraints on the duration and transition of voicing decisions.  The */
  156. /*   smoother takes these constraints into account. */
  157. /*        Finally, the energy estimates are updated along with the dither */
  158. /*   threshold used to calculate the zero crossing rate (ZC). */
  160. /* Inputs: */
  161. /*  VWIN      - Voicing window limits */
  162. /*              The indices read of arrays VWIN, INBUF, LPBUF, and BUFLIM */
  163. /*              are the same as those read by subroutine VPARMS. */
  164. /*  INBUF     - Input speech buffer */
  165. /*  LPBUF     - Low-pass filtered speech buffer */
  166. /*  BUFLIM    - INBUF and LPBUF limits */
  167. /*  HALF      - Present analysis half frame number */
  168. /*  MINAMD    - Minimum value of the AMDF */
  169. /*  MAXAMD    - Maximum value of the AMDF */
  170. /*  MINTAU    - Pointer to the lag of the minimum AMDF value */
  171. /*  IVRC(2)   - Inverse filter's RC's */
  172. /*              Only index 2 of array IVRC read under normal operation. */
  173. /*              (Index 1 is also read when debugging is turned on.) */
  174. /*  OBOUND    - Onset boundary descriptions */
  175. /*             Indices 1 through 3 read if (HALF .NE. 1), otherwise untouched.
  176. */
  177. /*  AF        - The analysis frame number */
  178. /* Output: */
  179. /*  VOIBUF(2,0:AF) - Buffer of voicing decisions */
  180. /*              Index (HALF,3) written. */
  181. /*              If (HALF .EQ. 1), skip down to "Read (HALF,3)" below. */
  182. /*              Indices (1,2), (2,1), (1,2), and (2,2) read. */
  183. /*              One of the following is then done: */
  184. /*                 read (1,3) and possibly write (1,2) */
  185. /*                 read (1,3) and write (1,2) or (2,2) */
  186. /*                 write (2,1) */
  187. /*                 write (2,1) or (1,2) */
  188. /*                 read (1,0) and (1,3) and then write (2,2) or (1,1) */
  189. /*                 no reads or writes on VOIBUF */
  190. /*              Finally, read (HALF,3) */
  191. /* Internal: */
  192. /*  QS        - Ratio of preemphasized to full-band energies */
  193. /*  RC1       - First reflection coefficient */
  194. /* AR_B      - Product of the causal forward and reverse pitch prediction gain
  195. s*/
  196. /* AR_F      - Product of the noncausal forward and rev. pitch prediction gain
  197. s*/
  198. /*  ZC        - Zero crossing rate */
  199. /*  DITHER    - Zero crossing threshold level */
  200. /*  MAXMIN    - AMDF's 1 octave windowed maximum-to-minimum ratio */
  201. /*  MINPTR    - Location  of minimum AMDF value */
  202. /*  NVDC      - Number of elements in each VDC vector */
  203. /*  NVDCL     - Number of VDC vectors */
  204. /*  VDCL      - SNR values corresponding to the set of VDC's */
  205. /*  VDC       - 2-D voicing decision coefficient vector */
  206. /*  VALUE(9)  - Voicing Parameters */
  207. /*  VOICE(2,3)- History of LDA results */
  208. /*              On every call when (HALF .EQ. 1), VOICE(*,I+1) is */
  209. /*              shifted back to VOICE(*,I), for I=1,2. */
  210. /*              VOICE(HALF,3) is written on every call. */
  211. /*              Depending on several conditions, one or more of */
  212. /*              (1,1), (1,2), (2,1), and (2,2) might then be read. */
  213. /*  LBE       - Ratio of low-band instantaneous to average energies */
  214. /*  FBE       - Ratio of full-band instantaneous to average energies */
  215. /*  LBVE      - Low band voiced energy */
  216. /*  LBUE      - Low band unvoiced energy */
  217. /*  FBVE      - Full band voiced energy */
  218. /*  FBUE      - Full band unvoiced energy */
  219. /*  OFBUE     - Previous full-band unvoiced energy */
  220. /*  OLBUE     - Previous low-band unvoiced energy */
  221. /*  REF       - Reference energy for initialization and DITHER threshold */
  222. /*  SNR       - Estimate of signal-to-noise ratio */
  223. /*  SNR2      - Estimate of low-band signal-to-noise ratio */
  224. /*  SNRL      - SNR level number */
  225. /*  OT        - Onset transition present */
  226. /*  VSTATE    - Decimal interpretation of binary voicing classifications */
  227. /*  FIRST     - First call flag */
  229. /* This subroutine maintains local state from one call to the next.  If */
  230. /* you want to switch to using a new audio stream for this filter, or */
  231. /* reinitialize its state for any other reason, call the ENTRY */
  232. /* INITVOICIN. */
  234. /* Subroutine */ int voicin_(integer *vwin, real *inbuf, real *
  235. lpbuf, integer *buflim, integer *half, real *minamd, real *maxamd, 
  236. integer *mintau, real *ivrc, integer *obound, integer *voibuf, 
  237. integer *af, struct lpc10_encoder_state *st)
  238. {
  239.     /* Initialized data */
  241.     real *dither;
  242.     static real vdc[100] /* was [10][10] */ = { 0.f,1714.f,-110.f,
  243.     334.f,-4096.f,-654.f,3752.f,3769.f,0.f,1181.f,0.f,874.f,-97.f,
  244.     300.f,-4096.f,-1021.f,2451.f,2527.f,0.f,-500.f,0.f,510.f,-70.f,
  245.     250.f,-4096.f,-1270.f,2194.f,2491.f,0.f,-1500.f,0.f,500.f,-10.f,
  246.     200.f,-4096.f,-1300.f,2e3f,2e3f,0.f,-2e3f,0.f,500.f,0.f,0.f,
  247.     -4096.f,-1300.f,2e3f,2e3f,0.f,-2500.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,
  248.     0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,
  249.     0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,
  250.     0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f,0.f };
  251.     static integer nvdcl = 5;
  252.     static real vdcl[10] = { 600.f,450.f,300.f,200.f,0.f,0.f,0.f,0.f,0.f,0.f }
  253.     ;
  255.     /* System generated locals */
  256.     integer inbuf_offset, lpbuf_offset, i__1, i__2;
  257.     real r__1, r__2;
  259.     /* Builtin functions */
  260.     integer i_nint(real *);
  261.     double sqrt(doublereal);
  263.     /* Local variables */
  264.     real ar_b__, ar_f__;
  265.     integer *lbve, *lbue, *fbve, *fbue;
  266.     integer snrl, i__;
  267.     integer *ofbue, *sfbue;
  268.     real *voice;
  269.     integer *olbue, *slbue;
  270.     real value[9];
  271.     integer zc;
  272.     logical ot;
  273.     real qs;
  274.     real *maxmin;
  275.     integer vstate;
  276.     real rc1;
  277.     extern /* Subroutine */ int vparms_(integer *, real *, real *, integer *, 
  278.     integer *, real *, integer *, integer *, integer *, integer *, 
  279.     real *, real *, real *, real *);
  280.     integer fbe, lbe;
  281.     real *snr;
  282.     real snr2;
  284. /*  Global Variables: */
  285. /*       Arguments */
  286. /* $Log: voicin.c,v $
  287. /* Revision 1.1.1.1  2002/07/28 05:23:51  freeman_yong
  288. /* lpc10 codec
  289. /*
  290.  * Revision 1.2  1996/08/20  20:45:00  jaf
  291.  * Removed all static local variables that were SAVE'd in the Fortran
  292.  * code, and put them in struct lpc10_encoder_state that is passed as an
  293.  * argument.
  294.  *
  295.  * Removed init function, since all initialization is now done in
  296.  * init_lpc10_encoder_state().
  297.  *
  298.  * Revision 1.1  1996/08/19  22:30:14  jaf
  299.  * Initial revision
  300.  * */
  301. /* Revision 1.3  1996/03/29  22:05:55  jaf */
  302. /* Commented out the common block variables that are not needed by the */
  303. /* embedded version. */
  305. /* Revision 1.2  1996/03/26  19:34:50  jaf */
  306. /* Added comments indicating which constants are not needed in an */
  307. /* application that uses the LPC-10 coder. */
  309. /* Revision 1.1  1996/02/07  14:44:09  jaf */
  310. /* Initial revision */
  312. /*   LPC Processing control variables: */
  314. /* *** Read-only: initialized in setup */
  316. /*  Files for Speech, Parameter, and Bitstream Input & Output, */
  317. /*    and message and debug outputs. */
  319. /* Here are the only files which use these variables: */
  321. /* lpcsim.f setup.f trans.f error.f vqsetup.f */
  323. /* Many files which use fdebug are not listed, since it is only used in */
  324. /* those other files conditionally, to print trace statements. */
  325. /*  integer fsi, fso, fpi, fpo, fbi, fbo, pbin, fmsg, fdebug */
  326. /*  LPC order, Frame size, Quantization rate, Bits per frame, */
  327. /*    Error correction */
  328. /* Subroutine SETUP is the only place where order is assigned a value, */
  329. /* and that value is 10.  It could increase efficiency 1% or so to */
  330. /* declare order as a constant (i.e., a Fortran PARAMETER) instead of as 
  331. */
  332. /* a variable in a COMMON block, since it is used in many places in the */
  333. /* core of the coding and decoding routines.  Actually, I take that back. 
  334. */
  335. /* At least when compiling with f2c, the upper bound of DO loops is */
  336. /* stored in a local variable before the DO loop begins, and then that is 
  337. */
  338. /* compared against on each iteration. */
  339. /* Similarly for lframe, which is given a value of MAXFRM in SETUP. */
  340. /* Similarly for quant, which is given a value of 2400 in SETUP.  quant */
  341. /* is used in only a few places, and never in the core coding and */
  342. /* decoding routines, so it could be eliminated entirely. */
  343. /* nbits is similar to quant, and is given a value of 54 in SETUP. */
  344. /* corrp is given a value of .TRUE. in SETUP, and is only used in the */
  345. /* subroutines ENCODE and DECODE.  It doesn't affect the speed of the */
  346. /* coder significantly whether it is .TRUE. or .FALSE., or whether it is 
  347. */
  348. /* a constant or a variable, since it is only examined once per frame. */
  349. /* Leaving it as a variable that is set to .TRUE.  seems like a good */
  350. /* idea, since it does enable some error-correction capability for */
  351. /* unvoiced frames, with no change in the coding rate, and no noticeable 
  352. */
  353. /* quality difference in the decoded speech. */
  354. /*  integer quant, nbits */
  355. /* *** Read/write: variables for debugging, not needed for LPC algorithm 
  356. */
  358. /*  Current frame, Unstable frames, Output clip count, Max onset buffer, 
  359. */
  360. /*    Debug listing detail level, Line count on listing page */
  362. /* nframe is not needed for an embedded LPC10 at all. */
  363. /* nunsfm is initialized to 0 in SETUP, and incremented in subroutine */
  364. /* ERROR, which is only called from RCCHK.  When LPC10 is embedded into */
  365. /* an application, I would recommend removing the call to ERROR in RCCHK, 
  366. */
  367. /* and remove ERROR and nunsfm completely. */
  368. /* iclip is initialized to 0 in SETUP, and incremented in entry SWRITE in 
  369. */
  370. /* sread.f.  When LPC10 is embedded into an application, one might want */
  371. /* to cause it to be incremented in a routine that takes the output of */
  372. /* SYNTHS and sends it to an audio device.  It could be optionally */
  373. /* displayed, for those that might want to know what it is. */
  374. /* maxosp is never initialized to 0 in SETUP, although it probably should 
  375. */
  376. /* be, and it is updated in subroutine ANALYS.  I doubt that its value */
  377. /* would be of much interest to an application in which LPC10 is */
  378. /* embedded. */
  379. /* listl and lincnt are not needed for an embedded LPC10 at all. */
  380. /*  integer nframe, nunsfm, iclip, maxosp, listl, lincnt */
  381. /*  common /contrl/ fsi, fso, fpi, fpo, fbi, fbo, pbin, fmsg, fdebug */
  382. /*  common /contrl/ quant, nbits */
  383. /*  common /contrl/ nframe, nunsfm, iclip, maxosp, listl, lincnt */
  384. /*  Parameters/constants */
  385. /*       Voicing coefficient and Linear Discriminant Analysis variables: 
  386. */
  387. /*       Max number of VDC's and VDC levels */
  388. /*       The following are not Fortran PARAMETER's, but they are */
  389. /*       initialized with DATA statements, and never modified. */
  390. /*       Actual number of VDC's and levels */
  391. /*       Local variables that need not be saved */
  392. /*       Note: */
  394. /*       VALUE(1) through VALUE(8) are assigned values, but VALUE(9) */
  395. /*       never is.  Yet VALUE(9) is read in the loop that begins "DO I = 
  396. */
  397. /*       1, 9" below.  I believe that this doesn't cause any problems in 
  398. */
  399. /*       this subroutine, because all VDC(9,*) array elements are 0, and 
  400. */
  401. /*       this is what is multiplied by VALUE(9) in all cases.  Still, it 
  402. */
  403. /*       would save a multiplication to change the loop to "DO I = 1, 8". 
  404. */
  405. /*       Local state */
  406. /*       WARNING! */
  408. /*       VOICE, SFBUE, and SLBUE should be saved from one invocation to */
  409. /*       the next, but they are never given an initial value. */
  411. /*       Does Fortran 77 specify some default initial value, like 0, or */
  412. /*       is it undefined?  If it is undefined, then this code should be */
  413. /*       corrected to specify an initial value. */
  415. /*       For VOICE, note that it is "shifted" in the statement that */
  416. /*       begins "IF (HALF .EQ. 1) THEN" below.  Also, uninitialized */
  417. /*       values in the VOICE array can only affect entries in the VOIBUF 
  418. */
  419. /*       array that are for the same frame, or for an older frame.  Thus 
  420. */
  421. /*       the effects of uninitialized values in VOICE cannot linger on */
  422. /*       for more than 2 or 3 frame times. */
  424. /*       For SFBUE and SLBUE, the effects of uninitialized values can */
  425. /*       linger on for many frame times, because their previous values */
  426. /*       are exponentially decayed.  Thus it is more important to choose 
  427. */
  428. /*       initial values for these variables.  I would guess that a */
  429. /*       reasonable initial value for SFBUE is REF/16, the same as used */
  430. /*       for FBUE and OFBUE.  Similarly, SLBUE can be initialized to */
  431. /*       REF/32, the same as for LBUE and OLBUE. */
  433. /*       These guessed initial values should be validated by re-running */
  434. /*       the modified program on some audio samples. */
  436. /*   Declare and initialize filters: */
  438.     dither = (&st->dither);
  439.     snr = (&st->snr);
  440.     maxmin = (&st->maxmin);
  441.     voice = (&st->voice[0]);
  442.     lbve = (&st->lbve);
  443.     lbue = (&st->lbue);
  444.     fbve = (&st->fbve);
  445.     fbue = (&st->fbue);
  446.     ofbue = (&st->ofbue);
  447.     olbue = (&st->olbue);
  448.     sfbue = (&st->sfbue);
  449.     slbue = (&st->slbue);
  451.     /* Parameter adjustments */
  452.     if (vwin) {
  453. --vwin;
  454. }
  455.     if (buflim) {
  456. --buflim;
  457. }
  458.     if (inbuf) {
  459. inbuf_offset = buflim[1];
  460. inbuf -= inbuf_offset;
  461. }
  462.     if (lpbuf) {
  463. lpbuf_offset = buflim[3];
  464. lpbuf -= lpbuf_offset;
  465. }
  466.     if (ivrc) {
  467. --ivrc;
  468. }
  469.     if (obound) {
  470. --obound;
  471. }
  472.     if (voibuf) {
  473. --voibuf;
  474. }
  476.     /* Function Body */
  478. /*       The following variables are saved from one invocation to the */
  479. /*       next, but are not initialized with DATA statements.  This is */
  480. /*       acceptable, because FIRST is initialized ot .TRUE., and the */
  481. /*       first time that this subroutine is then called, they are all */
  482. /*       given initial values. */
  484. /*       SNR */
  485. /*       LBVE, LBUE, FBVE, FBUE, OFBUE, OLBUE */
  487. /*       MAXMIN is initialized on the first call, assuming that HALF */
  488. /*       .EQ. 1 on first call.  This is how ANALYS calls this subroutine. 
  489. */
  491. /*   Voicing Decision Parameter vector (* denotes zero coefficient): */
  493. /*  * MAXMIN */
  494. /*    LBE/LBVE */
  495. /*    ZC */
  496. /*    RC1 */
  497. /*    QS */
  498. /*    IVRC2 */
  499. /*    aR_B */
  500. /*    aR_F */
  501. /*  * LOG(LBE/LBVE) */
  502. /*  Define 2-D voicing decision coefficient vector according to the voicin
  503. g*/
  504. /*  parameter order above.  Each row (VDC vector) is optimized for a speci
  505. fic*/
  506. /*   SNR.  The last element of the vector is the constant. */
  507. /*           E    ZC    RC1    Qs   IVRC2  aRb   aRf        c */
  509. /*  The VOICE array contains the result of the linear discriminant functio
  510. n*/
  511. /*   (analog values).  The VOIBUF array contains the hard-limited binary 
  512. */
  513. /*   voicing decisions.  The VOICE and VOIBUF arrays, according to FORTRAN
  514.  */
  515. /*   memory allocation, are addressed as: */
  517. /*     (half-frame number, future-frame number) */
  519. /*     |   Past    |  Present  |  Future1  |  Future2  | */
  520. /*     | 1,0 | 2,0 | 1,1 | 2,1 | 1,2 | 2,2 | 1,3 | 2,3 |  --->  time */
  522. /*   Update linear discriminant function history each frame: */
  523.     if (*half == 1) {
  524. voice[0] = voice[2];
  525. voice[1] = voice[3];
  526. voice[2] = voice[4];
  527. voice[3] = voice[5];
  528. *maxmin = *maxamd / max(*minamd,1.f);
  529.     }
  530. /*   Calculate voicing parameters twice per frame: */
  531.     vparms_(&vwin[1], &inbuf[inbuf_offset], &lpbuf[lpbuf_offset], &buflim[1], 
  532.     half, dither, mintau, &zc, &lbe, &fbe, &qs, &rc1, &ar_b__, &
  533.     ar_f__);
  534. /*   Estimate signal-to-noise ratio to select the appropriate VDC vector. 
  535. */
  536. /*   The SNR is estimated as the running average of the ratio of the */
  537. /*   running average full-band voiced energy to the running average */
  538. /*   full-band unvoiced energy. SNR filter has gain of 63. */
  539.     r__1 = (*snr + *fbve / (real) max(*fbue,1)) * 63 / 64.f;
  540.     *snr = (real) i_nint(&r__1);
  541.     snr2 = *snr * *fbue / max(*lbue,1);
  542. /*   Quantize SNR to SNRL according to VDCL thresholds. */
  543.     snrl = 1;
  544.     i__1 = nvdcl - 1;
  545.     for (snrl = 1; snrl <= i__1; ++snrl) {
  546. if (snr2 > vdcl[snrl - 1]) {
  547.     goto L69;
  548. }
  549.     }
  550. /*    (Note:  SNRL = NVDCL here) */
  551. L69:
  552. /*   Linear discriminant voicing parameters: */
  553.     value[0] = *maxmin;
  554.     value[1] = (real) lbe / max(*lbve,1);
  555.     value[2] = (real) zc;
  556.     value[3] = rc1;
  557.     value[4] = qs;
  558.     value[5] = ivrc[2];
  559.     value[6] = ar_b__;
  560.     value[7] = ar_f__;
  561. /*   Evaluation of linear discriminant function: */
  562.     voice[*half + 3] = vdc[snrl * 10 - 1];
  563.     for (i__ = 1; i__ <= 8; ++i__) {
  564. voice[*half + 3] += vdc[i__ + snrl * 10 - 11] * value[i__ - 1];
  565.     }
  566. /*   Classify as voiced if discriminant > 0, otherwise unvoiced */
  567. /*   Voicing decision for current half-frame:  1 = Voiced; 0 = Unvoiced */
  568.     if (voice[*half + 3] > 0.f) {
  569. voibuf[*half + 6] = 1;
  570.     } else {
  571. voibuf[*half + 6] = 0;
  572.     }
  573. /*   Skip voicing decision smoothing in first half-frame: */
  574. /*     Give a value to VSTATE, so that trace statements below will print 
  575. */
  576. /*     a consistent value from one call to the next when HALF .EQ. 1. */
  577. /*     The value of VSTATE is not used for any other purpose when this is 
  578. */
  579. /*     true. */
  580.     vstate = -1;
  581.     if (*half == 1) {
  582. goto L99;
  583.     }
  584. /*   Voicing decision smoothing rules (override of linear combination): */
  586. /*  Unvoiced half-frames:  At least two in a row. */
  587. /*  -------------------- */
  589. /*  Voiced half-frames:    At least two in a row in one frame. */
  590. /*  -------------------    Otherwise at least three in a row. */
  591. /*         (Due to the way transition frames are encoded) */
  593. /*  In many cases, the discriminant function determines how to smooth. */
  594. /* In the following chart, the decisions marked with a * may be overridden
  595. .*/
  597. /*   Voicing override of transitions at onsets: */
  598. /*  If a V/UV or UV/V voicing decision transition occurs within one-half 
  599. */
  600. /*  frame of an onset bounding a voicing window, then the transition is */
  601. /*  moved to occur at the onset. */
  603. /*  P 1F */
  604. /*  ----- ----- */
  605. /*  0   0   0   0 */
  606. /*  0   0   0*  1 (If there is an onset there) */
  607. /*  0   0   1*  0* (Based on 2F and discriminant distance) */
  608. /*  0   0   1   1 */
  609. /*  0   1*  0   0 (Always) */
  610. /*  0   1*  0*  1 (Based on discriminant distance) */
  611. /*  0*  1   1   0* (Based on past, 2F, and discriminant distance) */
  612. /*  0   1*  1   1 (If there is an onset there) */
  613. /*  1   0*  0   0 (If there is an onset there) */
  614. /*  1   0   0   1 */
  615. /*  1   0*  1*  0 (Based on discriminant distance) */
  616. /*  1   0*  1   1 (Always) */
  617. /*  1   1   0   0 */
  618. /*  1   1   0*  1* (Based on 2F and discriminant distance) */
  619. /*  1   1   1*  0 (If there is an onset there) */
  620. /*  1   1   1   1 */
  622. /*   Determine if there is an onset transition between P and 1F. */
  623. /*   OT (Onset Transition) is true if there is an onset between */
  624. /*   P and 1F but not after 1F. */
  625.     ot = ((obound[1] & 2) != 0 || obound[2] == 1) && (obound[3] & 1) == 0;
  626. /*   Multi-way dispatch on voicing decision history: */
  627.     vstate = (voibuf[3] << 3) + (voibuf[4] << 2) + (voibuf[5] << 1) + voibuf[
  628.     6];
  629.     switch (vstate + 1) {
  630. case 1:  goto L99;
  631. case 2:  goto L1;
  632. case 3:  goto L2;
  633. case 4:  goto L99;
  634. case 5:  goto L4;
  635. case 6:  goto L5;
  636. case 7:  goto L6;
  637. case 8:  goto L7;
  638. case 9:  goto L8;
  639. case 10:  goto L99;
  640. case 11:  goto L10;
  641. case 12:  goto L11;
  642. case 13:  goto L99;
  643. case 14:  goto L13;
  644. case 15:  goto L14;
  645. case 16:  goto L99;
  646.     }
  647. L1:
  648.     if (ot && voibuf[7] == 1) {
  649. voibuf[5] = 1;
  650.     }
  651.     goto L99;
  652. L2:
  653.     if (voibuf[7] == 0 || voice[2] < -voice[3]) {
  654. voibuf[5] = 0;
  655.     } else {
  656. voibuf[6] = 1;
  657.     }
  658.     goto L99;
  659. L4:
  660.     voibuf[4] = 0;
  661.     goto L99;
  662. L5:
  663.     if (voice[1] < -voice[2]) {
  664. voibuf[4] = 0;
  665.     } else {
  666. voibuf[5] = 1;
  667.     }
  668.     goto L99;
  669. /*   VOIBUF(2,0) must be 0 */
  670. L6:
  671.     if (voibuf[1] == 1 || voibuf[7] == 1 || voice[3] > voice[0]) {
  672. voibuf[6] = 1;
  673.     } else {
  674. voibuf[3] = 1;
  675.     }
  676.     goto L99;
  677. L7:
  678.     if (ot) {
  679. voibuf[4] = 0;
  680.     }
  681.     goto L99;
  682. L8:
  683.     if (ot) {
  684. voibuf[4] = 1;
  685.     }
  686.     goto L99;
  687. L10:
  688.     if (voice[2] < -voice[1]) {
  689. voibuf[5] = 0;
  690.     } else {
  691. voibuf[4] = 1;
  692.     }
  693.     goto L99;
  694. L11:
  695.     voibuf[4] = 1;
  696.     goto L99;
  697. L13:
  698.     if (voibuf[7] == 0 && voice[3] < -voice[2]) {
  699. voibuf[6] = 0;
  700.     } else {
  701. voibuf[5] = 1;
  702.     }
  703.     goto L99;
  704. L14:
  705.     if (ot && voibuf[7] == 0) {
  706. voibuf[5] = 0;
  707.     }
  708. /*  GOTO 99 */
  709. L99:
  710. /*   Now update parameters: */
  711. /*   ---------------------- */
  713. /*  During unvoiced half-frames, update the low band and full band unvoice
  714. d*/
  715. /*   energy estimates (LBUE and FBUE) and also the zero crossing */
  716. /*   threshold (DITHER).  (The input to the unvoiced energy filters is */
  717. /*   restricted to be less than 10dB above the previous inputs of the */
  718. /*   filters.) */
  719. /*   During voiced half-frames, update the low-pass (LBVE) and all-pass */
  720. /*   (FBVE) voiced energy estimates. */
  721.     if (voibuf[*half + 6] == 0) {
  722. /* Computing MIN */
  723. i__1 = fbe, i__2 = *ofbue * 3;
  724. r__1 = (*sfbue * 63 + (min(i__1,i__2) << 3)) / 64.f;
  725. *sfbue = i_nint(&r__1);
  726. *fbue = *sfbue / 8;
  727. *ofbue = fbe;
  728. /* Computing MIN */
  729. i__1 = lbe, i__2 = *olbue * 3;
  730. r__1 = (*slbue * 63 + (min(i__1,i__2) << 3)) / 64.f;
  731. *slbue = i_nint(&r__1);
  732. *lbue = *slbue / 8;
  733. *olbue = lbe;
  734.     } else {
  735. r__1 = (*lbve * 63 + lbe) / 64.f;
  736. *lbve = i_nint(&r__1);
  737. r__1 = (*fbve * 63 + fbe) / 64.f;
  738. *fbve = i_nint(&r__1);
  739.     }
  740. /*   Set dither threshold to yield proper zero crossing rates in the */
  741. /*   presence of low frequency noise and low level signal input. */
  742. /*   NOTE: The divisor is a function of REF, the expected energies. */
  743. /* Computing MIN */
  744. /* Computing MAX */
  745.     r__2 = sqrt((real) (*lbue * *lbve)) * 64 / 3000;
  746.     r__1 = max(r__2,1.f);
  747.     *dither = min(r__1,20.f);
  748. /*   Voicing decisions are returned in VOIBUF. */
  749.     return 0;
  750. } /* voicin_ */