Как нарисовать петлю гистерезиса с помощью визуализации данных TikZ?

Question

Извините за мой неправильный английский.

Я изменил код сздесьчтобы адаптировать его для TikZ datavisualization. Вот адаптированный код:

\documentclass{scrartcl}
\usepackage{tikz}
\usepackage[per-mode = fraction]{siunitx}
\usetikzlibrary{datavisualization.formats.functions}
\begin{document}
  \begin{tikzpicture}
    \datavisualization[
                       scientific axes = {clean, end labels},
                       all axes = {ticks and grid = {major at = 0}},
                       x axis = {label = $\frac{H}{\si{\A\per\m}}$},
                       y axis = {label = $\frac{B(H)}{\si{\tesla}}$},
                       data/format = function,
                       visualize as smooth line/.list = {left, right}
                      ]
    data[set = left] {
      var x : interval [-7 : 7];
      func y = 5 / (1 + exp(-1.7 * \value x - 1.5)) - 2.5;
    }
    data[set = right] {
      var x : interval [-7 : 7];
      func y = 5 / (1 + exp(-1.7 * \value x + 1.5)) - 2.5;
    }
    %accents important points
    info {
      \draw (visualization cs:x = {(-1.5 / 1.7)}, y = 0) circle [radius = 1pt]
        node [left, font = \footnotesize] {$H_C$};
    }
    info {
      \draw (visualization cs:x = 0, y = 1.59) circle [radius = 1pt]
        node [left, font=\footnotesize] {$B_R$};
    }
    info {
      \draw (visualization cs:x = 7, y = 2.5) circle [radius = 1pt]
        node [above, font=\footnotesize] {$B_S$};
    };
  \end{tikzpicture}
\end{document}

Но я получаю много таких ошибок:

! Размер слишком большой. \pgfmath@iterate...\pgfmath@xa \ifdim \pgfmath@x

Можно использовать TikZ fpuкак сделалздесь,здесьиздесь. Вот мой код:

\documentclass{scrartcl}
\usepackage{tikz}
\usepackage[per-mode = fraction]{siunitx}
\usetikzlibrary{datavisualization.formats.functions,fpu}
\begin{document}
  \begin{tikzpicture}
    \datavisualization[
                       scientific axes = {clean, end labels},
                       all axes = {ticks and grid = {major at = 0}},
                       x axis = {label = $\frac{H}{\si{\A\per\m}}$},
                       y axis = {label = $\frac{B(H)}{\si{\tesla}}$},
                       data/format = function,
                       visualize as smooth line/.list = {left, right}
                      ]
    data[set = left] {
      var x : interval [-7 : 7];
      \pgfkey{/pgf/fpu, /pgf/fpu/output format=fixed}
      func y = 5 / (1 + exp(-1.7 * \value x - 1.5)) - 2.5;
      \pgfkeys{/pgf/fpu=false}
    }
    data[set = right] {
      var x : interval [-7 : 7];
      \pgfkey{/pgf/fpu, /pgf/fpu/output format=fixed}
      func y = 5 / (1 + exp(-1.7 * \value x + 1.5)) - 2.5;
      \pgfkeys{/pgf/fpu=false}
    }
    %accents important points
    info {
      \draw (visualization cs:x = {(-1.5 / 1.7)}, y = 0) circle [radius = 1pt]
        node [left, font = \footnotesize] {$H_C$};
    }
    info {
      \draw (visualization cs:x = 0, y = 1.59) circle [radius = 1pt]
        node [left, font=\footnotesize] {$B_R$};
    }
    info {
      \draw (visualization cs:x = 7, y = 2.5) circle [radius = 1pt]
        node [above, font=\footnotesize] {$B_S$};
    };
  \end{tikzpicture}
\end{document}

Но тогда мой pdfLaTeXкомпилятор TeXnicCenterпопадает в бесконечный цикл. Что я делаю не так?

Кстати: линии оформления выше выглядят уже одинаково (меняются только координаты, выравнивание и текст). Есть ли способ объединить/оптимизировать это.

Заранее благодарю за помощь!

Answer 1

Кажется, когда значение x становится "большим" отрицательным числом, оно ломается. Если изменить нижний предел кривой leftна -6,5, а для rightкривой на -4,5, то здесь все отлично компилируется.

\documentclass{scrartcl}
\usepackage{tikz}
\usepackage[per-mode = fraction]{siunitx}
\usetikzlibrary{datavisualization.formats.functions}
\begin{document}
  \begin{tikzpicture}
    \datavisualization[
                       scientific axes = {clean, end labels},
                       all axes = {ticks and grid = {major at = 0}},
                       x axis = {label = $\frac{H}{\si{\A\per\m}}$},
                       y axis = {label = $\frac{B(H)}{\si{\tesla}}$},
                       data/format = function,
                       visualize as smooth line/.list = {left, right}
                      ]
    data[set = left] {
      var x : interval [-6.5 : 7];
      func y = 5 / (1 + exp(-1.7 * \value x - 1.5)) - 2.5;
    }
    data[set = right] {
      var x : interval [-4.5 : 7];
      func y = 5 / (1 + exp(-1.7 * \value x + 1.5)) - 2.5;
    }
    %accents important points
    info {
      \draw (visualization cs:x = {(-1.5 / 1.7)}, y = 0) circle [radius = 1pt]
        node [left, font = \footnotesize] {$H_C$};
    }
    info {
      \draw (visualization cs:x = 0, y = 1.59) circle [radius = 1pt]
        node [left, font=\footnotesize] {$B_R$};
    }
    info {
      \draw (visualization cs:x = 7, y = 2.5) circle [radius = 1pt]
        node [above, font=\footnotesize] {$B_S$};
    };
  \end{tikzpicture}
\end{document}