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);
}
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);
}
