微分演算子による連立微分方程式の解法①
微分演算子を使った連立微分方程式の解法
次のような連立微分方程式を考えます。
![](https://diff-eq.com/wp-content/uploads/2024/06/simultaneous_differential_equations_heavisides_operational_calculate_process_huge_img1-1024x337-1.png)
マセマティカで描画するとこんな感じです。
![](https://diff-eq.com/wp-content/uploads/2024/06/simultaneous_differential_equations_heavisides_operational_calculate_vector_ps2024img1-1.jpg)
ベクトル描写画像の枠が斜めって見えてますが、それは錯覚になります。
上記の2つの式を次のように置きます。
![](https://diff-eq.com/wp-content/uploads/2024/06/simultaneous_differential_equations_heavisides_operational_calculate_process_huge_img2-1024x148-1.png)
の式を変形させてそれを
の式に代入します。
![](https://diff-eq.com/wp-content/uploads/2024/06/simultaneous_differential_equations_heavisides_operational_calculate_process_huge_img3-1024x182-1.png)
![](https://diff-eq.com/wp-content/uploads/2024/06/simultaneous_differential_equations_heavisides_operational_calculate_process_huge_img4-1024x263-1.png)
このときの斉次解は、
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_heavisides_operator1_homogeneous_solution_.png)
となるのでこれを使って以下、
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_heavisides_operator1_homogeneous_general_solution_img1.png)
さらに特殊解においては、
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_heavisides_operator1_homogeneous_spcial_solution_img2.png)
となるので一般解は以下のようになります。
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_heavisides_operator1_homogeneous_general_solution_img-1024x68-1.png)
これを一階微分してまとめると以下、
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_heavisides_operator1_homogeneous_general_solution_differential_solution_img1.png)
次に上記式のを
の式に代入して計算していきます。
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_simultaneous_differential_equations_calculate_process_img1-1024x155-1.png)
これより、
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_simultaneous_differential_equations_calculate_process_img2-1024x103-1.png)
ここでこの一般解に対してという初期条件を加えると、
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_simultaneous_differential_equations_initial_condition_img1.png)
となるので、
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_simultaneous_differential_equations_calculate_result_img1.png)
さらに、
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_simultaneous_differential_equations_calculate_result_differential_img1.png)
となるので初期条件を加えると、
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_simultaneous_differential_equations_calculate_result_differential_img2.png)
よって初期条件を加えた場合の一般解は次のようになります。
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_simultaneous_differential_equations_calculate_result_differential_img3.png)
下の画像は同じ初期条件を満たす場合のこのときの同次解と非同次解との比較をグラフで表してあり、点線が同次になります。
![](https://diff-eq.com/wp-content/uploads/2024/06/heavisides_operational_calculus_simultaneous_differential_equations_calculate_result_graph_2024img.png)
![](https://diff-eq.com/wp-content/themes/x-t9/inc/patterns/images/sample-image-gray.png)
Title Text
Lorem ipsum dolor sit amet, consectetur adipiscing elit, sed do eiusmod tempor incididunt ut labore et dolore magna aliqua.
![](https://diff-eq.com/wp-content/themes/x-t9/inc/patterns/images/sample-image-gray.png)
Title Text
Lorem ipsum dolor sit amet, consectetur adipiscing elit, sed do eiusmod tempor incididunt ut labore et dolore magna aliqua.
![](https://diff-eq.com/wp-content/themes/x-t9/inc/patterns/images/sample-image-gray.png)
Title Text
Lorem ipsum dolor sit amet, consectetur adipiscing elit, sed do eiusmod tempor incididunt ut labore et dolore magna aliqua.
-
微分演算子による連立微分方程式の解法①
カテゴリー : 微分演算子による連立微分方程式の解法①オペレータ作用素をD=d/dtのようにおき、それを代数的に取り扱うことによって連立微分方程式の解をシステマティ…