MATLAB入門教程之數值分析
加入技術交流群
>>soln_2 = dsolve('Dy = 2*x*cos(y)^2','y(0) = pi/4')
ans= atan(x^2+1)
>>soln_3 = dsolve('Dy = 3*y + exp(2*x)',' y(0) = 3')
ans= -exp(2*x)+4*exp(3*x)
2．4非線性方程式的實根
要求任一方程式的根有三步驟：
先定義方程式。要注意必須將方程式安排成 f(x)=0 的形態，例如一方程式為sin(x)=3，
則該方程式應表示為 f(x)=sin(x)-3。可以 m-file 定義方程式。
代入適當范圍的 x, y(x) 值，將該函數的分布圖畫出，藉以了解該方程式的「長相」。
由圖中決定y(x)在何處附近(x0)與 x 軸相交，以fzero的語法fzero('function',x0) 即可求出在 x0附近的根，其中 function 是先前已定義的函數名稱。如果從函數分布圖看出根不只一個，則須再代入另一個在根附近的 x0，再求出下一個根。
以下分別介紹幾數個方程式，來說明如何求解它們的根。
例一、方程式為
sin(x)=0
我們知道上式的根有 ，求根方式如下：
>> r=fzero('sin',3) % 因為sin(x)是內建函數，其名稱為sin，因此無須定義它,選擇 x=3 附近求根
r=3.1416
評論