n1=1000
n2=2319
n3=44100
 
n1=2*2*2*3*5*7
n2=2*2*3*5*7
n3=2*2*3*3*7*41
 
a$=""
diviseur=2
 
print "recherche pgcd de :",n1,n2,n3
 
 
 
fin=0
tour=0
 
while fin=0
 
    if n1 mod diviseur =0 and n2 mod diviseur =0 and n3 mod diviseur =0 then
        'si la division des trois nombres par "diviseur" donnent trois entiers
        'alors "diviseur" est un facteur commun à ces nombres
        'on l'ajoute à a$
        a$=a$+str$(diviseur)+"*"
        'et on calcule les nombres pour le tour suivant
        n1=n1/diviseur
        n2=n2/diviseur
        n3=n3/diviseur
        print diviseur;" = facteur commun",n1,n2,n3
    else
        'sinon on ajoute 1 à divis pour le calcul suivant (diviseur n'a pas a être premier, puisqu'il est de + en + grand
        'print diviseur;" = pas commun"
        diviseur=diviseur+1
    end if
 
    'sortie de la boucle quand l'un au moins des nombres est devenu inférieur au diviseur
    'ou à partir d'un diviseur dont il est vraiscemblable qu'on ira jamais jusque là
    tour=tour+1
    if tour=100000 or (n1<diviseur or n2<diviseur or n3<diviseur)then fin=1
 
wend
 
print
print left$(a$,len(a$)-1)
 
 

 
'Géné BF 2 fréquences
'22 mars 2016 - 12h34

'ajout affichage forme d'onde : à revoir pour l'affichage du dernier point
'réglage d'amplitude manuel et automatique, balance gen1/gen2 : restent à créer

 
    'nomainwin

    NomWav$="pardefaut.wav"
    CheminWav$=DefaultDir$
 
    'dimensions de la fenêtre
    WindowWidth = 615
    WindowHeight = 495
 
    'formats d'affichage
    fonteGlobale$="font ms_sans_serif 0 16"
    fonteTitre$="!font ms_sans_serif 0 16 bold"
 
    'valeurs par défaut
    pi=3.1415926536
    freq1=880
    vol1=70
    freq2=440
    vol2=70
 
    duree=1
 
    'fréquences échantillonnage, Hz
    indexFech=1 'élément qui sera affiché dans la combobox
    FreqEch$(1)="11025"
    FreqEch$(2)="22050"
    FreqEch$(3)="44100"
 
    'résolution échantillonnage, bits
    indexRes=1
    WaveRes$(1)="8"
    WaveRes$(2)="16"
 
    'nombre de voies
    indexVoies=1
    WaveNbvoies$(1)="1-Mono"
    WaveNbvoies$(2)="2-Stereo"
 
    'création de l'interface utilisateur
    statictext #1.titreGen, "Générateurs :", 15, 15, 100, 20
 
    statictext #1.frequences, "Fréquence, Hz", 85, 40, 100, 20
    statictext #1.volumes, "Volume, 0 à 100", 195, 40, 100, 20
 
    statictext #1.gen1, "Gen. 1", 30, 65, 40, 20
    textbox #1.freq1, 80, 60, 100, 25
    textbox #1.vol1, 190, 60, 100, 25
 
    statictext #1.gen2, "Gen. 2", 30, 95, 40, 20
    textbox #1.freq2, 80, 90, 100, 25
    textbox #1.vol2, 190, 90, 100, 25
 
    button #1.genWav, "Générer", [genere_wav], UL, 300, 60, 80, 55
 
    largeurvisu=200
    hauteurvisu=120
    dim affichage(largeurvisu)
    graphicbox #1.visu, 390, 20, largeurvisu, hauteurvisu
 
    statictext #1.sortie, "Fichier de sortie :", 15, 140, 125, 20
    statictext #1.nomWav, NomWav$, 140, 140, 390, 20
    button #1.changeDest, "Emplacement et nom du fichier", [changeDestination], UL, 30, 165, 350, 25
    statictext #1.cheminWav, CheminWav$, 30, 195, 300, 40
 
    statictext #1.param, "Paramètres du fichier Wave :", 15, 245, 200, 20
 
    statictext #1.fech, "Freq. échantillonnage, Hz", 30, 270, 170, 20
    combobox #1.comboEch, FreqEch$(), [f_echantillonnage], 210, 268, 75, 20
 
    statictext #1.fech, "Résolution, bits", 30, 300, 170, 20
    combobox #1.comboRes, WaveRes$(), [Wav_resolution], 210, 297, 75, 20
 
    statictext #1.fech, "Nombre de voies", 30, 330, 170, 20
    combobox #1.comboVoi, WaveNbvoies$(), [Wav_nombre_voies], 210, 328, 75, 20
 
    statictext #1.dureeWav, "Durée, s", 30, 360, 170, 20
    textbox #1.bxDureeWav, 210, 358, 75, 25
    'button #1.okDuree, "OK", [valide_duree], UL, 290, 360, 25, 25

    statictext #1.avanti, "En avant !", 15, 400, 270, 20
 
    button #1.annulCalc, "Annuler", [annule_calcul], UL, 120, 425, 80, 25
    button #1.quit, "Quitter", [quitter], UL, 300, 425, 80, 25
 
    'ouverture fenêtre d'interface utilisateur
    open "Géné. BF" for window as #1
 
        'fonte global
        print #1, fonteGlobale$
 
        'fonte titres
        print #1.titreGen, fonteTitre$
        print #1.sortie, fonteTitre$
        print #1.param, fonteTitre$
        print #1.avanti, fonteTitre$
        print #1.titreGen, fonteTitre$
 
        'données pré-réglées
        print #1.freq1, freq1
        print #1.vol1, vol1
        print #1.freq2, freq2
        print #1.vol2, vol2
 
        print #1.visu, "down"
        print #1.visu, "fill darkgray"
        print #1.visu, "color lightgray"
        print #1.visu, "line 0 ";hauteurvisu/2;" ";largeurvisu;" ";hauteurvisu/2
        print #1.visu, "UP"
 
        print #1.comboEch,"selectindex ";indexFech
        fEchant=val(FreqEch$(indexFech))
 
        print #1.comboRes,"selectindex ";indexRes
        resolution=val(WaveRes$(indexRes))
 
        print #1.comboVoi,"selectindex ";indexVoies
        voies=val(WaveNbvoies$(indexVoies))
 
        print #1.bxDureeWav, duree
        print #1.annulCalc, "!disable"
        'peut-être une bonne idée de lancer d'ici une visu avec les paramètres
        'actuels
        print #1.freq1, "!setfocus"
        print #1, "trapclose [quitter]"
 
    wait
 
