LPP » INSTRU » Solar Orbiter LFR » LFR-FSW

Questions complémentaires.
1) Au niveau des intervalles des FFT, pour f2 par exemple, la bin 0 est la composant continue, la bin 1 va de la fréquence 0Hz à 1Hz et la bin 255 va de 254Hz à 255Hz je pense. Pourquoi ce 0.5Hz? Les limites tombent sur des Hz entiers non?
2) Il ne s'agit pas plutôt de LFR_BP2_F*2*/LFR_BP2_AUTO_F*2* au lieu de LFR_BP2_F*0*/LFR_BP2_AUTO_F*0*

1) Bah non. Quand on fait une FTT de $N$ points, donc avec $\Delta f = f_e/N$ de résolution, chacune de ses raie $i$ est centrée sur la fréquence $f_i = i \times \Delta f$. En d'autres termes, une raie $i$ est la raie correspondant à l'énergie comprise sur la bande de fréquence $[f_i - \Delta f/2 , f_i +\Delta f/2 ]$. Non ?

=> pour répondre à Gérald., le bin1 du canal F2 correspond à [0.5;1.5] Hz

Ensuite le choix d'un intervalle ouvert à droite ou à gauche est une convention à se fixer une fois pour toute.

Dans la FFT, le spectre est discrétisé, donc la question des intervalles de fréquences ne se pose peut-être pas. La question équivalente dans le domaine temporel reviendrait à se demander à quel interval appartiendrait un point qui ne serait pas pris à un instant d'échantillonnage.

Voici un début de réponse sur la façon dont le filtre est généré pour un mask (représenté par un char*, soit 16 octets), pour une fréquence (fonction définie dans tc_load_dump_parameter.c). A la relecture, j'ai utilisé deux fois la fonction floor, j'aurais peut-être dû utiliser ceil pour la borne sup de l'intervalle à rejeter.

void setFBinMask( unsigned char *fbins_mask, float rw_f, unsigned char deltaFreq, unsigned char flag )
{
    /** This function executes specific actions when a TC_LFR_UPDATE_INFO TeleCommand has been received.
     *
     * @param fbins_mask
     * @param rw_f is the reaction wheel frequency to filter
     * @param delta_f is the frequency step between the frequency bins, it depends on the frequency channel
     * @param flag [true] filtering enabled [false] filtering disabled
     *
     * @return void
     *
     */

    float fmin;
    float fMAX;
    int binBelow;
    int binAbove;
    unsigned int whichByte;
    unsigned char selectedByte;
    int bin;

    whichByte = 0;
    bin = 0;

    // compute the frequency range to filter [ rw_f - delta_f/2; rw_f + delta_f/2 ]
    fmin = rw_f - filterPar.sy_lfr_sc_rw_delta_f / 2.;
    fMAX = rw_f + filterPar.sy_lfr_sc_rw_delta_f / 2.;

    // compute the index of the frequency bin immediately below fmin
    binBelow = (int) ( floor( ((double) fmin) / ((double) deltaFreq)) );

    // compute the index of the frequency bin immediately above fMAX
    binAbove = (int) ( floor( ((double) fMAX) / ((double) deltaFreq)) );

    for (bin = binBelow; bin <= binAbove; bin++)
    {
        if ( (bin >= 0) && (bin<=127) )
        {
            if (flag == 1)
            {
                whichByte = bin >> 3;   // division by 8
                selectedByte = (unsigned char) ( 1 << (bin - (whichByte * 8)) );
                fbins_mask[whichByte] = fbins_mask[whichByte] & (~selectedByte);
            }
        }
    }
}

Je réalise qu'on est sur une mauvaise approche:
- la question n'est pas que vaux la largeur de la raie $i$ d'une FFT mais exactement l'inverse: lorsqu'on introduit un signal à une fréquence $f_{RW}$ qu'elles sont les raies de la FFT qui sont touchées ?
Et la réponse me semble bien plus drastique que ce que l'on a envisagé jusq'alors:
- même sans fenêrage de Hanning, si $f_{RW}$ tombe exactement entre 2 raies ces 2 raies vont être polluées à égalité => il faut virer ces 2 raies ...
- avec fenêtrage de Hanning, ce qui est notre cas toujours, même réponse si $f_{RW}$ tombe entre 2 raies mais lorsque $f_{RW}$ tombe pile sur une raie, les 2 raies adjacentes ne sont atténuées que par un facteur 0.25 (en énergie, les ASM sont des FFT au carré). Il sera donc peut-être utile d'envisager de virer aussi ces 2 raies adjacentes => 3 raies à virer et pas seulement 1 ...

Une proposition pourrait être la suivante:
- si l'intervalle $f_{RW} \pm \Delta f_{RW}$ est inclu dans l'intervalle $f_i \pm \Delta f/m$ (m valant a priori au moins 4) on vire la raie $f_i$ et les 2 raies adjacentes $f_{i-1}$ et $f_{i+1}$
- dans tous les autres cas en vire les 2 raies qui entourent $f_{RW}$

Ca me va!