Measurement Data: Add PT1000 measurements (current SW).

This commit is contained in:
Mario Hüttel 2020-02-05 23:11:36 +01:00
parent 825de053ce
commit a6fa35b1ea
3 changed files with 16152 additions and 9 deletions

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

View File

@ -344,10 +344,9 @@
"outputs": [],
"source": [
"# Read in data\n",
"kilo_ohm_sampling1 = pd.read_csv(r'1000OhmSampling1ms.csv')\n",
"kilo_ohm_sampling2 = pd.read_csv(r'1000OhmSampling1ms_later.csv')\n",
"kilo_ohm_sampling2_fast = pd.read_csv(r'1000OhmSamplingTimerTriggered800ARR.csv')\n",
"plot_data_tuples = [(kilo_ohm_sampling1, 'Day 1 Sampling'), (kilo_ohm_sampling2, 'Day 2 sampling'), (kilo_ohm_sampling2_fast,'Day 2 sampling faster (ARR=700)')]\n"
"kilo_ohm_sampling1 = pd.read_csv(r'1000OhmSamplingFixedStableCircuit.csv')\n",
"kilo_ohm_sampling1_ht = pd.read_csv(r'1000OhmSamplingFixedStableCircuitHT.csv')\n",
"plot_data_tuples = [(kilo_ohm_sampling1, 'Day 1 Sampling -- Stable circuit -- RT'), (kilo_ohm_sampling1_ht, 'Day 1 Sampling -- Stable circuit -- HT')]\n"
]
},
{
@ -360,14 +359,17 @@
"# def plot_histogram(ax, data, bin_count, title, signal_name):\n",
"# def calculate_temp_for_df(data_frame, resistance_col_name='ext_lf_corr', temp_col_name='temp_calculated'):\n",
"\n",
"fig, axes = plt.subplots(nrows=len(plot_data_tuples), ncols=2, sharex='col', figsize=(28, 16))\n",
"fig, axes = plt.subplots(nrows=len(plot_data_tuples), ncols=2, sharex='col', figsize=(28, 8*len(plot_data_tuples)))\n",
"\n",
"if len(plot_data_tuples) == 1:\n",
" axes = [axes]\n",
"\n",
"for df, title in plot_data_tuples:\n",
" calculate_temp_for_df(df, resistance_col_name='pt1000_res_raw_lf')\n",
" calculate_temp_for_df(df, resistance_col_name='pt1000_value')\n",
"\n",
"for (df,title),ax in zip(plot_data_tuples, axes):\n",
" plot_histogram(ax[0], df['pt1000_res_raw_lf'], 20, title, 'PT1000 Resistance')\n",
" plot_histogram(ax[1], df['temp_calculated'], 20, title, 'Calculated Temperature in °C')\n",
" plot_histogram(ax[0], df['pt1000_value'], 5, title, 'PT1000 Resistance')\n",
" plot_histogram(ax[1], df['temp_calculated'], 5, title, 'Calculated Temperature in °C')\n",
"\n",
"plt.tight_layout()\n",
"plt.show()"
@ -379,7 +381,7 @@
"metadata": {},
"outputs": [],
"source": [
"calc_temp(1097)"
"calc_temp(1100)"
]
}
],