home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / SMOOTHDE.C < prev    next >
Encoding:
C/C++ Source or Header  |  1989-04-11  |  1.5 KB  |  70 lines
  1.  
  2. # include <math.h>
  3. # include "smooth.h"
  4. # include <stdio.h>
  5. # include "miscio.h"
  6.  
  7. float xdata[100];
  8. float ydata[100];
  9. float smoothdata[100];
  10. int n;
  11. int i;
  12. int smoothnum;
  13. int derivnum;
  14. float weights[10];
  15. float wdivisor;
  16.  
  17.  
  18. void main()
  19. {
  20.  
  21.    printf("Data Smooth Weights \n" );
  22.    printf( "Smoothnum = ");
  23.    scanf("%d", &smoothnum);
  24.  
  25.    for ( i = 0; i <= smoothnum-1; ++i ) {
  26.       printf("weight # %3d =  ", i);
  27.       scanf("%f", &weights[i]);
  28.       printf("\n");
  29.    }
  30.    printf("Weight divisor =   ");
  31.    scanf("%f", &wdivisor);
  32.    printf("\n");
  33.    n = 100;
  34.    {
  35.       for ( i = 0; i <= n-1; ++i ) {
  36.          xdata[i] = i;
  37.      ydata[i] = cos(i / 10.0) + frandom();
  38.       }
  39.       DataSmoothWeights(ydata,n,smoothnum,weights,wdivisor,smoothdata);
  40.    }
  41.       for ( i = 0; i <= 99; ++i ) {
  42.          printf( "ydata = %7.2f   smoothdata =  %7.2f\n",
  43.           ydata[i],smoothdata[i]);
  44.       }
  45.  
  46.  
  47.    printf("\n Press Return to Continue\n" );
  48.    getch();
  49.    printf("Savitzky Golay Smoothing \n");
  50.    printf( "Smoothnum = ");
  51.    scanf("%d", &smoothnum);
  52.    printf( "Derivnum = ");
  53.    scanf("%d", &derivnum);
  54.  
  55.    n = 100;
  56.    {
  57.       for ( i = 0; i <= n-1; ++i ) {
  58.          xdata[i] = i;
  59.      ydata[i] = cos(i / 10.0) + frandom();
  60.       }
  61.       DataSmoothSg(ydata,n,smoothnum,derivnum,smoothdata);
  62.    }
  63.       for ( i = 0; i <=18; ++i ) {
  64.          printf( "ydata = %7.2f   smoothdata =  %7.2f\n",ydata[i],smoothdata[i]);
  65.       }
  66.    printf("\n Press Return to End\n" );
  67.    getch();
  68. }
  69.  
  70.