Exemple #1
0
        static void Main(string[] args)
        {
            string argument;

            // script python a executer pour ploter les resultats
            //-------------------------------------------------
            argument = "C:\\C#\\MathsLib\\Test_sysarg-v4.py" + " ";

            //test de plusieurs fonctionnalites
            //-------------------------------------------------

            int    nbechant;
            double Fe;

            Fe = 1000.0;   //1000 Hz



            //calcul du nombre de points a echantillonner
            //sur un crénau de 0.25 secondes (250 milisecondes) avec une fréquence d'échantillonnage de 1000 Hz
            // échantillonnage a 1 points toutes les 1/1000 = 0.001 secondes

            nbechant = SiTest.Nbechant(1000, 0.25);

            Console.WriteLine("nbechant = {0}", nbechant);

            decimal[] S2 = new decimal[nbechant];
            double[]  S1 = new double[nbechant];

            //Form a signal containing a 50 Hz sinusoid of amplitude 4
            //S2 = 4*sin(2*pi*50*t)

            S2 = SiTest.GenSin2piF(nbechant, 50.0, 4.0, Fe);

            // resultat de la DFT (Discret Fourier Transform) dans le tableau de complexe : Ts[]
            //----------------------------------------------------------------------------------
            Complex[] Ts = new Complex[nbechant];
            Ts = DFT.DFTv2(S2);

            //Recuperation du module de la DFT
            decimal[] Module = new decimal[nbechant];
            //Module = DFT.Dsdec(Ts);

            //Densite spectrale de puissance du signal
            //---------------------------------------------------------------------------
            decimal[] Module1 = new decimal[nbechant];
            Module1 = DSP.Dspdeci(Ts);

            //plot signal de test : Tableau Ti1[]
            //-------------------------------------------------------------------
            PlotPython.Plot1(S2, argument);


            //plot Densite spectrale de puissance du sgignal
            //------------------------------------------------------------------
            PlotPython.Plot1(Module1, argument);
        }
Exemple #2
0
        static void Main(string[] args)
        {
            string argument;

            // script python a executer pour ploter les resultats
            //-------------------------------------------------
            argument = "C:\\C#\\MathsLib\\Test_sysarg-v4.py" + " ";

            //test de plusieurs fonctionnalites
            //-------------------------------------------------

            int nbechant;

            //nombre de points a echantillonner
            nbechant = SiTest.Nbechant(1000, 0.25);  //crenau de 0.25 secondes (250 milisecondes) et frequence de 1000 Hz

            //Form a signal containing a 50 Hz sinusoid of amplitude 1 and a 120 Hz sinusoid of amplitude 1.
            //S = 1*sin(2*pi*50*t) + sin(2*pi*120*t);

            decimal[] S1 = new decimal[nbechant];
            decimal[] S2 = new decimal[nbechant];

            S1 = SiTest.GenSin(nbechant, 50, 1);
            S2 = SiTest.GenSin(nbechant, 120, 1);

            // Instantiate random number generator using system-supplied value as seed.
            // signal de bruit blanc ( random noise signal)
            //-----------------------------------------------------------------------
            decimal[] Rand = new decimal[nbechant];
            Rand = SiTest.Random(nbechant);


            // genere le signal de test : s1 + s2  + bruit blanc d Amplidtude 2
            //-----------------------------------------------------------------------
            decimal[] Ti = new decimal[nbechant];
            for (int i = 0; i < nbechant; i++)
            {
                Ti[i] = S1[i] + S2[i] + (2 * Rand[i]);
            }

            // resultat de la DFT (Discret Fourier Transform) dans le tableau de complexe : Ts[]
            //----------------------------------------------------------------------------------
            Complex[] Ts = new Complex[nbechant];
            Ts = DFT.DFTv2(Ti);

            //Recuperation du module de la DFT
            decimal[] Module = new decimal[nbechant];
            //Module = DFT.Dsdec(Ts);

            //Densite spectrale de puissance du signal
            //---------------------------------------------------------------------------
            decimal[] Module1 = new decimal[nbechant];
            Module1 = DSP.Dspdeci(Ts);

            // Densite spectrale de puissance du signal
            double[] Module3 = new double[nbechant];
            //Module3 = DSP.Dspdoub(Ts);

            //impression partie reelle (amplitude du signal transforme par la DFT
            //PlotPython.Plot1(Ti1, argument);

            //impression partie imaginaire (phase du signal transforme par la DFT
            //PlotPython.Plot1(Tc2, argument);

            //plot signal de test : Tableau Ti1[]
            //-------------------------------------------------------------------
            PlotPython.Plot1(Ti, argument);

            //plot magnitude du sgignal calcule par DFTv2
            //PlotPython.Plot1(Module, argument);

            //plot Densite spectrale de puissance du sgignal
            //------------------------------------------------------------------
            PlotPython.Plot1(Module1, argument);
        }