end
'***************sous-programmes********************

[quitter]
    'quitte le programme
    close #1
    end
 
[annule_calcul]
    'pas encore créé

    'on annule la création en cours (temps trop long, ou erreur de saisie
    'et vouloir relancer avant la fin réelle du traitement, ...)
    wait
 
[changeDestination]
    print "[changeDestination]"
    filedialog "Save to, entrer un nom pour le fichier :","*.wav",chemin$
    if chemin$<>"" then
        numCar=0
        antislash=0
 
        while antislash=0
            a$=mid$(chemin$,len(chemin$)-numCar,1)
            if a$="\" or numCar>len(chemin$) then antislash=1
            numCar=numCar+1
        wend
 
        NomWav$=right$(chemin$,numCar-1)
        CheminWav$=left$(chemin$,len(chemin$)-numCar+1)
 
        print #1.cheminWav, CheminWav$
        print #1.nomWav, NomWav$
 
    end if
    wait
 
 
[f_echantillonnage]
    print "[f_echantillonnage]"
 
    'lecture valeur combobox
    print #1.comboEch,"selectionindex? indexFech"
    fEchant=val(FreqEch$(indexFech))
 
    wait
 
 
[genere_wav]
    print "[genere_wav]"
    audio$="" 'la chaine qui contiendra les données audio
    'récupération fréquences et volumes des générateurs
    print #1.freq1, "!contents? a$":freq1=val(a$)
    print #1.vol1, "!contents? a$":vol1=val(a$)
    print #1.freq2, "!contents? a$":freq2=val(a$)
    print #1.vol2, "!contents? a$":vol2=val(a$)
 
    'calcul pulsations
    omega1=2*pi*freq1
    omega2=2*pi*freq2
 
    'récupération paramètres wav
    print #1.bxDureeWav, "!contents? a$":duree=val(a$)
 
    'déterminer la période décrivant complètement la somme des deux fréquences
    'calcul du nombre d'échantillons audio de la période
    tEchant=1/fEchant
    nbrEchant=fEchant/pgcd(freq1,freq2,fEchant)
 
    'et du nombre de fois où elle sera répétée pour obtenir la durée souhaitée
    repetitions=duree/(nbrEchant*tEchant)
    if duree mod nbrEchant*tEchant <> 0 then repetitions=int(repetitions)+1
    print "il y aura ";repetitions;" fois la période de la somme"
 
    'visualiser la période
    gosub [visu]
 
    'en fonction de la résolution et du nombre de voies on détermine le codage
    'pour la durée d'une "période somme"
    periodesomme$=""    'la chaine qui contiendra les échantillons de la période
    audio$=""           'la chaine qui contiendra les périodes avant transfert sur disque
    amplitudemax=400
    t=0
    if resolution=8 then
        if voies=1 then gosub [fichier_8bits_mono]
        if voies=2 then gosub [fichier_8bits_stereo]
    end if
 
    if resolution=16 then
        if voies=1 then gosub [fichier_16bits_mono]
        if voies=2 then gosub [fichier_16bits_stereo]
    end if
 
    if horslimite=0 then
        'on prépare les données pour la durée finale souhaitée
        for a=1 to repetitions
            audio$=audio$+periodesomme$
            'audio$=audio$+left$(periodesomme$,len(periodesomme$)-1)
        next
 
        'On prépare l'en-tête du fichier
        longueurAudio=len(audio$)
        nbOctet=resolution/8    'nombre d'octets par échantillon pour une voie
        fluxDonnees=nbOctet*voies*fEchant'débit en nombre d'octets par seconde

        'données du programme
        et$="RIFF"                              'toujours RIFF
        et$=et$+littleEndian$(36+longueurAudio,4)'taille fichier (devrait être 44+nbre d'octets audio ???)
        et$=et$+"WAVE"                          'toujours WAVE
        et$=et$+"fmt "                          'toujours "fmt "
        et$=et$+littleEndian$(16,4)              '16, 18 ou 40
        et$=et$+littleEndian$(1,2)               '1=PCM, Pulse Code Modulation, classique pour un wave ordinaire
        et$=et$+littleEndian$(voies,2)           'nombre de canaux interlacés (1=mono, 2=stéréo,...)
        et$=et$+littleEndian$(fEchant,4)         'fréquence échantillonnage
        et$=et$+littleEndian$(fluxDonnees,4)     'nombre moyen d'octets par seconde
        et$=et$+littleEndian$(nbOctet,2)         'taille du bloc de donnees
        et$=et$+littleEndian$(resolution,2)      'bits par échantillon (pour une voie)
        et$=et$+"data"                          'toujours "data", annonce les données audio
        et$=et$+littleEndian$(longueurAudio,4)   'longueur en octets des données audio

        'le futur fichier audio est complètement contenu dans la chaine wav$
        wav$=et$+audio$
 
        'on transfère wav$ dans le fichier
        Fichier$=CheminWav$+"\"+NomWav$ '
        dim info$(10, 10)
        files CheminWav$, NomWav$, info$()
 
        if val(info$(0, 0))<>0 then 'le fichier existe, le supprimer
            print "kill Fichier$ : ";Fichier$
            kill Fichier$
        end if
 
        print Fichier$
        open Fichier$ for output as #audio
            print #audio,wav$
        close #audio
 
        'contrôle (prévoir un arrêt prématuré
        playwave Fichier$
 
    else
        notice "Valeurs hors-limite trouvée, le fichier wave n'est pas généré"
    end if
 
    wait 'return provisoire 'ce serait-y pas un wait, des fois ?

 
