40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
  40m Log  Not logged in ELOG logo
Entry  Mon Sep 20 15:42:44 2021, Ian MacMillan, Summary, Computers, Quantization Code Summary 472px-Typical_State_Space_model.svg.pngBiquad_filter_DF-IIx.svg.pngSS-IIR-TF.pdf
    Reply  Wed Sep 22 14:22:35 2021, Ian MacMillan, Summary, Computers, Quantization Noise Calculation Summary IIR64-bitnoisespectrum.pdf
       Reply  Mon Sep 27 12:12:15 2021, Ian MacMillan, Summary, Computers, Quantization Noise Calculation Summary 
          Reply  Mon Sep 27 16:03:15 2021, Ian MacMillan, Summary, Computers, Quantization Noise Calculation Summary DeltaNoiseSpectrum.pdf
             Reply  Mon Sep 27 17:04:43 2021, rana, Summary, Computers, Quantization Noise Calculation Summary 
                Reply  Thu Sep 30 11:46:33 2021, Ian MacMillan, Summary, Computers, Quantization Noise Calculation Summary SS-IIR-Bode.pdfPSD_32bit.pdfPSD_64bit.pdf
                   Reply  Wed Nov 24 11:02:23 2021, Ian MacMillan, Summary, Computers, Quantization Noise Calculation Summary PSD_IIR_all.pdf
                      Reply  Wed Nov 24 13:44:19 2021, rana, Summary, Computers, Quantization Noise Calculation Summary 
Message ID: 16366     Entry time: Thu Sep 30 11:46:33 2021     In reply to: 16362     Reply to this: 16481
Author: Ian MacMillan 
Type: Summary 
Category: Computers 
Subject: Quantization Noise Calculation Summary 

First and foremost I have the updated bode plot with the mode moved to 10 Hz. See Attachment 1. Note that the comparison measurement is a % difference whereas in the previous bode plot it was just the difference. I also wrapped the phase so that jumps from -180 to 180 are moved down. This eliminates massive jumps in the % difference. 

Next, I have two comparison plots: 32 bit and 64bit. As mentioned above I moved the mode to 10 Hz and just excited both systems at 3.4283Hz with an amplitude of 1. As we can see on the plot the two models are practically the same when using 64bits. With the 32bit system, we can see that the noise in the IIR filter is much greater than in the State-space model at frequencies greater than our mode.

Note about windowing and averaging: I used a Hanning window with averaging over 4 neighbor points. I came to this number after looking at the results with less averaging and more averaging. In the code, this can be seen as nperseg=num_samples/4 which is then fed into signal.welch

Attachment 1: SS-IIR-Bode.pdf  94 kB  Uploaded Thu Sep 30 12:56:21 2021  | Hide | Hide all
SS-IIR-Bode.pdf
Attachment 2: PSD_32bit.pdf  26 kB  Uploaded Thu Sep 30 13:13:39 2021  | Hide | Hide all
PSD_32bit.pdf
Attachment 3: PSD_64bit.pdf  18 kB  Uploaded Thu Sep 30 13:14:11 2021  | Hide | Hide all
PSD_64bit.pdf
ELOG V3.1.3-