matlab中变.,Matlab变系数微分方程

附件给出的是齿轮动力学振动模型和微分方程

其中x为要求的位移响应,Zeta为齿轮副相对阻尼系数,Zeta =0.07,m为等效质量m=3;Fn为法向载荷Fn=100N。

k(t)为齿轮副的时变啮合刚度

对k(t)我做了关于时间t的傅里叶变换:

k=750527161.016+1.2466e+008*sin(2*pi*f0*t-0.86045)+8.3314e+007*sin(4*pi*f0*t-0.2292)+4.9789e+007*sin(6*pi*f0*t+0.41062)+1.7784e+007*sin(8*pi*f0*t+1.0182)+7.633e+006*sin(10*pi*f0*t-1.2305);

这个方程是变系数微分方程,其中就k(t)一个时变系数,k(t)该如何表示?是否可把k(t)作为t的函数直接代入方程?

下面是我编的M文件和求解程序,各位大神帮我看看处理的对不对?谢谢

M文件

function Yd= fun7( t,x,k,zeta,m,Fn)

% FUN7.M: 时变微分方程例子

Fn=100;

zeta=0.07;

m=3;

n1=3000;

z1=45;

Tz=0.00044444;

f0=1/Tz;

k=750527161.016+1.2466e+008*sin(2*pi*f0*t-0.86045)+8.3314e+007*sin(4*pi*f0*t-0.2292)+4.9789e+007*sin(6*pi*f0*t+0.41062)+1.7784e+007*sin(8*pi*f0*t+1.0182)+7.633e+006*sin(10*pi*f0*t-1.2305);

dx=zeros(2,1);

dx(1)=x(2);

Yd=[x(2);Fn/m-2*zeta*sqrt(k/m)*x(2)-k/m*x(1)];

求解程序:

clear

m=3;

Tz=0.00044444;

f0=1/Tz;

tspan=linspace(0,0.00044444,120);

Y0=[1.5040e-007;0];

[t,XX]=ode45('fun7',tspan,Y0);

X=XX(:,1);

Y=XX(:,2);

plot(t,X),xlabel('经历时间/(s)'),ylabel('位移x'),figure;

plot(t,Y),xlabel('经历时间/(s)'),ylabel('速度v'),figure;

plot(X,Y)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/566931.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

ztree 标准得json数据格式_酷站推荐 - json-c.github.io/json-c - json-c API

json-c.github.io/json-c/http://json-c.github.io/json-c/JSON:JavaScript 对象表示法(JavaScript Object Notation)。 JSON 是存储和交换文本信息的语法。类似 XML。 JSON 比 XML 更小、更快,更易解析。JSON(JavaScript Object …

matlab randn state 2,MATLAB?中的randn函数

matlab函数randn:产生正态分布的随机数或矩阵的函数randn产生均值为0,方差 σ^2 1,标准差σ 1的正态分布的随机数或矩阵的函数。用法:Y randn(n)返回一个n*n的随机项的矩阵。如果n不是个数量,将返回错误信息。Y ra…

qt弹框输入密码_Android仿支付宝密码输入框

实现效果图:实现流程:1、定义6位密码输入View思路:要绘制边框矩形,绘制分割线,绘制圆点。绘制圆的数目要与字符串的长度有关,添加或者删除都要修改字符串,输入6位后就是要关闭弹框,拿…

php 前置匹配,浅析PHP正则表达式匹配的特定实现

PHP正则表达式匹配的特定实现是什么呢?其实我们知道,在实际的匹配操作的过程中我们操作的不仅仅是单个的字母或是数字,那么我们面对像单词的或是一组数字的时候该如何处理呢?PHP正则表达式匹配的特定实现会使用到PHP正则表达式内置…

esc指令检查打印状态_Z.115 胶片自助打印设备

一、概述病人在完成影像检查后,需要快速获取结果,传统的方式是向病人发放胶片和诊断报告,随着大型医院病人检查量的快速增加,这种方式有诸多弊病。自助打印机将胶片和诊断报告打印集成起来,消除时空障碍,实…

java执行class找不到main函数_你所不知道的HelloWorld背后的执行原理

专注于Java领域优质技术,欢迎关注作者:饭谈编程【今日最佳】对于程序员而言,所谓的二八定律指的是 花百分之八十的时间去学习日常研发中不常见的那百分之二十的原理。据说阿里某程序员对书法十分感兴趣,退休后决定在这方面有所建树…

matlab窗函数 响应,matlab窗函数设计方案.doc

matlab窗函数设计方案PAGE \* MERGEFORMATPAGE \* MERGEFORMAT 13目录TOC \o "1-3" \h \u HYPERLINK \l _Toc18947 目录 PAGEREF _Toc18947 1HYPERLINK \l _Toc5619 摘要 PAGEREF _Toc5619 2HYPERLINK \l _Toc16857 1 MATLAB的概况 PAGEREF _Toc16857 2HYPERLINK \l …

