しままるの雑記帳

しままるの思ったことや考えたことをアウトプットする場。「PC/プログラミング」「写真/旅行」「科学(勉強)」「心の問題」が主なテーマ。

プログラミング

gnuplotでハートの形を描く

投稿日:2017年5月6日 更新日:

gnuplotにハートの関数を入れて描画させよう

今回は、対話型グラフ描画ソフトgnuplotで「ハートの形を描く」ということをやってみようと思います。

ハート関数

ハートの形を描画できる関数は、以下に示すような陰関数です。
参考サイト:プログラムモグモグ|美しすぎるハートの数式(こちら

gnuplotで陰関数を疑似的に描く

残念ながらgnuplotで陰関数を描くことはできません。そこで、上記の式を

として、3次元プロットでz=0における等高線を描けば、ハートの形が描けそうです。つまり、3次元プロットを使って、疑似的に陰関数を描いています。やってみましょう。

※サンプルプログラム ここから※
f(x,y)= (x**2+y**2-1)**3-x**2*y**3
len = 3 #座標軸長さ定義
set xrange [-len:len]
#x軸の長さは-3~3です。(len = 3)
set yrange [-len:len]
#y軸の長さは-3~3です。(len = 3)
set zrange [-len:len]
set isosample 200,200
#立体のメッシュの細かさはこのくらいで
set ticslevel len/(-len-len)
#z軸におけるx-y平面の位置はこのくらいで(詳細はこちら)
set contour base
#等高線を表示してください
set cntrparam levels discrete 0.0
#z=0における等高線を表示してください
set nosurface
#表面の表示は要りません
set view equal xyz
#xy軸の比率は同じにしてください
unset border
#グラフの枠は要りません
set xtics axis
#x軸の目盛りはx軸上に書いてください。
set ytics axis
#y軸の目盛りはx軸上に書いてください。
set ztics axis
#y軸の目盛りはx軸上に書いてください。
unset xtics
#x軸の目盛りは要りません。
unset ytics
#y軸の目盛りは要りません。
unset ztics
#y軸の目盛りは要りません。
set zeroaxis ls -1
set view 0,0
#x軸方向に0°、z軸方向に0°傾けた形で描画してください。(xy平面を描画してください)

splot f(x,y) lt 2 lw 2 lc rgb ‘red’

set nokey
set title font “Times New Roman,25”
set title “heart-function”

#pause -1
※サンプルプログラム ここまで※

上記サンプルプログラムを実行します。

描けました。アルゴリズム的には

で表される平面を一旦描いて、そのz=0における等高線を表示させて、平面の表示を消しています。

詳細説明

はこんな形をしています。縦軸の長さは3で、x,y,z軸の交点は原点(0,0,0)です。

ここで例えばz=2における等高線を描くならば

こんな形になりましたz=0の等高線を描けば

このようにハートの形が描けるのです。

できなかったこと

ハートの色を赤色にしたかったのですが、残念ながらできませんでした。等高線の色を変えることができれば実現できるかと思っています。gnuplotに詳しい方、やり方を知っている方がいらっしゃれば、ご教示のほどお願いいたします。

-プログラミング

執筆者:


comment

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です

関連記事

初心者向けLinux入門:ディストリビューションとは?

目次1 はじめに2 ディストリビューションとは=頒布3 Linuxの種類≒ディストリビューション4 どのディストリビューションがいいか?5 まとめ Linuxディストリビューションの意味について解説 …

初心者向けLinux入門:ubuntuのダウンロード

簡単図解!Ubuntuをダウンロードしよう! ここでは、初心者向けLinux入門(下リンク)の一環として、ubuntuのダウンロード方法について説明したいと思います。 初心者向けLinux入門:Lin …

イントロクイズの作り方 ~soxの使い方を解説~

目次1 soxとは?1.1 対話型音声ファイル編集ソフト sox1.2 コマンドで編集を進める意味2 soxのコマンド2.1 切り取り2.2 フェードイン・フェードアウト2.3 その他3 イントロクイ …

[備忘録] gnuplotのダウンロード・インストール方法と使い方

目次1 インストール2 対話的使用法3 スクリプトを使った方法4 gnuplot実践4.1 描きたい式4.2 gnuplotにおける極座標4.3 描画5 まとめ 対話型グラフ描画ソフト gnuplot …

[備忘録] htmlにおけるベージ内ジャンプの方法

目次1 ページ内リンクのプログラム2 サンプルプログラム2.1 プログラムの説明3 動作確認4 使用例5 応用例5.1 応用例 動作確認 ページ内リンクのプログラム  ブログを書いていて、とりあえず結 …