%
clc;clear;close all;
hold on;box on;xlim([0.96,1.52]);% 生成空的坐标轴
N=80; % 设置迭代次数
F1=ones(1,N); % 生成空的数组
for F=0.96:.001:1.52;
for n=1:N;
[t,x]=ode45('pendulum',[0,66],[0,n/N*2],[],F); % 求解微分方程
T(n)=x(end,2); % 提取关键点
end
F % 实时显示F的取值
plot(F1*F,T,'k.','markersize',1); % 绘图
pause(0.01); % 暂停一下
end
xlabel('\itF','Fontsize',16,'Fontname','Times new roman'); % x轴标注
ylabel('{\it\theta} \prime','Fontsize',16,'Fontname','Times new roman','Rotation',0); % y轴标注
set(gca,'Xtick',0.96:0.07:1.52); % 设置坐标轴刻度
xlim([0.96,F]);
% clc;
% [t,x]=ode45('pendulum',[0,29],[0,rand],[],0.96);
% x(end,2),plot(t,x);hold on;
% tn=1:12;
% plot(tn*2*pi,zeros(1,12),'r*')
set(gcf,'Color','w');
