Masassiah Blog

Masassiah のブログです。主に読書で得た気づきをまとめています。

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 の設定

2020年5月14日更新 GNU Octave を起動したとき,Octave 環境に個人的な変更を加えるファイルについて説明する。 www.gnu.org Windows10 の場合,「.octaverc」という名前のファイルに,Octave 環境に個人的な変更を加えるコマンドを記載し,ホームディレクト…

文字抽出

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

音響FDTD2次元

以前に書いた音響FDTD1次元バージョンを拡張して、同じくMATLABで音響FDTD(Finite Difference Time Domain)2次元バージョンのプログラムを作ってみた。 音圧と粒子速度に関する運動方程式と連続の式を、空間、時間領域で差分方程式に展開して、逐次計算して…

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

2021年3月7日更新 はじめに 流体力学分野において,東工大の矢部らによって提案された高次精度移流方程式解法である Constrained Interpolation Profile Method(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

俺のMATLAB環境を変更するコマンドを記述したスタートアップファイル(startup.m)です. startup.mをユーザのMATLABパス上に保存することで,MATLAB環境をカスタマイズすることができます.MATLABは起動時にstartup.mファイルが存在すれば,その設定を読み込…