布料仿真先导2-带阻尼的单个小球单摆下的拉格朗日方程列些和matlab仿真
在先导1中,我们对保守系统的双球系统进行了仿真,现在我们先从单个小球摆动入手,考虑加入了阻尼后,如果通过拉个朗日分析来得到拉格朗日函数进而得到小球的运动方程。
首先回顾一下,先导1中的拉个朗日方程使用的是保守场中的形式,也就是说所有主动力都是保守力,即保守力的做工与路径无关,通过保守力的性质我们可以得到拉格朗日量L = T - V;但是现在的情况是,出现了一个与速度相关的阻尼力,这个阻尼力不是保守力,所以不能通过拉个拉格朗日的保守场形式,需要做一下变动。
考虑单个小球单摆下,受到重力和阻尼力,重力任然可以归结为势能处理,阻尼力需要通过一般形式拉格朗日方程中的广义力来计算求得。公式推导如下所示:
这里有一个要注意的地方,就是重力我还是将其归结为势能,归类到保守场的拉格朗日方程中,如果没有阻尼力的话,保守场拉格朗日方程右边应该为0,但是现在还有一个非保守力阻尼力,所以现在拉格朗日方程的右边就是这个阻尼力对应的广义力,求解出来带入就得到系统的运行方程了。matlab代码如下所示:
%基于ode45的单小球摆动带阻尼的运动方程
%此处是基于拉个朗日函数分析法给出的微分方程求解close all;clearvars;clc;
dbstop if error;global m l g k;
m = 5;
l = 1;
g = 9.8;
k = 2;[t , y] = ode45(@vdp1,[0 25],[30*pi/180 0]);x1 = y(: , 1) .* 180/pi;
h1 = figure;
for i = 1 : 1 : length(t)set(0,'CurrentFigure',h1);title(num2str(i));scatter(0 , 0 , '*');hold on;x = l*sind(x1(i));y = -1*l*cosd(x1(i));scatter(x , y);hold on;plot([0 x(1)] , [0 y(1)] , 'r');hold on;xlim([-3 3]);ylim([-3 3]);pause(0.05);hold off;
endfunction dydt = vdp1(t,y) %y = [phi1 phi1']输入为弧度global m l g k;dydt = zeros(2 , 1);dydt(1) = y(2);dydt(2) = (-2*k*l/m)*y(2) - (g/l)*sin(y(1));
end
布料仿真先导2-带阻尼的单个小球单摆下的拉格朗日方程列些和matlab仿真相关推荐
- 通信使用matlab三元码实验报告,Matlab仿真DSBSC信号与双极性不归零码通信原理实验报告二实验...
Matlab仿真DSBSC信号与双极性不归零码通信原理实验报告二实验 Matlab 仿真 DSB-SC 信号与双极性不归零码通信原理实验报告二实验题目:Matlab 仿真 DSB-SC 信号与双极性不 ...
- 数字信号载波传输的Matlab仿真
摘要:数字信号与信息科学的迅猛发展,计算机和通信系统的迅速发展的理论和实践,通信系统是将信息从信源发送到一个或者多个目的地,用来达到各种信息传输过程的一般技术体系.数字信号调制与解调过程中起着重要的信 ...
- 单边带调制(SSB调制)的理论基础和MATLAB仿真
文章目录 前言 一.SSB调制的起源 1. 信号的幅度调制: 2. AM和DSB-SC: 3. SSB调制: 二.MATLAB仿真 1.Hilbert变换的仿真: 2."sin(t)&quo ...
- zigbee抗干扰matlab仿真,基于Zigbee的扩频通信MATLAB仿真.doc
基于Zigbee的扩频通信MATLAB仿真 天津大学仁爱学院本科生毕业设计(论文)开题报告 课题名称基于Zigbee的扩频通信MATLAB仿真系名称信息工程系专业名称通信工程学生姓名指导教师 扩频通信 ...
- 先进pid控制matlab仿真程序,先进PID控制Matlab仿真第4版-仿真程序-上交
文件名大小更新时间 先进PID控制Matlab仿真第4版-仿真程序-上交02016-03-07 先进PID控制Matlab仿真第4版-仿真程序-上交\chap102016-03-07 先进PID控制M ...
- matlab采集网口数据,使用Signal Tap II采集到的数据进行Matlab仿真(转载)
在使用FPGA进行无线通信或者进行信号处理时,一般按照这样的步骤进行: (1)利用matlab进行浮点算法仿真 (2)将matlab浮点算法转换为定点算法,进行仿真,并与浮点算法的性能进行比较 (3) ...
- dc dc变换器的建模及matlab仿真,基于Matlab的AC/DC变换器的系统建模和仿真.pdf
基于Matlab的AC/DC变换器的系统建模和仿真 Research.Developmentl 基于 Matlab的AC/DC变换器的 系统建模和仿真 1'lleM odelingand Simula ...
- matlab帽檐图像,图像分割技术与MATLAB仿真
中南民族大学 毕业论文(设计) 学院: 计算机科学学院 专业: 自动化 年级:2012 题目: 图像分割技术与MATLAB仿真 学生姓名: 高宇成 学号:2012213353 指导教师姓名: 王黎 职 ...
- 同步电机matlab设计,同步发电机matlab仿真 同步电机模型的MATLAB仿真毕业设计.doc...
同步发电机matlab仿真 同步电机模型的MATLAB仿真毕业设计 PAGE PAGE 1 同步发电机matlab仿真 同步电机模型的MATLAB仿真毕业设计 导读:就爱阅读网友为您分享以下" ...
- 最短路径算法【matlab仿真】
最短路径问题因为其问题的普遍性,以及应用的实际性,不仅是数据结构的热点问题,也是数学信息学科.计算机学科.地理信息学科等学科的一个研究热点.由于科学技术的不断进步,使得应用数学中的图论与计算机算法与结 ...
最新文章
- 融资 2000 万美元后,他竟将核心代码全开源,这……能行吗?
- Sitecore 9有什么新功能
- 不同操作系统之间的网络配置
- 萤石的云存储到期怎么启动内存卡_3款萤石智能入户单品 让普通入户门秒变智能...
- 华为动态NAT小实验演示
- 详细介绍如何在Eclipse中使用SVN
- 2020牛客多校第1场I-1 or 2一般图最大匹配带花树
- mysql三种引擎_MySQL常见的三种存储引擎
- leetcode最长递增子序列问题
- (很全面)SpringBoot 使用 Caffeine 本地缓存
- Apache Storm技术实战之3 -- TridentWordCount
- 计算机网络的各层及其协议,计算机网络的体系结构 (architecture) 是计算机网络的各层及其协议的集合...
- 计算机网络每日一题,【计算机考研】每日一题-数网操-35
- 魔域进游戏老是显示服务器繁忙,魔域2014年春节-温情卡诺萨
- 我做产品的三大思维:发散思维、纵横思维和表里思维(下篇)
- swagger2报错Illegal DefaultValue null for parameter type integer
- session与cookie有效时间
- 蛋白质二级结构预测Linux,蛋白质稳定性预测工具-Rosetta ddg_monomer
- 国外数据挖掘方面的经典博客
- 物流计算机流程图,物流管理流程图.doc
热门文章
- python识别图像里文字_使用 Python 识别并提取图像中的文字
- C#实现对即插U盘文件的自动拷贝
- putty拷贝服务器文件,Windows下拷贝Linux的文件到本地(Putty)
- java登录验证_java实现登录验证码
- 如何将NCM转换成MP3 - NCM转MP3转换软件
- 网络安全kali渗透学习 web渗透入门 Google搜索引擎的使用技巧
- 公众号管理系统 html,Java SpringBoot+Mybatis Layui+JQuery+html微信公众号后台管理系统...
- PR转场预设 放大特效带有重影效果的PR视频转场预设
- ENVI监督分类后背景值也被分成一种地物,解决方案和转移矩阵制作方法
- Au入门系列之五:轨道与轨道控制​