ブログ

Python embeddableにはpython**.zipなるファイルが付いている．
なんだこれと思っていたらこれにPythonのモジュールなどのフォルダをそのまま突っ込むと使用できるようになる．
もちろん展開に時間がかかるので多少レスポンスがおちる．

小さなモジュール，Pure Pythonモジュール類は問題なく追加できるが，NumpyやScipyはそれができない．
実はpydの実態がWindows dllであるということ．
Windowsではzipファイルからdllをロードできなくなったのでpydをロードできない．したがって，NumpyやScipyはzipファイルに入れては使用できない（厳密にはpydを使っていない部分だけなら使用できるが，内部でどのように依存しあっているかは外からはわからない）．

まず周波数解析の例を示す．周波数軸を作る部分は
np.linspace(0, サンプリング周波数, データ長さ)
で大体あっているのだが，微妙にずれる．
このずれの分が気になったのでよーく考えると以下の例の場合には
freq = np.linspace(0, fsmp - 1/(t_len + t[1] - t[0]), num=len(t));
ということになる．

ピったしの周波数軸が作れる．

import numpy as np
import matplotlib.pyplot as plt

# 例の信号を作成
t = np.linspace(0.001, 4, 4000);
z = 0.1 + 0.2 * np.sin(t * 10 * 2 * np.pi) + 0.2 * np.sin(t * 33 * 2 * np.pi);

# サンプリング周波数
fsmp = 1 / (t[1] - t[0]);
# 解析時間
t_len = t.max() - t.min();

Fz = np.fft.fft(z) / z.shape[0] * 2; # 折り返すのでパワーが２分の１になっている．
Fz[0] = Fz[0] / 2; # 平均成分は折り返さない．
Fz_abs = np.abs(Fz);

freq = np.linspace(0, fsmp, num=len(t));

fname = 'text.csv';
np.savetxt(fname, np.array([freq, np.real(Fz), np.imag(Fz)]).transpose(), fmt='%.18g', delimiter=',');

plt.figure(1)
plt.subplot(311)
plt.plot(t, z)

plt.subplot(312)
plt.plot(freq, Fz_abs)
plt.xlim([0, fsmp/2]);

#136/4
plt.subplot(313)
plt.plot(range(len(t)), Fz_abs)
#plt.xlim([0, 500]); # 後半(高周波側)は前半の折り返し

plt.show();

例えば以下のような演算ができる．
以下の例では虚部と実部を分けて指定する方法と極座標表示から得る方法の例．

import cmath

# 虚部と実部を分けて指定
a1 = complex(1, 0);
a2 = complex(2, 0);

# 偏角を使った指定
mag1 = 1.0;
arg1 = 0.0;
a1 = complex(mag1 * cmath.cos(arg1), mag1 * cmath.sin(arg1))

mag2 = 2.0;
arg2 = 0.0;
a2 = complex(mag2 * cmath.cos(arg2), mag2 * cmath.sin(arg2))


trial = complex(100, 0);

a21 = a2 - a1;
b = a21 / trial

print(b)

いったん複素数になってしまえば掛け算，割り算は可能．
（例の偏角はラジアンなので注意）

IDLEを起動するとPython ShellというWindowが現れる．このとき，ファイルを開く方法は何ら難しいことはないと思う．
ファイルを開くとエディタの中にファイルが開く．
よくIDEを使った開発になれている方であればF5キーでプログラムが実行されるもの推測しやすいと思う．

ただ，ブレイクポイントはどのように設定するのか？

じつは，エディタの側ではなくPython Shellの側のメニューバーのDebug->Debuggerを選択してデバッグをオンにしておかなければならない．

これをオンにしておくとエディタで行を右クリックした時に，Set Breakpoint / Clear Breakpointという選択ができるようになっている．
毎回Visual Studio Codeをインストールしていた．もっと早く教えてよ．

参考
・PythonのIDLEでデバッグを行う方法（GAMMSOFT: https://gammasoft.jp/python/debugging-python-idle/）

