6.matlab在求解方程中的应用.pptVIP

在求解方程中的应用 许海洋 * * 6.1 代数方程 6.2 微分方程 6.3 非线性方程求解 求多项式零点 6.4 常微分方程数值解法 6.1 代数方程 solve (f,t) 功能：对变量t 解方程f=0，t 缺省时默认为x 或最接近字母x 的符号变量。 例：求解一元二次方程f=a*x^2+b*x+c的实根 syms a b c x f=a*x^2+b*x+c; solve (f,x) ans= [1/2/a*(-b+(b^2-4*a*c)^ (1/2))] [1/2/a*(-b-(b^2-4*a*c)^ (1/2))] g1=x+y+z-1; g2=x-y+z-2; g3=2*x-y-z-1; [x,y,z]=solve(g1,g2,g3) x = 2/3 y = -1/2 z = 5/6 求解方程组 6.2 微分方程 dsolve(‘eq’, ’s1’, ’s2’,…, ’x’) 其中eq为方程；s1,s2,……为初始条件，缺省时给出含任意常数c1,c2,……的通解；x为自变量，缺省时默认为t 。 例：求微分方程的通解 dsolve(Dy=1+y^2) ans= tan(t+c1) y的一阶导数—— Dy, y的二阶导数—— D2y 微分方程符号解 y = dsolve(Dy=1/(1+x^2)-2*y^2,y(0) = 0,x) y = 2*x/(2*x^2+2) 符号解： pretty(y) y(x)= x / (1 + x 2) ezplot(y,-6,6) syms P t P=dsolve(DP=0.02*P*(1-P/500),P(0)=76) P = 500/(1+106/19*exp(-1/50*t)) ezplot(P,0,200) pretty(P) 例 P’ =0.02P ( 1 – P/500 ) 初始条件: P(0)=76 求一元函数零点 fzero(fun, x0) X = fzero(fun, x0)寻找fun在x0附近的零点 6.3 非线性方程求解 已知函数 f(x)，求 x 使 f(x)=0 x y (k = 1,2,·······) 例 x = tan x f=inline(x-tan(x)) x=fzero(f,3.5) x = 4.4934 求多项式零点: roots(p) 例1 求方程 x3 – 4x + 5 = 0 的解 P=[1 0 -4 5]; roots(P) ans = -2.4567 1.2283 + 0.7256i 1.2283 - 0.7256i 例2 求 x3 –8x2 +6x– 30=0的解 P=[1 -8 6 -30]; r=roots(P) r = 7.7260 0.1370 + 1.9658i 0.1370 - 1.9658i 6.4 常微分方程数值解法 求解微分方程初值问题 在区间[t0,tf]的数值解的函数为 [t,x]=ode23(‘fun’,[t0,tf],x0,’选项’) ----2/3阶RKF(Runge-Kutta-Felhberg)法 [t,x]=ode45(‘fun’,[t0,tf],x0,’选项’) ----4/5阶RKF(Runge-Kutta-Felhberg)法 其中函数fun的编写格式如下: function xdot=fun(t,x) 编写函数文件verderpol.m： function xprime = verderpol(t,x) global MU xprime = [x(2);MU*(1-x(1)^2)*x(2)-x(1)]; 再在命令窗口中执行： global MU MU = 7; Y0=[1;0] [t,x] = ode45(verderpol,0,40,Y0); x1=x(:,1);x2=x(:,2); plot(t,x1,t,x2) 图形结果 例2 食饵-捕食者模型 其中x(t),y(t)分别表示食饵和捕食者t时刻的数量。 求r=1,d=0.5,a=0.1,b=0.02,x0=25,y0=2时的数值解，并画出x(