The programme named “psd_ECG-COMON.exe” calculates the potential PSD, the current PSD, and the noise impedance Zn
(square root of the ratio PSD V/PSD I) of a EN data file consisting of three columns (time in s, potential
noise in V, and current noise in A) separated by tabulators. The programme can also divide the time record in one
or several sections over which it averages the PSDs. It uses the configuration file “config_psd.txt”.
These two files must be in the directory where the data (noise ASCII file) is. Only the configuration file has to
be modified according to your data. The structure of the configuration file is the following:

fichier_in hanning (0 = no, 1 = yes) nb_section nbpoint_section fichier_out

(fichier_in: name of the EN data file (ASCII file); hanning: Hann window will be used (1) or not (0);
nb_section: number of sections in which the EN data file is divided; nbpoint_section: number of points/samples
in each section; fichier_out: name of the file with the PSD and Zn data)


For example:

	Huet_RR1Mo_10Hz.txt 1 10 2048 Huet_RR1Mo_10HzPSD.txt

for a data file “Huet_RR1Mo_10Hz.txt” containing 20480 pairs of V and I values. In that case, the data file is
divided in 10 sections of 2048 pairs of V and I values and the Hann window is used in each section. The result
is stored in the file “Huet_RR1Mo_10HzPSD.txt” as follows:

data file = Huet_RR1Mo_10Hz.txt    hanning = 1    nb_section = 10    nbpoint_section = 2048
  frequency	    PSD_v	     PSD_i	      Zn
4.882813e-003	1.901942e-011	2.974846e-024	2.528520e+006
9.765625e-003	1.913229e-012	9.364153e-025	1.429385e+006
1.464844e-002	5.789599e-013	5.085475e-025	1.066985e+006
...


Please note:

• no character "blank" in the names of fichier_in and fichier_out
• nbpoint_section must be a power of 2
• in the noise ASCII file the values of potential and current noise have to be in V and A (if not, ask F. Huet for a modified programme)



For modifications of the programme, please feel free to ask F. Huet (francois.huet@upmc.fr).