simulation/interf2.BAS

   1 DIM SHARED freq(1 TO 100)\r
   2 DIM SHARED phase(1 TO 100)\r
   3 \r
   4 SCREEN 12\r
   5 \r
   6 PSET (0, 0)\r
   7 FOR x = 0 TO 639\r
   8 \r
   9  fs = (x - 320) / 5000 + 1\r
  10 ' fs = (x - 320) / 50 + 1\r
  11 \r
  12   y = 0\r
  13   FOR t = 1 TO 5000 STEP 5\r
  14     y1 = SIN(t)\r
  15     y2 = SIN(t * fs)\r
  16     y = y + ABS(y1 + y2)\r
  17   NEXT t\r
  18 \r
  19   y = y / 5\r
  20   IF y > 470 THEN y = 470\r
  21   IF y < 0 THEN y = 0\r
  22 \r
  23   LINE -(x, 479 - y), 15\r
  24 \r
  25 NEXT x\r
  26 \r
  27 SUB getfreq\r
  28 \r
  29 FOR a = 1 TO 100\r
  30   freq(a) = RND / 7 + 1\r
  31 NEXT a\r
  32 \r
  33 FOR a = 1 TO 100\r
  34   phase(a) = RND * 100\r
  35 NEXT a\r
  36 \r
  37 \r
  38 END SUB\r
  39 \r
  40 SUB gety (t)\r
  41 \r
  42 y = 0\r
  43 FOR a = 1 TO 100\r
  44   y = y + SIN(t * freq(a) + phase(a))\r
  45 NEXT a\r
  46 \r
  47 \r
  48 END SUB\r
  49 \r