MATLAB入門(mén)教程之數值分析
加入技術(shù)交流群
>>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非線(xiàn)性方程式的實(shí)根
要求任一方程式的根有三步驟：
先定義方程式。要注意必須將方程式安排成 f(x)=0 的形態(tài)，例如一方程式為sin(x)=3，
則該方程式應表示為 f(x)=sin(x)-3?？梢?m-file 定義方程式。
代入適當范圍的 x, y(x) 值，將該函數的分布圖畫(huà)出，藉以了解該方程式的「長(cháng)相」。
由圖中決定y(x)在何處附近(x0)與 x 軸相交，以fzero的語(yǔ)法fzero('function',x0) 即可求出在 x0附近的根，其中 function 是先前已定義的函數名稱(chēng)。如果從函數分布圖看出根不只一個(gè)，則須再代入另一個(gè)在根附近的 x0，再求出下一個(gè)根。
以下分別介紹幾數個(gè)方程式，來(lái)說(shuō)明如何求解它們的根。
例一、方程式為
sin(x)=0
我們知道上式的根有 ，求根方式如下：
>> r=fzero('sin',3) % 因為sin(x)是內建函數，其名稱(chēng)為sin，因此無(wú)須定義它,選擇 x=3 附近求根
r=3.1416
評論