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
评论