// Q 1 // pour effacer les variables clear all; //*********************************************************** // génération d'une fenêtre triangulaire de longueur 32 : tri1 = window('tr',32); // domaine fréquentiel : TFD de rect1 TRI1=fft(tri1,-1); TRI1=fftshift(TRI1); f1 = -0.5:1/32:0.5-1/32; // fréquences normalisées échantillonnées (32 points) // affichage xset("window",0);clf(); // selectionne (ou crée) la fenetre 0 et l'efface subplot(211);plot2d(0:31,tri1,-9) xtitle("Fenetre triangulaire de longueur N=32") subplot(212);plot2d(f1,20*log10(abs(TRI1/length(TRI1)+%eps)),-9,"011","",[-0.5 -100 0.5 0]) // NB : normalisation / NFFT xtitle("TFD sur 32 points (en dB)","f (kHz)") //********************************************** // Complétion de tri1 par des zéros : //tri2 = [zeros(1,240),tri1,zeros(1,240)]; //// TFD de tri2 //TRI2=fft(tri2,-1); //TRI2=fftshift(TRI2); //f2 = -0.5:1/512:0.5-1/512; // fréquences normalisées échantillonnées (512 points) // //// affichage //xset("window",1);xbasc(); //subplot(211);plot2d(0:511,tri2,1) //xtitle("Fenetre triangulaire de longueur N=32") //subplot(212);plot2d(f2,20*log10(abs(TRI2/length(TRI2)+%eps)),1,"011","",[-0.5 -150 0.5 0]) // NB : normalisation / NFFT //xtitle("TFD sur 512 points (en dB)","f (kHz)")