物理MATLAB可视化——弹簧摆

越月学习空间

于 2022-11-18 14:45:11 发布

阅读量1.9k

点赞数 8

文章标签： matlab python 开发语言

本文链接：https://blog.csdn.net/pmclinziyouzuo/article/details/127922245

版权

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档

物理MATLAB可视化——弹簧摆

一、模型建立
二、代码实现
三、程序运行结果

一、模型建立

在这里插入图片描述

二、代码实现

%% 
clear,clc,close all
%% 参数初始化
theta0=pi/10;
m=1;k=80;g=9.8;
L0=1;
L=L0+m*g/k;
%% 求常微分方程
F=@(t,u)[u(2);u(1)*u(4)^2+g*cos(u(3))-k/m*(u(1)-L+m*g/k);
       u(4);-2*u(2)*u(4)/u(1)-g*sin(u(3))/u(1)];
[t,u1]=ode45(F,[0:0.02:15],[L0 0 theta0 0]);
[y1,x1]=pol2cart(u1(:,3),u1(:,1));%将极坐标换为直角坐标
y1=-y1;
%% 设置坐标和其他图形参数
ymax=max(abs(y1));
axis([-1.2 1.2 -1.2*ymax 0.2]);%设置坐标范围
axis off
title('弹簧摆','FontSize',14)
hold on
R=0.055;
yy=-L0:0.01:0;
xx=R*sin(yy./L0*30*pi);%用正弦曲线表示垂直位置的弹簧
[a,r]=cart2pol(xx,yy);%将弹簧的直角坐标换成极坐标
a=a+theta0;%通过极角的转动将垂直弹簧的极坐标换回直角坐标
%% 画x轴，摆球，弹簧
line([-1,1],[0,0],'color','b','linewidth',2);
ball=line(xx(1),yy(1),'color','r','marker','.','markersize',70);
%track=line(xx(1),yy(1),'color','g','linestyle','-','linewidth',1.3);
spring=line(xx,yy,'color','k','linewidth',2);
pause(0.5)
%% 
aviobj=VideoWriter('TanHuangBai.avi');%新建叫example.avi的文件
open(aviobj); %打开
h = animatedline;
for i=1:length(t)
    yy=-u1(i,1):0.01:0;
    xx=R*sin(yy./u1(i,1)*30*pi);
    [a,r]=cart2pol(xx,yy);
    a=a+u1(i,3);
    [xx,yy]=pol2cart(a,r);
    %plot(x1,y1,'r');
    addpoints(h,x1(i),y1(i));%添加点
    set(ball,'XData',x1(i),'YData',y1(i));
    set(spring,'XData',xx,'YData',yy);
    drawnow;
    currFrame = getframe;
    writeVideo(aviobj,currFrame);
end
close(aviobj); %关闭