last modified: 2018/11/03 18:19:40 UTC+0900
2018.10.31(水) 前日 当日 翌日 コメントする 雑記帳
- 晴れの朝。朝のホームでの電車待ちが肌寒く感じる。
-
[授業(秋1, 水曜, 第8週)]
- 定例卒研報告会
いろいろと指示を出した。
- 定例卒研報告会
- ついに8週目
- 論文P投稿
- 昨日の計算の最初の立式にミスを見つける......orz...... HMHDの\(\langle L_\xi V|L^\dagger_\xi V\rangle \)書けた。 MHD極限も一致した。 \(C\int\Psi\cdot(\nabla\times\Psi)dV-\int\Phi\cdot(\nabla\times\Phi)dV\) の形なので正定値性は望むべくもないような気がする。
- とある書類．
- 自分用の追尾ハードウェアを作成した。 かなり適当なコードで動いた。 OpenCVの顔検出の誤検出をどうやり過ごすか。
- ORCIDに自分の過去の論文のリストを作成した。 クリックだけなので簡単だった。
- 快晴の朝
- [授業(秋1, 火曜, 第7週)]
- 長崎大学に電話
- 第2変分の計算。MHD極限と合わない。
- pigpio.spi_xfer, pigpio::spiXfer の動作の再確認。 5ビットのMOSIのコマンドを1バイトの中に詰め込むことができると分かる。 それに合わせてはてなブログのコードも書き換えた。
- 勉強会。 追尾できるようになった。
- 晴れの朝。島々の姿は霞んでいる。 駅のホームの向かいの土手の桜の葉がほとんど散ってしまっていることに気付く。
-
[授業(秋1, 月曜, 第7週)]
- 知能機械工学実験２(第7回)引張試験
- 組合のメールボックスが一杯だと連絡があった
- 対応1件．
- \(\langle \boldsymbol{\Psi} | \nabla\times\boldsymbol{\Psi} \rangle \)の符号次第？
- 晴れの一日。2次元非粘性流の復習。 Laplacian の固有値問題の境界条件を設定し忘れていた。 2次元流なので流れ関数の等高線が流線であり、 関数値の定数分のずれは物理的には意味が無いことを思い起こすのに時間を食った。 ITC27/APPTC2018の原稿の手入れ…というか、 夏休みの暮れに書き殴ったアイディアを整理したりとか、 文献参照を整理したりとか。 Noether の文献は Nachrichten von der Gesellschaft der Wissenschaften zu Goettingen, Mathematisch-Physikalische Klasse (1918), Volume: 1918, page 235-257 である。どう省略するんだ？ 2次の摂動は Lagrangian だと \( \langle L_\xi V|L_\xi V\rangle + \langle L_\xi V|L^\dagger_\xi V\rangle \) の形で、 Hamiltonian だと \( \langle L^\dagger_\xi V|L^\dagger_\xi V\rangle + \langle L_\xi V|L^\dagger_\xi V\rangle \) の形。 いじっている空間の微分がもろに出る。 \(L^\dagger_\xi V\) の方が固有値問題の組み込みが簡単。
- 快晴の朝。雨上がりで地面が濡れている。大学に上がる■ 昨日、動作を確認した pigpio/Python ライブラリの spi_open(), spi_xfer() のパラメータの解読を進めた。pigpio/C の spiXfer() のコードも書いた。 いずれも pigpio のドキュメントにある auxiliary SPI の CH2 を動作させて確認した。 これはクックブックが無いとパラメータの渡し方と、 戻ってきたデータの組み替え方が分からない！■ ミニマルコードはブログにまとめた： pigpio ライブラリの spi_xfer() 関数を用いて Raspberry Pi 3B から SPI 接続のA/Dコンバータ MCP3208 を利用する■ 以前に Python/pigpio::write(), read() で作った readadc() が任意のポートで動作することを今日になってようやく確認した。 このコードを C に移植しようとしたが、できなかった■ 勉強会。 サーボモータのブラケットをブックエンドに仮止めして動作をやらせる。 向きを正しく取り付けていたので、コードの改良が捗った。
- 快晴の朝。天気予報は今夜、天気が崩れると告げる。 ようやく議事録のリリース
- 書類．
- 柏木君にコードの ver.2 を渡して、動作を見る。 予定通りの動作をしたので嬉しい．
- pigpioによるI2CとSPIインタフェースの実装に掲載されているコードを用いて、 pigpio の spi_open(), spi_xfer() を用いた MCP3208 からのデータ取得がようやくできた。
- make コマンドの存在を教える
- ちょっと． 広汎性発達障害？
- 快晴の朝。昨夜は道端にタヌキの姿を見かける。 今朝、「今月今夜」と書かれた地区の神社の御神燈が立っていることに気付く。
- [授業(県立大3Q,第5週)]
複素関数論(第5回)
第4回
第6回
5.複素関数の微分(教科書 第3章 p.42-57)
微分の形式的な定義と\(\Delta z\to0\)の意味。実1変数と実2変数のTaylor展開の図形的意味を Microsoft Math を用いて説明。全微分可能の定義。原点で偏導関数の値が方向によって異なる関数の例; \(\frac{x}{\sqrt{x^2+y^2}}\)。導関数の定義からの計算とCauchy-Riemannの関係式の導出。Cauchy-Riemannの関係式から調和関数の導出。調和関数の例：\(x^2−y^2\), \(x^3−3xy^2\). 正則の定義; 目的:点の情報をそれを囲むループの情報から計算するので, Greenの定理が領域全体で使えないといけない; 内点、外点、境界点とは。\(\mathbb{Q}\)は境界点しかない集合
来週は定義から導出できる公式: 線形性, Leibniz rule, chain rule - [授業(秋1,木曜,第7週)]
- アドバンストセミナー(第4回)
先週休んでいた[ぴー]君に「ラズパイでLチカ」ができるように、他のゼミ生に指導をさせる。 PRi3B+/nano で Python コードを書く練習。
- アドバンストセミナー(第4回)
- 先日の案件の始末書．
- 急激な腹下しで4限目のFD講演会をパスする
- 線形代数の質問(3×3行列の対角化)に答える。 答えの検算は RPi3B+/Mathematica を用いた。 ミスの頻度が小さいので、計算のうまい学生だと思った。
- Linuxでkbhit関数作成 (C言語) を参考にして、 RPi3B+/C言語 でキー入力によるイベント駆動型プログラムのひな型を作成した。
-
RPi3B+/pigpio/C言語 のコードのコンパイル(cf. pigpio C Interface)
gcc test.c -pthread -lpigpio -lrt
実行は(cf. pigpio Can't lock /var/run/pigpio.pid)sudo killall pigpiod sudo ./a.out
- 晴れの一日。ここしばらく出勤時にいつもの東向き区間で、 正面から日が差して前が見づらい。
-
[授業(秋1, 水曜, 第7週)]
- 定例卒研報告会
- 午前は卒研指導のためのコード書きと 「ロボット知能化演習」で使っているコードを書き換えたものの動作確認。 エラーを吐いたところが、去年もひっかかっていたところと後で気付いてメゲる。
- 勉強会。 サーボモータのブラケットの組み上げ。
- 昼過ぎより雨
- [授業(秋1, 火曜, 第6週)]
- 10時よりヒアリング1件。14時よりヒアリング1件。
- 体調不良との連絡あり。
- Tic-tac-toe にQ学習を組み込んでみた。 パラメータが悪いのか、評価が悪いのか、全然、強くならない。
- 快晴の朝
-
[授業(秋1, 月曜, 第6週)]
- 知能機械工学実験２(第6回)荒木担当分なし
- ほぼ一日、卒研関連のRPi3B+環境整備とコード書きなど。 cv2.VideoCapture(0)がUSBカメラの認識にコケる場合がある。
- 勉強会。RPiに新しいOSのmicroSDカードを入れて復活
- 晴れの一日。惰眠。RPi3B+の育成。新品の NOOBS microSD からの起動。 Qiita のページ『ラズパイ3にOpenCV3を簡単に導入』に従って実行したら使えるようになった。 これで RPi3B+/OpenCV/pigpio の環境を一つ増やした。 卒研用のバックアップが一つ増えてほっとする。 日本語フォルダ名を英語にする方法。 乱数で手を打つ Tic-Tac-Toe の Python コードを書く。 夕方より舜と二人でテラスの下で焼肉をする。 炭火を起こすのに時間を取られる。
- 快晴の朝。舜、参観日（「おまつり」らしい）。 卒論のコード書きで大学に上がる。 RPi3B を燃やしたかもしれないとの報告が。 tkinter で event_generate("<Button-1>") が期待と異なるウィジェットへのイベントを投げるバグで悩む (帰りのバスでバグの原因に気付いて、解決した)。 RPi3B+ の起動。 sudo apt-get install xrdp しただけで windows からのリモートデスクトップ接続ができるようになった。 emacs と samba も入れて Windows からの SSH 接続経由のファイル編集と共有が自在に。
- 快晴の朝。
- [授業(県立大3Q,第4週)]
複素関数論(第4回)
第3回
第5回
逆三角関数の定義と解。
4. 1, 2, 3のまとめの演習
複素関数の計算: \(\{0,\pi,\pi+i\pi,i\pi\}\)を頂点とする矩形を\(z^2\), \(e^z\)で写像。 円\(|z|=e\)を位相を\(0\leq\theta<2\pi\)に制限して\(\ln z\)で写像。 \(\sin z\) を加法定理を用いて実関数で書く。 \(\sin z=2\) の解を求める。 - [授業(秋1,木曜,第6週)]
- 県立大の非常勤講師控室のプリンタを利用してみた
- 県立大からの帰り道、 備中国分寺と作山古墳の間に広がる田んぼの稲刈りが意外に進んでいなかったので、 少し驚く。
- 学生より気になる連絡．
- 面談1件．
- 会議があった
- 会議があった
- 気がかりな連絡．
- 昨夜の雨が嘘のような快晴の朝。
-
[授業(秋1, 水曜, 第6週)]
- 定例卒研報告会
- アドバンストセミナー(第3回)
ラズパイの再起動。 LANのIPアドレススキャンツールの導入。 Linuxのコマンド ls cd cat .. [TAB]補完 等の指導。
- とある案件をお断りする
- 勉強会。RPiが起動しなくなる。
- なんか立て続けにRPiの不具合が続く
- 曇りがちの朝
-
[授業(秋1, 月曜, 第5週)]
- 知能機械工学実験２(第5回)荒木担当分なし
-
Raspberry Pi 3B + OpenCV 機に
KP-CAM32 を2台接続して動作した。ミニマルPythonコード：
# -*- coding: utf-8 -*- import cv2 cap0 = cv2.VideoCapture(0) cap1 = cv2.VideoCapture(1) while(True): ret0, frame0 = cap0.read() ret1, frame1 = cap1.read() cv2.imshow('frame0', frame0a) cv2.imshow('frame1', frame1a) if cv2.waitKey(1) == 27: break cap0.release() cap1.release() cv2.destroyAllWindows()VGA(640x480), 30fps らしい。 - Raspberry Pi 3B+ が届いていた
- 晴れの一日。午前、惰眠。 午後、だらだらと「知能化演習」用のコードを自分流に書き換えて過ごす。
- 快晴の朝。大学に上がる。 2年生に「数値計算法」関連の Gnuplot インストールの面倒を見てあげたりとか。 4年生の Django の書籍を1章分丁寧に読む。 勉強会。
- 雨の朝。駅のいつものお祭りの旗が片付けられていることに気付く。
- [授業(県立大3Q,第3週)]
複素関数論(第3回)
第2回
第4回
3.複素関数、写像
[教科書,第2章,p.22-] 複素関数の作り方：整数べきの級数(収束性)；実2変数関数2個(可微分性)。 多項式、指数関数、三角関数(指数関数で定義)と実関数での表現 \(\cos(z)=\cos(x)\cosh(y)-i \sin(x)\sinh(y)\)、 双曲線関数(指数関数による定義、双曲線関数の読み)、 n乗根、対数関数と位相に由来する多価性、 複素数の一般べき \(a^z=e^{z\ln a}\), \(2^i=(\exp(\ln(2)+2\pi i k))^i=\exp(-2\pi k+i\ln(2)) ={\underbrace{(e^{2\pi})}_{535.4...}}^{-k} [\underbrace{\cos(\overbrace{\ln(2)}^{0.6931...})}_{0.7692...} +i\underbrace{\sin(\ln(2))}_{0.6389...}] \) \(i^i=(\exp(i\frac\pi2+2k\pi i ))^i=\exp(-\frac\pi2-2k\pi) \)
逆三角関数の定義と解は次回へ。 - [授業(秋1,木曜,第5週)]
- アドバンストセミナー(第2回)
ラズパイの起動と終了。いわゆる「Lチカ」のデモ。
- アドバンストセミナー(第2回)
- 県立大から清音を抜け、高梁川を渡ると真備である。 二万口の方面へと向かう道すがら、道端の家々を見ると、 1階と言わず2階と言わずサッシの窓がみな外されて、 そこから部屋の中の板を外された天井が垣間見えた。
- Threadripperマシンの調子がおかしいので修理に出した
- 会議
- 曇りがちの朝。昨日の事、蛇の脱け殻が落ちてたと言ってきたので、 「金運のおまじないになるね」と教えたら、長さ1.5m以上あろうかという アオダイショウか何かの脱け殻まるまる一匹分を 拾ってきていた。これは立派だ。ここまで成長したということは、 生命力にあやかれるかもしれないね。
-
[授業(秋1, 水曜, 第5週)]
- 定例卒研報告会
- 晴れの朝。 いつもの峠道を下りて盆地側に入ると稲刈りが済んで稲木を作ってある田が何枚か。 ホームの電車待ちで向かいの土手を見ると、 桜の葉が色づき始めて、時折はらはらと落ちはじめていることに気付く。
- [授業(秋1, 火曜, 第4週)]
- pip install django でインストールできたのかな？
- 流力中四九＠長崎のアナウンスを流す
- 快晴の朝
-
[授業(秋1, 月曜, 第4週)]
- 知能機械工学実験２(第4回)データ整理
- 知能機械工学実験２(第4回)データ整理
- 2019年度学科時間割 1st draft を投げる
- 留学生の受講指導1件
-
\(y=\sin x\)の\((x,y)\)座標データを作成して
テキストファイルに出力するコードを書きだして、
プログラミングの授業で「Cの文法」は教えても
「Cのコーディングの初歩的なノウハウ」はロクに教えていないと気づく。
例えば：
- ソースコードを自動整形してくれるソフトウェアを使う
-
数値には意味に合わせた「名前」をつける。すなわち変数を作る。
[目的1] 2か所以上で同じデータ使うとき、数値の書き換えが1か所で済む
[目的2] マジックナンバー(意味が即座に理解できない数値)を書かない - 定数(プログラム中で意味と値が固定された数)は宣言文で「名前」と「値」を定義する。
#include <math.h> #include <stdio.h> int main (void) { FILE *file= fopen("test.txt","w"); if ( file ) { const double xMin= -3.0 * M_PI, xMax= 3.0 * M_PI, nDiv= 200, dx= ( xMax - xMin )/nDiv; double x; for ( x = xMin; x < xMax + dx; x += dx ) { double y= sin(x); fprintf(file,"%f %f\n",x,y); } fclose(file); } }
- 曇り空の朝。業務のため出勤。 毎日、書き溜めていたこの日記をほぼ1か月ぶりにアップする。 駅に帰り着くとホームにお祭りの旗が例年のように飾られていることに気付いた。
- 晴れの一日。台風25号の強風圏のせいで、風が強い。 岡山のAMeDASのデータは10時から17時までの気温が30度C以上であったことを示す。
- 英文校正の書き写し完了
- 勉強会で大学に出る
- 対数法則 \(f(xy)=f(x)+f(y)\)から出発して、 eの極限の式や論点先取な指数関数\((1+1/t)^t\)の導入を避けて、 対数関数の導関数を求める穴埋め問題を作成するのが、意外に難しい。
- 曇り空の朝。
- [授業(秋1, 金曜, 第4週)]
- 2017年に卒業したKさんが学校の「休み」を利用して訪ねてきてくれた。 ご褒美の腕時計をいつも使っていると見せてくれた。 学校での苦労話を聞かせてもらえた。 「いつものK」は健在であったｗｗｗ
- 実数\(a\)(ただし\(a>1\)とする)と十分に大きな自然数\(N\)に対し \(a\approx\left(1+\frac1N\right)^p\) となる自然数\(p\)を探すと、\(p/N\)が\(a\)の対数の近似値になる： \(a\approx\left(1+\frac1N\right)^p = \underbrace{ \left(1+\frac{p/N}{p}\right)^p }_{整数回の四則演算} = \underbrace{ \left[\left(1+\frac1N\right)^N\right]^{p/N} }_{こっちは指数が有理数} \approx e^{p/N}.\) \(p/N\)をいきなり実数\(x\)に置き換え、指数関数 \(\displaystyle \left(1+\frac{x}{p}\right)^p \approx e^x:= \underbrace{ \lim_{p\to\infty}\left(1+\frac{x}{p}\right)^p }_{整数回の四則演算の極限} =\sum_{k=0}^{\infty}\frac{x^k}{k!} \) が芽生えてくる。
- 曇り空の朝。高梁川はいまだに泥水が流れている。
- [授業(県立大3Q,第2週)]
複素関数論(第2回)
第1回
第3回
2.Taylor級数と微分方程式、複素平面上の集合
2階同次線形微分方程式\(f''-(\alpha+\beta)f'+\alpha\beta f=0\)を、3項間漸化式として解いて、結果の\(f^{(n)}\)をTaylor級数の式に代入して、微分方程式の解を求める。漸化式の解法が左固有ベクトルの作用になっていることを示す■ 線分、円のパラメーター表示；PCでCコードを使って表示する。 - [授業(秋1,木曜,第4週)]
- AO入試実施方法説明会
- 会議
- 晴れの朝。
-
[授業(秋1, 水曜, 第4週)]
- 定例卒研報告会
- 「専門英語」のアンケートに回答した
-
(800,600)ピクセルのウィンドウに -4.0≦x≦4.0, -3.0≦y≦3.0
の範囲の座標をマップする C コードのサンプル。
gint(), pset(), gend() は『Windowsプログラムを10倍**に作る』のAPI.
#include <math.h> #include "../xlib/gint.f" #include "../xlib/pset.f" const int WHITE= 7; const int Width= 800, Height= 600; const double xMin= -4.0, xMax= 4.0; const double yMin= -3.0, yMax= 3.0; int i (double x) { return ( x - xMin )/( xMax - xMin ) * Width; } int j (double y) { return ( yMax - y )/( yMax - yMin ) * Height; } int main ( void ) { if ( ! gint(Width,Height) ) { double x, dx= ( xMax - xMin )/(double)Width; for ( x = xMin; x < xMax + dx; x+= dx ) { double y= cos(x); pset(i(x), j(y), WHITE ); } gend(); } } - 勉強会
- 晴れの朝。いつもの峠道では台風の雨水が今朝もまだ側溝や法面を流れていた。
- [授業(秋1, 火曜, 第3週)]
- 英文校正が戻ってきた。真っ赤。
- 台風一過の晴れの朝。体がだるい。昨夜は夜中に目が覚める。 ここしばらく左膝が痛い。 いつもの峠道には木の葉があちこちに散っていて、倒木、倒竹もちらほら。 JRの運行情報では、 今朝は里庄から岡山側で始発から運行していたが(里庄から福山は運転見合わせ)、 里庄からでは車両のやりくりをどうするんだろうと思っていたら、 案の定、ダイヤが乱れていて岡山方面のホームは人が溢れかえっていた。 午前に授業がないので、某所で仕事を進める。
-
[授業(秋1, 月曜, 第3週)]
- 知能機械工学実験２(第3回)引張試験
山田担当分(第2回)を代理で実施(あらき担当分第1回は松浦先生に)。 壁回避のQ-学習を実装
- 知能機械工学実験２(第3回)引張試験
- 進め方についてご指導いただく
- 書類を作成
- 次年度のご担当伺い
- 推薦書を1通作成中
- スマートホンの動作が安定しない。 突然、機内モードになったり(通話中にタッチパネルのどこかを触ったのか？)、 相手側に音声が届かなくなっ(ているらしい反応だっ)たり。
- 64bit Windows 10 + gcc -mwindows なのだが、『Windowsプログラムを10倍**に作る』の掲載プログラムが意外に使える。正方形のフィールド内の十字形の壁避けの Q-learning のプロセスを図にして眺めてみる。 何か蟻の動きを見ているような気分にさせられる。
9月 11月 最新月 あらきトップページ 知能機械工学科 岡山理科大学
(C) あらきけいすけ 2018：
