*応用数学2B 2015年度 [#ff95d280] ** &color(green){[第1回]}; 2015.9.30 概要(ラプラス変換とフーリエ変換の関係)、ラプラス変換の定義と例 [#i8a1301a] #ref(2015.09.30-1.jpg,left,noimg,板書 #1); #ref(2015.09.30-2.jpg,left,noimg,板書 #2); #ref(2015.09.30-3.jpg,left,noimg,板書 #3); #ref(2015.09.30-4.jpg,left,noimg,板書 #4); #ref(2015.09.30-5.jpg,left,noimg,板書 #5); -Q: シラバスは配布しないのか? -A: シラバスのホームページからダウンロードしてください。[[ここ>https://www.nagaoka-ct.ac.jp/wp-content/themes/ngaoka_kousen2009/syllabus/h27/31266.pdf]]にあります。 -Q: 声をもっと張ってほしい -A: すみません。改善します。 -Q: 広義積分の極限を取るときに値が収束するか等の証明はテストで回答する際に必要か? -A: 不要です。 ** &color(green){[第2回]}; 2015.10.7 ラプラス変換の性質 [#ac8b48bf] #ref(2015.10.07-1.jpg,left,noimg,板書 #1); #ref(2015.10.07-2.jpg,left,noimg,板書 #2); #ref(2015.10.07-3.jpg,left,noimg,板書 #3); #ref(2015.10.07-4.jpg,left,noimg,板書 #4); #ref(2015.10.07-5.jpg,left,noimg,板書 #5); -Q: L[sin t]を \[\sin t = \frac{e^{it} - e^{-it}}{2i}\] で求めると s > i になるが、部分積分法で求めると s > 0 になるのは何故か? -A: どちらも、s の実部が正、となります。s > i は正しくは、Re s > Re i、つまり Re s > 0 です。「Re」は実部の意味です。s > 0 も同様。\[e^{it}, e^{-it}\] は t によらず常に振幅が 1 の単位円上の複素数であり、無限区間の積分値は収束しません。これを収束させるため、s の実部が正となる必要があります。sin t と同じ状況です。 -Q: \[e^{\alpha t} \sin \omega t \] の負に注意? -A: 話がこの部分に進んだとき、その前の部分で話すつもりだったことを思い出し、話を戻して説明を加えたため、混乱を与えたと思います。また、その説明が中途半端でした。すみません。伝えたかったことは、三角関数について知っている結果 \[ \sin (-\omega t) = - \sin \omega t, \quad \cos(-\omega t) = \cos \omega t \] が、そのラプラス変換 \[ {\cal L}\left[\sin ((-\omega) t)\right] = \frac{(-\omega)}{s^2 + (-\omega)^2} = - \frac{\omega}{s^2 + \omega^2} = {\cal L}\left[ -\sin\omega t \right] \] \[ {\cal L}\left[\cos ((-\omega) t)\right] = \frac{s}{s^2 + (-\omega)^2} =\frac{s}{s^2 + \omega^2} = {\cal L}\left[ \cos\omega t \right] \] を介しても理解できる、ということです。 ** &color(green){[第3回]}; 2015.10.21 ラプラス変換の性質(つづき)〜逆ラプラス変換 [#df2567e8] #ref(2015.10.21-1.jpg,left,noimg,板書 #1); #ref(2015.10.21-2.jpg,left,noimg,板書 #2); #ref(2015.10.21-3.jpg,left,noimg,板書 #3); #ref(2015.10.21-4.jpg,left,noimg,板書 #4); #ref(2015.10.21-5.jpg,left,noimg,板書 #5); -Q: 逆ラプラス変換時に、留数和による方法は使っても大丈夫でしょうか? -A: 大丈夫です。 ** &color(green){[第4回]}; 2015.10.28 逆ラプラス変換(つづき:部分分数展開を用いた一般的な方法)〜微分方程式への応用 [#f5bb7eff] #ref(2015.10.28-1.jpg,left,noimg,板書 #1); #ref(2015.10.28-2.jpg,left,noimg,板書 #2); #ref(2015.10.28-3.jpg,left,noimg,板書 #3); #ref(2015.10.28-4.jpg,left,noimg,板書 #4); -Q: 部分分数展開の分母からどのようにおいていいかわからない。自由度についてよくわからない。 -A: 説明だけして具体的に式を書かなかったため、理解しにくかったかもしれません。 以下、はしょった部分を補足します。足りない場合はまた聞いてください。 -- まず、例題15(1)でなぜ \[ \frac{s}{(s+2)(s+3)(s+4)} = \frac{A}{s+2} + \frac{B}{s+3} + \frac{C}{s+4} \] とおくか。そもそも、上式左辺はかけ算で、右辺は足し算です。ここでは逆ラプラス変換をすることが目的のため、簡単な1〜2次の有理関数の足し算で表したい要望があります(もし、上式を満足する定数 A, B, C が求まれば、右辺は直ちに逆ラプラス変換できる)。&br; 次に、左辺から右辺をどう構成するかについてですが、左辺の極(分母多項式=0)が右辺の極でもあるようにします。つまり、左辺に(s+a)という因子があったら、右辺に1/(s+a)を含める(ただし、未知係数をかける)。上の例では、s を -2 とすると左辺は発散するので、右辺も発散させる必要があります。つまり、共通の極を持つということです。 -- では、例題15(2)の場合はどうか。 \[ \frac{1}{(s+1)^2(s-3)} = \frac{A}{(s+1)^2} + \frac{B}{s-3} \] とおけば十分ではないか?左辺も右辺も、s = -1, 3 を共通の極として持つのだから、と思うかもしれません。しかし、自由度が足りません。これは、以下の説明とも重複しますが、左辺の分子多項式が一般の二次多項式の場合、3つの係数が含まれます。ところが右辺の分子多項式には A と B の二つの自由度しかないため、上式を満たすように A, B を定めることが一般にはできません(たまたま左辺の二次多項式が \[ a(s-3)+b(s+1)^2 = bs^2 + (2b+a)s + (b -3a) \] という特殊な形をしていれば、 \[ \frac{bs^2+(2b+a)s + (b-3a)}{(s+1)^2(s-3)} = \frac{a(s-3)+b(s+1)^2}{(s+1)^2(s-3)} \] \[ \frac{A}{(s+1)^2} + \frac{B}{s-3} = \frac{A(s-3)+B(s+1)^2}{(s+1)^2(s-3)} \] より、A = a, B = b と定まりますが)。実際、今の例では分子多項式が 1 なので、s^2 の係数を 0 とするために B = 0 でなくてはなりません。さらに、s^1 の係数を 0 とするために A も 0 となってしまい破綻します。 -- 例題15(1) で分子が s^3 の場合:定数項が必要です。これを D とすると \[ \frac{s^3}{(s+2)(s+3)(s+4)} = \frac{A}{s+2} + \frac{B}{s+3} + \frac{C}{s+4} + D \] とおけます。実際、 \begin{eqnarray} (上式右辺)&=& \frac{A(s+3)(s+4) + B(s+2)(s+4) + C(s+2)(s+3) + D(s+2)(s+3)(s+4)}{(s+2)(s+3)(s+4)} \\ &=& \frac{A(s^2 + 7s + 12) + B(s^2 + 6s + 8) + C(s^2 + 5s + 6) + D(s^3 + 9s^2 + 26s + 24)}{(s+2)(s+3)(s+4)} \\ &=& \frac{D s^3 + (A+B+C+9D)s^2 + (7A+6B+5C+26D)s + (12A+8B+6C+24D)}{(s+2)(s+3)(s+4)} \end{eqnarray} となって、D=1 となります(逆に言えば、D=0では、s^3 を表せない)。 残りの A, B, C は、s^2, s^1, s^0 の係数が全て 0 となるように定める必要がありますが、今、未知数が3で三本の(一次独立な)式があるので、一意に定まります。 今の例のように分子多項式が 3 次の場合には、その係数を表すためにA,B,C,Dの四つの自由度が必要となります(分子多項式が2次の場合は三つの自由度が必要)。 より一般的に、分子が定数 a_i を持つ任意の三次多項式 \[ a_3 s^3 + a_2 s^2 + a_1 s + a_0 \] の場合には、 \[ \frac{a_3 s^3 + a_2 s^2 + a_1 s + a_0}{(s+2)(s+3)(s+4)} = \frac{D s^3 + (A+B+C+9D)s^2 + (7A+6B+5C+26D)s + (12A+8B+6C+24D)}{(s+2)(s+3)(s+4)} \] となり、s^3, s^2, s^1, s^0 の各係数が両辺で等しくなるように A, B, C, D が a_i で表せます。 ** &color(green){[第5回]}; 2015.11.4 微分方程式への応用(つづき:初期値の取り扱い)〜たたみこみ [#r3db4c31] #ref(2015.11.04-1.jpg,left,noimg,板書 #1); #ref(2015.11.04-2.jpg,left,noimg,板書 #2); #ref(2015.11.04-3.jpg,left,noimg,板書 #3); #ref(2015.11.04-4.jpg,left,noimg,板書 #4); ** &color(green){[第6回]}; 2014.11.14 線形システムの伝達関数とデルタ関数 [#b7f2f953] #ref(2015.11.11-1.jpg,left,noimg,板書 #1); #ref(2015.11.11-2.jpg,left,noimg,板書 #2); #ref(2015.11.11-3.jpg,left,noimg,板書 #3); #ref(2015.11.11-4.jpg,left,noimg,板書 #4); □1 の黒板にδ(t)とe^{-t}のグラフと、以下追記(写真取り忘れ、すみません): インパルス応答 ... t=0の瞬間の入力による、時刻tにおける正味の出力 -Q: 丸4の重みが分からなかった -A: 説明がイマイチでした。説明したかったことは次のことです: -- 積分の中身において、\[ e^{-(t - \tau)} \] は、\[ x(\tau) \] に対する重みとして働く。たとえば x(0) と x(t) に対する重みはそれぞれ \[ e^{-t}, \quad 1 \] となり、 y(t)を定めるためには x(t) の方が x(0) よりも重視される(x(0) は時間 t の経過とともに忘却される)。ただし、これは H(s) = 1/(s+1) の場合の話。 -- 一般の H(s) の場合、x(0) が時間 t の経過とともにどうなるか? ... これを表しているのがインパルス応答 h(t) である(時刻 0 でしか値を持たないデルタ関数の影響が、時間の経過とともにどうなるかを表現しているのだから)。 なお、これは説明しませんでしたが、H(s) = 1/s の場合には、「忘却」されずに一定に保持されつづけます。 //■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ &color(black,red){&size(20){!!!!!!!!!!!!!!!!!!!!!!!!!!!!!以下は過去の情報です!!!!!!!!!!!!!!!!!!!!!!!!!!};}; ** &color(green){[第7回]}; 2014.11.19 中間テスト [#ha213f87] ** &color(green){[第8回]}; 2012.11.21 フーリエ解析の概要、周期2πのフーリエ級数 [#kef8fc2b] #ref(2014.12.03-1.jpg,left,noimg,板書 #1); #ref(2014.12.03-2.jpg,left,noimg,板書 #2); #ref(2014.12.03-3.jpg,left,noimg,板書 #3); #ref(2014.12.03-4.jpg,left,noimg,板書 #4); #ref(2014.12.03-5.jpg,left,noimg,板書 #5); -Q: フーリエ級数展開とフーリエ変換の違いは、周期が有限か無限かだけで、やっていることは同じ? -A: おおざっぱに言えば、同じです。もう少し正確に言えば、 とても長〜い周期 2l の関数 f(x) が与えられたとき、これに対する複素フーリエ係数 C_n が、f(x) の一周期分だけ取り出した関数(g(x)とします)に対するフーリエ変換 G(u) に対応します。もっと厳密に言うと次のようになります。 \[ \Delta u = \frac{\pi}{l} \] とおくと、l が十分大きいのでこれは十分小さくなります。 次に、フーリエ変換の引数 u に対して \[ n = \frac{u}{\Delta u} \] とおきます。(フーリエ変換の引数 u が変わることが、複素フーリエ係数 C_n の n が変わることに対応します) 以上のもとで、 \[ \frac{1}{2\pi} G(u) \Delta u \simeq C_n \] が近似的に成り立ちます。 f(x) を具体的な関数の例(たとえば、-1<= x <= 1 の範囲で 1、それ以外で 0 という関数が周期 2l で繰り返す、とか)で考えるとイメージし易いかもしれません。 複素フーリエ級数が、周期を∞にもっていくとフーリエ変換の反転公式になる、という証明にも関連します。 -Q: 最終成績の評価方法について -A: 最終成績は、中間試験と期末試験の平均値です。授業の最初の方で伝えたつもりでしたが、伝えていなかったかもしれません。次回、口頭でアナウンスします。 -Q: 中間試験結果の分布を知りたい -A: &ref(hist.jpg); ** &color(green){[第9回]}; 2014.12.10 一般の周期のフーリエ級数 [#t3e94843] #ref(2014.12.10-1.jpg,left,noimg,板書 #1); #ref(2014.12.10-2.jpg,left,noimg,板書 #2); #ref(2014.12.10-3.jpg,left,noimg,板書 #3); #ref(2014.12.10-4.jpg,left,noimg,板書 #4(二箇所誤りあり)); #ref(2014.12.10-5.jpg,left,noimg,板書 #4(赤字で誤り訂正)); ...授業後、2箇所の誤りを指摘してもらいました。申し訳ありません。 -Q:anまたはbnが0となることが図から明らかな場合、計算せずにan=0などと書いて良いか? -A:良いです。 ** &color(green){[第10回]}; 2014.12.24 複素フーリエ級数 [#e7b9b7b6] #ref(2014.12.24-1.jpg,left,noimg,板書 #1); #ref(2014.12.24-2.jpg,left,noimg,板書 #2); #ref(2014.12.24-3.jpg,left,noimg,板書 #3); #ref(2014.12.24-4.jpg,left,noimg,板書 #4); #ref(2014.12.24-5.jpg,left,noimg,板書 #5); -Q: 期末テストの過去問を配るか? -A: 配ります。 -Q: 昨年度の電子制御工学実験のレポート課題で既知の内容かもしれない -A: ご指摘ありがとうございます。他科目との重複については以前から聞いていましたが、レポートの内容まで見たことはありませんでした。取り寄せてみます。 ** &color(green){[第11回]}; 2015.1.7 フーリエ級数の偏微分方程式への応用 [#d1994ab7] #ref(2015.01.07-1.jpg,left,noimg,板書 #1); #ref(2015.01.07-2.jpg,left,noimg,板書 #2); #ref(2015.01.07-3.jpg,left,noimg,板書 #3); #ref(2015.01.07-4.jpg,left,noimg,板書 #4); -Q:ただk,l,f(x)をu(x,t)に代入するだけでなく、そこまでの過程も知りたい -A:もっともな指摘です。今日は説明する時間がなくすみません。解の導出の概要がわかれば、あとは教科書の説明をフォローできるのではないかと思います。以下のpdfファイルは両端の温度が0でない場合への拡張を扱っていますが、温度0の場合についても参考になる部分があると思います。 両端の温度が0でない場合について、一昨年度のQAを以下に記載します(pdfファイルは一部誤りの修正版): -Q:熱伝導問題の解について、両端の温度が 0 でない場合、sin が cos になったりするのか? -A:講義の際、答えがあやふやだったため補足します。 結論は、cos になることはありません。両端の温度が 0 でない場合は、その両端の温度を直線でつないだ温度分布が、最終的にいきつく温度分布となります。 その最終的な温度分布からの差分が、フーリエ級数で表現されることになるので、フーリエ級数の部分は、両端の温度が0のときと同じ形(t->∞で0に収束する指数関数と、sin関数を持つ)になります。詳細は、以下のpdfファイルを見てください。 #ref(ex_fixed.pdf); ** &color(green){[第12回]}; 2015.1.14 フーリエ変換と積分定理 [#k0d14fec] #ref(2015.01.14-1.jpg,left,noimg,板書 #1); #ref(2015.01.14-2.jpg,left,noimg,板書 #2); #ref(2015.01.14-3.jpg,left,noimg,板書 #3); #ref(2015.01.14-4.jpg,left,noimg,板書 #4); #ref(2015.01.14-5.jpg,left,noimg,板書 #5); -Q:卒研が忙しい -A:がんばってください。 ** &color(green){[第13回]}; 2015.1.21 フーリエ変換の性質 [#p602c46e] #ref(2015.01.21-1.jpg,left,noimg,板書 #1); #ref(2015.01.21-2.jpg,left,noimg,板書 #2); #ref(2015.01.21-3.jpg,left,noimg,板書 #3); #ref(2015.01.21-4.jpg,left,noimg,板書 #4); #ref(2015.01.21-5.jpg,left,noimg,板書 #5); #ref(2015.01.21-6.jpg,left,noimg,板書 #6); -Q:基本的な関数のフーリエ変換が5つ出て来たが、これらは実用的な場面が多くあるのか? -A:5つのうち、(iii)は次回、熱伝導方程式の解を構成する際に出てきます。熱伝導問題に限定して考えると、フーリエ変換される関数f(x)は棒の初期温度分布に対応します。すると質問は、初期温度分布としてどのようなものが実用上多くありそうか、ということになります。たとえば棒のある一部の区間だけ100度で外が0度のような初期状態がステップ状の関数Paに対応します。このような単純な状況は非現実的であり、実用的でない、と言えなくもありませんが、複雑な現象を検討する前にまず単純な状況を検討するという進め方はよくあります。教科書p.110の練習問題2の2には、(ダクト中を伝播する音波など)一次元の波動方程式が載っています。この場合、フーリエ変換される対象は、初期圧力分布です(ただし二つの進行方向に対応して二つの関数があります)。やはりこの場合も、ステップ状の初期圧力分布は、他の複雑な現象を考える上で基礎を与えるため、教科書に良く載っています。 答えになっていない場合は、また聞いてください。 ** &color(green){[第14回]}; 2015.2.4 偏微分方程式への応用 [#u52ddb3a] #ref(2015.02.04-1.jpg,left,noimg,板書 #1); ... □1を飛ばして□2 から書き始めてしまいました。 #ref(2015.02.04-2.jpg,left,noimg,板書 #2); #ref(2015.02.04-3.jpg,left,noimg,板書 #3); #ref(2015.02.04-4.jpg,left,noimg,板書 #4); -Q:中間試験と期末試験の再試は別々に行うのか? -A:別々に行います(それぞれの難易度は本試験と同じです) -Q:試験で u(x,t) を求めた後、初期条件を満たすことを示す必要はあるか? -A:ありません。 -Q:偏微分方程式が変わったら、解も変わるのか? -A:変わります。たとえば正の実数 k で一般化した偏微分方程式 \[ \frac{\partial u}{\partial t} = k \frac{\partial^2 u}{\partial x^2} \] に対する解は \[ u(x,t) = \frac{1}{2\sqrt{\pi k t}} \phi(x) \ast e^{-\frac{x^2}{4kt}} \] となります。