FrontPage

2011/06/16からのアクセス回数 2913

ここで紹介したSageワークシートは、以下のURLからダウンロードできます。

http://www15191ue.sakura.ne.jp:8000/home/pub/4/

また、Sageのサーバを公開しているサイト(http://www.sagenb.org/, http://www15191ue.sakura.ne.jp:8000/)にユーザIDを作成することで、ダウンロードしたワークシートを アップロードし、実行したり、変更していろいろ動きを試すことができます。

グラフの使い方

Sageでのグラフの使い方について、説明します。 レファレンスマニュアル を参考にしながら見てください。

3次元多項式をもう一度

sage/text/Sageを使ってみよう でプロットした3次元多項式について、plot関数の使い方も含めて詳しくみてみましょう。

plotの呼び出しは、以下の形式で覚えると便利です。

plot(関数, [変数名, 最小, 最大], オプション)

オプションは、省略可能です。plotのオプションは、plot.optionsで知ることができます。 それ以外にもGraphicsのオプションも使えます。よく使うオプションを以下にしめします。

  • 描画範囲指定のxmin, xmax, ymin, ymax
  • グラフの比率を指定するaspect_ratio
  • 線の色指定のrgbcolor

3次元多項式の場合には、最初にプロットで使用する変数xをvar関数で定義します。 次に関数fを定義します。今回はplotの結果をf_pltに代入していますので、最後に show関数でf_pltを表示します。

f_plt変数に代入することによって、後で他のグラフと重ね合わせて表示することが できます。

sageへの入力:

x = var('x')
f = x^3 - x^2 - 2*x
f_plt = plot(f, [x, -2.5, 2.5]); show(f_plt)

out1.png

plot.optionsを知るには、plot?でヘルプを表示したり、plot_optionsを表示するとよいでしょう。

sageへの入力:

# plotのオプションを知る
plot.options
{'fillalpha': 0.5, 'detect_poles': False, 'plot_points': 200, 
'thickness': 1, 'alpha': 1, 'adaptive_tolerance': 0.01, 'fillcolor': 
'automatic', 'adaptive_recursion': 5, 'exclude': None, 'legend_label': None, 
'rgbcolor': (0, 0, 1), 'fill': False}

接線を求める

点\((x_0, y_0)\)での接線の傾きは、関数fの微分の\(x_0\)での値から得ることができますので、 接線の式は、以下のように求めることができます。 $$ (y - y_0) = f'(x_0) (x - x_0) $$ $y_0$の値は、$f(x_0)$ですから、接線の式は以下のようになります。 $$ y = f'(x_0) (x - x_0) + f(x_0) $$

この計算をSageを使って計算してみましょう。まずSageの変数x0を定義し、y0にf(x=x0)の値を代入します。 f1にfの微分をセットします。

sageへの入力:

# (x0, y0)での接線を求める
x0 = var('x0')
y0 = f(x=x0)
f1 = diff(f, x); view(f1)

out2.png

材料がそろったので、接線の式を定義し表示します。確認のために\(x_0 = 0\)での式も表示します。

\(x_0 = 0\)での接線は、原点を通り傾き-2と求まりました。

sageへの入力:

# x0 = 0での接線の式
y = f1(x=x0)*(x - x0) + y0; view(y)
view(y(x0 = 0))

out3.png

この直線をプロットしてみましょう。 x0=0での接線の式なので、plotの引数にはyではなく、y(x0=0)が渡されていることに注意してください。

重ね合わせたときに区別できるように線の色を緑にします。線の色は、rgbcolorオプションで指定します。

sageへの入力:

y_plt = plot(y(x0=0), [x, -2.5, 2.5], rgbcolor='green')
show(y_plt)

out4.png

接線を計算したポイントを示すために、point関数を使って点を表示します。

point関数の使い方は、以下の通りです。座標は、リストまたはタプル形式で与えます。

point(座標, オプション属性(pointsize, rgccolor, faceted等))

sageへの入力:

# x = 0のポイントをセット
pt = (0, f(x=0))
pt_plt = point(pt, rgbcolor='red', pointsize=30)

グラフの重ね合わせ

Sageのグラフ表示機能の最大の特徴は、重ね合わせです。 これまで計算した以下の結果を同時に表示してみましょう。

  • f_plt: 3次元多項式(青)
  • y_plt: 接線(緑)
  • pt_plt: 接点(赤)

重ね合わせは至って簡単で上記の3つの変数を足し合わせて、showメソッドで表示するだけです。

sageへの入力:

(f_plt + y_plt + pt_plt).show()

out5.png

ユーザ指定(インタラクティブ)

Sageのノートブックでのもう一つの特徴にユーザが逐次値を指定することができる、 インタラクティブ機能があります。

接点の座標x0をユーザ任意に指定できるように、_MyGraph関数を定義して、 与えられたx1でのグラフの重ね合わせを返すようにします。

sageへの入力:

# 与えられたx座標でのfと接線のグラフを返す関数を定義
def _MyGraph(x1):
    y_plt = plot(y(x0=x1), [x, -2.5, 2.5], rgbcolor='green')
    pt = (x1, f(x=x1))    # x0, y0の座標をセット
    pt_plt = point(pt, rgbcolor='red', pointsize=30)
    return (f_plt + y_plt + pt_plt) # f_pltは先の結果を再利用

インタラクティブ機能を使う場合、

@interact	

の後に関数を定義します。ここでは、_(アンダーバー)という特殊な関数名を指定します。

関数の引数で与えられた変数をユーザが入力することができるようになります。 ここでは、slider関数を使って-2.5から2.5の範囲で1/10の刻みのスライダー バーを表示してx1を指定できるようにします。

関数の内容は、_MyGraph関数でx1で指定された座標の接線とポイントの重ね合わせ グラフを返し、showメソッドでグラフを表示します。グラフ表示範囲がx1の値でかわらない ようにymin, ymaxで範囲を指定します。

sageへの入力:

@interact
def _(x1 = slider(-2.5, 2.5, 1/10)):
    graph = _MyGraph(x1)
    graph.show(ymin=-20, ymax=10)

out6.png

アニメーション

グラフの変化を分かりやすく表示する方法として、アニメーション機能も提供しています。

アニメーションでは、グラフのリストを作成し、それをanimate関数に渡し、描画範囲を xmin, xmax, ymin, ymaxで指定するだけで簡単にアニメーションが作成できます。 作成されたアニメーションの表示もshowメソッドを使います。

sageへの入力:

g_list = [_MyGraph(x1) for x1 in srange(-2.5,2.5,0.2)]
a = animate(g_list, xmin=-2.5, xmax=2.5, ymin=-20, ymax=10)
a.show()

out7.gif

並べて表示

Sageのプロット結果は、縦に表示されてしまうため、結果を比較したときには不便です。

このような場合には、htmlコマンドにプロットしたいグラフのリストを渡すと、 表形式に変換して表示してくれます。表示するグラフのオプションにfigsizeでグラフの 大きさを小さくしておくとよいでしょう。

以下にsin曲線とcos曲線の表示例を示します。

sageへの入力:

# グラフを並べて表示する
sin_plt = plot(sin(x), [x, -pi, pi], rgbcolor= 'green', figsize=4)
cos_plt = plot(cos(x), [x, -pi, pi], figsize=4)
html.table([['sin曲線', 'cos曲線'], [sin_plt, cos_plt]])

out8.png

リストプロット

データを表示するのに便利なのが、list_plot関数です。

list_plot関数は、以下のように使用します。

list_plot(プロットするリスト, オプション)	

リストの要素に[x, y]形式で座標を指定すると分布図が表示されます。 また、plot関数でエラーで表示できない場合の代替手段としてlist_plot関数を使うこともあります。

以下にデータ点のプロットの場合とデータを結ぶ場合の例を示します。

sageへの入力:

# データを点でプロット
list_plot([1, 2, 4, 3, 6])

out9.png

sageへの入力:

# データを結ぶ
list_plot([1, 2, 4, 3, 6], plotjoined=True)

out10.png

異なるデータを同時に表示する場合には、rgbcolorオプションで色をセットしたり、 legen_labelでデータの種類をセットすると便利です。

sageへの入力:

# データ(座標)のプロット
c1 = [[1,2],[1,4],[2,4]]
c2 = [[2,1],[5,1],[4,2]]
# プロットして分布を確認
pl1 = list_plot(c1, rgbcolor='red', legend_label='c1')
pl2 = list_plot(c2, rgbcolor ='blue', legend_label='c2')
(pl1+pl2).show(xmin=0, xmax=5, ymin=0)

out11.png

関係を表示

プロットしたい関数ではなく、条件などの関係式を表示したいときがあります。 このような場合には、implicit_plotを使うと便利です。

implicit_plotの例として、単位円を表示してみます。 単位円の条件は、以下のように与えられます。 $$ x^2 + y^2 = 1 $$

これをimplicit_plotで表示すると以下のようになります。

sageへの入力:

# 関係を満たすグラフを表示する
x, y = var('x y')
i_plt = implicit_plot(x^2 + y^2 == 1, [x, -1.5, 1.5], [y, -1.5, 1.5])
i_plt.show(aspect_ratio=1)

out12.png

パラメトリックプロット(媒介変数表示)

x, yが時間変数t(パラメータ)によって表される場合、parametric_plot関数 を使って表示すると簡単です。

パラメトリックプロットの例として、以下のサイクロイド曲線をparametric_plot関数 で表示してみます。

$$ \begin{eqnarray} x = 2 (t - sin(t)) \\ y = 2 (1 - cos(t)) \end{eqnarray} $$

sageへの入力:

# 媒介変数表示
t = var('t')
x = 2*(t-sin(t))
y = 2*(1-cos(t))
# サイクロイドのプロット(パラメトリックプロットの例)
parametric_plot([x, y], (t, 0, 2*pi))

out13.png

3次元グラフ

3次元のグラフには、plot3d関数を使用します。重ね合わせは2次元のグラフと同じように使えます。

表示された3次元図形はマウスで自由に拡大、回転することができます。(驚きました)

sageへの入力:

x, y = var('x y')
plot3d(sin(x*y),(x,-pi,pi),(y,-pi,pi), mesh=True)

out14.png

等高線図

3次元情報を表示する場合によく使用する方法に等高線図があります。

等高線図は、contour_plot関数を使って表示します。 前のグラフと同じものをcontour_plot関数で表示すると以下のようになります。

sageへの入力:

var('x y')
contour_plot(sin(x*y), [x, -pi, pi], [y, -pi, pi], aspect_ratio=1)

out15.png

基本図形

以下によく使う基本図形を以下に示します。

  • 円: circle
  • 文字列: text
  • 線: line
  • 点: point
  • ポリゴン: polygon

円は以下のように表示します。

circle((座標), 半径)

circleの例を以下に示します。座標は原点(0, 0)、半径は1です。

sageへの入力:

circle((0,0), 1)

out16.png

文字列

次に文字列textです。表示する文字列には、$で囲んでlatexの数式を表示することもできます。 残念ながら、日本語が含まれているとエラーまたは文字化けします。

textは以下の形式で使用します。文字列の中心が指定した座標になるようにプロットされます。

text(文字列, (座標))

textの例を以下に示します。

sageへの入力:

test_plt = text('test', (1, 1))
eq_plt = text('$f(x) = x^2 + 1$', (0.5,0.5))
(test_plt + eq_plt).show(xmin=0, ymin=0)

out17.png

線(line)は、指定された座標のリストを線で結びます。 lineの使い方は、簡単です。

line([(開始座標), (終了座標)])

以下に例を示します。日本語が使えないため、タイトルをhtml関数で先に表示し、その後に図形を表示します。 show関数で表示領域を指定し、すべての図形が表示されるようにします。

sageへの入力:

html('<center>テスト</center>')
c = circle((0.5,0.5), 1)
l = line([(0,0), (1, 1)])
pt = point((0.5, 0.5), rgbcolor='white', pointsize=30, faceted=True)
(c + l + pt).show(xmin=-1, xmax=2)

out18.png

ポリゴンの塗りつぶし

polygon関数を使うとリストで指定した座標の図形を塗りつぶします。

sageへの入力:

polygon([(0,0), (1,1), (0,1)])

out19.png

コメント

選択肢 投票
おもしろかった 3  
そうでもない 0  
わかりずらい 0  

皆様のご意見、ご希望をお待ちしております。

  • p=plot(sin,(-2*pi,2*pi); p.save('demo.png') で保存した図をノートブックに読み込む方法を教えてください。よろしくお願いします。 -- [[ysato ]] 2014-09-24 (水) 23:33:48
  • p=plot(sin,(-2*pi,2*pi)); 最後の括弧が抜けてしまったので訂正します。 -- [[ysato ]] 2014-09-24 (水) 23:36:30
  • ysatoさま、私はp.save(DATA + 'demo.png')で保存し、html('<img src="demo.png"/>)で表示しています。 -- 竹本 浩? 2014-09-26 (金) 23:02:54
  • 'sage/text/Rとの連携' で定義された関数 postGraph に使われていることを確認しました。ありがとうございました。 -- ysato? 2014-09-28 (日) 22:44:59
  • h4>等高線図</h4>、h4>線</h4> の表示があります。小さなことですみません。 -- ysato? 2014-09-29 (月) 12:11:09
  • ysatoさま、ご指摘ありがとうございます。修正しました。 -- 竹本 浩? 2014-10-02 (木) 00:17:10

(Input image string)


添付ファイル: fileout19.png 454件 [詳細] fileout18.png 429件 [詳細] fileout17.png 454件 [詳細] fileout16.png 447件 [詳細] fileout15.png 484件 [詳細] fileout14.png 455件 [詳細] fileout13.png 471件 [詳細] fileout12.png 489件 [詳細] fileout11.png 459件 [詳細] fileout10.png 449件 [詳細] fileout9.png 501件 [詳細] fileout8.png 481件 [詳細] fileout7.gif 477件 [詳細] fileout6.png 501件 [詳細] fileout5.png 454件 [詳細] fileout4.png 505件 [詳細] fileout3.png 458件 [詳細] fileout2.png 513件 [詳細] fileout1.png 490件 [詳細]

トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2014-10-07 (火) 09:27:22 (902d)
SmartDoc