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

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

Linux/Unix编程

开发平台：

Visual C++

mload.c：源码内容
  1. /*
  3. $Log: mload.c,v $
  4. Revision 1.1.1.1  2002/07/28 05:23:36  freeman_yong
  5. lpc10 codec
  7.  * Revision 1.1  1996/08/19  22:31:25  jaf
  8.  * Initial revision
  9.  *
  11. */
  13. #ifdef P_R_O_T_O_T_Y_P_E_S
  14. extern int mload_(integer *order, integer *awins, integer *awinf, real *speech, real *phi, real *psi);
  15. #endif
  17. /*  -- translated by f2c (version 19951025).
  18.    You must link the resulting object file with the libraries:
  19. -lf2c -lm   (in that order)
  20. */
  22. #include "f2c.h"
  24. /* ***************************************************************** */
  26. /*  MLOAD Version 48 */
  28. /* $Log: mload.c,v $
  29. /* Revision 1.1.1.1  2002/07/28 05:23:36  freeman_yong
  30. /* lpc10 codec
  31. /*
  32.  * Revision 1.1  1996/08/19  22:31:25  jaf
  33.  * Initial revision
  34.  * */
  35. /* Revision 1.5  1996/03/27  23:59:51  jaf */
  36. /* Added some more accurate comments about which indices of the argument */
  37. /* array SPEECH are read.  I thought that this might be the cause of a */
  38. /* problem I've been having, but it isn't. */
  40. /* Revision 1.4  1996/03/26  19:16:53  jaf */
  41. /* Commented out the code at the end that copied the lower triangular */
  42. /* half of PHI into the upper triangular half (making the resulting */
  43. /* matrix symmetric).  The upper triangular half was never used by later */
  44. /* code in subroutine ANALYS. */
  46. /* Revision 1.3  1996/03/18  21:16:00  jaf */
  47. /* Just added a few comments about which array indices of the arguments */
  48. /* are used, and mentioning that this subroutine has no local state. */
  50. /* Revision 1.2  1996/03/13  16:47:41  jaf */
  51. /* Comments added explaining that none of the local variables of this */
  52. /* subroutine need to be saved from one invocation to the next. */
  54. /* Revision 1.1  1996/02/07 14:48:01  jaf */
  55. /* Initial revision */
  58. /* ***************************************************************** */
  60. /* Load a covariance matrix. */
  62. /* Input: */
  63. /*  ORDER            - Analysis order */
  64. /*  AWINS            - Analysis window start */
  65. /*  AWINF            - Analysis window finish */
  66. /*  SPEECH(AWINF)    - Speech buffer */
  67. /*                     Indices MIN(AWINS, AWINF-(ORDER-1)) through */
  68. /*                             MAX(AWINF, AWINS+(ORDER-1)) read. */
  69. /*                     As long as (AWINF-AWINS) .GE. (ORDER-1), */
  70. /*                     this is just indices AWINS through AWINF. */
  71. /* Output: */
  72. /*  PHI(ORDER,ORDER) - Covariance matrix */
  73. /*                    Lower triangular half and diagonal written, and read.*/
  74. /*                     Upper triangular half untouched. */
  75. /*  PSI(ORDER)       - Prediction vector */
  76. /*                     Indices 1 through ORDER written, */
  77. /*                     and most are read after that. */
  79. /* This subroutine has no local state. */
  81. /* Subroutine */ int mload_(integer *order, integer *awins, integer *awinf, 
  82. real *speech, real *phi, real *psi)
  83. {
  84.     /* System generated locals */
  85.     integer phi_dim1, phi_offset, i__1, i__2;
  87.     /* Local variables */
  88.     integer c__, i__, r__, start;
  90. /*       Arguments */
  91. /*       Local variables that need not be saved */
  92. /*   Load first column of triangular covariance matrix PHI */
  93.     /* Parameter adjustments */
  94.     --psi;
  95.     phi_dim1 = *order;
  96.     phi_offset = phi_dim1 + 1;
  97.     phi -= phi_offset;
  98.     --speech;
  100.     /* Function Body */
  101.     start = *awins + *order;
  102.     i__1 = *order;
  103.     for (r__ = 1; r__ <= i__1; ++r__) {
  104. phi[r__ + phi_dim1] = 0.f;
  105. i__2 = *awinf;
  106. for (i__ = start; i__ <= i__2; ++i__) {
  107.     phi[r__ + phi_dim1] += speech[i__ - 1] * speech[i__ - r__];
  108. }
  109.     }
  110. /*   Load last element of vector PSI */
  111.     psi[*order] = 0.f;
  112.     i__1 = *awinf;
  113.     for (i__ = start; i__ <= i__1; ++i__) {
  114. psi[*order] += speech[i__] * speech[i__ - *order];
  115.     }
  116. /*   End correct to get additional columns of PHI */
  117.     i__1 = *order;
  118.     for (r__ = 2; r__ <= i__1; ++r__) {
  119. i__2 = r__;
  120. for (c__ = 2; c__ <= i__2; ++c__) {
  121.     phi[r__ + c__ * phi_dim1] = phi[r__ - 1 + (c__ - 1) * phi_dim1] - 
  122.     speech[*awinf + 1 - r__] * speech[*awinf + 1 - c__] + 
  123.     speech[start - r__] * speech[start - c__];
  124. }
  125.     }
  126. /*   End correct to get additional elements of PSI */
  127.     i__1 = *order - 1;
  128.     for (c__ = 1; c__ <= i__1; ++c__) {
  129. psi[c__] = phi[c__ + 1 + phi_dim1] - speech[start - 1] * speech[start 
  130. - 1 - c__] + speech[*awinf] * speech[*awinf - c__];
  131.     }
  132. /*   Copy lower triangular section into upper (why bother?) */
  133. /*       I'm commenting this out, since the upper triangular half of PHI 
  134. */
  135. /*       is never used by later code, unless a sufficiently high level of 
  136. */
  137. /*       tracing is turned on. */
  138. /*  DO R = 1,ORDER */
  139. /*     DO C = 1,R-1 */
  140. /*        PHI(C,R) = PHI(R,C) */
  141. /*     END DO */
  142. /*  END DO */
  143.     return 0;
  144. } /* mload_ */