pgfplots 3D с логарифмической осью Y

pgfplots 3D с логарифмической осью Y

Мне нужна логарифмическая ось Y для моего 3D-графика, но почему-то кажется, что сгенерированный график также находится в отрицательном диапазоне Y (не метки осей, а область Y). Область Y, похоже, равна -1:1 (на самом деле это так, но я думал, что это должно работать из-за "10^y" -> область Y=0.1:10 (см. комментарий в коде)).

Я использовал этот пост для своего кода: 3D-график поверхности с логарифмическими осями x и y

Участок 1 - линейная ось Y, Участок 2 - логарифмическая ось Y (попытка)

введите описание изображения здесь

МВЭ:

\documentclass{article}
\usepackage[utf8]{inputenc}
\usepackage{pgfplots}
\usepackage{tikz}

\begin{document}

\begin{figure}[h!]
\begin{tikzpicture}
\begin{axis} 
[scale = 1,
xlabel = {$\sigma$},
xmin = -10, xmax = 10,
ylabel = {$\omega$},
ymin = 0.1, ymax = 10,
zlabel={$A$},
zmin = -40, zmax = 20,
colormap/viridis]

\addplot3[
surf,
samples=40,
domain=-1:10,
domain y=0:10]
(x, y, {20*log10(1/sqrt(((1.40845*y)+(2*x*y))^2+(1+(1.40845*x)+x^2-y^2)^2))});
\end{axis}
\end{tikzpicture}
\caption{Plot 1}
\end{figure}


\begin{figure}[h!]
\begin{tikzpicture}
\begin{axis} 
[scale = 1,
xlabel = {$\sigma$},
xmin = -10, xmax = 10,
ylabel = {$\omega$},
ymin = 1e-1, ymax = 1e1,
zlabel={$A$},
zmin = -40, zmax = 20,
ymode=log,        %added ymode
colormap/viridis]

\addplot3[
surf,
samples=40,
domain=-1:10,
domain y=-1:1]    %new y-domain (10^y in the next line) 
(x, 10^y, {20*log10(1/sqrt(((1.40845*y)+(2*x*y))^2+(1+(1.40845*x)+x^2-y^2)^2))});
\end{axis}
\end{tikzpicture}
\caption{Plot 2}
\end{figure}

\end{document}

Я ожидаю чего-то вроде этого (график создан в Mathematica — с логарифмической осью Y):

введите описание изображения здесь

решение1

Я думаю, вы думали, что это слишком сложно. Просто адаптируйтесь domain yв соответствии с ymin/ ymaxи вы получите желаемый результат ...

% used PGFPlots v1.16
\documentclass[border=5pt]{standalone}
\usepackage{pgfplots}
    \pgfplotsset{compat=1.16}
\begin{document}
\begin{tikzpicture}
    \begin{axis}[
        view={50}{50},
        width=\axisdefaultwidth,
        height=5cm,
        xlabel={$\sigma$},
        xmin=-1, xmax=10,
        ylabel={$\omega$},
        ymin=1e-1, ymax=1e1,
        zlabel={$A$},
        zmin=-40, zmax=20,
        ymode=log,
        colormap/viridis,
]
        \addplot3[
            surf,
            samples=40,
            domain=-1:10,
            domain y=0.1:10,    % <-- adapted
        % (no need to use a parametric plot)
        ] {20*log10(1/sqrt(((1.40845*y)+(2*x*y))^2+(1+(1.40845*x)+x^2-y^2)^2))};
    \end{axis}
\end{tikzpicture}
\end{document}

изображение, показывающее результат кода выше

решение2

Я думаю, что если вы подставите каждый y на 10^y в уравнение, вы можете получить более похожий результат, посмотрите это:

\documentclass{article}
\usepackage[utf8]{inputenc}
\usepackage{pgfplots}
\usepackage{tikz}

\begin{document}

\begin{tikzpicture}
\begin{axis} 
[scale = 1,
xlabel = {$\sigma$},
xmin = -10, xmax = 10,
ylabel = {$\omega$},
ymin = 1e-1, ymax = 1e1,
zlabel={$A$},
zmin = -40, zmax = 20,
ymode=log,        %added ymode
view={60}{45},
colormap/viridis]

\addplot3[
surf,
samples=100,
domain=-1:9,
domain y=-1:1]    %new y-domain (10^y in the next line) 
(x, 10^y, {20*log10(1/sqrt(((1.40845*(10^y))+(2*x*(10^y)))^2+(1+(1.40845*x)+x^2- 
(10^y)^2)^2))});
\end{axis}
\end{tikzpicture}

\end{document}

введите описание изображения здесь

Закругленная форма в нижней части оранжевой поверхности может быть способностью Mathematica обрезать поверхность и делать более непрерывные графики. В качестве альтернативы, используя tikz, вы можете достичь желаемой формы поверхности, выполняя вычитания других математических функций, таких как F(x,y,z)= G(x,y,z)-h(x,y,z)

Связанный контент