如何画伯努利双纽线?

发布于 2021-09-08 13:08:58

我的参考代码如下:采用的极坐标的思路

\documentclass[tikz]{standalone}
\usetikzlibrary{arrows.meta}
\begin{document}
\begin{tikzpicture}[samples=200,scale=1.2]
\draw[-Stealth] (-2.5,0)--(2.5,0) node[below] {$x$};
\draw[-Stealth] (0,-1)--(0,1) node[left] {$y$};
\node at (-3pt,-3pt) {$o$};
\draw[domain=0:2*pi] plot (xy polar cs:angle=\x r,radius={(2*cos(2*\x r))^(1/2)}); 
\node[below] at (1.414,0) {$a$};
\node[below] at (-1.414,0) {$-a$};
\node[align=center] at (0,-1.5) {$\rho^2=a^2\cos2\theta$};
\end{tikzpicture}
\end{document}

报错是这样的

Package PGF Math Error: I cannot calculate the logarithm of `-0.01544' (in '(2*cos(2*0.78926r))^(1/2)'). ...:angle=\x r,radius={(2*cos(2*\x r))^(1/2)})

TikZ 能不能画隐函数的方程?

$(x^2+y^2)^2=a^2(x^2-y^2)$

查看更多

关注者
0
被浏览
3.8k
寄存器
寄存器 2021-09-08
这家伙很懒,什么也没写!

可以用tkz-fct宏包进行函数曲线绘制:

\documentclass[margin=5pt]{standalone}

\usepackage{tkz-fct}

\begin{document}

\begin{tikzpicture}
  \tkzInit [xmin=-8,xmax=8,ymin=-8,ymax=8,xstep=1,ystep=1]
  \tkzFctPolar[domain=0:2*pi,samples=1000]%
  { sqrt(2*cos(2*t)) }
\end{tikzpicture}

\end{document}

基结果为:
2021-09-08_15-55.png


TikZ能可以画隐函数的方程,做法是构造一个f(x,y)=(x^2+y^2)^2-a^2(x^2-y^2)函数,然后进行绘制,如:

\documentclass{standalone}

\usepackage{tikz}

\begin{document}
  \begin{tikzpicture}[xscale=1,yscale=.25]
    \draw plot[id=curve, raw gnuplot] function{
      f(x,y) = ((x**2 + y**2)**2) - 2*(x**2 - y**2);
      set xrange [-8:8];
      set yrange [-8:8];
      set view 0,0;
      set isosample 2000,2000;
      %set size square;
      set cont base;
      set cntrparam levels incre 0,0.1,0;
      unset surface;
      splot f(x,y)
    };
  \end{tikzpicture}
\end{document}

其结果为:
2021-09-08_15-54.png


好像画出来的不一样,不知对与不对?


注意,由于使用了外部gnuplot进行计算,所以需要安装gnuplot工具,并在编译时使用--shell-escape参数。

5 个回答
其实有点卡
其实有点卡 2021-09-08
这家伙很懒,什么也没写!

利用伯努利双纽线的极坐标表达式,弄了一个简单例子。

\documentclass{standalone}

\usepackage{tikz}

\begin{document}
  \begin{tikzpicture}% 注意定义域
    \draw[->] (-3,0) -- (3,0);
    \draw[->] (0,-2) -- (0,2);
    \draw[domain=-pi/4:pi/4,samples=400] plot ({deg(\x)}:{2*sqrt(cos(2*deg(\x)))});
    \draw[domain=-pi/4:pi/4,samples=400] plot ({deg(\x)}:{-2*sqrt(cos(2*deg(\x)))});
  \end{tikzpicture}
\end{document}


\begin{tikzpicture}
  \coordinate (O) at (0,0);
  \draw[->] (-3,0) -- (3,0);
  \draw[->] (0,-2) -- (0,2);
  \draw[domain=-pi/4:pi/4,samples=400]
    plot ({1.5*sqrt(2*cos(2*deg(\x)))*cos(deg(\x))},
          {1.5*sqrt(2*cos(2*deg(\x)))*sin(deg(\x))});
  \draw[domain=-pi/4:pi/4,samples=400]
    plot ({-1.5*sqrt(2*cos(2*deg(\x)))*cos(deg(\x))},
          {-1.5*sqrt(2*cos(2*deg(\x)))*sin(deg(\x))});
\end{tikzpicture}%%% 双纽线参数方程,这个画出来有一些问题

image.png

画曲线的代码有两行,感觉不太完美,主要是为了避免给负数开方。

伯努利双纽线的极坐标方程很奇怪,感觉有某些约定。

-XM-
-XM- 2021-09-08
这家伙很懒,什么也没写!

分左右两段去画即可(抄自https://tex.stackexchange.com/questions/299218/cassini-oval-curve-in-tikz/299225#299225

\documentclass{standalone}
\usepackage{tikz}
\begin{document}
\begin{tikzpicture}[
        x=2cm,
        y=2cm,
    ]
    \def\samples{100}
    \def\c{1}
    \pgfmathsetmacro\cc{\c*\c}
    \pgfmathsetmacro\cccc{\cc*\cc}
    \def\a{1.01}
    \draw (0, 0) --
    plot[
            variable=\t,
            domain=-45+1/\samples:45-1/\samples,
            samples=\samples,
            smooth,
        ]
    (\t:{\cc*sqrt(\cc*cos(2*\t) + \cc*cos(2*\t))})
    -- cycle;
    \draw (0, 0) --
    plot[
            variable=\t,
            domain=180-45+1/\samples:180+45-1/\samples,
            samples=\samples,
            smooth,
        ]
    (\t:{\cc*sqrt(\cc*cos(2*\t) + \cc*cos(2*\t))})
    -- cycle;
    %\path plot[mark=x] coordinates { (-\c, 0) (\c, 0) };
\end{tikzpicture}
\end{document}
零蛋大
零蛋大 2021-09-08
知乎主页:https://www.zhihu.com/people/tangshaodong

向老师的 tikz 绘图笔记二 里看到了方案, 也是分段解决, 代码如下

\documentclass[margin=5pt]{standalone}
\usepackage{tikz}
\usetikzlibrary{arrows.meta}
\begin{document}
\begin{tikzpicture}[samples=200,scale=1.2]
\draw[-Stealth] (-2,0)--(2,0) node[below] {$x$};
\draw[-Stealth] (0,-0.6)--(0,0.8) node[left] {$y$};
\draw[domain=-pi/4:pi/4] plot ({1.414*sqrt(cos(2*\x r))*cos(\x r)},{1.414*sqrt(cos(2*\x r))*sin(\x r)})--(0,0);
\draw[domain=-pi/4:pi/4,rotate=180] plot ({1.414*sqrt(cos(2*\x r))*cos(\x r)},{1.414*sqrt(cos(2*\x r))*sin(\x r)})--(0,0);
\node at (-4pt,-9pt) {$o$};
\node[below] at (-1.65,0) {$-a$};
\node[below] at (1.55,0) {$a$};
\node[align=center] at (0,-1.1) {$(x^2+y^2)^2=a^2(x^2-y^2)$\\
$\rho^2=a^2\cos2\theta$};
\end{tikzpicture}
\end{document}
越来越少
越来越少 2021-09-09
这家伙很懒,什么也没写!
\tikz{%r^2=a^2*cos(2*i)
  \def\a{2}
  \def\tempsave{}
  \foreach \i in {-45,-44.5,...,45,135,135.5,...,225}
  {
    \pgfmathparse{sqrt(\a^2*cos(2*\i))}
    \xdef\tempsave{\tempsave(\i:\pgfmathresult)}
  }
  \draw plot[smooth] coordinates {\tempsave};
}

cb3e54f03f9a696c85cc4ab5325dbe78.png

撰写答案

请登录后再发布答案,点击登录

发布
问题

分享
好友

手机
浏览

扫码手机浏览