创意平板折叠桌动态展开Matlab源码

发表于   |   分类于

我在参加2014年大学生数学建模竞赛时，花两天时间写出来的创意平板折叠桌展开动画程序。

  • 文件 mymovie1.m 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
% 1.加上钢筋木条动画
num = 1;
for alpha = 0:0.015:asin(5/6)
    plot3([25 25 -25 -25 25],[-60 60 60 -60 -60],zeros(1,5),'-')
    hold on
    axis equal
    tt = 0:0.01:2*pi;asin(5/6);
    plot(25.*sin(tt),25.*cos(tt))
    x = linspace(-25,25,20);                    % 平面分析x坐标
    y = sqrt(625 - x.^2);                       % 平面分析y坐标
    t = sqrt(y.^2+900-60.*cos(alpha).*y);       % 轴到铰链点距离
    cosBelta =( y.^2+t.^2-30^2)./(2*t.*y);
    x1 = x;                                     % 桌脚边缘线x1
    Belta = acos(cosBelta);
    y1 =y - (60-y).*cos(Belta);                 % 桌脚边缘线y1
    z1 = - (60-y).*sin(Belta);                  % 桌脚边缘线z1
    plot3(x1,y1,z1,'r')                         % 做出桌腿边缘线
    handle = plot3(x1,y1,z1,'r');
    set(handle,'linewidth',3)
    xx = reshape([x1;x],1,[]);
    yy = reshape([y1;y],1,[]);
    zz = reshape([z1;zeros(1,20)],1,[]);
    for i = 1:2:40                              % 绘出木条
        line([xx(i),xx(i+1)],[yy(i),yy(i+1)],[zz(i),zz(i+1)])
    end
    x2 = x;                                     % 钢筋x2
    y2 = y - t.*cos(Belta);                     % 钢筋y2
    z2 = - t.*sin(Belta);                       % 钢筋z2
    plot3(x2,y2,z2);                            % 画出钢筋
    xlabel('x')
    ylabel('y')
    zlabel('z')
    view(-131,32)
    hold off
    mov(num) = getframe(gcf);
    num = num+1;
end
movie2avi(mov, 'mymovie1.avi', 'compression', 'None')
  • 文件mymovie2.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
% 2.桌腿边缘线延时动画
plot3([25 25 -25 -25 25],[-60 60 60 -60 -60],zeros(1,5),'-')
hold on
axis equal
num = 1;
tt = 0:0.02:2*pi;asin(5/6)
plot(25.*sin(tt),25.*cos(tt))
for alpha = 0:0.02:asin(5/6)
    x = linspace(-25,25,20);                    % 平面分析x坐标
    y = sqrt(625 - x.^2);                       % 平面分析y坐标
    t = sqrt(y.^2+900-60.*cos(alpha).*y);       % 轴到铰链点距离
    cosBelta =( y.^2+t.^2-30^2)./(2*t.*y);
    x1 = x;                                     % 桌脚边缘线x1
    Belta = acos(cosBelta);
    y1 =y - (60-y).*cos(Belta);                 % 桌脚边缘线y1
    z1 = - (60-y).*sin(Belta);                  % 桌脚边缘线z1
    plot3(x1,y1,z1,'r')
    xlabel('x')
    ylabel('y')
    zlabel('z')
    view(-56,22)
    mov(num) = getframe(gcf);
    num = num+1;
end
movie2avi(mov, 'mymovie2.avi', 'compression', 'None')

  • 效果图如下