急求MATLAB解常微分方程组数值解的解决方案

[急求物业公司薪酬方案]去文库，查看完整内容> 内容来自用户：6361xueyu 物业薪酬管理制度前言鉴于物业公司目前运行的薪酬制度结构单一、缺乏激励机制，在吸引优秀人才、激励员工业绩和成长、留住优秀...+阅读

急求MATLAB解常微分方程组数值解的解决方案

调用方法：[t,x] = ode45(eqx,[0 0.2] , [0 1 1.1 1.25])x的初始值写在一个矩阵里. 另外你的函数定义时K少一维K

(4)，我给你加了一个： function xdot = eqx(t,x); m = 0.5; alf = 0.15; K=[0.0265,0.1,0.15,0.12]; Q=20; a=0.8; xdot = zeros(4,1); xdot

(1) = Q - alf * x

(1) * x

(2) + a*m*(x

(1)+x

(2)+x

(3)+x

(4)); xdot

(2) = x

(2) * (-m + K

(2) * alf * x

(1) - alf * x

(3)); xdot

(3) = x

(3) * (-m + K

(3) * alf * x

(2) - alf * x

(4)); xdot

(4) = x

(4) * (-m + K

(4) * alf * x

(3));

matlab常微分方程组求解

上面两题用dsolve都解不出，用ode45数值解法如下：

1题中x1=i,x2=s，则：

syms x1(t) x2(t)

lamda=0.5;u=0.1;

V = odeToVectorField(diff(x1) ==lamda*x1*x2-u*x2,diff(x2) ==-lamda*x1*x2+0.1)

M = matlabFunction(V,'vars', {'t','Y'})

%计算区间[0 1]，初值[0 0]

options = odeset('OutputFcn',odeplot);

sol = ode45(M,[0 1],[0 0],options)

2题中x=x1,y=x2,z=x3，则：

syms x1(t) x2(t) x3(t)

V = odeToVectorField(diff(x1) ==-8/3*x1+x2*x3,diff(x2) ==-10*x2+10*x3,diff(x3)==-x1*x2+28*x2-x3)

M = matlabFunction(V,'vars', {'t','Y'})

%计算区间[0 1]，初值[0 0 1.0e-8]

options = odeset('OutputFcn',odeplot);

sol = ode45(M,[0 4],[0 0 1.0e-8],options)

figure

x = linspace(0,4,100);

y = deval(sol,x,1:3);

plot(y(1,:),y(2,:))

title('x-y相图')

xlabel('x'),ylabel('y')

grid on

figure

plot(y(1,:),y(3,:))

title('x-z相图')

xlabel('x'),ylabel('z')

grid on

figure

plot3(y(1,:),y(2,:),y(3,:))

title('x-y-z相图')

xlabel('x'),ylabel('y'),zlabel('z')

grid on

编写Matlab程序求常微分方程组的解

先在一个空白的M文件里创建一个M函数，文件名要与函数名相同，例如：文件名：differential.m函数内容：function [dxdt] = differential(t,x)% x：二维列向量，x

(1)代表f(t)的值，x

(2)代表g(t)的值% dxdt：与x对应的对时间导数项% 这里假设a=b=1a = 1;b = 1;dxdt = zeros(2,1);dxdt

(1) = a*x

(1)+b*x

(2);dxdt

(2) = -a*x

(1)+b*x

(2);endM函数编写完成后，在相同文件目录下，再建一个空的M文件，文件名任意，写入如下内容：t0 = 0； % 起始时刻tf = 10； % 终止时刻f0 = 10; % f的初始值g0 = 10; % g的初始值Hfunc = differential； % 创建函数句柄[T,Y] = ode45(Hfunc,[t0 tf],[f0 g0])； % 解算%运行完后，结果以列向量形式存储，时间存在T中，f存在Y的第一列中（Y(:,1)),g存在Y的第 %二列中（Y(:,2)），可以执行如下绘图指令查看结果figure

(1)grid onplot(T,Y(:,1)） % 绘制t-f曲线figure

(2)grid onplot(T,Y(:,2)） % 绘制t-g曲线...