水平有摩擦力K的变化对运动曲线的影响

%%  水平有摩擦力K的变化对运动曲线的影响

clc;clear;global K R 

figure('numbertitle','off','name','水平有摩擦力K的变化对运动曲线的影响');

for K=0:0.5:3

    R=0.25;

    t1=[0,80];

    H1=[0;0.1];

    [t,H] = ode45(@dH,t1,H1);

    plot(t,H(:,1),'LineWidth',1.5); hold on

    axis([0 60 -0.3 1.9])

end

set(gca,'linewidth',1.2);legend('有摩擦力K=0','有摩擦力K=0.5','有摩擦力K=1','有摩擦力K=1.5','有摩擦力K=2','有摩擦力K=2.5','有摩擦力K=3');

    xlabel('t');ylabel('角度');

    title(',角度与时间关系图');

function dHdt=dH(t,H)

global K R 

dHdt=zeros(2,1);

dHdt(1)=H(2);

dHdt(2)=-sin(H(1))*(1-K*cos(H(1)))-R*H(2);

end