1)-k2*s(2)];s0 = [0; 0]; % 初始条件，注意自行设置tf = 10; % 仿真时间，可根据需要自行修改[t,s] = ode45(ds,[0 tf],s0); % 使用ode求解器求解微分方程% 绘图plot(t,s)xlabel \itt; ylabel '{\its}_1, {\its}_2'legend('{\its}_1(\itt)',...
该二阶微分方程可以ode（）函数求得，其方程的数值解，然后用plot（）函数绘制其x—y、x—dy/dx函数图形，以及y—dy/dx函数的相平面图形。主要执行代码
用matlab求解微分方程组可以用ode（）函数求得其数值解。求解方法：1、首先，自定义微分方程组函数odefun(t,x)，其主要内容 dy(1)=x(1)*(1-x(1))-x(2)+1/x(1);dy(2)=x(2)*(-1+x(1));2、再用ode（）函数求得其数值解。使用方法 [t,x]=ode45(@odefun,[0 20],y0);3、用...
[t0,Y0]=ode45('dif_dsdedi',tspan,y0)[t1,Y1]=ode45('dif_dsdedi',tspan,y1)figure subplot(2,1,1)plot(t0,Y0(:,1),'-b',t0,Y0(:,2),'-g',t0,Y0(:,3),'-k')grid on legend('s','e','i')xlabel('t-time')title('y0=[0.5 0 0.5 ]')subplot(2,1,2)plot(...
dy = @(x,y) -1303.5*(1-y./(y+4.2442)).*y;y0 = 200; % 初始条件 [X,Y] = ode45(dy,[0 0.05],y0);Y2 = Y./(Y+4.2442);figure plot(X,Y)ylabel('y1');figure plot(X,Y2)ylabel('y2');
>> y=dsolve('Dy=y-1','y(2)=10','x')y =9*exp(-2)*exp(x) + 1
y(1)=x1*(1-x1/150000-0.5*x2/30000);y(2)=x2*(-1-x2/30000+0.5*x1/150000-2*x3/2200);y(3)=x3*(-1-x3/2200+2*x2/30000);2解方程并画图 [t,x]=ode23(@fun,[0,1],[121000,27000,2000]);hold on plot(t,x(:,1))plot(t,x(:,2))plot(t,x(:,3))3结果...
H(jw)的虚部与实部之比的反正切是与系统的相频特性有关。具体画法当给定c,k,a数值时即可进行。对于二阶系统的情况与(1)类似，先要给出二阶系统的微分方程：my''+cy'+ky=ax (2) 假定初值为0。以下按一阶系统完全一样的讨论方法做下去即可，不再细述。不知道这是不是您要求的解答。
题主的问题，可以这样解决：第一步：自定义微分方程组函数 func=@(t,x)[-5*x(1)+x(2)+1;x(1)-0.6*x(2)+1];第二步：确定x0的初值及时间范围 x0=[1;1];tspan=[0,10];第三步：使用ode45函数，求解t、x、y [t,x]=ode45(@(t,x)func(t,x),tspan,x0);第四步：使用plot...
)];[t,x]=ode45(dx,[0 1],[0.5 0.5]);plot(t,x(:,2))你现在的问题有很多不清楚的地方：1、第一个表达式很乱，我不确定是不是正确理解了你的意图；2、没给初始条件；3、常数n是多大？这个微分方程似乎属于刚性的，不同的n可能会有影响，而使用不同的ode求解器结果可能也不一样。