2018年12月23日日曜日

開発環境

気になった時はとりあえずグラフの描画してみる。ということで描画してみた。

コード(Emacs)

Python 3

#!/usr/bin/env python3
from sympy import pprint,  symbols, I, sqrt, pi, plot, sin

x = symbols('x', real=True)
f = I / 2 * (((1 - sqrt(3) * I) / 2) ** (3 * x / pi) -
             ((1 + sqrt(3) * I) / 2) ** (3 * x / pi))


real, imag = f.as_real_imag()

for i, t in [('関数', f), ('実部', real), ('虚部', imag)]:
    print(i)
    for s in [t, t.expand(), t.simplify(), t.factor()]:
        pprint(s)
        print()
    print()

p = plot(f, legend=True, show=False)

p.save('sample.png')

p = plot(sin(x), legend=True, show=False)
p.save('sample1.png')

入出力結果(Terminal, cmd(コマンドプロンプト), Jupyter(IPython))

$ ./sample.py
関数
  ⎛          3⋅x             3⋅x⎞
  ⎜          ───             ───⎟
  ⎜           π               π ⎟
  ⎜⎛1   √3⋅ⅈ⎞      ⎛1   √3⋅ⅈ⎞   ⎟
ⅈ⋅⎜⎜─ - ────⎟    - ⎜─ + ────⎟   ⎟
  ⎝⎝2    2  ⎠      ⎝2    2  ⎠   ⎠
─────────────────────────────────
                2                

 -3⋅x              3⋅x    -3⋅x              3⋅x
 ─────             ───    ─────             ───
   π                π       π                π 
2     ⋅ⅈ⋅(1 - √3⋅ⅈ)      2     ⋅ⅈ⋅(1 + √3⋅ⅈ)   
────────────────────── - ──────────────────────
          2                        2           

 -3⋅x    ⎛          3⋅x             3⋅x⎞
 ─────   ⎜          ───             ───⎟
   π     ⎜           π               π ⎟
2     ⋅ⅈ⋅⎝(1 - √3⋅ⅈ)    - (1 + √3⋅ⅈ)   ⎠
────────────────────────────────────────
                   2                    

   ⎛            x             x⎞ ⎛          2⋅x             x           x     
   ⎜            ─             ─⎟ ⎜          ───             ─           ─     
   ⎜            π             π⎟ ⎜           π              π           π     
   ⎜  ⎛1   √3⋅ⅈ⎞    ⎛1   √3⋅ⅈ⎞ ⎟ ⎜⎛1   √3⋅ⅈ⎞      ⎛1   √3⋅ⅈ⎞  ⎛1   √3⋅ⅈ⎞    ⎛1
-ⅈ⋅⎜- ⎜─ - ────⎟  + ⎜─ + ────⎟ ⎟⋅⎜⎜─ - ────⎟    + ⎜─ - ────⎟ ⋅⎜─ + ────⎟  + ⎜─
   ⎝  ⎝2    2  ⎠    ⎝2    2  ⎠ ⎠ ⎝⎝2    2  ⎠      ⎝2    2  ⎠  ⎝2    2  ⎠    ⎝2
──────────────────────────────────────────────────────────────────────────────
                                             2                                

        2⋅x⎞ 
        ───⎟ 
         π ⎟ 
   √3⋅ⅈ⎞   ⎟ 
 + ────⎟   ⎟ 
    2  ⎠   ⎠ 
─────────────
             


実部
    ⎛          3⋅x⎞     ⎛          3⋅x⎞
    ⎜          ───⎟     ⎜          ───⎟
    ⎜           π ⎟     ⎜           π ⎟
    ⎜⎛1   √3⋅ⅈ⎞   ⎟     ⎜⎛1   √3⋅ⅈ⎞   ⎟
  im⎜⎜─ - ────⎟   ⎟   im⎜⎜─ + ────⎟   ⎟
    ⎝⎝2    2  ⎠   ⎠     ⎝⎝2    2  ⎠   ⎠
