基于MATLAB的飞机襟副翼运动建模分析【GUI开发】

%初始化参数
l1=80;l2=250;l3=180;l4=300;%各杆长
omega1=180;
alpha1=0;
%弧度制与角度制转换
hu_du=pi/180;
jiao_du=180/pi;
%计算四杆机构位移，角速度，角加速度
for n=0:360
    theta1=n*hu_du;
    [theta,omega,alpha]=crank_rocker(theta1,omega1,alpha1,l1,l2,l3,l4);
    theta2(n+1)=theta(1);
    theta3(n+1)=theta(2);
    omega2(n+1)=omega(1);
    omega3(n+1)=omega(2);
    alpha2(n+1)=alpha(1);
    alpha3(n+1)=alpha(2);
end

% %%
n=1:361;
%%
if get(handles.radiobutton1,'value')
    plot(n,theta2*jiao_du,n,theta3*jiao_du,'k');
    title('角位移曲线图');
    xlabel('曲柄转角\theta_1/\circ')
    ylabel('角位移/\circ')
    grid on;
    hold on;
    text(140,100,'\theta_3')
    text(140,30,'\theta_2')
elseif get(handles.radiobutton2,'value')
    plot(n,omega2,n,omega3,'k');
    title('角速度曲线图');
    xlabel('曲柄转角\theta_1/\circ');
    ylabel('角速度/rad\cdots^{-1}');
    grid on;
    hold on;
    text(250,50,'\omega_2')
    text(130,50,'\omega_3')
elseif get(handles.radiobutton3,'value')
    plot(n,alpha2,n,alpha3,'k')
    title('角加速度曲线图');
    xlabel('曲柄转角\theta_1/\circ')
    ylabel('角加速度/rad\cdots^{-2}')
    grid on;
    hold on;
    text(230,0.1e4,'\alpha_2')
    text(150,-0.7e4,'\alpha_3')
elseif get(handles.radiobutton4,'value')
    x(1)=0;
    y(1)=0;
    x(2)=l1*cos(70*hu_du);
    y(2)=l1*sin(70*hu_du);
    x(3)=l4+l3*cos(theta3(70));
    y(3)=l3*sin(theta3(70));
    x(4)=l4;
    y(4)=0;
    x(5)=0;
    y(5)=0;
    plot(x,y);
    grid on;
    hold on;
    plot(x(1),y(1),'o');
    plot(x(2),y(2),'o');
    plot(x(3),y(3),'o');
    plot(x(4),y(4),'o');
    title('平面四杆机构');
    xlabel('mm')
    ylabel('mm')
    
elseif get(handles.radiobutton5,'value')
    %四杆机构运动仿真
    fig = figure('NumberTitle','off',...
        'Color','w',...
        'menubar','none',...
        'toolbar','none',...
        'Resize','off');
    title('平面四杆机构');
    xlabel('mm');
    ylabel('mm');
    grid on;
    box on;
    axis([-100 330 -100  200]);
    m = moviein(20);
    j=0;
    hold on;
    for n=1:5:360
        x(1)=0;
        y(1)=0;
        x(2)=l1*cos((n-1)*hu_du);
        y(2)=l1*sin((n-1)*hu_du);
        x(3)=l4+l3*cos(theta3(n));
        
        y(3)=l3*sin(theta3(n));
        x(4)=l4;
        y(4)=0;
        x(5)=0;
        y(5)=0;
        if n == 1
            line0 = plot(x,y);
            line1 = plot(x(1),y(1),'o');
            line2 = plot(x(2),y(2),'o');
            line3 = plot(x(3),y(3),'o');
            line4 = plot(x(4),y(4),'o');
        else
            set(line0,'XData',x,'YData',y);
            set(line1,'XData',x(1),'YData',y(1));
            set(line2,'XData',x(2),'YData',y(2));
            set(line3,'XData',x(3),'YData',y(3));
            set(line4,'XData',x(4),'YData',y(4));
        end
        j=j+1;
        m(j)=getframe;
    end
    movie(m,100);
    hold off;
end