Axis2

本日のセミナーの課題は，２次自己回帰過程 ]のサンプル時系列から自己共分散関数とパワースペクトルを推定することです．]は，平均0，分散の白色ノイズです．とりあえず，，でやってみて，他のパラメタでも調べてみてください．これまでの知識を使えば簡単に実行できると思います．追加の課題として，白色ノイズ]の部分を，単位インパルス に置き換えた場合の振舞も数値的に調べてみてください．インパルスをいれるタイミングは0でなくても構いません。例えば，初期値をy[1] = 0，y[2] = 0として，ｗ[n]を{0, 0, 1, 0, 0, 0, 0, 0, 0, 0,... ,0}にしてください． 2次自己回帰…

今日のセミナーの課題は， 問題 数値的に生成した1次自己回帰過程のサンプル時系列の，自己共分散関数とパワースペクトル（ピリオドグラム）を推定し，理論値と比較せよ．推定は，複数のサンプル時系列を使って平均として求めよ．解答としては，以下のような図を想定しています．ただし，この図では，ラグと周波数の領域を欲張りすぎで，実用上は，結果として参照するラグの上限はN/10（最大でもN/4）程度，周波数の上限は10/N（最小でも4/N）程度にしてください． 実行例（時系列の長さN=2^13=8192，サンプル数N.samps=2^6=64）上の図のような数値実験の結果を作るＲスクリプトは，以下です（今回は…

今日のセミナーの課題は，以下です． 問題1 1次自己回帰過程の時系列サンプルを生成せよ． ここで，は，の定数，は，平均0，分散の白色ノイズである． ，とせよ． 問題2 1次自己回帰過程のパワースペクトルを推定せよ．100個の時系列サンプルのパワースペクトルを計算して，周波数毎に平均をとること． パワースペクトルの推定量は，と，教科書の定義（自己共分散関数のフーリエ変換）の2つの方法を使い，結果を比較すること． 問題3 1次自己回帰過程のパワースペクトルの理論曲線（解析的に求めた結果）と，上で計算したパワースペクトル（数値的に求めた結果）のグラフを重ねて描け． 時系列の生成 Rで時系列を再生する…

今回はFreeCADで断面性能を表示させるPythonスクリプトを説明しますAutoCADにはリージョンの断面性能を表示させるマスプロパティ（MASSPROP）というコマンドがあります 複雑な断面形状の図心位置や断面2次モーメント等を知りたいとき大変重宝する機能ですね これと同じ様なことをFreeCADでもやりたいというお話です(*´ω｀*) 前提条件として FreeCAD上で定義された単一の面（Face）を対象とします XY平面，XZ平面, YZ平面のいずれかに平行に配置されていることを想定しています FreeCADのv0.20でしか動作検証してませんPythonスクリプトを以下に示します【…

現在の状況 右スティック側のAXIS2, 3が0.01超過である。 実際はもっとスティックの反応がビクビクしていて痙攣しているような感じ。 絶対値で0.01以下だとよいらしい。 gamepad-tester.com dualsense 右スティック不調

Raspberry Pi PicoでDVI出力（HDMI）を試してみました。 前回の「VGA出力」はアナログRGBなのに対して、今回はデジタルRGBです。上の写真はAdafruitが公開しているサンプル「hello_dvi」を実行中の様子です。 使用した部品、「Adafruit DVIブレークアウトボード」です。スイッチサイエンスで330円（+送料）でした。 この他、「Raspberry Pi Pico用DVIアダプタ基板（990円）」という部品も使えます。商品説明が「HDMIっぽい」となってますが、HDMIです。 説明文を元に書いた回路図。 VGA出力にくらべてシンプルです。公式もこっちを採…

桜が見頃になってきましたね！ ついこの間咲き始めたと思ったら、いつの間にやらもう満開。季節の巡りは早いものです。 今回は天体観測をします。 春の夜には早々と消えていきますが... 私の大好きなオリオン座の星雲を観てみます。 赤外線天文衛星「あかり」などの望遠鏡によって撮影された写真をいじってみたり - Qiita こちらの記事を参考にさせていただきました。ありがとうございました。 Astroqueryというモジュールを使うことで、SkyViewの天文写真をPythonで取得することができます。望遠鏡を覗き込まずに、プログラムから天体観測ができるなんてとても素敵です。 では、いってみましょう。 …

みなさまこんにちは。プロダクト開発本部の岸田です。 以前に「ハードウェア側に機械学習環境を立てて推論を行い、クラウドに結果を収集して分析状況を確認する」ユースケースをこちらの記事で考察しました。手軽にローカルデバイスとクラウドを連携するサービスとして有名なものが「AWS IoT Greengrass」ですよね。公式のドキュメントでもエッジ推論システムをAWS IoT Greengrassで構築するユースケースが紹介されています。 docs.aws.amazon.com そこで上記を参考に、弊社の製品であるJetson TX2搭載のEDGEPLANT T1を利用して、AWSサービスを活用したエッ…

今回の電子工作は「PC接続したゲームパッドのスティック入力をBluetoothでラズパイに飛ばしてDCモーターおよびサーボモーターを駆動させる」というものです。

torch.maxとtorch.argmaxがゴタゴタしていたのでメモ torch.maxは1つ目に最大値、2つ目に最大値の入っているindexを返す。 torch.argmaxはtorch.maxの2つ目のみを返す。 axisを指定しないと1行にreshapeされた後にargmaxが行われる。 import torch a = torch.rand((2,3,4,5)) print(a.shape) # torch.Size([2, 3, 4, 5]) print(a.argmax(axis=0).shape) # torch.Size([3, 4, 5]) print(a.argmax(a…

薬剤の投与量を考えるときに、クレアチニンクリアランスというものを考慮しないといけないのだが、クレアチニンクリアランスにはクレアチニンと男女と体重が必要である。 そんなものは電子カルテで自動計算してくれればいいのだが、デフォルトではしてくれない仕様である。 クレアチニンクリアランスで、だいたい50以下、30以下、10以下、くらいがわかればいいのでとりあえずやってみる。 クレアチニンクリアランスは、クレアチニン値、年齢、体重 として で、女性は0.85倍される。体重60kgとすると、1.5 であるので、そんなに高齢でなければCre 2 くらいまではなんにも腎調節のことは考えずに初回投与してもたぶん…