MATLAB

どうもこんにちは。 今回は仕事上でAWSのLinuxサーバー上にMatlabRuntimeをインストールして画像の解析アプリケーションを作動させることになりいろいろ準備していたのですが Gstreamerがインストールされてないよ！ という旨のエラーを吐いたので一時中断。 要するにGstreamerというソフトウェアをインストールすればよいのだがバージョン管理などの問題が発生したため備忘録的に記録する。 まずエラーの根本原因 ・Gstreamer1.0以上をインストールしてない これに尽きる。デフォルトではGstreamer0.1などが入ってたりするらしくそれでエラーが出やすい。（今回ははいっ…

Simulinkを用いたDCモータシミュレーションを、トランジスタ技術2022年2月号のとおりに行いました。 モータモデルやフィードバック制御は把握していても、どのようにシミュレーションを行うか疑問に思っている方は私以外にも多いと思います。「トランジスタ技術2022年2月号 モータ物理シミュレーション」が非常に参考になったため紹介します。 DCモータシミュレーションとは DCモータシミュレーションの方法 シミュレーション結果 シミュレーションの内容 機械のモデル化 モータのモデル化 トルク目標値への制御 DCモータのモデルの概要 個人的に作業に詰まった点について まとめ DCモータシミュレーシ…

CNN用にtrainとtestにデータを分ける CNNでMRIのシーケンスの疾患の分類精度を比較する 各シーケンス毎に同じ患者のデータが入っているので、各シーケンスでtrainとtestが同じ患者データになるように分けるフォルダ構造はこんな感じ 手順 フォルダ構造を得る 保存用のフォルダを作成する 分割用のインデックスを乱数を用いて発生させる movefileでファイルを移動 各シーケンスで繰り返す フォルダ構造を得る main_fd = uigetdir() [a,b,c]= fileparts(main_fd); sub_fd =dir(main_fd); %remove waste fi…

