Télécharger multirec.procedur

Retour à la liste

Numérotation des lignes :

  1. * MULTIREC PROCEDUR PICA 92/07/03 21:27:30 619
  2. *----------------------------------------------------------------------*
  3. * *
  4. * PROCEDURE MULTIREC *
  5. * *
  6. * CALCUL DE LA RECOMPOSITION D'UN SIGNAL A PARTIR D' UNE SEQUENCE DE *
  7. * FILTRE ONDELETTE ORTHOGONALE *
  8. * *
  9. *----------------------------------------------------------------------*
  10. * *
  11. * ENTREE *
  12. * *
  13. * COARSE X/Y CONTENANT LE SIGNAL DE BASE (COARSE SIGNAL). DX DOIT ETRE *
  14. * CONSTANT DANS LE PROCESSUS DE DECOMPOSITION DIADIC. *
  15. * *
  16. * ONDELL X/Ys CONTENANT LA DECOMPOSITION EN ONDELETTE DANS L'ORDRE DES *
  17. * BANDES DE FREQUENCE CROISSANTE *
  18. * *
  19. * H Y CONTENANT LE FILTRE (NON FORCEMENT SYMETRIQUE). *
  20. * MMH POINTS NEGATIFS (MMH<0--->SYMETRIQUE) *
  21. * *
  22. * G Y CONTENANT LE FILTRE ORTHOGONAL (NON FORCEMENT SYMETRIQUE). *
  23. * MMG POINTS NEGATIFS (MMG<0--->SYMETRIQUE) *
  24. * *
  25. * *
  26. * OPTIO MOT CONTENANT LE MODE DE CORRELATION (DEFAUT='SYME') *
  27. * *
  28. * N NOMBRE DE NIVEAU DE RECOMPOSITION DESIRE *
  29. * *
  30. *----------------------------------------------------------------------*
  31. * *
  32. * SORTIE *
  33. * *
  34. * M NOMBRE DE NIVEAU EFFECTIVEMENT ATTEINT *
  35. * *
  36. * SIGNAL SIGNAL RECOMPOSE *
  37. *----------------------------------------------------------------------*
  38. 'DEBPROC' MULTIREC COARSE*'EVOLUTION' ONDELL*'EVOLUTION'
  39. H*'LISTREEL' MMH*'ENTIER'
  40. G*'LISTREEL' MMG*'ENTIER'
  41. OPTIO/'MOT' N/'ENTIER';
  42. *----------------------------------------------------------------------*
  43. YPHI='COPIER' ('EXTRAIRE' COARSE 'ORDO' 1);
  44. *
  45. NPT='DIMENSION' YPHI;
  46. NPTT='DIMENSION' ('EXTRAIRE' ONDELL 'ABSC' 1);
  47. 'SI' (NPT 'NEG' NPTT);
  48. 'MESSAGE' '*** nb de point COARSE neg nb de point ONDELL(1)';
  49. 'MESSAGE' '--> pas de calcul';
  50. M=0;
  51. 'SINON';
  52. M='DIMENSION' ONDELL;
  53. 'FINSI';
  54. *
  55. 'SI' (MMH '<' 0); MMHT='TEXT' ' ';
  56. MMH=('DIMENSION' H) - 1;
  57. MPH= MMH;
  58. 'SINON'; MMHT='TEXT' NPNE MMH;
  59. MPH=('DIMENSION' H) - (MMH+1);
  60. 'FINSI';
  61. *
  62. 'SI' (MMG '<' 0); MMGT='TEXT' ' ';
  63. MMG=('DIMENSION' G) - 1;
  64. MPG=MMG;
  65. 'SINON'; MMGT='TEXT' NPNE MMG;
  66. MPG=('DIMENSION' G) - (MMG+1);
  67. 'FINSI';
  68. *
  69. NNPT=2*NPT;
  70. 'SI' ((MMG '>' NNPT) 'OU' (MPG '>' NNPT));
  71. 'MESSAGE' '*** nb de point insuffisant pour Phi au niveau 1';
  72. 'MESSAGE' '--> pas de calcul';
  73. M=0;
  74. 'FINSI';
  75. 'SI' ((MMH '>' NNPT) 'OU' (MPH '>' NNPT));
  76. 'MESSAGE' '*** nb de point insuffisant pour Psi au niveau 1';
  77. 'MESSAGE' '--> pas de calcul';
  78. M=0;
  79. 'FINSI';
  80. *
  81. 'SI' (M 'NEG' 0);
  82. 'SI' (('NON' ('EXISTE' OPTIO)));
  83. OPTIO='TEXT' SYME; 'FINSI';
  84. 'SI' (('NON' ('EXISTE' N))); N=1000; 'FINSI';
  85. 'SI' (M > N); M=N; 'FINSI';
  86. *
  87. I=0;
  88. 'REPETER' LAB1 M;
  89. *
  90. I=I+1;
  91. YPSI='EXTRAIRE' ONDELL 'ORDO' I;
  92. *
  93. YYPSI='DIADIQUE' YPSI 'INVE';
  94. YYPHI='DIADIQUE' YPHI 'INVE';
  95. YPSI ='CVOLUTION' YYPSI G MMGT OPTIO;
  96. YPHI ='CVOLUTION' YYPHI H MMHT OPTIO;
  97. YYPHI=2 * (YPSI + YPHI);
  98. YPHI=YYPHI;
  99. *
  100. 'FIN' LAB1;
  101. 'FINSI';
  102. *
  103. XXXX='EXTRAIRE' ONDELL 'ABSC' M;
  104. XX ='DIADIQUE' XXXX 'IVIN';
  105. 'TITRE' 'Signal recompose sur ' M ' niveaux';
  106. SIGNAL='EVOLUTION' 'MANUEL' 'temps' XX 'Y' YPHI;
  107. 'FINPROC' M SIGNAL;
  108.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales