%spring2.m-用解析法求二自由度阻尼系统的自由振动(先求系统状态方程组)
m1=1;m2=9;k1=4;k2=2; 
c1=input('c1=');c2=input('c2=');  %输入阻尼系数 
x0=[1;0];xd0=[0;-1];  %初始条件
tf=50;dt=0.1;  %时间向量
M=[m1,0;0,m2];K=[k1+k2,-k2;-k2,k2];  %二阶参数矩阵
C=[c1+c2,-c2;-c2,c2];
A=[zeros(2,2),eye(2);-M\K,-M\C];  %四阶参数矩阵
y0=[x0;xd0];  %四元变量的初始条件
for  i=1:round(tf/dt)+1   %设定计算点,作循环计算
    t(i)=dt*(i-1);
    y(:,i)=expm(A*t(i))*y0;   %循环计算矩阵指数
end
subplot(2,1,1),plot(t,y(1,:)),grid
subplot(2,1,2),plot(t,y(2,:)),grid

matlab用解析法求二自由度阻尼系统的自由振动(先求系统状态方程相关推荐

  1. matlab用卷积积分求单自由度阻尼系统响应(零极点模型)

    %spring1.m-用卷积积分求单自由度阻尼系统的强迫振动(先求系统的零极点模型) m=1;c=4;K=100;dt=0.015;  %dt为卷积积分时的采样间隔 w0=sqrt(K/m);  %固 ...

  2. matlab单自由度系统,【Matlab代码】外力单自由度阻尼系统的数值解

    function dx=odefun2(t,x) dx=zeros(2,1);%存储x1,x2的导数 dx(1)=x(2);%第一个方程 dx(2)=-20*x(2)-100*x(1);%第二个方程 ...

  3. 线性二自由度汽车模型对前轮角阶跃输入的稳态响应分析(matlab/simulink)

    线性二自由度汽车模型对前轮角阶跃输入的稳态响应分析(matlab/simulink) 一.线性二自由度汽车模型 为了便于掌握操纵稳定性的基本特性,我们将多自由度的整车模型简化成为二自由度模型.主要有以 ...

  4. 嵌入式 ARM 汇编编程例题(二维数组按规律求和,求两数 gcd / lcm,求数组 min / max,字符串复制,排序)

    文章目录 0x00. 整数加减乘除 0x01. 一维数组按某种规律求和 0x01.1 求 1~100 之和 0x01.2. 求一维数组的和 0x01.2. 求一维数组的所有奇数的和 0x02. 二维数 ...

  5. matlab 二自由度振动仿真,基于Adams和Matlab二自由度系统振动的仿真分析研究

    化学工程与装备 2017 年 第 11 期 28 Chemical Engineering & Equipment 2017 年 11 月 基于 Adams 和 Matlab 二自由度 系统振 ...

  6. 多自由度有阻尼matlab,多自由度阻尼系统固有振型的MATLAB求解程序

    多自由度阻尼系统固有振型的MATLAB求解程序function [v,w,zeta]=vbr_sf(m,d,k) %vbr_sf  vbr_sf(m,d,k) %        [v,w,zeta]= ...

  7. matlab 自由度,matlab二自由度系统振动

    利用Adams 和Matlab 对二自由度系统振动进行仿真与分析 一.实验思想 Adams 是一种可以对一些典型运动进行高效仿真的软件,本实验是利用Adams 对二自由度系统振动进行仿真及分析,再和理 ...

  8. 二自由度系统强迫振动matlab,清华大学出版社-图书详情-《振动力学基础与MATLAB应用》...

    前言 振动是客观世界最普遍的运动形式之一,在自然世界.工程领域.社会活动和日常生活中,普遍存在着物体往复运动或空间状态往复变化的振动现象.振动力学已经成为机械.航空.土木.水利.动力和交通运输等工程领 ...

  9. 二自由度振动仿真:matlab直接解微分方程virtual.lab motion仿真

    一.二自由度振动解析法与多种数值算法见:https://blog.csdn.net/lijil168/article/details/67647924 二.用matlab符号建立微分方程,并用龙格库塔 ...

最新文章

  1. A Grain of Sand 一粒沙子
  2. c#退出窗口跳转_关于winform如何如在关闭一个窗口时打开另外一个窗口
  3. 每次digital painting 之后,都可以把作品放到这里,比较好看,也和nft相关度比较大
  4. RecyclerView复用item导致数据混乱
  5. Opserver配置Redis、SqlServer监控
  6. Python学习之路——Python基础之运算符
  7. 2.4. myisamchk — MyISAM Table-Maintenance Utility
  8. 聚类算法 距离矩阵_论文阅读9——AP聚类算法
  9. 停就删稿!这个“可怕”的 GitHub 开源项目值得一试
  10. eclips调试linux内核,使用Eclipse调试Qemu及Linux Kernel
  11. 在oracle表中增加字段,并调整字段的顺序
  12. Android心电图动画效果,手把手教你打造一个心电图效果View Android自定义View
  13. 新华书店牵手阿里云 要用云计算打造“智慧悦读”
  14. LambdaQueryWrapper构建查询条件、模糊查询、范围查询、排序
  15. linux cd是什么目录,linux cd是什么意思?
  16. 大华平台linux密码,大华DSS平台低权限账户越权直接修改system密码
  17. python数据类型的性能
  18. JavaCV的配置及使用
  19. fetion发送短信介绍(比如天气预报)
  20. 2023,创客永不停歇

热门文章

  1. ES6知识点汇总(全)
  2. 微服务教程--什么是 Nacos
  3. 这部日本「神作」彻底拉低了我入门AI的门槛
  4. Homebrew命令详解
  5. BGP——本地优先级选路+BGP路由水平分割机制(讲解+配置命令)
  6. Hibernate 中出现 xxx表 is not mapped xxx的问题
  7. C#LeetCode刷题之#643-子数组最大平均数 I( Maximum Average Subarray I)
  8. 给你安利一款国产良心软件uTools
  9. javascript语法_JavaScript传播语法简介
  10. JavaScript状态2018