[fichier_8bits_mono]
    print "[fichier_8bits_mono]"
 
    maxoct=0:minoct=256 '(c'est inversé, c'est normal ici)

    for a=0 to nbrEchant-1
        'calcul échantillon
        somme=vol1*sin(omega1*t)+vol2*sin(omega2*t)
        echantillon=int(128+somme*256/amplitudemax)   'échantillons

        'arrondir
        entiere=int(echantillon)
        decimale=somme-int(echantillon)
        if decimales>=0.5 then echantillon=entiere+1 else echantillon=entiere
 
        'stoker l'échantillon en chaine
        periodesomme$=periodesomme$+chr$(echantillon)
        'print "somme = ";somme;", échantillon = ";echantillon

        'prépare pour échantillon suivant
        t=t+tEchant
    next
 
    return
 
 
[fichier_8bits_stereo]
    print "[fichier_8bits_stereo]"
 
    for a=0 to nbrEchant-1
 
        'voie gauche
        'calcul échantillon
        somme=vol1*sin(omega1*t)+vol2*sin(omega2*t)
        echantillon=int(128+somme*256/amplitudemax)   'échantillons

        'arrondir
        entiere=int(echantillon)
        decimale=somme-int(echantillon)
        if decimales>=0.5 then echantillon=entiere+1 else echantillon=entiere
 
        'stoker l'échantillon en chaine
        periodesomme$=periodesomme$+chr$(echantillon)
 
 
        'voie droite
        'calcul échantillon
        somme=vol1*sin(omega1*t)+vol2*sin(omega2*t)
        echantillon=int(128+somme*256/amplitudemax)   'échantillons
        'arrondir
        entiere=int(echantillon)
        decimale=somme-int(echantillon)
        if decimales>=0.5 then echantillon=entiere+1 else echantillon=entiere
 
        'stoker l'échantillon en chaine
        periodesomme$=periodesomme$+chr$(echantillon)
 
        t=t+tEchant
    next
 
