Skip to content
Snippets Groups Projects
Commit afdc10e5 authored by Thomas Bock's avatar Thomas Bock :speech_balloon:
Browse files

T se3

parent b8e5bc0f
No related branches found
No related tags found
No related merge requests found
No preview for this file type
No preview for this file type
No preview for this file type
......@@ -168,12 +168,12 @@ gibt den vom Temperaturnormal angezeigten Wert (zitiert als $T_{PTB}$) wieder.
Zur Auswertung bzw. zur Ermittlung der Sensorkorrekturen $k_i$ wurde der
arithmetische Mittelwert der Differenzen $T_{PTB} - T_{i}$
gebildet. $k_i$ muss zur Korrektur eines Anzeigewertes addiert
gebildet. Um nach der Kalibrierung eine gemessene Temperatur $T_i'$ zu korrigieren muss $k_i$ addiert
werden.
\begin{equation}
\label{eq:k}
T_{korr,i} = T_i + k_i \rightarrow T_i +\frac{1}{N}\sum_{j=1}^N( T_{PTB} - T_{i})_j \rightarrow T_{PTB}
T_{korr,i} = T_i' + k_i = T_i' + (T_{PTB} - T_{i})
\end{equation}
In nachfolgender Rechnung werden beide Messreichen (
......@@ -181,7 +181,7 @@ In nachfolgender Rechnung werden beide Messreichen (
und die Ergebnisse in den |Result|-Abschnitt der jeweiligen Dokumente geschrieben.
%% begin.rcode
%% begin.rcode echo=FALSE
%% library(R4CouchDB)
%% library(vacLabTools)
%% cdb <- cdbIni()
......@@ -290,11 +290,12 @@ berücksichtigen:
% # Unsicherheit des Temperaturnormals (KS-Angabe ist 5mK)
% # hizu kommt Langzeitstabilität Abschätzung aus
% # Veränderung des Widerstandswertes
% # (2016:Ohm, 2018: 100.021Ohm):
% # (2016: 100.019Ohm, 2018: 100.021Ohm):
% R.2016 <- 100.019
% R.2018 <- 100.021
% u.lt <-(R.2016 / R.2018 -1 )*23.0
% u1 <- 0.0025 + u.lt # K
% u.lt <-(R.2016 / R.2018 -1 )*269 # in K
% u.tr <- 0.005 # K
% u1 <- (0.0025^2 + u.lt^2 + u.tr^2)^0.5 # K
% u1
% # Digitalisierung
% u2 <- 2.9e-5 # K
......@@ -310,9 +311,9 @@ berücksichtigen:
% u5 <- abs(mean(corr.mean - corr.17)) * 2 # K
% u5 # K
% # experimentelle Streuung wird aus der Standardabw.
% # der Einzelwerte eines Messpunktes abgeschätzt
% # (doppelter Mittelwert)
% u6 <- mean(corr.sd) * 2 # K
% # der Korrekturen abgeschätzt
% u6 <- mean(corr.sd) *2 # K
% u6
% plot(x, corr.sd,
% ylab="mittlere Standardabw. der Einzelwerte",
% xlab="Kanal",
......
No preview for this file type
......@@ -218,12 +218,12 @@ gibt den vom Temperaturnormal angezeigten Wert (zitiert als $T_{PTB}$) wieder.
Zur Auswertung bzw. zur Ermittlung der Sensorkorrekturen $k_i$ wurde der
arithmetische Mittelwert der Differenzen $T_{PTB} - T_{i}$
gebildet. $k_i$ muss zur Korrektur eines Anzeigewertes addiert
gebildet. Um nach der Kalibrierung eine gemessene Temperatur $T_i'$ zu korrigieren muss $k_i$ addiert
werden.
\begin{equation}
\label{eq:k}
T_{korr,i} = T_i + k_i \rightarrow T_i +\frac{1}{N}\sum_{j=1}^N( T_{PTB} - T_{i})_j \rightarrow T_{PTB}
T_{korr,i} = T_i' + k_i = T_i' + (T_{PTB} - T_{i})
\end{equation}
In nachfolgender Rechnung werden beide Messreichen (
......@@ -232,67 +232,9 @@ und die Ergebnisse in den |Result|-Abschnitt der jeweiligen Dokumente geschriebe
\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlkwd{library}\hlstd{(R4CouchDB)}
\hlkwd{library}\hlstd{(vacLabTools)}
\hlstd{cdb} \hlkwb{<-} \hlkwd{cdbIni}\hlstd{()}
\hlstd{cdb}\hlopt{$}\hlstd{DBName} \hlkwb{<-} \hlstr{"vl_db"}
\hlstd{cdb}\hlopt{$}\hlstd{id} \hlkwb{<-} \hlstr{"cal-2018-mkt50-ik-4040_0001"}
\hlstd{doc} \hlkwb{<-} \hlkwd{cdbGetDoc}\hlstd{(cdb)}\hlopt{$}\hlstd{res}
\hlstd{a} \hlkwb{<-} \hlkwd{abbrevList}\hlstd{(doc)}
\hlstd{doc} \hlkwb{<-} \hlkwd{refreshResult}\hlstd{(doc)}
\hlstd{t.ms.b} \hlkwb{<-} \hlkwd{getConstVal}\hlstd{(a}\hlopt{$}\hlstd{cmv}\hlopt{$}\hlstd{Time,} \hlstr{"begin"}\hlstd{)}
\hlstd{t.ms.e} \hlkwb{<-} \hlkwd{getConstVal}\hlstd{(a}\hlopt{$}\hlstd{cmv}\hlopt{$}\hlstd{Time,} \hlstr{"end"}\hlstd{)}
\hlstd{t.h} \hlkwb{<-} \hlstd{(t.ms.b} \hlopt{-} \hlstd{t.ms.b[}\hlnum{1}\hlstd{])}\hlopt{/}\hlnum{1000}\hlopt{/}\hlnum{3600}
\hlstd{t.ptb} \hlkwb{<-} \hlkwd{getConstVal}\hlstd{(a}\hlopt{$}\hlstd{cmv}\hlopt{$}\hlstd{Temperature,} \hlstr{"mkt50"}\hlstd{)}
\hlkwd{par}\hlstd{(}\hlkwc{mfrow}\hlstd{=}\hlkwd{c}\hlstd{(}\hlnum{1}\hlstd{,}\hlnum{1}\hlstd{))}
\hlstd{k}\hlkwb{=}\hlnum{1}
\hlstd{ch} \hlkwb{<-} \hlkwd{c}\hlstd{(}\hlnum{1001}\hlopt{:}\hlnum{1030}
\hlstd{,} \hlnum{2001}\hlopt{:}\hlnum{2030}
\hlstd{,} \hlnum{3001}\hlopt{:}\hlnum{3030}\hlstd{)}
\hlstd{ch.N} \hlkwb{<-} \hlkwd{length}\hlstd{(ch)}
\hlstd{corr.sd} \hlkwb{<-} \hlkwd{rep}\hlstd{(}\hlnum{NA}\hlstd{, ch.N)}
\hlstd{corr.mean} \hlkwb{<-} \hlkwd{rep}\hlstd{(}\hlnum{NA}\hlstd{, ch.N)}
\hlkwa{for}\hlstd{( i} \hlkwa{in} \hlnum{1}\hlopt{:}\hlstd{ch.N)\{}
\hlstd{ch.name} \hlkwb{<-} \hlkwd{paste0}\hlstd{(}\hlstr{"ch_"}\hlstd{,ch[i])}
\hlstd{t.ch} \hlkwb{<-} \hlkwd{getConstVal}\hlstd{(a}\hlopt{$}\hlstd{cmv}\hlopt{$}\hlstd{Temperature, ch.name)}
\hlstd{corr} \hlkwb{<-} \hlstd{t.ptb} \hlopt{-} \hlstd{t.ch}
\hlstd{corr.sd[i]} \hlkwb{<-} \hlkwd{sd}\hlstd{(corr)}
\hlstd{corr.mean[i]} \hlkwb{<-} \hlkwd{mean}\hlstd{(corr)}
\hlstd{doc}\hlopt{$}\hlstd{Calibration}\hlopt{$}\hlstd{Result}\hlopt{$}\hlstd{Values}\hlopt{$}\hlstd{Temperature} \hlkwb{<-}
\hlkwd{setCcl}\hlstd{(doc}\hlopt{$}\hlstd{Calibration}\hlopt{$}\hlstd{Result}\hlopt{$}\hlstd{Values}\hlopt{$}\hlstd{Temperature,}
\hlkwd{paste0}\hlstd{(}\hlstr{"corr_"}\hlstd{, ch.name),}
\hlstr{"K"}\hlstd{,}
\hlstd{corr.mean[i])}
\hlkwa{if}\hlstd{(k}\hlopt{==}\hlnum{1}\hlstd{)\{}
\hlkwd{plot}\hlstd{(corr}
\hlstd{,} \hlkwc{ylim}\hlstd{=}\hlkwd{c}\hlstd{(}\hlopt{-}\hlnum{0.3}\hlstd{,}\hlnum{0.2}\hlstd{)}
\hlstd{,} \hlkwc{main}\hlstd{=cdb}\hlopt{$}\hlstd{id}
\hlstd{,} \hlkwc{xlab} \hlstd{=} \hlstr{"Messpkt."}
\hlstd{,} \hlkwc{ylab}\hlstd{=}\hlstr{"k in K"}
\hlstd{)}
\hlstd{\}}\hlkwa{else}\hlstd{\{}
\hlkwd{points}\hlstd{(corr,} \hlkwc{col} \hlstd{= k)}
\hlstd{\}}
\hlstd{k} \hlkwb{=} \hlstd{k}\hlopt{+}\hlnum{1}
\hlstd{\}}
\end{alltt}
\end{kframe}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}
\includegraphics[width=\maxwidth]{figure/unnamed-chunk-1-1}
\begin{kframe}\begin{alltt}
\hlkwa{if}\hlstd{(}\hlopt{!}\hlnum{FALSE}\hlstd{)\{}
\hlstd{cdb}\hlopt{$}\hlstd{dataList} \hlkwb{<-} \hlstd{doc}
\hlkwd{cdbUpdateDoc}\hlstd{(cdb)}\hlopt{$}\hlstd{res}
\hlstd{\}}
\end{alltt}
\begin{verbatim}
\begin{kframe}\begin{verbatim}
## $ok
## [1] TRUE
##
......@@ -300,7 +242,7 @@ und die Ergebnisse in den |Result|-Abschnitt der jeweiligen Dokumente geschriebe
## [1] "cal-2018-mkt50-ik-4040_0001"
##
## $rev
## [1] "147-ba39188bcd24a6fb2f81195b8a1e6749"
## [1] "158-352f0089de7a62497a7fe5e6e4d905cf"
\end{verbatim}
\end{kframe}
\end{knitrout}
......@@ -343,15 +285,16 @@ berücksichtigen:
\hlcom{# Unsicherheit des Temperaturnormals (KS-Angabe ist 5mK)}
\hlcom{# hizu kommt Langzeitstabilität Abschätzung aus}
\hlcom{# Veränderung des Widerstandswertes}
\hlcom{# (2016:Ohm, 2018: 100.021Ohm):}
\hlcom{# (2016: 100.019Ohm, 2018: 100.021Ohm):}
\hlstd{R.2016} \hlkwb{<-} \hlnum{100.019}
\hlstd{R.2018} \hlkwb{<-} \hlnum{100.021}
\hlstd{u.lt} \hlkwb{<-}\hlstd{(R.2016} \hlopt{/} \hlstd{R.2018} \hlopt{-}\hlnum{1} \hlstd{)}\hlopt{*}\hlnum{23.0}
\hlstd{u1} \hlkwb{<-} \hlnum{0.0025} \hlopt{+} \hlstd{u.lt} \hlcom{# K}
\hlstd{u.lt} \hlkwb{<-}\hlstd{(R.2016} \hlopt{/} \hlstd{R.2018} \hlopt{-}\hlnum{1} \hlstd{)}\hlopt{*}\hlnum{269} \hlcom{# in K}
\hlstd{u.tr} \hlkwb{<-} \hlnum{0.005} \hlcom{# K}
\hlstd{u1} \hlkwb{<-} \hlstd{(}\hlnum{0.0025}\hlopt{^}\hlnum{2} \hlopt{+} \hlstd{u.lt}\hlopt{^}\hlnum{2} \hlopt{+} \hlstd{u.tr}\hlopt{^}\hlnum{2}\hlstd{)}\hlopt{^}\hlnum{0.5} \hlcom{# K}
\hlstd{u1}
\end{alltt}
\begin{verbatim}
## [1] 0.00204009657972
## [1] 0.00775772177769
\end{verbatim}
\begin{alltt}
\hlcom{# Digitalisierung}
......@@ -378,9 +321,14 @@ berücksichtigen:
\end{verbatim}
\begin{alltt}
\hlcom{# experimentelle Streuung wird aus der Standardabw.}
\hlcom{# der Einzelwerte eines Messpunktes abgeschätzt}
\hlcom{# (doppelter Mittelwert)}
\hlstd{u6} \hlkwb{<-} \hlkwd{mean}\hlstd{(corr.sd)} \hlopt{*} \hlnum{2} \hlcom{# K}
\hlcom{# der Korrekturen abgeschätzt}
\hlstd{u6} \hlkwb{<-} \hlkwd{mean}\hlstd{(corr.sd)} \hlopt{*}\hlnum{2} \hlcom{# K}
\hlstd{u6}
\end{alltt}
\begin{verbatim}
## [1] 0.00156246202678
\end{verbatim}
\begin{alltt}
\hlkwd{plot}\hlstd{(x, corr.sd,}
\hlkwc{ylab}\hlstd{=}\hlstr{"mittlere Standardabw. der Einzelwerte"}\hlstd{,}
\hlkwc{xlab}\hlstd{=}\hlstr{"Kanal"}\hlstd{,}
......@@ -408,7 +356,7 @@ Man erhält eine erweiterte Gesamtunsicherheit von:
\hlstd{k} \hlopt{*} \hlkwd{sqrt}\hlstd{(u1}\hlopt{^}\hlnum{2} \hlopt{+} \hlstd{u2}\hlopt{^}\hlnum{2} \hlopt{+} \hlstd{u3}\hlopt{^}\hlnum{2} \hlopt{+}\hlstd{u4}\hlopt{^}\hlnum{2} \hlopt{+}\hlstd{u5}\hlopt{^}\hlnum{2} \hlopt{+}\hlstd{u6}\hlopt{^}\hlnum{2}\hlstd{)}
\end{alltt}
\begin{verbatim}
## [1] 0.0122605597411
## [1] 0.0193494790025
\end{verbatim}
\end{kframe}
\end{knitrout}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment