如何定义函数文件 数学软件和 与数学实验 .pptVIP

如何定义函数文件 [T,Y] = ode45(odefun,tspan,y0) [T,Y] = ode23(odefun,tspan,y0) odefun？ ode45、ode23 能解什么样的 ODE？ 如何定义函数文件 [T,Y] = ode45(odefun,tspan,y0) [T,Y] = ode23(odefun,tspan,y0) ode45、ode23 等函数可用于求解显式常微分方程 当 是向量函数时，所对应的方程即为微分方程组 odefun 举例说明 fun=inline(-2*y+2*x^2+2*x,x,y); [x,y]=ode23(fun,[0,0.5],1); 解法一：使用 inline 定义微分方程 odefun odefun 为方程右端项 f(t,y) 可以用 inline 定义（只适合于单个方程的情形） 通过函数文件定义，然后用函数句柄调用（适合所有情形） 注：自变量必须在前面，因变量在后面！ 举例说明（单个方程） function dy = myfun1(x,y) dy = -2*y+2*x^2+2*x; 解法二：通过函数文件定义微分方程 odefun 1、先编写函数文件 myfun1.m clear; [x,y]=ode23(@myfun1,[0,0.5],1); 2、编写主文件 main1.m 或直接在 Matlab 命令窗口输入上面的语句。 举例说明（方程组） 解：此时只能通过函数文件定义微分方程 odefun 例：求 , , 的数值解。 function dy = myfun2(t,y) dy = zeros(3,1); % dy must be a column vector! dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2); 1、先编写函数文件 myfun2.m clear; [T,Y]=ode45(@myfun2,[0,12],[0,1,1]); 2、编写主文件 main2.m dy = [ y(2)*y(3); -y(1)*y(3); ... -0.51*y(1)*y(2)]; 例 function dy = myfun3(t,y) dy = zeros(2,1); dy(1) = y(2) + t; dy(2) = t - 2; 例：解初值问题： ， ， function out = myfun3(t,y) out = [y(2) + t; t – 2]; function yprime = myfun3(x,y) yprime = [y(2) + x; x – 2]; clear; [T,Y]=ode45(@myfun3,[0,10],[1,1]); 2、主文件 main3.m 1、函数文件 myfun3.m 高阶常微分方程 高阶常微分方程 例： Van der Pol 初值问题 一阶常微分方程组 变量替换 化为 参数怎么处理？ 用全局变量传递 参数传递 1、函数文件 verderpol.m function xprime = verderpol(t,x) global mu; xprime = [x(2); mu*(1-x(1)^2)*x(2)-x(1)]; 2、主文件 vdp1.m clear; global mu; y0=[1;0]; mu=7; [t,x]=ode45(verderpol,[0,40],y0); % [t,x]=ode45(@verderpol,[0,40],y0); plot(t,x(:,1),r-, t,x(:,2),b-);