MATLAB

この記事ではモデル誤差抑制補償器（Model Error Compensator, MEC）についてまとめます。モデル誤差抑制補償器に関する動画・関連記事・関連論文・MATLABリンクは最下部に置いています。モデル誤差抑制補償器は，制御系のロバスト性を高める方法です。既存の制御系のロバスト性を高めたり，ノミナルの制御手法と併用してロバスト性を向上したりする目的で使用できます。本記事は以下の総合論文をベースにしてモデル誤差抑制補償器を紹介したものです。式展開や理論の詳細は原稿や最下部の関連論文をご覧ください。 執筆者：制御工学の研究を20年行っている国立大学教員 >>（総合論文，計測と制御PDF…

MRI画像をRescale slope, Scale factorで補正する。 ↑はPHILIPSの場合ちなみにこれはSignal intensityをそのまま使うときようのものであって、SNR、CNRなど”比”をとるときは必要ない。 〜 dcm_path = "DICOMデータのpathを書く" dcm_data = dicomread(dcm_path); dcm_data = double(dcm_data); info = dicominfo(dcm_path); RescaleSlope = info.RescaleSlope; Scale_factor = info.Private…

DICOM のタグを任意に編集して保存する。今回は、すべてのタグを一度からにして、任意のタグの値を編集する。 % DICOMファイルからinfoを読み込む dcm_path = "パスを書く"; info = dicominfo(dcm_path) % タグの数を得る names=fieldnames(info) length(names) ％タグを一旦空にする emp_info = info for x= 1:length(names) emp_info = setfield(emp_info,names{x},""); end % タグを編集する ed_info = emp_info ed…

DICOMデータをCSVへ一括変換する。 MRIのT1 mapなどの処理に有用。 % dcm 2 csv main_fd = "DICOMデータが入っているフォルダを指定"; imgs = dir(main_fd); imgs = imgs(3:end); csv_fd = "CSVファイルを保存するフォルダを指定"; for x = 1:length(imgs) img_path = fullfile(main_fd, imgs(x).name); img = dicomread(img_path); nm = split(imgs(x).name,"."); csv_path = full…

はじめに MATLABとは MATLAB EXPO 2023 JAPANの概要 MATLAB EXPOの雰囲気 講演会の内容 まとめ はじめに MathWorksが主催するMATLAB EXPO 2023 JAPANに参加してきました。このイベントでは、MATLABの最新の技術や活用事例が紹介され、MATLABの可能性を感じることができました。今回は、会場の雰囲気や聴講した基調講演の内容についてお伝えします。 www.matlabexpo.com MATLAB EXPO MATLABとは MATLABは、MathWorksというアメリカの企業が提供している技術計算言語環境です。MATLABを利…