調べてみるとQ-360はHIPMには対応していないようだ．

そこで，「PCの動作が突然止まってしまう！SSDの謎の不具合「プチフリ」の原因と対策は？」「Windows10のスリープ不安定問題 その1」を参考に設定をいじってみた．
これで設定がうまくいっているのかわからない．．．．

レジストリエディタを使わないと設定できない項目だったのか

20190118 追記：　ちなみにどうも私の場合にはプチフリーズの原因は違ったらしい．（遅いネットワーク経由でWindows Updateの準備が始まっていたらしい）

1. WinPythonをダウンロードする．

2. ダウンロードしたファイルを実行して適当なフォルダにインストールする．

3. WinPython Command Prompt.exeを実行して，
pip install scipy
を実行する．すると，numpyもインストールされる．
pip install matplotlib
も実行しておく．
4. IDLE (Python GUI).exeを実行して，参考ページからコードをコピーしてくる．
そしてF5キーでこれを実行してみる．
おそらく，
「TypeError: slice indices must be integers or None or have an __index__ method」とエラーになるはず．

5. y1 = signal.savgol_filter(y, n/4+1, 5)などとなっている部分の「/」が問題で「//」に置き換える必要がある．
（y1 = signal.savgol_filter(y, n//4+1, 5)とするということ）

最近のPythonでは「/」が実数を返してしまうということであった．

レイトバインディングならdllへの依存性が減らせるが，プログラムの作成中に面倒なことにオートコンプリートが効かない．Microsoft.Office.Interop.Excelを参照に追加するアーリーバインディングが普通です．

using Excel = Microsoft.Office.Interop.Excel;

これが参考に載っていたやり方．ただし，using Microsoft.Office.Interopでもよいかも．

namespace app
{
    public partial class Form1 : Form
    {
        private Excel.Application objApp;
        private Excel.Workbook objWorkbook;
        private Excel.Worksheet objWorksheet;


        public Form1()
        {
            InitializeComponent();
        }

        private void button1_Click(object sender, EventArgs e)
        {
            objApp = new Excel.Application();
            objWorkbook = objApp.Workbooks.Add();
            objWorksheet = (Excel.Worksheet)objWorkbook.Sheets[1];
            Excel.Range objRange = objWorksheet.Range["A1..B2"];//objApp.get_Range("A1");
            object[,] tbl = new object[2, 2];
            tbl[0, 0] = "TEST";
            tbl[0, 1] = "TEST2";
            tbl[1, 0] = 0;
            tbl[1, 1] = 0.1;
            objRange.Value = tbl;
            objApp.Visible = true;
        }
    }
}

いったんここまで出来上がってからレイトバインディングに直すのが宜しい．
昔は，VBとC#では大きな違いがあったが，最近はC#でもdynamicなる属性が追加されており，InvokeMember関数を最小限にして置き換えができるようになった．

dynamic属性に付け替えExcel.Application()の部分を汎用のType.GetTypeFromProgID()に変更するだけでよい

    public partial class Form1 : Form
    {
        private dynamic objApp;
        private dynamic objWorkbook;
        private dynamic objWorksheet;

        public Form1()
        {
            InitializeComponent();
        }

        private void button1_Click(object sender, EventArgs e)
        {
            Type t = Type.GetTypeFromProgID("Excel.Application");
            objApp = Activator.CreateInstance(t);

            objWorkbook = objApp.Workbooks.Add();
            objWorksheet = objWorkbook.Sheets[1];
            dynamic objRange = objWorksheet.Range["A1..B2"];
            object[,] tbl = new object[2, 2];
            tbl[0, 0] = "TEST";
            tbl[0, 1] = "TEST2";
            tbl[1, 0] = 0;
            tbl[1, 1] = 0.1;
            objRange.Value = tbl;
            objApp.Visible = true;
        }
    }

参考

【Excelの起動から】 C#でExcelを操作する【保存まで】　（ Milk's Memo Note: https://www.milkmemo.com/entry/csharp_excel）