Função definida por partes mal colorida

Question

Eu gostaria de traçar isso

\documentclass{standalone}
\usepackage{tikz}
\usepackage{pgfplots}
\usetikzlibrary{calc}
\pgfplotsset{compat=newest}% <- set a compat!! (current version is 1.14)
%%%%%%%
\pgfkeys{/pgf/declare function={argsinh(\x) = ln(\x + sqrt(\x^2+1));}}
\pgfkeys{/pgf/declare function={argcosh(\x) = ln(x + sqrt(-1 + x)*sqrt(1 + x));}}
%%%%%%
\begin{document}
\begin{tikzpicture}
\begin{axis}[   
    cycle list name= color list,    
    tick align=center,
    axis lines = left,
    xmin = 0,
    xmax = 2,
    ymin = 0,
    ymax = 1.1,   
    ylabel = {$\mathopen|H(j\Omega)\mathclose|$},    
    xtick = {1,1.2},
    xticklabels = {1, $\Omega_s$},
    ytick = {0,0.1, 0.5, 1},
    yticklabels = {0, $\frac{1}{1+\varepsilon^2/k_1^2}$,$\frac{1}{1+\varepsilon^2}$,1},  
    ]

    \foreach \n in {5,6} {
    \addplot+[domain=0:1,samples=201]{1/(1+1^2*(cos(\n*(acos(\x))))^(2))};

   \addplot+[domain=1:2,samples=201]{1/(1+1^2*(cosh(\n*(argcosh(\x))))^(2))};
    }

    \addplot[black,dotted] coordinates {
        (1.2,1)
        (1.2,.1)
        (2,.1)
    };
    \addplot[black,dotted] coordinates {
        (0,.5)
        (1,.5)
        (1,0)
    };
\end{axis}
\end{tikzpicture}
\end{document}

Mas tenho que evitar essas quatro cores e usar apenas duas para ambas \n. Então percebi que isso pode ser feito definindo uma nova função por partes única

%%%%%%
\begin{tikzpicture}[
\declare function={
chebyshev(\x)= (\x<=1) * (1/(1+1^2*(cos(\n*(acos(\x))))^(2)))+
     and(\x>1) * (1/(1+1^2*(cosh(\n*(argcosh(\x))))^(2)));
     }]
%%%%%% ...

\foreach \n in {5,6} {

    \addplot+[domain=0:2,samples=400]{chebyshev(\x)};

}

Mas sempre dá erro e simplesmente não consigo ver onde está o problema. Pode me ajudar?

Answer 1

O problema é que acos(x)só está definido para x \in [0,1]e você está tentando plotar para x \in [0,2]. A coisa mais fácil a fazer seria manter a variante de duas funções e fazer backup do índice de cores em um antes de plotar o segundo usando \pgfplotsset{cycle list shift=-1}.

\documentclass{standalone}
\usepackage{tikz}
\usepackage{pgfplots}
\usetikzlibrary{calc}
\pgfplotsset{compat=newest}% <- set a compat!! (current version is 1.14)
%%%%%%% 
\pgfkeys{/pgf/declare function={argsinh(\x) = ln(\x + sqrt(\x^2+1));}}
\pgfkeys{/pgf/declare function={argcosh(\x) = ln(x + sqrt(-1 + x)*sqrt(1 + x));}}
%%%%%% 
\begin{document}
\begin{tikzpicture}
  \begin{axis}[   
    cycle list name= color list,    
    tick align=center,
    axis lines = left,
    xmin = 0,
    xmax = 2,
    ymin = 0,
    ymax = 1.1,   
    ylabel = {$\mathopen|H(j\Omega)\mathclose|$},    
    xtick = {1,1.2},
    xticklabels = {1, $\Omega_s$},
    ytick = {0,0.1, 0.5, 1},
    yticklabels = {0, $\frac{1}{1+\varepsilon^2/k_1^2}$,$\frac{1}{1+\varepsilon^2}$,1},  
    ]

    \foreach \n in {5,6} {
      \addplot+[domain=0:1,samples=201]{1/(1+1^2*(cos(\n*(acos(\x))))^(2))};
      \pgfplotsset{cycle list shift=-1}
      \addplot+[domain=1:2,samples=201]{1/(1+1^2*(cosh(\n*(argcosh(\x))))^(2))};
    }

    \addplot[black,dotted] coordinates {
      (1.2,1)
      (1.2,.1)
      (2,.1)
    };
    \addplot[black,dotted] coordinates {
      (0,.5)
      (1,.5)
      (1,0)
    };
  \end{axis}
\end{tikzpicture}
\end{document}