基于MATLAB实现的一个仿真单摆非线性振动的程序，可以计算单摆的分岔图相图+使用说明文档.zip

% 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');