python  ·  28日 4月 2021

    Python+VSCode(StoreのPythonで)(2)

    以前作ったスクリプトも動くかどうか試してみたい。

    で、「numpyでCPFSK再び」でやったやつにしてみる。

    いきなり

      File "d:\Documents\python3-vscode\study2.py", line 2, in <module>

        from IPython import get_ipython

    ModuleNotFoundError: No module named 'IPython'

    と言われるので、

     

    from IPython import get_ipython

    get_ipython().magic('reset -sf')

     

    だったところを、

     

    import sys

    sys.modules[__name__].__dict__.clear()

     

    ってする。

     

    で、

     

      File "d:\Documents\python3-vscode\study2.py", line 5, in <module>

        import matplotlib.pyplot as plt

    ModuleNotFoundError: No module named 'matplotlib'

     

    って言われるので、matplotlibとscipyを入れる。matplotlibを入れたところでnumpyも入っちゃう。Anacondaと違ってpipっていうコマンドで入れる。Storeで入手したPythonは勝手にPATHが通してあるので、コマンドでpipってするだけ。

     

    次は

     

      File "d:\Documents\python3-vscode\study2.py", line 68, in <module>

        np.save('output/fsk_o_t',t)

      File "<__array_function__ internals>", line 5, in save                                         ocalCache\local-packages\Python39\site-packages\numpy\lib\npyio.py", line 525, in save

      File "C:\Users\hogehoge\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\numpy\lib\npyio.py", line 525, in save

        file_ctx = open(file, "wb")

    FileNotFoundError: [Errno 2] No such file or directory: 'output/fsk_o_t.npy'

     

    って言われる...フォルダ作ればいいのかな?

    コマンド画面でmkdirで作れるらしい。

    PS D:\Documents\python3-vscode> mkdir output

     

     

        ディレクトリ: D:\Documents\python3-vscode

     

     

    Mode                LastWriteTime         Length Name

    ----                -------------         ------ ----

    d-----       2021/04/27     23:16                output

     

    ってフォルダができる。

     

    で、実行するといろんなwarningが出るけど(たぶんSpyderでも出てたんだと思う)、終わったっぽい...グラフが出ない。

     

    Spyderでは要らなかったけど、こちらの環境ではplt.show()が必要らしい。

     

    うまくいきました。

    変更したのは2，3行目で、80行目を追加です。

    1. # 変数クリア
    2. import sys
    3. sys.modules[__name__].__dict__.clear()
    4.  
    5. import matplotlib.pyplot as plt
    6. import numpy as np
    7. import scipy.signal as signal
    8.  
    9. # シンボル列を引き延ばす
    10. def extend_symbol_to_simulation_sampling_rate(symbols,t_symbol,t_simulation_sampling):
    11.     # シンボル数
    12.     size_of_symbols=symbols.size
    13.     # 最後のシンボルが終わる時間
    14.     end_time=size_of_symbols*t_symbol
    15.     # 計算タイミング
    16.     t=np.arange(0,end_time,t_simulation_sampling)
    17.     # 各シンボルの開始（終了）時間
    18.     t_periods_of_symbols=np.append(0,np.arange(1,size_of_symbols,1)*t_symbol)
    19.     # 結果データ領域確保
    20.     ex_symbols=np.zeros_like(t)
    21.     # 先頭データは入れておく
    22.     ex_symbols[0]=symbols[0]
    23.     # 各シンボルについて、その期間のインデックスを取得し、データを入れ込む
    24.     for i in range(0,size_of_symbols):
    25.         ex_symbols[np.where(t>t_periods_of_symbols[i])]=symbols[i]
    26.     return ex_symbols,end_time
    27.          
    28. # よく使う変数
    29. pi=np.pi
    30. deg2rad=pi/180.0
    31. twopi=2*pi
    32.  
    33. # 設定
    34. f_carrier=433.92e6            # 搬送波周波数 in Hz
    35. f_deviation=1e6         # 周波数偏差 in Hz
    36. f_sym=1e6                 # シンボルレート in Hz
    37. f_samp_simulation=f_carrier*10     # 計算のサンプリングレート in Hz
    38. phase_initial=45*deg2rad  # 搬送波の初期位相(何でもいい) in rad
    39. symbols_in_base=signal.max_len_seq(7)[0] # シンボル
    40. #symbols_in_base[0::2]=1
    41. #symbols_in_base[1::2]=0
    42.  
    43. # 後で使う変数
    44. t_sym=1/f_sym
    45. t_samp_simulation=1/f_samp_simulation
    46. omega_carrier=twopi*f_carrier
    47. omega_deviation=twopi*f_deviation
    48.  
    49. # シンボルを計算用に拡張
    50. [symbols,t_end]=extend_symbol_to_simulation_sampling_rate(symbols_in_base,t_sym,t_samp_simulation)
    51. symbols=symbols*2-1 # -1 or +1 にする
    52.  
    53. # 計算するタイミング
    54. t=np.arange(0,t_end,t_samp_simulation)
    55.  
    56. # 全計算タイミングにおける瞬間的な角速度
    57. omega_in_moment=omega_carrier+omega_deviation*symbols
    58.  
    59. # 全計算タイミングにおける瞬間的な位相変化量
    60. phase_change_in_moment=omega_in_moment*t_samp_simulation
    61.  
    62. # 全計算タイミングでの位相
    63. phase_t=np.cumsum(phase_change_in_moment)+phase_initial
    64.  
    65. # 波形生成
    66. amp_t=np.exp(1j*phase_t)
    67.  
    68. np.save('output/fsk_o_t',t)
    69. np.save('output/fsk_o_amp_t',amp_t)
    70.  
    71. # プロット
    72. fig=plt.figure()
    73. fig.add_subplot(3,1,1)
    74. plt.plot(t,symbols)
    75. fig.add_subplot(3,1,2)
    76. plt.plot(t,amp_t)
    77. fig.add_subplot(3,1,3)
    78. plt.plot(t,amp_t)
    79. plt.xlim(t_sym*7.95,t_sym*8.05) # ちょうどシンボルの変わり目を表示
    80. plt.show()

     

    環境を整えるのがAnacondaと違って一発で完了ではないけど(VSCodeの設定が手間)。まぁまぁいい感じ。

    カテゴリ： 非現実