題目：設一系統的常系數差分方程為y(n)-ay(n-1)=x(n),若系統的輸入序列為x(n)=δ(n),試求系統的輸出序列y(n)【代碼不是很懂】 ...

2.2 差分格式 列出幾個常用的數值微分公式。 引理 2.2.1 設 \(h>0\) 和 \(c\) 為常數 如果 \(g(x) \in C^2[c-h, c+h]\)，則有 \[g(c) = \frac{1}{2} [g(c-h) + g ...

Kuiken 利用相似變換，得到如下非線性微分方程 滿足如下邊界條件 其中， 表示對 求導，為普朗特數. 此方程是耦合的非線性邊值問題，在無窮遠點具有奇性. 當 時，使用Matlab的bvp4c求解如下： 將原方程轉化為一階方程組 % kuikenode.m ...

　　譜方法(Spectral Method)是配點法(Collocation Method)的一種。一般來說，配點法包括有限元方法(Finite Element)和譜方法(Spectral Method)。配點法的一般思路是：選取合適的函數基底，這些函數基底的導數都是已知的，求得疊加系數 ...

根據線性代數中求解方程組的基本知識，首先應判斷系數矩陣的秩是否和增廣矩陣的秩相等，若不等，則無解；若有解，根據秩和未知量個數的關系，判斷是唯一解還是無窮多解；若為無窮多解，其通解為齊次方程組的通解加非齊次方程組的特解。 求非齊次線性方程組Ax=b的特解，可直接使用命令A\b，求解齊次 ...

例如方程組： 法1：左除法 >> A=[3 1 -1;1 2 4;-1 4 5];b=[3.6;2.1;-1.4]; >> x=A\b x = 1.4818 -0.4606 0.3848 法2：求逆法 ...

3.2 五點差分格式 3.2.1 五點差分格式的建立 (1) 建立差分格式 將區間 \([a,b]\) 做 \(m\) 等分，記 \[h_1 = \frac{b-a}{m}, \quad x_i = a + ih_1, \quad i=0,\cdots,m \] 將區間 \([c ...

簡介 van der Pol 方程 code ...