return
 
 
[fichier_16bits_mono]
    print "[fichier_16bits_mono]"
    horslimite=0
 
 
 
    'on va placer le bit de signe sur l'octet de poids fort, bit de poids fort
    'les valeurs ne vont plus de 0 à 65536 mais de -32767 à +32768
    '
    print "échantillon définitif, octetFort, octetFaible, contrôle"
 
    for a=0 to nbrEchant-1
 
        'calculer l'échantillon
        somme=vol1*sin(omega1*t)+vol2*sin(omega2*t)
 
        echantillon=somme*65535/amplitudemax   'échantillons

        'arrondir l'échantillon
        entiere=int(echantillon)
        decimale=echantillon-int(echantillon)
 
        echantillon=entiere
        if decimale>=0.5 then echantillon=entiere+1
        if decimale<=-0.5 then echantillon=entiere-1
 
        'convertir l'échantillon
        if echantillon<0 then
            'en dessous de zéro
            echantillon=2^15-abs(echantillon)
            octetFort = abs(int(echantillon/256)+128)
            octetFaible = echantillon mod 256
        else
            'zéro et plus
            octetFort=int(echantillon/256)
            octetFaible=echantillon mod 256
            verif=octetFort*256+octetFaible
        end if
 
        'verifs
        'print echantillon,octetFort,octetFaible, verif

        'enregistrement
        periodesomme$=periodesomme$+chr$(octetFaible)+chr$(octetFort)
 
    t=t+tEchant
    next
 
return
 
 
[fichier_16bits_stereo]
        print "[fichier_16bits_mono]"
    for a=0 to nbrEchant
 
        'voie gauche
        'calcul échantillon
        somme=vol1*sin(omega1*t)+vol2*sin(omega2*t)
 
        echantillon=somme*65535/amplitudemax   'échantillons

        'arrondir l'échantillon
        entiere=int(echantillon)
        decimale=echantillon-int(echantillon)
 
        echantillon=entiere
        if decimale>=0.5 then echantillon=entiere+1
        if decimale<=-0.5 then echantillon=entiere-1
 
        'convertir l'échantillon
        if echantillon<0 then
            'en dessous de zéro
            echantillon=2^15-abs(echantillon)
            octetFort = abs(int(echantillon/256)+128)
            octetFaible = echantillon mod 256
        else
            'zéro et plus
            octetFort=int(echantillon/256)
            octetFaible=echantillon mod 256
            verif=octetFort*256+octetFaible
        end if
 
        'verifs
        'print echantillon,octetFort,octetFaible, verif

        'enregistrement
        periodesomme$=periodesomme$+chr$(octetFaible)+chr$(octetFort)
 
        'voie droite
        'calcul échantillon
        somme=vol1*sin(omega1*t)+vol2*sin(omega2*t)
 
        echantillon=somme*65535/amplitudemax   'échantillons

        'arrondir l'échantillon
        entiere=int(echantillon)
        decimale=echantillon-int(echantillon)
 
        echantillon=entiere
        if decimale>=0.5 then echantillon=entiere+1
        if decimale<=-0.5 then echantillon=entiere-1
 
        'convertir l'échantillon
        if echantillon<0 then
            'en dessous de zéro
            echantillon=2^15-abs(echantillon)
            octetFort = abs(int(echantillon/256)+128)
            octetFaible = echantillon mod 256
        else
            'zéro et plus
            octetFort=int(echantillon/256)
            octetFaible=echantillon mod 256
            verif=octetFort*256+octetFaible
        end if
 
        'verifs
        'print echantillon,octetFort,octetFaible, verif

        'enregistrement
        periodesomme$=periodesomme$+chr$(octetFaible)+chr$(octetFort)
 
        t=t+tEchant
    next
 
    return
 
 