GANや超解像で作成した画像データ（JPEG, PNG）を元のDICOMデータの情報を付加してDICOM変換する。 % 元のDICOMデータのlist作成 dcm_fd = フォルダ名; dcm_fls = dir(dcm_fd); dcm_fls = dcm_fls(3:end); dcm_fls = {dcm_fls.name}; % JPEGデータのlist作成 （今回は超解像をやったことにてHR (high reso) としている） HR_fd = フォルダ名; HR_fls = dir(HR_fd); HR_fls = HR_fls(3:end); HR_fls = {HR_fls.…

ファルダ内に隠しファイルがあるので、３つ目以降を得ることで隠しファイルを削除するのがコツ 〜 fd_path = "OOOOOOO"; fls = dir(fd_path); fls = fls(3:end); for x = 1:length(fls) fls_nm{x} = {fls(x).name} ; end 〜

昨日は飲み会行ってたのでポピーは枯れました。 ゴールデンウィーク実質2日目。昨日教えてもらったアイドルユニットの曲聞いてたら1日終わった。 みんな会社の人に向かってチクチクしてる俺もするか。まじで開発職なのにドット（微分）の意味分からんのおかしいだろ。高校レベルの物理も分からんやつ結構おるしどうやって生きてたんだよマジで。こんなのと同じ扱いされるのはさすがに恥ずかしいわ。 ただ、俺より仕事できるんだよね。ホウレンソウを始めとするコミュ力はあるしやる気もあるし。結局、評価されるのは俺じゃないんだろうね。俺はホウレンソウもせず、勝手に間違って方向に突き進むだけだし。初心者が使うランス？ こんなしょ…

熊本大学情報電気工学科は３つの教育プログラムに分かれています。この３つについて説明されているページが無かったことから新設したのが以下のサイトになります。 熊本大学情報電気工学科・情報電気工学専攻 熊本大学工学部情報電気工学科info infoページでは、概要や学生向けページ、企業向けページなどをまとめています。 電気工学教育プログラム 研究や教育についてまとめています。 情報電気工学科 - 電気工学教育プログラム (kumamoto-u.ac.jp) 電子工学教育プログラム 研究や教育についてまとめています。 情報電気工学科 - 電子工学教育プログラム (kumamoto-u.ac.jp) 情…

本ブログでは、外れ値がセンサ信号に含まれる系での状態推定について説明します。まず、外れ値を含む系の状態推定問題を定式化します。ここでは、外れ値の問題に対処するために複数の候補を並列に定め、その中から外れ値の影響を受けていない推定値を選択する新たな状態推定法を提案しています。 [link] H. Okajima, Y. Kaneda and N. Matsunaga: Full article: State estimation method using median of multiple candidates for observation signals including outlier…

この記事は雑記記事です。熊本大学工学部情報電気工学科学生向けの内容になっています。 クリアファイルデータの発掘 情報電気工学科チートシート クリアファイルデータの発掘 制御工学チャンネルのクリアファイルがほとんどなくなったので新しいものを作ります。それにあたって、過去のファイルを漁っていたところ、なぜかどの立場でかわかりませんが、情報電気工学科学生向けのチートシートを作っていたのが見つかりました。 自分が阪大に入学したころ、履修登録や種々の手続きは、大学に入ってからの友達に頼り切りだったように記憶しています。大学の手続きは複雑であり、これで躓くこともあるかもしれません。割とうちの学科は、ガイダ…

プログラミングっても概ねSBC(ほぼraspi)でなんかシンプルな制御か計測など単機能サーバのための簡単なプログラミングだけしかやってない。 だいぶ前はほぼperlしか書いてなかったが、gpio関係は比較的pythonの方がやりやすい状況だったので、近年は pythonとperlが半々位になってきていた。 たしかほんの一ヶ月前程度に perl を廃して python 一本にしようとか記事にしてたのだが、いまはちょっと julia 面白そうという気分になってきた。 私は高度な要求条件を駆使しなければならないことはほぼないのだが、こんな言語だと言われるとついついかじってみたくなるほどには興味心は持…

この記事ではシステム同定・動的システムのモデリングについてまとめます。関連記事リンクや解説記事へのリンクは最下部に置いています。 システム同定の概要 モデル表現 システム同定のシミュレーション例 モデルの相互変換について システム同定の動画・関連記事 自己紹介 システム同定の概要 それではモデリングシステム同定について説明をしていきたいと思います。まず制御工学において数式モデルは重要な役割を果たします。モデルベースの制御を行う場合には、まず制御対象の入出力特性に基づいてモデルを導出します。 モデルに基づいた制御設計 そしてその得られた数理モデルに対して必要な性能が得られるような制御機器を設計す…

特定の応用を持っているわけではないですが、画像処理や画像認識が比較的見栄えもよくて処理することが多いので、OpenCV の本を買ってみましたが、実際のところはほとんど matlab で書いていました。 matlab は R や Python で定番のパッケージを後追いで実装して toolbox にしているイメージがあります。もちろん商用ソフトなので、その実装にはそれなりの責任が伴うので単純な話ではないですがちょっと微妙だなと思うこともあります。 あと、産総研だと matlab はアカデミック価格では買えません。そもそもアカデミック価格というのは学生に撒き餌をするようなものですから産総研が対象外…