画像のResizeをparforで行う。parforは並列で処理を行なってくれるので、高速で作業が終わります。手順はこんな感じ フォルダの構造を得る 保存用フォルダを作成する for文でコードを作成 最後にresiezeに関わる部分をparforに変更する フォルダ構造はこんな感じ ・まずはメインフォルダを選択して内部構造を得ていく main_fd = uigetdir(); [a,b,c]= fileparts(main_fd); sub_fd =dir(main_fd); %remove waste files sub_fd =sub_fd (~ismember({sub_fd.name}…

MATLABで配列xとyをプロットする時はplot(x,y)ですよね。plotのオプションを指定する時は、 Name,Value の引数ペアをコンマ区切りで加えていきます。例えば、点線で線の太さを2ポイントにしたければ、plot(x, y, ’LineStyle', ':', 'LineWidth', 2)です。 plotの線を透明化するときは、この手法は使えません。plot関数はchart line オブジェクトとよばれるオブジェクトを出力することができるので、そのオブジェクトを使って線を透明化します。例えば、次のようにやると透過度が0.2の線（0で透明、1で非透明）となります。 p=plo…

この記事は、MATLAB/Simulink Advent Calendar 2021の7日目の記事として書かれています。 qiita.com また、8日目も記事を投稿していて、もしよろしければこちらもご覧いただけると幸いです。 kentapt.hatenablog.com また、ここで用いるコードなどはこちらからダウンロード可能です。 もし、本記事に誤りや誤字が含まれていれば、教えていただけますと幸いです。 github.com はじめに Spatial CNNとよばれる手法を用いて、運転動画から車線検出を行います。ここで得られる結果の例を以下に示します。自分のいるレーンの車線とその横の車線の…

年末調整の給与所得を算定するための『訳の分からない表』をプロットしてみた。 『訳の分からない表』はこれ↓ プロットした結果がこれ↓。給与収入が百万円（106）から一億円（108）までをプロットしている。赤線が控除後の所得で、黒点線は控除無しの所得（給与収入の全額が給与所得）である。したがって、給与が低い人ほど多く控除してもらっているのが分かる。あの『訳の分からない表』をプロットすると結構美しいというか、滑らかに接続していると思う。 これをプロットしたmatlabのコードは以下のとおり。 close all clear variables %%%%%%%%%%%%%%%%%%%%%%%%%%%%…

インパルス応答 OATSPの数式 Common Lispで試す 時間領域にもっていく 関数にした Juliaで同じものを書いてみると…… インパルス応答 頭部伝達関数を測定したり、スピーカーの応答を測ったり、室内音響指標の計算に使ったり……インパルス応答が測定できると色々と便利です。すでにMatlab、R、Juliaなどでインパルス応答測定のプログラムは書いてあるのですが、勉強中のCommon Lispで行うのが目下の目標です。 まずはインパルス応答を測定するための信号の生成を行います。今日は生成だけで、ファイルへの保存や実際にオーディオデバイスから出力するところなどは後々やります。 代表的な…

McCabe とは 循環的複雑度とは 試してみる インストール 実行方法 flake8 のプラグインとして まとめ McCabe とは github.com Python コードの複雑さをチェックしてくれるモジュールです。複雑さをチェックすることで、保守性の悪いコードをリリース前に検知することが可能です。もちろん、コードレビューで検知することも出来ますが、複雑と感じるかは個人差がありますので、予めチーム内で指標を決めて機械的にチェックする方が生産性が高そうです。 尚、複雑さ は具体的には 循環的複雑度 と呼ばれる指標を用います。 循環的複雑度とは 循環的複雑度（英: Cyclomatic co…

こんにちは。モフモフスキー(rab (@movmovsky) / Twitter)です。 今日は私（記事執筆時で学部3年）が所属している計数工学科の授業を紹介したいと思います。新しく計数に内定された方や進学を考えている方にも、学科HPには載っていない授業の所感などをお届けできたらなと考えています。 （底点等の進学するまでの情報って結構潤沢ですけど、入った後の情報って比較的少ないですよね..） なお、現在私が学部3年なこともあり、本記事では計数進学後の2Aセメスターと3Sセメスターのみに絞って話していくので、その点ご承知おきください。 3Aセメスター以降の授業に関してはこちらの記事などに詳しく書…

前回は、eeglabを起動して、脳波データを読み込むところまでやりました。 今日は、土曜出勤しているので、もっとどんどん進めますよ～♡ 前回の続きのイベントファイルを読み込むところから書きますね。 １－２ イベントファイルの読み込み file > Import event info > From Matlab array or ASCII fileをクリックします。 Import event infoのウィンドウが出てくるので、右上のBrowseをクリックし、事前に作っておいたイベントファイルを選択する。 ↓ 読み込んだイベントファイルを開いてみると、上にlatencyとtypeというラベルが…

先月末「Common LispでFFTライブラリを使ってみる - 丸井綜研」を書きました。このエントリの最後に「関数化した上でWAVファイルを読み込んで分析してみたい」と書いていたので、今日はそれをやってみました。 WAVファイルを読み込むのには「Common Lispで音ファイルを読み込む - 丸井綜研」でも使ったcl-wavを使います。以下のようにすれば、音を読み込んだあとでWAVに関する情報が得られます。ここではWindowsに最初から入っている44.1 kHz / 16 bit / 2 chの音ファイルを使ってみました。 (ql:quickload :cl-wav) (defun wa…

前置き(いらない) 最近気温がちょっとずつ下がってきて秋の訪れを感じます。 何もしてないのに時間が過ぎていてちょっと怖いです…。 CHMPWMとは 大容量インバータなどでPWM周波数をどうしても下げなければいけない場合、どうしても高調波含有量による電動機損失が問題になります。(テンプレ) この問題を最小限にするため、電流の高調波成分が最小になるようにノッチ角を設定するPWMを高調波電流最小化PWM(CHM: Current Harmonic Minimum)と言います。自己満足です。電流の高調波含有量を最小にすることで、電流リプルによるトルクリプルあるいは電動機の温度上昇を最小限に抑えます。 …

はじめに OpenCVとはインテルが開発・公開したオープンソースのコンピュータビジョン向けライブラリで、画像処理・画像解析および機械学習等の機能を持つC/C++、Java、Python、MATLAB用ライブラリとなっています。 OpenCVには以下のような機能があります。 グレースケール化 2値化 輪郭の検出 顔検出 ノイズ除去 テンプレートマッチング OpenCVは上記などの画像処理だけではなく動画の表示や保存をすることができます。動画を取得して1フレームごとに画像処理を行うこともできます。例えば1フレームごとに人物検出を行い、検出された場合にはその部分に四角形を描画して新たな動画として保存…

各確率分布に従う確率変数の4次までのモーメント 二項分布 B(n,p) 次数 モーメントE[Xn] 中心モーメントE[(X-E[X])n] 1次 2次 3次 4次 ポアソン分布 次数 モーメントE[Xn] 中心モーメントE[(X-E[X])n] 1次 2次 3次 4次 正規分布 次数 モーメントE[Xn] 中心モーメントE[(X-E[X])n] 1次 2次 3次 4次 ガンマ分布 次数 モーメントE[Xn] 中心モーメントE[(X-E[X])n] 1次 2次 3次 4次 n次モーメントの計算方法 次モーメントは，確率変数 の積率母関数 の 階微分に を代入した値である（参考：積率母関数）． ま…

先日Common LispでFFTが使えるようになったので、それを使って畳み込みをします。畳み込みには環状畳み込み（円状畳み込みや循環畳み込みとも呼びます）と線形畳み込みの2種類がありますので、両方とも作ってみようと思います。 ちょっとした道具 実際の畳み込みのプログラムに入る前に、Matlabを参考にしてnextpow2を作ります。これは与えられた値より大きい最小の2のべき乗の指数を計算するものです。たとえば5の次の2のべき乗は8=23なので3、12の次の2のべき乗は16=24なので4が返ります。（ceilingは多値を返すので実際には余りも出てきますが、ここでは1つ目の返り値だけを使います…

お久しぶりです．せかPです．さて，最近は朝晩の風が涼しくなり，理工学基礎実験 IIB が脳内をちらつく季節かと思います．冗談は置いておいて，きそじのレポートを Word でやりたくないな， TeX でやってもいいのかな？と思う人もいるはずです． そんな方に，ちょっとだけアドバイスを！という記事です． それ以外の方は読んでも特にメリットが無いと思います． 結論 先に結論を述べておきます． TeX で書いても，レポートの評価の観点では特に問題はない． 数式+文章 TeX ，グラフはエクセルで作成して貼り付けなら割と楽 TikZ でグラフを指定の書式に従って書くのが面倒 そんなことより考察をしっかり…

2つのブログを初投稿させていただいてから早いもので数週間。 自己紹介 - AIソリューションアーキテクトのおしゃべり データ活用人材とはどんな人々か - AIソリューションアーキテクトのおしゃべり 思っていた以上に反響があり、 「面白くて、あっという間に読めた」 「ブログとリアルのキャラが一致しなくて困惑した」 「怪しい壺が買いたくなった」*1 など、数々のご好評をいただきまして非常に感謝してます。 壺は売らないという決意を揺るがせないでください。 時系列解析って何か色々あってわかりません！(憤怒) さて、本日取り上げたいのは時系列解析。 最近では新型コロナウイルスの感染者数という大々的に報じ…

はじめに OpenCVとはインテルが開発・公開したオープンソースのコンピュータビジョン向けライブラリで、画像処理・画像解析および機械学習等の機能を持つC/C++、Java、Python、MATLAB用ライブラリとなっています。 OpenCVには以下のような機能があります。 グレースケール化 2値化 輪郭の検出 顔検出 ノイズ除去 テンプレートマッチング 画像内での物体の位置を把握する方法の一つに輪郭検出があります。OpenCVを用いて輪郭検出を行うには findContours() およびcv2.polylines()などが使用されます。 本記事では、OpenCVのcv2.threshold(…