matlab在求解方程中的.pptxVIP

会计学;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);微分方程符号解;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);求一元函数零点 fzero(fun, x0) X = fzero(fun, x0)寻找fun在x0附近的零点;;求多项式零点: roots(p);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);;;;例2 食饵-捕食者模型 其中x(t),y(t)分别表示食饵和捕食者t时刻的数量。 求r=1,d=0.5,a=0.1,b=0.02,x0=25,y0=2时的数值解，并画出x(t),y(t)的图形以及相图(x,y)。 function xdot=shier(t,x) r=1;d=0.5;a=0.1;b=0.02; xdot=diag([r-a*x(2),-d+b*x(1)])*x;;s_b.m ts=0:0.1:15; x0=[25,2]; [t,x]=ode45(shier,ts,x0); [t,x] plot(t,x),grid,gtext(x1(t)),gtext(x2(t)), plot(x(:,1),x(:,2)),grid,xlabel(x1),ylabel(x2);第16页/共22页;马尔萨斯(Malthus)人口模型; 生物学家研究发现，战争期间捕鱼量减少,食肉鱼总数会急增，数学家将鱼分成捕食者和被捕食者，分别以y2(t) , y1(t)表示它们在时刻t 的总数。假定被捕食者(食用鱼)有充分的食物。;t0=0;tf=15; y0=[20 20]; [t,y]=ode23(lotka,[t0,tf],y0); plot(t,y(:,1),t,y(:,2),r) ;洛伦兹模型与“蝴蝶效应”;Are you ready?