java接口文档生成工具_【分享】接口文档生成工具apipost

一、为什么要写接口文档?正规的团队合作或者是项目对接,接口文档是非常重要的,一般接口文档都是通过开发人员写的。一个工整的文档显得是非重要。项目开发过程中前后端工程师有一个统一的文件进行沟通交流开发,项目维护中或者项目…

matlab boll源码,【每日一策】Matlab量化交易策略之 布林带突破+头寸管理

function Strategyn(freq)%targetList traderGetTargetList();%获取目标资产信息HandleList traderGetHandleList();%获取账户句柄global record;global state;for k1:length(targetList);%--------------------仓位、K线、当前bar的提取-----------------------------%%获取…

jacobi matlab程序,jacobi迭代法实验MATLAB程序数值分析

jacobi迭代法实验MATLAB程序数值分析 例1. 求线性方程组 得近似解。精确解为x*[3,2,1]’。 解:对方程进行移项就得 记为Axb,或写为xB0 xf,其中 取初始值,代入原方程组可得再将把它代入可得.反复利用这个计算过程,得到一向量序列和…

docker部署python web应用_安装docker并部署web项目

一.docker简介1、docker定义:docker是一个用来装应用的容器,就像杯子可以装水,笔筒可以装笔,书包可以放书一样。你可以把“hello world!”放到docker中,也可以把网站放到docker中,你可以把任何你…

java 水印乱码,linux java程序加水印及中文乱码方案

在linux,centos环境下的,生成的带文字的水印图片在显示为方框乱码。img.setFont(new Font("宋体", Font.BOLD, 20));写中文进入图片之前必须设置字体,而且这个字体必须支持中文,否则就会出现乱码或者方框、问号等等。Bu…

resultmap拿不到数据_阿里巴巴国际站每日电商运营工作数据表格

一、日常运营工作表1.数据日报统计每天的流量数据:2.数据周报统计每周的流量数据:3.数据月报统计每月流量数据:前半部分后半部分4.P4P月数据总览统计P4P的流量数据:大图包含上面所有的方案5.P4P日数据总览6.爆款数据统计每月统计爆…

php oracle新增数据类型,oracle数据类型

数据类型字符串类型字符串类型用来存储字符串数据char用于描述定长的字符型数据0~2000nchar用来存储Unicode字符集的定长字符型数据0~1000varchar2用于描述可变长的字符型数据0~4000nvarchar2用来存储unicode字符集的可变长字符型数据0~1000long用来存储边长的字符串0~2GBvarch…

jenkins ssh 远程部署_Jenkins部署jar到远程服务器

首选得确定Jenkins的安装完整,还有插件的安装,除默认插件,此处必备插件:Publish Over SSH,SSH plugin,因为我源码在gitlab所以关于gitlab的插件我也安装了,看各自的情况选择安装。然后到系统管理…

jieba分词_Jieba.el – 在Emacs中使用jieba中文分词

jieba.el在Emacs中使用jieba中文分词众所周知, Emacs并没有内置中文分词系统, 以至于 forward-word 和 backward-word 以及 kill-word 等以单词为单位的操作只能粗暴的标点符号进行确定所谓”词汇”, (其实是中文里的句子).jieba.el 利用nodejieba对buffer中的中文句子进行分割…

python的request请求401_Python模拟HTTPS请求返回HTTP 401 unauthorized错误

Python模拟HTTPS请求返回HTTP 401 unauthorized错误开始是使用的 httplib模块,代码如下:header {"Content-type": "application/json", "Accept": "*/*" }params { ‘source‘:‘en‘, ‘target‘:‘es‘, ‘…

java php 女生数量,萌妹子告诉你php和java如何选

原标题:萌妹子告诉你php和java如何选phpPHP 独特的语法混合了C、Java、Perl 以及 PHP 自创新的语法。它可以比 CGI 或者 Perl 更快速的执行动态网页PHP具有非常强大的功能,所有的CGI的功能PHP都能实现,而且支持几乎所有流行的数据库以及操作系…

new 结构体指针_Go:我应该用指针替代结构体的副本吗?

logo对于许多 golang 开发者来说,考虑到性能,最佳实践是系统地使用指针而非结构体副本。我们将回顾两个用例,来理解使用指针而非结构体副本的影响。1. 数据分配密集型让我们举一个简单的例子,说明何时要为使用值而共享结构体&…

Oracle 同音字查询,汉字的演变过程100字,汉字的演变图片

中国的汉字博大精深,历史悠久。早期甲骨文主要是线条图,统一而优美。后来,它逐渐演变成一幅由笔画组成的图画,这是当代汉字的基础商周时期,甲骨文逐渐成为金文。西周时期,青铜器铸造,铭文用于记…