[visu]
    print "[visu]"
 
    'calcul durée de la période reproductible
    affperiode=1/pgcd(freq1,freq2,freq1)
    pas=affperiode/largeurvisu
 
    'calcul échantillons et limites pour visu
    minaff=0:maxaff=0
 
    for a=0 to largeurvisu-1
        ta=a*pas
        affichage(a)=vol1*sin(omega1*ta)+vol2*sin(omega2*ta)
        if affichage(a)>maxaff then maxaff=affichage(a)
        if affichage(a)<minaff then minaff=affichage(a)
    next a
 
    amplitude=maxaff-minaff
    amplitudegraf=hauteurvisu-30
    pasy=amplitudegraf/amplitude
    yzero=25+amplitudegraf/2
 
    'affichage de la période
    print #1.visu, "color cyan"
    print #1.visu, "down"
    print #1.visu, "cls"
    print #1.visu, "fill darkgray"
 
    print #1.visu, "place 10 15"
    print #1.visu, "backcolor darkgray"
    duree$=using("#####.###",affperiode*1000)
    print #1.visu, "\durée : ";duree$;" ms"
 
    print #1.visu, "place 0 ";yzero;" "
 
    for a=0 to largeurvisu-1
        ech=affichage(a)*pasy
        'print "visu, élément ";a;"  : ";ech;" affiché en ";yzero-ech
        print #1.visu, "goto ";a;" ";yzero - ech
        'print "bhem goto ";a;" ";yzero - ech
    next
 
    'dessin axe y
    print #1.visu, "color lightgray"
    print #1.visu, "line 0 ";yzero;" ";largeurvisu;" ";yzero
    print #1.visu, "up"
 
    return
 
 
[Wav_nombre_voies]
    print "[Wav_nombre_voies]"
    print #1.comboVoi,"selectionindex? indexVoies"
    voies=val(WaveNbvoies$(indexVoies))
    wait
 
 
[Wav_resolution]
    print "[Wav_resolution]"
    print #1.comboRes,"selectionindex? indexRes"
    resolution=val(WaveRes$(indexRes))
    wait
 
 
'************** fonctions ***************

function littleEndian$(nombre,longueur)
    print
    print "function littleEndian$ ";
    print "nombre=";nombre;"  longueur=";longueur
    code$="":reste=nombre
 
    for a=longueur-1 to 1 step -1
 
        octet=int(nombre/256^a)
        code$=chr$(octet)+code$
 
        reste=nombre-octet*256^a
 
        nombre=reste
 
        print "boucle=";a;" octet=";octet;"  nombre=";nombre;"  code$=";code$
    next
 
    code$=chr$(reste)+code$
    littleEndian$=code$
 
    'vérification
    for a=1 to longueur
        print asc(mid$(code$,a,1));" ";
    next
    print
end function
 
 
function pgcd(n1,n2,n3)
 
    diviseur=2
    pgcd=1
    fin=0
    tour=0
 
    while fin=0
 
        if n1 mod diviseur =0 and n2 mod diviseur =0 and n3 mod diviseur =0 then
            pgcd=pgcd*diviseur
            n1=n1/diviseur
            n2=n2/diviseur
            n3=n3/diviseur
            'print diviseur;" = facteur commun",n1,n2,n3
        else
            'print diviseur;" = pas commun"
            diviseur=diviseur+1
        end if
 
        'sortie de la boucle quand l'un au moins des nombres est devenu inférieur au diviseur
        'ou à partir d'un diviseur dont il est vraiscemblable qu'on ira jamais jusque là
        tour=tour+1
        if tour=100000 or (n1<diviseur or n2<diviseur or n3<diviseur)then fin=1
 
    wend
 
end function