- ───────────────── + ─────────────────
          2                   2        

    ⎛          3⋅x⎞     ⎛          3⋅x⎞
    ⎜          ───⎟     ⎜          ───⎟
    ⎜           π ⎟     ⎜           π ⎟
    ⎜⎛1   √3⋅ⅈ⎞   ⎟     ⎜⎛1   √3⋅ⅈ⎞   ⎟
  im⎜⎜─ - ────⎟   ⎟   im⎜⎜─ + ────⎟   ⎟
    ⎝⎝2    2  ⎠   ⎠     ⎝⎝2    2  ⎠   ⎠
- ───────────────── + ─────────────────
          2                   2        

 -3⋅x  ⎛    ⎛          3⋅x⎞     ⎛          3⋅x⎞⎞
 ───── ⎜    ⎜          ───⎟     ⎜          ───⎟⎟
   π   ⎜    ⎜           π ⎟     ⎜           π ⎟⎟
2     ⋅⎝- im⎝(1 - √3⋅ⅈ)   ⎠ + im⎝(1 + √3⋅ⅈ)   ⎠⎠
────────────────────────────────────────────────
                       2                        

    ⎛          3⋅x⎞     ⎛          3⋅x⎞
    ⎜          ───⎟     ⎜          ───⎟
    ⎜           π ⎟     ⎜           π ⎟
    ⎜⎛1   √3⋅ⅈ⎞   ⎟     ⎜⎛1   √3⋅ⅈ⎞   ⎟
- im⎜⎜─ - ────⎟   ⎟ + im⎜⎜─ + ────⎟   ⎟
    ⎝⎝2    2  ⎠   ⎠     ⎝⎝2    2  ⎠   ⎠
───────────────────────────────────────
                   2                   


虚部
  ⎛          3⋅x⎞     ⎛          3⋅x⎞
  ⎜          ───⎟     ⎜          ───⎟
  ⎜           π ⎟     ⎜           π ⎟
  ⎜⎛1   √3⋅ⅈ⎞   ⎟     ⎜⎛1   √3⋅ⅈ⎞   ⎟
re⎜⎜─ - ────⎟   ⎟   re⎜⎜─ + ────⎟   ⎟
  ⎝⎝2    2  ⎠   ⎠     ⎝⎝2    2  ⎠   ⎠
───────────────── - ─────────────────
        2                   2        

  ⎛          3⋅x⎞     ⎛          3⋅x⎞
  ⎜          ───⎟     ⎜          ───⎟
  ⎜           π ⎟     ⎜           π ⎟
  ⎜⎛1   √3⋅ⅈ⎞   ⎟     ⎜⎛1   √3⋅ⅈ⎞   ⎟
re⎜⎜─ - ────⎟   ⎟   re⎜⎜─ + ────⎟   ⎟
  ⎝⎝2    2  ⎠   ⎠     ⎝⎝2    2  ⎠   ⎠
───────────────── - ─────────────────
        2                   2        

 -3⋅x  ⎛  ⎛          3⋅x⎞     ⎛          3⋅x⎞⎞
 ───── ⎜  ⎜          ───⎟     ⎜          ───⎟⎟
   π   ⎜  ⎜           π ⎟     ⎜           π ⎟⎟
2     ⋅⎝re⎝(1 - √3⋅ⅈ)   ⎠ - re⎝(1 + √3⋅ⅈ)   ⎠⎠
──────────────────────────────────────────────
                      2                       

 ⎛    ⎛          3⋅x⎞     ⎛          3⋅x⎞⎞ 
 ⎜    ⎜          ───⎟     ⎜          ───⎟⎟ 
 ⎜    ⎜           π ⎟     ⎜           π ⎟⎟ 
 ⎜    ⎜⎛1   √3⋅ⅈ⎞   ⎟     ⎜⎛1   √3⋅ⅈ⎞   ⎟⎟ 
-⎜- re⎜⎜─ - ────⎟   ⎟ + re⎜⎜─ + ────⎟   ⎟⎟ 
 ⎝    ⎝⎝2    2  ⎠   ⎠     ⎝⎝2    2  ⎠   ⎠⎠ 
───────────────────────────────────────────
                     2                     


$

ということで、グラフは正弦と一致。

0 コメント:

コメントを投稿