当前位置: 首页 > news >正文

MATLAB-振动问题:两自由度耦合系统自由振动

一、基本理论

 

二、MATLAB实现

        以下是两自由度耦合系统自由振动质量块振动过程动画显示的MATLAB程序。

clear;
clc;
close allx0 = 1;
D1 = 40;
D12 = 8;
D2 = D1;
m1 = 1;
omega0 = sqrt(D1/m1);
k1 = D12 / D1;
k2 = D12 / D2;
k = sqrt(k1 * k2);
omegazh = omega0 * sqrt(1 + k);
omegaaf = omega0 * sqrt(1 - k);
for t = 0:0.1:10%%%T = 0:0.1:t;X1 = x0 * cos((omegazh - omegaaf) * T/2) .* cos((omegazh + omegaaf) * T/2) ;X2 = x0 * sqrt(k2/k1) * sin((omegazh - omegaaf) * T/2) .* sin((omegazh + omegaaf) * T/2) ;subplot(2,2,1)plot(T, X1,'b','LineWidth',1.5);xlabel('t');ylabel('Amplitude');axis([0 10 -1 1]);subplot(2,2,3)plot(T, X2,'r','LineWidth',1.5);xlabel('t');ylabel('Amplitude');axis([0 10 -1 1]);%%%X1 = x0 * cos((omegazh - omegaaf) * t/2) .* cos((omegazh + omegaaf) * t/2) ;X2 = x0 * sqrt(k2/k1) * sin((omegazh - omegaaf) * t/2) .* sin((omegazh + omegaaf) * t/2) ;subplot(2,2,2)plot(X1, 0, '.', 'MarkerSize', 30);xlabel('x1');axis([-1 1 -1 1]);subplot(2,2,4)plot(X2, 0, '.', 'MarkerSize', 30);xlabel('x2');axis([-1 1 -1 1]);pause(0.1)
endT = 0:0.1:10;
X1 = x0 * cos((omegazh - omegaaf) * T/2) ;
X2 = x0 * sqrt(k2/k1) * sin((omegazh - omegaaf) * T/2)  ;figure
subplot(2,1,1)
plot(T, X1,'b','LineWidth',1.2);
xlabel('t');
ylabel('Amplitude');
axis([0 10 -1 1]);
title('包络');
set(gca, 'FontSize', 16);
set(0,'defaultfigurecolor','w');subplot(2,1,2)
plot(T, X2,'r','LineWidth',1.2);
xlabel('t');
ylabel('Amplitude');
axis([0 10 -1 1]);
set(gca, 'FontSize', 16);
set(0,'defaultfigurecolor','w');

        程序执行结果包括两个振子的位置随时间变换的轨迹(图1中左图)以及两个振子振动情况的动画显示(图1中右图)。 

图1  代码运行时动图显示

        图2为两个振子振动信号的包络。

图2
http://www.lryc.cn/news/389917.html

相关文章:

  • 人工智能-NLP简单知识汇总01
  • Spring Boot中的异步编程技巧
  • 深度解密Spark性能优化之道
  • 在U盘/移动硬盘上安装热插拔式Ubuntu系统,并将Docker目录挂载到NTFS硬盘
  • 商城小程序论文(设计)开题报告
  • 15. Java的 CAS 操作原理
  • 修改element-ui日期下拉框datetimePicker的背景色样式
  • Linux—— 逻辑运算符,压缩和解压缩
  • 音视频入门基础:H.264专题(6)——FFmpeg源码:从H.264码流中提取NALU Header、EBSP、RBSP和SODB
  • STM32实现按键单击、双击、长按、连按功能,使用状态机,无延时,不阻塞
  • C#之Delta并联机械手的视觉同步分拣
  • 01:Linux的基本命令
  • GNSS 载波、测距码和导航电文的关系简介
  • deepE 定位系统卡顿问题实战(一) ----------- 锁造成的阻塞问题
  • YOLOv5改进 | 主干网络 | ODConv + ConvNeXt 增强目标特征提取能力
  • TIA博途WinCC通过VB脚本从 Excel中读取数据的具体方法介绍
  • 第5篇 区块链的技术架构:节点、网络和数据结构
  • vue长列表,虚拟滚动
  • 【实战场景】记一次UAT jvm故障排查经历
  • 线性代数--行列式1
  • tensorflow神经网络
  • Python基础001
  • 【udp报文】udp报文未自动分片,报文过长被拦截问题定位
  • 某网页gpt的JS逆向
  • 【python脚本】批量检测sql延时注入
  • 在C++中如何理解const关键字的不同用法(如const变量、const成员函数、const对象等)
  • JavaSEJava8 时间日期API + 使用心得
  • 【亲测解决】Python时间问题
  • Linux屏幕驱动开发调试笔记
  • Nginx Http缓存的必要性!启发式缓存有什么弊端?