音響振動制御特論2018 の履歴(No.6) - Kobayashi Laboratory Homepage

履歴一覧
差分を表示
現在との差分を表示
ソースを表示
音響振動制御特論2018 へ行く。
- 1 (2018-04-26 (木) 12:16:49)
- 2 (2018-04-26 (木) 13:21:29)
- 3 (2018-04-26 (木) 17:06:56)
- 4 (2018-04-27 (金) 11:46:39)
- 5 (2018-05-09 (水) 17:09:47)
- 6 (2018-05-14 (月) 18:08:32)
- 7 (2018-05-14 (月) 18:15:00)
- 8 (2018-05-31 (木) 14:22:54)
- 9 (2018-06-12 (火) 09:27:36)
- 10 (2018-07-03 (火) 12:35:43)

授業

第1回(2018.4.26)†

参考資料：B2.pdf の 3.1節、3.2節、付録A, B
複素振幅による調和振動の表現 ... 3.1節、式(2)
圧力と粒子（体積）速度に関する基本の式 ... 付録A、式(27)と(31)
特性音響インピーダンス
波動方程式の解（ヘルムホルツ方程式からも導出できる）

第2回(2018.4.27)†

前回のつづき
- 粒子（体積）速度に関する基本の式(31)
- 特性音響インピーダンス
周波数応答実験結果とモデルの比較
- 周波数応答実験の結果（3年学生実験の装置）: spk1.dat spk2.dat
- mファイル: freqresp_distx_p.m p_tilde_down.m p_tilde_up.m

板書1

板書2

↑

第3回(2018.5.9)†

管内の圧力分布（前回の続き）:freqresp_distx_p_fixed.m
- 周波数が高くなると、（圧力の）腹と節の数が増える
- マイクの位置が節に重なる周波数で、ゲイン特性が下向きのピークを示す
- wからz,yまでの周波数応答における違い：リファレンスマイクは何度も節を横切るが、エラーセンサは周波数が十分高くなって初めて節を横切る
- ...
むだ時間要素（導波管、2センサ法、指向性センサ、etc. の基本）の伝達関数と周波数応答

2センサ法による進行波圧力成分の計算

load_data_raw_AB_fixed_rows_changed_PCB.m

pc2_jw = frdata(:,4).*exp(i*pi/180.0*frdata(:,5))/pcb_coeff_12638; % 
ps1_jw = frdata(:,8).*exp(i*pi/180.0*frdata(:,9))/pcb_coeff_12992; % dummy
ps2_jw = frdata(:,2).*exp(i*pi/180.0*frdata(:,3))/pcb_coeff_12637; % 
pc1_jw = frdata(:,6).*exp(i*pi/180.0*frdata(:,7))/pcb_coeff_12991; % dummy

A1_jw = ( exp(-i*k*dl).*ps1_jw - exp(-i*k*(ls+dl)).*pc1_jw)./( 2*i*sin(k*ls));
B1_jw = (-exp( i*k*dl).*ps1_jw + exp( i*k*(ls+dl)).*pc1_jw)./( 2*i*sin(k*ls));
A2_jw = (-exp( i*k*dl).*ps2_jw + exp( i*k*(ls+dl)).*pc2_jw)./( 2*i*sin(k*ls));
B2_jw = ( exp(-i*k*dl).*ps2_jw - exp(-i*k*(ls+dl)).*pc2_jw)./( 2*i*sin(k*ls));

板書1

板書2

↑

第4回(2018.5.14)†

前回のつづき：2センサ法による進行波圧力成分の計算 ... https://doi.org/10.5687/iscie.28.392 (22)式の下の式
指向性センサ

/home/baba/spkcont/spkcont_app.c (参考：spkcont_module.c spkcont.h):

#ifdef PCB
    As1 = sqrt(msg->as1*msg->as1 + msg->bs1*msg->bs1) / PCB_COEFF_12638;
    Ac1 = sqrt(msg->ac1*msg->ac1 + msg->bc1*msg->bc1) / PCB_COEFF_12637;
    As2 = sqrt(msg->as2*msg->as2 + msg->bs2*msg->bs2) / PCB_COEFF_12992;
    Ac2 = sqrt(msg->ac2*msg->ac2 + msg->bc2*msg->bc2) / PCB_COEFF_12991;
#else // NAGANO_KEIKI
    As1 = sqrt(msg->as1*msg->as1 + msg->bs1*msg->bs1) * NAGANO_KEIKI_COEFF;
    Ac1 = sqrt(msg->ac1*msg->ac1 + msg->bc1*msg->bc1) * NAGANO_KEIKI_COEFF;
    As2 = sqrt(msg->as2*msg->as2 + msg->bs2*msg->bs2) * NAGANO_KEIKI_COEFF;
    Ac2 = sqrt(msg->ac2*msg->ac2 + msg->bc2*msg->bc2) * NAGANO_KEIKI_COEFF;
#endif
               
// Phase
Ps1 = atan(msg->bs1/msg->as1) + M_PI * (msg->as1 < 0);
Pc1 = atan(msg->bc1/msg->ac1) + M_PI * (msg->ac1 < 0);
Ps2 = atan(msg->bs2/msg->as2) + M_PI * (msg->as2 < 0);
Pc2 = atan(msg->bc2/msg->ac2) + M_PI * (msg->ac2 < 0);

A1 = sqrt(As1*As1 - 2*As1*Ac1*cos(k*L + Ps1 - Pc1) + Ac1*Ac1) / (2*fabs(sin(k*L)));
B1 = sqrt(As1*As1 - 2*As1*Ac1*cos(k*L + Pc1 - Ps1) + Ac1*Ac1) / (2*fabs(sin(k*L)));
A2 = sqrt(Ac2*Ac2 - 2*Ac2*As2*cos(k*L + Pc2 - Ps2) + As2*As2) / (2*fabs(sin(k*L)));
B2 = sqrt(Ac2*Ac2 - 2*Ac2*As2*cos(k*L + Ps2 - Pc2) + As2*As2) / (2*fabs(sin(k*L)));

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!以下は過去の情報です!!!!!!!!!!!!!!!!!!!!!!!!!!

長さlの管の伝達マトリクス（入出力関係が非因果的）
進行波圧力振幅を入出力とするシステム（因果的）

#ref(): File not found: "2015.04.30-1.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.04.30-2.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.04.30-3.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.04.30-4.jpg" at page "音響振動制御特論2018"

↑

第3回(2015.5.29)†

ナイキストの安定判別（開ループ系が安定な場合）
むだ時間系の安定判別
MIMO版ナイキストの安定判別
ループ管の安定判別

#ref(): File not found: "2015.05.29-1.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.05.29-2.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.05.29-3.jpg" at page "音響振動制御特論2018"

↑

第4回(2015.8.26)†

2センサ法
指向性センサ

#ref(): File not found: "2015.08.26-1.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.08.26-2.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.08.26-3.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.08.26-4.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.08.26-5.jpg" at page "音響振動制御特論2018"

第5回(2015.8.27)

リニアモータの電気-音響特性に基づく管路の計算手順
CHAIN とその逆変換

#ref(): File not found: "2015.08.27-1.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.08.27-2.jpg" at page "音響振動制御特論2018"

#ref(): File not found: "2015.08.27-3.jpg" at page "音響振動制御特論2018"