[Mathematics][MIT 18.03] Detailed Explanation of the Frequency Problems in Second-Order Differential Equation of Oscillation System
Well, to begin with, I'd like to say thank you to MIT open courses twice. It's their generosity that gives me the priviledge to enjoy the most outstanding education resources.
Okay, come to the point -- The Frequency in Oscillation System.
In general, there are two cases, undamped and damped. Let's begin with the first and easier one.
Case 1. Undamped Oscillation System.
In this system, the general form of the equation is $y'' + \omega_{0}^{2}y=f(x)$, in which $f(x)$ refers to the input or, more practically, driving or forcing term.
In this case, there is only one frequency: The Natural Frequency of this system $\omega_{0}$, AND if $f(x)$ is a trigonometric function, to maximize the amplitude or give rise to authentic resonance, we just need to let the frequency of $f(x)$ equals to $\omega_{0}$.
Okay, let's do the inferences. First of all, if we want to get the natural frequency, we have to come back to the reduced or associated homogeneous equation(just let $f(x)$ equal to $0$).
So, we get $y'' + \omega_{0}^{2}y=0$. It's easy to solve by set $y = e^{rt}$, in which $t$ is the usually independent variable, specifically time in practical usage, and $r$ is the unknown constant.
This process is quite simple. So I just omit it and present the result. We finally get $y = c_{1}\cos\omega_{0} t + c_{2}\sin\omega_{0} t$. Let's choose a simple particular solution.
$y = \cos\omega_{0} t$
So, coming to this step, it's quite clear to see that this system possess a natural frequency $\omega_{0}$.
But, further, with the driving term $\cos\omega t$, the solution will be the real part of $\displaystyle \tilde{y} = \frac{e^{i\omega t}}{\omega_{0}^{2}-\omega^{2}}$. And, by some analysis, we will get the conclusion that the resonance will occur if $\omega \to \omega_{0}$.
Case 2. Damped Oscillation System.
In the same way, the general form of the equation is $y'' + 2py' + \omega_{0}^{2}y = f(x)$, in which $f(x)$ refers to the input or, more practically, driving or forcing term and $p \ne 0$.
In this case, there exists three frequency: The Natural Frequency of the associated undamped oscillation system $\omega_{0}$, The Damping Frequency, which is a pseudo frequency because actually the solution is not a pure trigonometric function, $\omega_{1}$, which equals to $\sqrt{\omega_{0}^{2}-p^{2}}$, The Frequency which driving force needs to maximize its pseudo amplitude or give rise to pseudo resonance, $\omega_{r}$, which equals to $\sqrt{\omega_{0}^{2}-2p^{2}}$ or $\sqrt{\omega_{1}^{2}-p^{2}}$.
Okay, so let's begin with the simplest one, the natural frequency. Just delect the first-order term and remove the $f(x)$ and according to the case 1, we will get the result. That's simple, so it does not deserve detailed explanation.
Secondly, the damping frequency. Above all, I have to clarify that in order to generate oscillation, this system has to be underdamped, which means that $p^{2}-\omega_{0}^{2}<0$, and we need to go back to the associated homogeneous equation too. So, by using the characteristic function, we get the roots are $-p \pm \sqrt{\omega_{0}^{2}-p^{2}}i$. That's why we set $\omega_{1}=\sqrt{\omega_{0}^{2}-p^{2}}$. As a result, we will finally get one simple particular solution $y = e^{-pt}\cos\omega_{1}t$. And, that's pseudo frequency because despite the behaviors it possesses, like frequency, period, oscillation, its amplitude actually fades with time. So the oscillation actually decays with time.
Thirdly, we come to the most difficult part, which frequency maximizes its pseudo amplitude or gives rise to pseudo resonance. Above all, by same practice, we set $f(x) =\cos\omega t$ and assume that $p(\alpha) \ne 0$. (NOTE: 1. $p(\alpha)$ means we substitude every $D$(the note of differentiation) in the polynomial with $\alpha$, in this specific case, the polynomial is $D^{2}+2pD+\omega_{0}^{2}$ 2. $\alpha$ is the complex index number of the exponential form of driving term after we complexify the $f(x)$.) So we get the complex solution $\displaystyle \tilde{y} = \frac {e^{i\omega t}}{-\omega^{2}+\omega_{0}^{2}+2p\omega i}$.
Above all, we need to simplify the problem.
Step 1. We need to do some fixes on the denominator. $\omega_{0}^{2}+2p\omega i - \omega^{2} = \omega_{1}^{2}+p^{2}+2p\omega i-\omega^{2}$. Let's say $\bar{\omega}^{2} = \omega_{1}^{2}+p^{2}-\omega^{2}$(NOTE: bar here doesn't mean conjugation). So we have $\displaystyle \tilde{y} = \frac {e^{i\omega t}}{\bar{\omega}^{2}+2p\omega i}$.
Step 2. Extract the real part. $\displaystyle \tilde{y} = \frac {e^{i\omega t}}{\bar{\omega}^{2}+2p\omega i}=\frac {e^{i\omega t}(\bar{\omega}^{2}-2p\omega i)}{\bar{\omega}^{4}+4p^{2}\omega^{2}}$, so $\displaystyle y = \frac {\bar{\omega}^{2}}{\bar{\omega}^{4}+4p^{2}\omega^{2}}\cos\omega t + \frac{2p\omega}{\bar{\omega}^{4}+4p^{2}\omega^{2}}\sin\omega t$.
Step 3. So the amplitude is $\displaystyle \sqrt{(\frac {\bar{\omega}^{2}}{\bar{\omega}^{4}+4p^{2}\omega^{2}})^{2} + (\frac{2p\omega}{\bar{\omega}^{4}+4p^{2}\omega^{2}})^{2}}$, which is $\displaystyle \frac{1}{\sqrt{\bar{\omega}^{4} + 4p^{2}\omega^{2}}} = \frac{1}{\sqrt{(\omega_{1}^{2}+p^{2}-\omega^{2})^{2} + 4p^{2}\omega^{2}}}$.
Step 4. In order to get the maximize amplitude, we need to minimize the denominator. And using some senior high school's knowledge, this would be $\omega^{2} = \omega_{1}^{2}-p^{2}$. And that is $\omega_{r}$.
And there remains an assumption to deal with. $p(\alpha) = 0$, is this possible? Let us see the denominator $\omega_{1}^{2}+p^{2}+2p\omega i-\omega^{2}$. If it's zero, then $p = 0$ or $\omega = 0$. Based on the condition, we know that $p \ne 0$ and $\omega = 0$ is meaningless. So we can safely exclude this case.
[Mathematics][MIT 18.03] Detailed Explanation of the Frequency Problems in Second-Order Differential Equation of Oscillation System的更多相关文章
- [Mathematics][MIT 18.03] Proof of a Theory about the Solution to Second-order Linear Homogeneous Differential Equation
At first, I'd like to say thank you to MIT open courses which give me the privilege to enjoy the mos ...
- [Mathematics][MIT 18.02]Detailed discussions about 2-D and 3-D integral and their connections
Since it is just a sort of discussion, I will just give the formula and condition without proving th ...
- PYTHON替代MATLAB在线性代数学习中的应用(使用Python辅助MIT 18.06 Linear Algebra学习)
前言 MATLAB一向是理工科学生的必备神器,但随着中美贸易冲突的一再升级,禁售与禁用的阴云也持续笼罩在高等学院的头顶.也许我们都应当考虑更多的途径,来辅助我们的学习和研究工作. 虽然PYTHON和众 ...
- Docker 18.03 Centos7.6 安装 内网
首先访问https://download.docker.com/linux/centos/7/x86_64/stable/Packages/获取对应版本的rpm包docker包docker-ce-18 ...
- 18/03/18 04:53:44 WARN TaskSchedulerImpl: Initial job has not accepted any resources; check your cluster UI to ensure that workers are registered and have sufficient resources
1:遇到这个问题是在启动bin/spark-shell以后,然后呢,执行spark实现wordcount的例子的时候出现错误了,如: scala> sc.textFile()).reduceBy ...
- windows的docker开始支持linux的镜像 ,Version 18.03.0-ce-win59 (16762)
LCOW containers can now be run next to Windows containers.Use '--platform=linux' in Windows containe ...
- [MIT 18.06 线性代数]Intordution to Vectors向量初体验
目录 1.1. Vectors and Linear Combinations向量和线性组合 REVIEW OF THE KEY IDEAS 1.2 Lengths and Dot Products向 ...
- Docker 18.03导入导出
docker中分容器和镜像,简单可以理解为容器是运行中的实例,镜像是运行实例所需的静态文件. 导入导出既可以对容器做操作,也可以对镜像做操作.区别在于镜像可以随时导出,容器必须要停止之后才可以导出,否 ...
- ApacheCN 编程/大数据/数据科学/人工智能学习资源 2019.12
公告 我们的所有非技术内容和活动,从现在开始会使用 iBooker 这个名字. "开源互助联盟"已终止,我们对此表示抱歉和遗憾.除非特地邀请,我们不再推广他人的任何项目. 公众号自 ...
随机推荐
- eclipse中maven web项目部署时缺少classes文件或者resources文件
写这篇博客的原因 问题描述 昨天发现eclipse中maven web项目部署时缺少classes文件或者resources文件 本来以为是很常见的原因, 依次检查"Java Build P ...
- 十六 StudentManagerSystem的一些业务实现
1 删除学生的JSP实现: <%@ page language="java" contentType="text/html; charset=UTF-8" ...
- Censoring「USACO 2015 Feb」
题目描述 有一个S串和一个T串,长度均小于1,000,000,设当前串为U串,然后从前往后枚举S串一个字符一个字符往U串里添加,若U串后缀为T,则去掉这个后缀继续流程. 输入格式 包含两行,第一行为S ...
- [noip2014]P2312 解方程
P2312 解方程 其实这道题就是求一个1元n次方程在区间[1, m]上的整数解. 我们枚举[1, m]上的所有整数,带进多项式中看看结果是不是0即可. 这里有一个技巧就是秦九韶算法,请读者自行查看学 ...
- JuJu团队1月8号工作汇报
JuJu团队1月8号工作汇报 JuJu Scrum 团队成员 今日工作 剩余任务 困难 飞飞 实现三维Dense 将crossentrophy和softmax连接起来 无 婷婷 完善main.jl ...
- windows驱动不要签名
BCDEDIT -SET LOADOPTIONS DISABLE_INTEGRITY_CHECKSBCDEDIT -SET TESTSIGNING ON
- [转载]@Component 和 @Bean 的区别
@Component 和 @Bean 的区别 @Component 和 @Bean 的区别 Spring帮助我们管理Bean分为两个部分,一个是注册Bean,一个装配Bean. 完成这两个动作有三种方 ...
- 使用PYaudio录制音频和视频(自己)
参考:https://blog.csdn.net/zhaoyun_zzz/article/details/84341801 音频录制:简洁版 import pyaudioimport waveimpo ...
- CentOS 7安装/卸载Redis,配置service服务管理
Redis简介 Redis功能简介 Redis 是一个开源(BSD许可)的,内存中的数据结构存储系统,它可以用作数据库.缓存和消息中间件. 相比于传统的关系型数据库,Redis的存储方式是key-va ...
- PHPstudy2018 后门简单使用
首先声明,仅为记录使用. 测试用例php-5.4.45 + Apache index.php 使用Firefox 浏览器,可以编辑并且重发http请求 打印出“net user” base64 加密后 ...