Clauser+plot 源代码

[T,Y]=ode45('rigid',[0 12],[0 1 1]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')

clear x=0:0.01:1; y=-5*(1-x).^(4/5)/8+5*(1-x).^(6/5)/12+5/24; plot(x,y,'*')

t=0:0.1:2*pi; x=t.*sin(3*t); y=t.*sin(t).*sin(t); plot(x,y);

[T,Y]=ode45('rigid',[0 12],[0 1 1]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')

%radbas satlins n = -5:0.1:5; a = radbas(n); plot(n,a) %------ n = -5:0.1:5; a = satlins(n); plot(n,a)

[t,x]=ode15s('c7impode',[0,2],[1,0,0,1]); plot(t,x)

[t,x]=sim('c7mdde2',[0,10]); plot(t,x)

a=[1,8,26,73,90]; n=[3.5,3.1,2.3,1.2,0.5]; t=0:0.002:10; u=90*sin(t.^2); y=glfode(u,t,a,n); subplot(211), plot(t,y); subplot(212…

x=-2:0.01:2; y=tansig(x); plot(x,y) x=-2:0.01:2; y=logsig(x); plot(x,y)

[t,y]=fdiff('c7fun5',[0,1],[1,4],50); plot(t,y)