代码搜索结果
找到约 12,326 项符合
M 的代码
shier2.m
function dy=shier2(t,y)
dy=zeros(2,1);
dy(1)=y(1)*(0.9-0.1*y(2));
dy(2)=y(2)*(-0.6+0.02*y(1));
eq1.m
function dy=eq1(x,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=1/5*sqrt(1+y(1)^2)/(1-x);
ff6.m
x0=0;xf=0.99999;
[x,y]=ode23('eq1',[x0 xf],[0 0]);
Y=0:0.01:2;
plot(1,Y,'g.')
hold on
plot(x,y(:,1),'b*')
shier1.m
function dx=shier1(t,x)
dx=zeros(2,1);
dx(1)=x(1)*(0.7-0.1*x(2));
dx(2)=x(2)*(-0.8+0.02*x(1));
ff2.m
y=dsolve('D2y+4*Dy+29*y=0','y(0)=0,Dy(0)=15','x')
chase3.m
t0=0;tf=10;
[t,y]=ode45('eq3',[t0 tf],[0 0]);
T=0:0.1:2*pi;
X=10+20*cos(T);
Y=20+15*sin(T);
plot(X,Y,'-')
hold on
plot(y(:,1),y(:,2),'*')
ff1.m
dsolve('Du=1+u^2','t')
vdp1000.m
function dy=vdp1000(t,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=1000*(1-y(1)^2)*y(2)-y(1);
eq2.m
function dy=eq2(t,y)
dy=zeros(2,1);
dy(1)=5*(1-y(1))/sqrt((1-y(1))^2+(t-y(2))^2);
dy(2)=5*(t-y(2))/sqrt((1-y(1))^2+(t-y(2))^2);
chase4.m
t0=0;tf=100;
[t,y]=ode45('eq4',[t0 tf],[0 0]);
T=0:0.1:2*pi;
X=10+20*cos(T);
Y=20+15*sin(T);
plot(X,Y,'-')
hold on
plot(y(:,1),y(:,2),'*')