Masassiah Blog

現役サラリーマンのスキルアップのための読書まとめ

MATLAB・Octave

RL 直列回路の過渡現象

自己インダクタンス $L$ と抵抗 $R$ との直列回路に直流電源電圧 $V$ が印加されているときに回路に流れる電流 $i$ の時間的変化とインダクタンス $L$ 両端の電圧 $v_\text{L}$ を MATLAB・Octave で描画する。 ただし,直流電源電圧 $V$ の大きさは 1 [V],…

交流 RL 直列回路の電圧・電流波形

2020年6月13日作成,2020年6月15日更新 交流 RL 直列回路の電圧・電流波形を MATLAB・Octave で描いてみた。 RL 直列回路 抵抗 R,インダクタンス L を直列に接続した回路に交流電源電圧 e を加えたとき,電流 i が流れたとする。 交流電源電圧の波形と RL …

三角形状分布荷重を受ける片持ちはりの SFD,BMD,たわみ曲線の計算・グラフ化

2020年3月20日更新 はじめに 概要 プログラムのソースコード 計算条件 プログラムの説明 グラフの説明 せん断力図(SFD) 曲げモーメント図(BMD) たわみ曲線 はじめに 本稿では,三角形状分布荷重を受ける片持ちはりのせん断力図(SFD),曲げモーメント図…

集中荷重を受ける片持ちはりの SFD,BMD,たわみ曲線の計算・グラフ化

2020年3月20日更新 はじめに 概要 プログラムのソースコード 設定条件 グラフの説明 せん断力図(SFD) 曲げモーメント図(BMD) たわみ曲線 はじめに 本稿では,集中荷重を受ける片持ちはりのせん断力図(SFD),曲げモーメント図(BMD),たわみ曲線を MAT…

等分布荷重を受ける片持ちはりの SFD,BMD,たわみ曲線の計算・グラフ化

2020年8月11日更新 はじめに 概要 検討モデル プログラムのソースコード 設定条件 グラフの説明 せん断力図(SFD) 曲げモーメント図(BMD) たわみ曲線 はじめに 本稿では,等分布荷重を受ける片持ちはりのせん断力図(SFD),曲げモーメント図(BMD),た…

Octave の精義【第二版】

2019年11月25日更新 数値計算の結果を簡単に可視化(visualize)できることは,Octave の魅力の一つ。その魅力を存分に活用するためには,色々なノウハウが必要です。 Octave を使ったグラフィックについてノウハウを得るため,2019年10月25日に第二版第1冊…

集中荷重を受ける単純支持はりの SFD,BMD,たわみ曲線の計算・グラフ化

集中荷重を受ける単純支持はりのせん断力,曲げモーメント,たわみを MATLAB・Octave により計算し,SFD,BMD,たわみ曲線をグラフ化する方法について解説する。

等分布荷重を受ける単純支持はりの SFD,BMD,たわみ曲線の計算・グラフ化

2020年3月20日更新 はじめに 概要 プログラムのソースコード 設定条件 グラフの説明 せん断力図(SFD) 曲げモーメント図(BMD) たわみ曲線 はじめに 本稿では,等分布荷重を受ける単純支持はりのせん断力図(SFD),曲げモーメント(BMD),たわみ曲線を M…

三角形状分布荷重を受ける単純支持はりの SFD,BMD,たわみ曲線を計算・グラフ化

2020年3月20日更新 はじめに 概要 プログラムのソースコード 設定条件 グラフの説明 せん断力図(SFD) 曲げモーメント図(BMD) たわみ曲線 はじめに 本稿では,三角形状分布荷重を受ける単純支持はりのせん断力図(SFD),曲げモーメント図(BMD),たわみ…

GNU Octave 起動時に読み込む .octaverc の設定

数値解析を目的としたプログラミング言語である GNU Octave を起動したとき,Octave 環境に個人的な変更を加えるファイル「.octaverc」について説明する。 www.gnu.org GNU Octave は,MATLAB とほぼ互換性のある数値実験用プログラミング言語として使用でき…

文字抽出

MATLABを使って文字抽出を行ってみました。文字抽出というより、文字とその他の部分で画像を2値化するプログラムです。 以前のボリュームデータの可視化で使ったisosurfaceを使って実現しました。isosurfaceは3次元でしか使えない(?)ので画像データ(もちろん…

MATLAB を使って音波の伝搬を FDTD 法でシミュレーション

MATLAB で音波(sound wave)の伝搬を FDTD(Finite Difference Time Domain)法でシミュレーションしてみた。 はじめに FDTD 法による音波の伝搬解析 伝搬解析に用いたパラメータ MATLABプログラム 参考文献 更新履歴 はじめに 音波は空気その他の物質,す…

高次精度移流方程式解法 CIP 法

はじめに 流体力学分野において,東京工業大学名誉教授である矢部 孝らによって提案された高次精度移流方程式解法である Constrained Interpolation Profile Method(CIP 法)*1を試してみた。 CIP 法は多次元の問題にも適用できるが,とりあえず,一次元の…

音響FDTD1次元バージョン

MATLABで音響FDTDプログラムを作成しました.とりあえず1次元バージョンです. *自由音場における音圧p,粒子速度uに関する偏微分方程式を差分化して,逐一音圧と粒子速度を交互に計算していきます. *密度と体積弾性率は室温20度の値を用いています. *Mur…

HRIR を wav ファイルに畳み込み

HRIR*1を wav ファイル*2に畳み込む MATLAB プログラムを作成した。 load LHRIR load RHRIR [wn,fs,bits]=wavread('hoge1.wav'); voiceLR(:,1)=conv(LHRIR,wn); voiceLR(:,2)=conv(RHRIR,wn); wavwrite(voiceLR,fs,bits,'hoge2.wav'); 解説 LHRIR が左耳の H…

ボリュームデータの可視化

2020年2月11日更新 MATLAB を使って三次元のボリュームデータを可視化する方法を紹介します。 %===parameter===% th=0.5; SIGMA=10; x=0:30;y=0:30;z=0:30; [X,Y,Z]=meshgrid(x,y,z); Vol=exp(-((X-15).^2+(Y-15).^2+(Z-15).^2)/2/SIGMA/SIGMA); Vol=Vol+0.2…

俺の startup.m

俺のフレンチ,俺のイタリアンなど,「俺の○○」が流行っていた当時(2008年2月),俺の startup.m というタイトルで記事を書いてみた。 startup.m は,MATLAB 環境を変更するコマンドを記述したスタートアップファイルであり,startup.m をユーザの MATLAB …