1

陽に y = f(x) の形で表されている関数のグラフを (x, y) 平面上に描くことを考えます。ただし x_min ≦ x ≦ x_max の範囲に限定してプロットするとします。

素朴には、x_min から x_max まで少しずつ x を変更しながら (x, f(x)) を計算して、それらの点を線分で繋ぐとグラフが描けそうです。

しかしこの方法だと特異点があるグラフに対して上手く動きません。たとえば反比例のグラフ y = 1 / x は x = 0 において定義されていません。上の素朴なアルゴリズムを Scratch で実装して動かしてみると、y 軸付近に本来は存在しない縦線が引かれてしまいます (緑の太線がグラフです)。

y = 1000 / x のグラフを描こうとした実行結果です。y 軸付近にも縦の線が引かれてしまっています。

(以下はこのグラフを描くためのスクリプトです。クリックすると拡大します)

y = 1000 / x のグラフを描くための Scratch のソースコードです。

対して WolframAlpha にグラフを描かせてみると、上手く特異点付近を除いてプロットしています。

WolframAlpha に y = 1 / x のグラフをプロットさせた結果です

実際 Wolfram Mathematica の機能を見ると、特異点を検知して除去している模様です。

また、y = 1 / x 程度であれば手計算で x = 0 が特異点と分かるので手動でそこだけ除去すれば良いのですが、y = tan(x) など特異点が複数現れると面倒です。

質問

関数 f の定義式は入力として与えられるとして、y = f(x) のグラフをプロットするときに上手く特異点の周りをプロットしないようにするアルゴリズムはありますか? 数値的に解析するものでも、定義式から解析するものでも構いません。必要であれば f の形を制限して構いません。個人的には y = 1 / x や y = tan(x) 程度が描ければ満足です。ある程度汎用性を持たせたいので、特異点を手計算して除く以外の方法が知りたいです。

※書きやすかったので上では Scratch で書きましたが、アルゴリズムが知りたいだけなので他のプログラミング言語でも疑似コードでも方針だけでも構いません。

1 件の回答 1

1

描画するときに、1つ前のステップのf(x)の値を持っておき、前ステップとのf(x)の差の絶対値が一定値以上になったら線で結ばないようにしたらどうでしょうか?

1
  • 1
    +1: 実用上はこれが簡単そうに思います。ただこの方法は完璧ではなくて、閾値以上になったものすべてを消してしまうので、傾きの大きな直線や、指数関数の端の方、振幅が大きく周期の短い三角関数でも消えてしまいます。(そのようなグラフが必要となる場合があるかは別問題ですが……。)
    – nekketsuuu
    Commented 2019年3月27日 18:05

この質問に回答するには、ログインする必要があります。

求めていた回答ではありませんか? のタグが付いた他の質問を参照する。