[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von MaikG:
Hab schon überlegt ob man ein Sample mit dem Basic Programm
in die Ram Disk legt und dann mit einem C Programm das Sample
auswertet.

---

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von MaikG:

Ich bin zu blöd um Irgendein Source hier umzuwandeln, das ist
zuviel und vorallem Funktionen womit ich noch nie zu tun habe.

---

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von Holger:

Der Beispielcode muss nur 8000 Mal pro Sekunde mit einem Samplewert gefüttert werden, und spuckt die decodierten Zeichen aus. Den kann man umsetzen, ohne mit .fd-Files rumhantieren zu müssen...

---

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von Mad_Dog:
Sollte man nicht immer mit der doppelten Frequenz des Eingangssignals abtasten (so habe ich das jedenfalls gelernt)?
Wenn Du mit 8kHz abtastest, kannst Du maximal eine Frequenz von 4kHz sauber abtasten...

---

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von Holger:

Zitat:

---

Original von Mad_Dog:
Sollte man nicht immer mit der doppelten Frequenz des Eingangssignals abtasten (so habe ich das jedenfalls gelernt)?
Wenn Du mit 8kHz abtastest, kannst Du maximal eine Frequenz von 4kHz sauber abtasten...

---

Zum Erkennen von Signaltönen nimmt man noch ein bisschen mehr mit, um Störungen leichter erkennen zu können. Aber bei einem Maximum von einer Signalkombination von 941 Hz und 1633 Hz reicht das locker aus...

---

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von Mad_Dog:
BTW.: Ich kann mich vage an ein Programm von einer Demo CD erinnern, welches in Echtzeit von dem gesampelten Sound vom Parallelportsampler eine Spektralanalyse gemacht hat (also braucht's wohl doch nicht soviel Rechenpower...). Mal schauen, vielleicht find ich's ja mal wieder. ;-)

---

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von MaikG:
Library aufrufe sind normalerweise nicht schlimm. Kompliziertes
C umzuwandeln schon.

---

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

code:

---

REM aktualisiert 04.09.2006 - 11:53

REM $Nolibrary
REM $Nolines
REM $NoOverflow
REM $NOEVENT

DEFINT A-Z

'$INCLUDE Exec.bh
'$INCLUDE Timer.bh

REM $INCLUDE BLib/ExecSupport.bas

LIBRARY OPEN "exec.library", LIBRARY_MINIMUM&

CONST SECSPERMIN& = 60
CONST SECSPERHOUR& = 3600
CONST SECSPERDAY& = 86400

DIM currentval(timeval_sizeof  2)
DIM SHARED junk&,dummy$

FUNCTION create_timer&(BYVAL unit&)
    STATIC r&, timerport&, timerIO&
    create_timer& = NULL&
    timerport& = CreatePort&(NULL&, 0)
    IF timerport& <> NULL& THEN
        timerIO& = CreateExtIO&(timerport&, timerequest_sizeof)
        IF timerIO& <> NULL& THEN
            r& = OpenDevice&(SADD("timer.device" + CHR$(0)), unit&, timerIO&, 0)
            IF r& = 0 THEN
                create_timer& = timerIO&
            ELSE
                delete_timer timerIO&
            END IF
        ELSE
            DeletePort timerport&
        END IF
    END IF
END FUNCTION

SUB delete_timer(BYVAL tr&)
    STATIC tp&
    IF tr& <> NULL& THEN
        tp& = PEEKL(tr& + tr_node + IORequestio_Message + mn_ReplyPort)
        IF tp& <> 0 THEN DeletePort tp&
        CloseDevice tr&
        DeleteExtIO tr&
    END IF
END SUB

FUNCTION time_delay(BYVAL tv&, BYVAL unit&)
    STATIC tr&
    time_delay = FALSE&
    tr& = create_timer&(unit&)
    IF tr& <> NULL& THEN
        wait_for_timer tr&, tv&
        delete_timer tr&
        time_delay = TRUE&
    END IF
END FUNCTION

SUB wait_for_timer(BYVAL tr&, BYVAL tv&)
    POKEW tr& + tr_node + IORequestio_Command, TR_ADDREQUEST&
    CopyMem tv&, tr& + tr_time, timeval_sizeof    ' structure assignment
    junk& = DoIO&(tr&)
END SUB

CONST SAMPLING_RATE%=8000
CONST MAX_BINS%=8
CONST GOERTZEL_N%=92

DIM SHARED q1#(MAX_BINS%) 'double      q1[ MAX_BINS ];
DIM SHARED q2#(MAX_BINS%) 'double      q2[ MAX_BINS ];
DIM SHARED r#(MAX_BINS%)  'double      r[ MAX_BINS ];
DIM SHARED freqs#(MAX_BINS%)
DIM SHARED sample_count% 'INT         sample_count;
DIM SHARED coefs#(MAX_BINS%)
DIM SHARED row_col_ascii_codes$(4,4)


' double      freqs[ MAX_BINS] = 
'    697,
'    770,
'    852,
'    941,
'    1209,
'   1336,
'    1477,
'    1633
'};

freqs#(0)=697
freqs#(1)=770
freqs#(2)=852
freqs#(3)=941
freqs#(4)=1209
freqs#(5)=1336
freqs#(6)=1477
freqs#(7)=1633


row_col_ascii_codes$(0,0)="1"
row_col_ascii_codes$(0,1)="2"
row_col_ascii_codes$(0,2)="3"
row_col_ascii_codes$(0,3)="A"
row_col_ascii_codes$(1,0)="4"
row_col_ascii_codes$(1,1)="5"
row_col_ascii_codes$(1,2)="6"
row_col_ascii_codes$(1,3)="B"
row_col_ascii_codes$(2,0)="7"
row_col_ascii_codes$(2,1)="8"
row_col_ascii_codes$(2,2)="9"
row_col_ascii_codes$(2,3)="C"
row_col_ascii_codes$(3,0)="*"
row_col_ascii_codes$(3,1)="0"
row_col_ascii_codes$(3,2)="#"
row_col_ascii_codes$(3,3)="D"


'double      coefs[ MAX_BINS ] ;
' *  calc_coeffs

'void
'calc_coeffs()
'{
'int n;
'
'  for(n = 0; n < MAX_BINS; n++)
'  {
'    coefs[n] = 2.0 * cos(2.0 * 3.141592654 * freqs[n] / SAMPLING_RATE);
'  }
'
'}

REM SUB calc_coeffs
 FOR n%=0 TO 7 REM MAX_BINS%-1
  coefs#(n%)=2.0*COS(2.0*3.141592654 * freqs#(n%)/SAMPLING_RATE%)
 NEXT n%
REM END SUB

' *  post_testing
' * This is where we look at the bins and decide if we have a valid signal.

'void
'post_testing()

'{
'int         row, col, see_digit;
'int         peak_count, max_index;
'double      maxval, t;
'int         i;

SUB post_testing
STATIC row%, col%, see_digit%
STATIC peak_count%, max_index%, i%
STATIC maxval#,t#


 'char *  row_col_ascii_codes[4][4] = {
 '       {"1", "2", "3", "A"},
 '       {"4", "5", "6", "B"},
 '       {"7", "8", "9", "C"},
 '       {"*", "0", "#", "D"}};


'    /* Find the largest in the row group. */
'    row = 0;
'     maxval = 0.0;
'     for ( i=0; i<4; i++ )
'     {
'         if ( r[i] > maxval )
'         {
'             maxval = r[i];
'             row = i;
'         }
'     }

row%=0
maxval#=0.0
FOR i%=0 TO 3
 IF r#(i%)>maxval# THEN maxval#=r#(i%):row%=i%
NEXT i%



'     /* Find the largest in the column group. */
'     col = 4;
'     maxval = 0.0;
'     for ( i=4; i<8; i++ )
'     {
'         if ( r[i] > maxval )
'         {
'             maxval = r[i];
'             col = i;
'         }
'     }

col%=4
maxval#=0.0
FOR i%=4 TO 7
 IF r#(i%)>maxval# THEN maxval#=r#(i%):col%=i%
NEXT i%



'   /* Check FOR minimum energy */
'   if ( r[row] < 4.0e5 )   /* 2.0e5 ... 1.0e8 no change */
'    {
'        /* energy not high enough */
'    }
'    else if ( r[col] < 4.0e5 )
'    {
'        /* energy not high enough */
'    }
'    else
'    {
'        see_digit = TRUE;


    IF r#(row%)<4*10^5 THEN
     dummy$="energy not high enough"
     ELSEIF r#(col%)<4*10^5 THEN dummy$="energy NOT high enough"
     ELSE see_digit% = TRUE&
    END IF

'        if ( r[col] > r[row] )
'        {
'            /* Normal twist */
'            max_index = col;
'            if ( r[row] < (r[col] * 0.398) )    /* twist > 4dB, error */
'                see_digit = FALSE;
'        }
'        else /* if ( r[row] > r[col] ) */
'        {
'            /* Reverse twist */
'            max_index = row;
'            if ( r[col] < (r[row] * 0.158) )    /* twist > 8db, error */
'                see_digit = FALSE;
'        }

IF r#(col%)>r#(row%) THEN
 max_index%=col%
 IF r#(row%)<(r#(col%)*0.398) THEN see_digit% = FALSE&
 ELSE
 max_index%=row%
 IF r#(col%)<(r#(row%)*0.158) THEN see_digit% = FALSE&
END IF

'        if ( r[max_index] > 1.0e9 )
'            t = r[max_index] * 0.158;
'        else
'            t = r[max_index] * 0.010;
'
'        peak_count = 0;
'        for ( i=0; i<8; i++ )
'        {
'            if ( r[i] > t )
'                peak_count++;
'        }
'        if ( peak_count > 2 )
'            see_digit = FALSE;
'
'        if ( see_digit )
'        {
'            printf( "%s", row_col_ascii_codes[row][col-4] );
'            fflush(stdout);
'        }
'    }
'}

IF r#(max_index%)>1*10^9 THEN t#=r#(max_index%)*0.158 ELSE t#=r#(max_index%)*0.010

peak_count%=0

REM FOR i%=0 TO 7
REM IF r#(i%)>t# THEN INCR peak_count%
REM IF peak_count%>2 THEN see_digit%=FALSE&
REM IF see_digit% THEN LOCATE 2,2:PRINT row_col_ascii_codes$(row%,col%-4)
REM NEXT i%
 
FOR i%=0 TO 7
 IF r#(i%)>t# THEN INCR peak_count%
NEXT i%
 IF peak_count%>2 THEN see_digit%=FALSE&
 IF see_digit% THEN LOCATE 2,2:PRINT row_col_ascii_codes$(row%,col%-4)





' *  goertzel
'void
'goertzel( int sample )
'{
'double      q0;
'ui32        i;
'
'    if ( sample_count < GOERTZEL_N )
'    {
'        sample_count++;
'        for ( i=0; i<MAX_BINS; i++ )
'        {
'            q0 = coefs[i] * q1[i] - q2[i] + sample;
'            q2[i] = q1[i];
'            q1[i] = q0;
'        }
'    }
'    else
'    {
'       for ( i=0; i<MAX_BINS; i++ )
'        {
'            r[i] = (q1[i] * q1[i]) + (q2[i] * q2[i]) - (coefs[i] * q1[i] * q2[i]);
'            q1[i] = 0.0;
'            q2[i] = 0.0;
'        }
'        post_testing();
'        sample_count = 0;
'    }
'}

END SUB


FUNCTION goertzek%(sample%)
STATIC q0#, i%
 IF sample_count%<GOERTZEL_N% THEN
  INCR sample_count%
  FOR i%=0 TO 7 REM MAX_BINS%-1
   q0#=coefs#(i%)*q1#(i%)-q2#(i%)+sample%
   q2#(i%)=q1#(i%)
   q1#(i%)=q0#
  NEXT i%
 ELSE
  FOR i%=0 TO 7 REM MAX_BINS%-1
   r#(i%)=(q1#(i%)*q1#(i%)) + (q2#(i%)*q2#(i%)) - (coefs#(i%)*q1#(i%)*q2#(i%))
   q1#(i%)=0.0
   q2#(i%)=0.0
  NEXT i%
  post_testing
  sample_count%=0
 END IF
END FUNCTION


POKEB &hBFE301,0

1
IF PEEKB(&hBFE101)>180 THEN
 PRINT "jetzt"
 POKEL VARPTR(currentval(0)) + tv_secs, 0
 POKEL VARPTR(currentval(0)) + tv_micro, 000125 REM 500000=0.5s 000125=8000HZ
 2
junk = time_delay(VARPTR(currentval(0)), UNIT_MICROHZ&)
INCR mycount%:IF mycount%=8000 THEN LOCATE 1,1:PRINT TIME$:mycount%=0 
sa%=PEEKB(&hBFE101)
IF sa%>127 THEN
 sa%=sa%*128
 ELSEIF sa%<127 THEN sa%=sa%*(-128)
 ELSE sa%=0
 END IF
 bla%=goertzek%(sa%)
REM bla%=goertzek%(PEEKB(&hBFE101))
 IF INKEY$=CHR$(27) THEN 5
 GOTO 2
END IF
 IF INKEY$=CHR$(27) THEN 5
GOTO 1

5 END

---

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

1. Die Variablen die oben bei C definiert werden sind für alle Unterprogramme zugänglich?

---

Zitat:

---

2. In freqs werden nur 8 Werte abgelegt bei Basic geht fängt ein Array von 0 an, ist es jetzt richtig
die werte von 1-8 zu schreiben oder sollten diese von 0-7 gehen?

---

Zitat:

---

3. Bei FOR ist das ähnlich es geht von 0 bis 8 ist das korrekt oder müsste es von 0 bis 8-1 heissen?

---

Zitat:

---

4. z.B. 4.0e5 entspricht das 4.0^5 ?

---

Zitat:

---

5. Wie gebe ich da überhaupt meinen Sample hinein? sample% kann ja kein Zeiger auf ein Speicherbereich sein.

---

Nein, ein einzelne Sample-Wert. Deshalb schrieb ich ja, Du musst die Funktion 8000 Mal pro Sekunde aufrufen (z.B. per Timer). Damit kannst Du zwischen den Aufrufen auch andere Dinge machen...

Übrigens gehört letztere Frage zu denen, die ich am Anfang auch gegrübelt habe. Der code ist wirklich unzureichen dokumentiert, aber dafür bist Du anscheinend schon weit gekommen, was meine These unterstützt, das er doch nicht so kompliziert ist Wenn Du noch weitere Fragen hast, immer her damit.

mfg
--
Good coders do not comment. What was hard to write should be hard to read too.

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

So, ich hab jetzt das erwähnte Programm (bzw. mehrere Programme) gefunden, welche in Echtzeit FTT machen. Die Teile Stammen von der Demo-Gruppe Infect und waren auf der CD "Demos are forever" mit drauf.

Ich hab die Teile mal in ein Archiv gepackt und auf meinen Webspace geschmissen: http://w3.norman-interactive.com/files/Infect_Scope.lha

Besonders interssant ist in diesem Zusammenhang das Programm scopext2.
Es macht in Echtzeit eine Spektralanalyse der Daten vom Parallelportsampler. Mit meinem "Stereo Master" funktioniert das prima - dürfte vermutlich auch mit anderen Samplern funktionieren.

Doch vorsicht: Die Programme öffnen einen Pal-Screen. Und beim Beenden schmieren sie leider (zumindest bei mir) ab.

Wie genau die Spektralanalyse hier ist, kann ich nicht sagen.
Ist auch leider kein Sourcecode dabei, aber dafür ausführliche Textdateien.
--
http://www.norman-interactive.com

[ Dieser Beitrag wurde von Mad_Dog am 02.09.2006 um 18:45 Uhr geändert. ]

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von Mad_Dog:
Besonders interssant ist in diesem Zusammenhang das Programm scopext2.
Es macht in Echtzeit eine Spektralanalyse der Daten vom Parallelportsampler. Mit meinem "Stereo Master" funktioniert das prima - dürfte vermutlich auch mit anderen Samplern funktionieren.

Doch vorsicht: Die Programme öffnen einen Pal-Screen. Und beim Beenden schmieren sie leider (zumindest bei mir) ab.

Wie genau die Spektralanalyse hier ist, kann ich nicht sagen.
Ist auch leider kein Sourcecode dabei, aber dafür ausführliche Textdateien.

---

Scopext4 habe ich hier, incl. Sourcecode! Author: Smack/Infect! (Michael Henke). Da hier auch ein MuForce Hit beim beenden auftritt, hatt eich Michael Henke einen Bugreport geschickt. Leider (so weit ich mich noch erinnern kann) hat er keine Hardware oder Interesse.

Wer den Sourcecode haben will soll eben bescheid geben. Wäre auch schön, wenn jemand dann gleich den Bug beheben könnte.
Scopext4 rulet einfach nur. Hier mit dem TechnoSoundTurbo2 Sampler.

Sorry for being offtopic.

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von MaikG:
Was sind eigentlich 8000 HZ in Microsekunden?

---

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von MaikG:
>In C fangen Indizes immer mit 0 an, wenn die Größe 8 ist, gehen die
>gültigen Indizes von 0-7.

#define MAX_BINS 8

double q1[ MAX_BINS ];


Heiss doch aber dim q1#(8), also hat man ein 9er Array erstellt
q1#(0) bis q1#(8) besitzen gültigkeit. Also 1 ist da irgendwie
überflüssig.

---

Zitat:

---

Also ich soll 8000 Mal(ob das Basic schafft?) einen Integer
wert übergeben also 0-255?

---

Zitat:

---

Wie soll man aus einer einzelnen Amplitude dann irgendeine Frequenz rausbekommen oder wird die irgendwo zwischengespeichert?

---

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von MaikG:
IF sa%>127 THEN
sa%=sa%*128
ELSEIF sa%<127 THEN sa%=sa%*(-128)
ELSE sa%=0
END IF

---

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von MaikG:
sa% ist Integer, lege ich PEEKB(xyz) da rein ist sa% eine
zahl von 0-255, 127 ist dabei die 0.
Die Konventerierung erzeugt eine reguläre 0, und das Vorzeichen
und 16 Bit werte in akteptabler Amplitude.

---

[ - Antworten - Zitieren - Direktlink - ]

Zitat:

---

Original von MaikG:
Was sind eigentlich 8000 HZ in Microsekunden?

---

Im Physikunterricht gepennt?

1 Hz = 1/s

T = 1/n (T = Periodendauer, n = Frequenz)

1 Hz bedeutet "ein mal pro Sekunde".
D.h. wenn Du was 8000 mal pro Sekunde machen willst und dich fragst, wie lange "ein Durchgang" dauert, dann benutzt Du die Formel

T = 1/n

hier:

T = 1/8000
T = 0,000125 s

Siehe auch http://de.wikipedia.org/wiki/Frequenz

Zum Thema "milli, mikro, nano usw." : http://de.wikipedia.org/wiki/Vors%C3%A4tze_f%C3%BCr_Ma%C3%9Feinheiten

1 Millisekunde = 1 * 10 ^ -3 s
1 Mikrosekunde = 1 * 10 ^ -6 s

Also

0,000125 s
= 125 * 10 ^ -6 s
= 125 Mikrosekunden


--
http://www.norman-interactive.com

[ Dieser Beitrag wurde von Mad_Dog am 04.09.2006 um 11:37 Uhr geändert. ]

[ - Antworten - Zitieren - Direktlink - ]

[ - Antworten - Zitieren - Direktlink - ]