modules/devices/control_system/filters/biquad/test/results/get_telemetry.m

   1 \r
   2 \r
   3 \r
   4 filename = 'output.log';\r
   5 \r
   6 display('Lecture telemetrie ...')\r
   7 tableau = telemetry(filename);\r
   8 display('Fin telemetrie')\r
   9 clear filename;\r
  10 \r
  11 affich_telemetry(tableau);\r
  12 \r
  13 % titre\r
  14 subplot(2, 1, 1);\r
  15 title('input : pseudorandom noise');\r
  16 \r
  17 subplot(2, 1, 2);\r
  18 title('filtered output');\r
  19 \r
  20 \r
  21 tableau = tableau(2:end, :);\r
  22 \r
  23 \r
  24 % display spectrums\r
  25 \r
  26 \r
  27 figure;\r
  28 Y1 = fft(tableau(:, 1),512);\r
  29 Pyy1 = Y1.* conj(Y1) / 512;\r
  30 f = 1000*(0:256)/512;\r
  31 semilogy(f,Pyy1(1:257))\r
  32 \r
  33 hold on;\r
  34 \r
  35 Y1 = fft(tableau(:, 2),512);\r
  36 Pyy1 = Y1.* conj(Y1) / 512;\r
  37 f = 1000*(0:256)/512;\r
  38 semilogy(f,Pyy1(1:257), '-r')\r
  39 \r
  40 xlabel('frequency (Hz)')\r
  41 legend('in', 'out');