单摆的动画表示源码程序 - matlab其它源码 - 谷速源码
下载频道> 资源分类> matlab源码> 其它源码> 单摆的动画表示源码程序

标题：单摆的动画表示源码程序
分享到：

所属分类: 其它源码 资源类型： 文件大小: 31.76 KB 上传时间: 2016-01-24 22:12:31 下载次数: 3 资源积分：1分 提 供 者: matlab源码 单摆的动画表示源码程序
内容：
单摆的动画表示源码程序，程序员在编程的过程中可以参考学习使用，希望对IT程序员有用，此源码程序简单易懂、方便阅读，有很好的学习价值！
部分代码如下：
% Parameters and initial conditions
%   mass, link length, initial angles, simulation time
m = 1;
L = 1;
theta1 = 3*pi/4;
theta2 = 3*pi/8;
t = linspace(0, 10, 200);
 
% Solving ODE of a double pendulum
[T,Y] = ode45(@(t, x) double_pendulum(t, x, m, L), ...
    t, [theta1, theta2, 0, 0]);
 
% Calculating joint coordinates for animation purposes
x = [L*sin(Y(:,1)),  L*sin(Y(:,1))+L*sin(Y(:,2))];
y = [-L*cos(Y(:,1)), -L*cos(Y(:,1))-L*cos(Y(:,2))];
 
% Convert radians to degrees
ang = Y(:,1:2)*180/pi;
 
% Set up first frame
figure('Color', 'white')
subplot(2,1,1)
plot(T, ang, 'LineWidth', 2)
hh1(1) = line(T(1), ang(1,1), 'Marker', '.', 'MarkerSize', 20, ...
    'Color', 'b');
hh1(2) = line(T(1), ang(1,2), 'Marker', '.', 'MarkerSize', 20, ...
    'Color', 'r');
xlabel('time (sec)')
ylabel('angle (deg)')
 
subplot(2,1,2)
hh2 = plot([0, x(1,1);x(1,1), x(1,2)], [0, y(1,1);y(1,1), y(1,2)], ...
    '.-', 'MarkerSize', 20, 'LineWidth', 2);
axis equal
axis([-2*L 2*L -2*L 2*L])
ht = title(sprintf('Time: %0.2f sec', T(1)));
 
% Get figure size
pos = get(gcf, 'Position');
width = pos(3);
height = pos(4);
 
% Preallocate data (for storing frame data)
mov = zeros(height, width, 1, length(T), 'uint8');
 
% Loop through by changing XData and YData
for id = 1:length(T)
    % Update graphics data. This is more efficient than recreating plots.
    set(hh1(1), 'XData', T(id), 'YData', ang(id, 1))
    set(hh1(2), 'XData', T(id), 'YData', ang(id, 2))
    set(hh2(1), 'XData', [0, x(id, 1)], 'YData', [0, y(id, 1)])
    set(hh2(2), 'XData', x(id, :), 'YData', y(id, :))
    set(ht, 'String', sprintf('Time: %0.2f sec', T(id)))
    
    % Get frame as an image
    f = getframe(gcf);
    
    % Create a colormap for the first frame. For the rest of the frames,
    % use the same colormap
    if id == 1
        [mov(:,:,1,id), map] = rgb2ind(f.cdata, 256, 'nodither');
    else
        mov(:,:,1,id) = rgb2ind(f.cdata, map, 'nodither');
    end
end
 
% Create animated GIF
imwrite(mov, map, 'animation.gif', 'DelayTime', 0, 'LoopCount', inf)
 
%%
% *Animated GIF*
%
% <<../animation.gif>>
 

文件列表（点击上边下载按钮，如果是垃圾文件请在下面评价差评或者投诉）:

单摆的动画表示源码程序/
单摆的动画表示源码程序/Animation/
单摆的动画表示源码程序/Animation/html/
单摆的动画表示源码程序/Animation/Animation.m
单摆的动画表示源码程序/Animation/double_pendulum.m
单摆的动画表示源码程序/Animation/html/Animation.html
单摆的动画表示源码程序/Animation/html/Animation.png
单摆的动画表示源码程序/Animation/html/Animation_01.png

关键词： 单摆 源码 程序

相关推荐

matlab源码下载排行

Top_arrow
回到顶部

友情链接

matlab源码 网站源码 编程语言 网络编程 手机通讯编程 KUNAG 官网 KUNAG 电商 更多
联系方式| 版权声明| 招聘信息| 广告服务| 银行汇款| 法律顾问| 兼职技术| 付款方式| 关于我们|
网站客服网站客服 程序员兼职招聘 程序员兼职招聘
沪ICP备19040327号-3
公安备案号：沪公网安备 31011802003874号
库纳格流体控制系统（上海）有限公司 版权所有
Copyright © 1999-2014, GUSUCODE.COM, All Rights Reserved