课程笔记
文章平均质量分 81
忧郁奔向冷的天
力学,信号处理,数值计算,机器学习,编程。
I come from China and I'm interested in mathematics and physics.
I love programming.
I love to listen to and sing Cantonese songs.
I love keeping fit.
展开
-
张量学习(一)
零阶张量,标量:密度,温度,压力一阶张量,矢量:速度,力二阶张量:物理量有大小和两个方向。应力,应变原创 2022-02-23 20:53:01 · 2776 阅读 · 0 评论 -
力学—质点运动学
矢量微积分切向量 t^=r⃗˙/∣r⃗∣\hat{t} = \dot{\vec{r}}/|\vec{r}|t^=r˙/∣r∣点积可以交换位置,叉积不可以交换位置A⃗⋅A⃗=A2\vec{A}\cdot \vec{A} = A^2A⋅A=A2常矢量求导为0,如直角坐标系三个常矢量一元矢量函数求导A⃗(t)=Ax(t)x^+Ay(t)y^+Az(t)z^\vec{A}(t) = A_x(t) \hat{x} + A_y(t) \hat{y} + A_z(t) \hat{z}A(t)=Ax(t)x^原创 2022-02-16 09:44:14 · 407 阅读 · 0 评论 -
Condition Number
矩阵条件数条件数cond(A)=∣∣A−1∣∣∣∣A∣∣cond(A)=||A^{-1}|| ||A||cond(A)=∣∣A−1∣∣∣∣A∣∣可作为矩阵求逆的病态程度的度量,条件数越大,则认为病态越严重。条件数与所取的范数有关,常用的有:谱条件数 cond2(A)=∣∣A−1∣∣2∣∣A∣∣2cond_2(A) = ||A^{-1}||_2||A||_2cond2(A)=∣∣A−1∣∣2∣∣A∣∣2F条件数 condF(A)=∣∣A−1∣∣F∣∣A∣∣Fcond_F(A) = ||A^{-1}原创 2021-09-30 14:01:44 · 252 阅读 · 0 评论 -
施密特正交化
https://zh.wikipedia.org/wiki/格拉姆-施密特正交化施密特正交化在线性代数中,如果内积空间上的一组向量能够组成一个子空间,那么这一组向量(linearly independent)就称为这个子空间的一个基。Gram-Schmidt正交化提供了一种方法,能够通过这一子空间上的一个基得出子空间的一个正交基,并可进一步求出对应的标准正交基。以施密特命名,其实Laplce与Cauchy比施密特更早发现这一方法。向量u在向量v上的投影 projvu=<u,v><原创 2021-09-23 11:27:34 · 3509 阅读 · 0 评论 -
杜哈梅积分
振动力学 倪震华https://zh.wikipedia.org/wiki/杜哈梅积分二阶常系数齐次微分方程的通解是如何推导出来的? https://wenku.baidu.com/view/d8b279e0b84ae45c3a358c30.html贝塞尔方程 Bessel Differential Equation1797年-1872年Jean-Marie Constant Duhamel 是一名法国数学家与物理学家。单自由度欠阻尼系统振动力学方程 mx′′+cx′+kx=p(t)mx'原创 2021-09-16 23:07:19 · 5750 阅读 · 1 评论 -
涡量与涡旋
流体力学中为什么要引入涡量? - 更高的天空的回答 - 知乎https://www.zhihu.com/question/31159018/answer/377602533流体力学中为什么要引入涡量? - 朱辉的回答 - 知乎https://www.zhihu.com/question/31159018/answer/51047540液滴击打水面时产生的像卡门涡街的涡量分布 http://blog.sciencenet.cn/blog-739225-889268.html为什么说刚体线速度的旋度等于角原创 2021-09-13 12:16:51 · 6771 阅读 · 0 评论 -
实验设计方法绪论
电子科技大学 实验设计方法绪论实验:去做实验验证某个已知结论;试验:去做实验探索某个未知结论。实验设计方法:把数学上优化理论、技术应用于试验设计中,科学的安排试验、处理试验结果的方法。或者是,采用科学的方法去安排试验,处理试验结果,以最少的人力和物力消费,在最短的时间内取得更多、更好的生产和科研成果的最有效的技术方法。田口正交设计法,Taguchi 田口玄一。华罗庚优选法,华罗庚是中国的实验设计方法之父。方开泰均匀设计法。绪论实例为提高某化工产品的转化率,选取三个有关的因素(A反应温度、B反原创 2021-08-21 11:28:58 · 616 阅读 · 0 评论 -
进化计算、群体优化、搜索求解的笔记
进化计算进化计算是以达尔文的进化论、孟德尔的遗传学作为基础的智能方法。现实中生物进化包括:繁殖、基因重组、基因变异、个体竞争、自然选择,进化计算主要包括:选择、重组、变异。进化计算方法:遗传算法、协同进化、进化策略、进化规划、遗传编程遗传算法一种随机搜索方法,不需要数学的求导交叉:概率选择交叉点变异:概率选择变异点A*算法遗传算法、模拟退火算法粒子群算法、蚁群算法...原创 2021-08-09 18:05:31 · 404 阅读 · 0 评论 -
主动隔振与被动隔振
主动隔振主动隔振:将作为振源的机器与地基隔离。用单自由度系统举例,隔振材料位k与c隔振前机器传到地基的力F0eiwtF_0e^{iwt}F0eiwt隔振后机器传到地基的响应x=F0kβei(wt−θ1)x=\frac{F_0}{k}\beta e^{i(wt-\theta_1)}x=kF0βei(wt−θ1),其中F0k\frac{F_0}{k}kF0是静变形,β=1(1−s2)2+(2ζs)2\beta=\frac{1}{\sqrt{(1-s^2)^2+(2\zeta s)^2}}β原创 2021-08-05 22:41:42 · 4593 阅读 · 0 评论 -
量纲分析
针对某一问题,要确定R与或者是等式右边左左除第一种形式 确定出a,b,c...后,第二种形式 确定出a,,b,c...后,Pi定理:设影响某现象的物理量数为n个,这些物理量的基本量纲为m个,则该物理现象可用n-m个独立的无量纲数群(准数)关系式表示量纲法则:1.只有量纲相同的物理量,才能彼此相加、相减和相等; 2.指数函数、对数函数和三角函数的宗量应当是量纲1的。零的量纲解析 一词怎么解释?...原创 2021-03-08 15:06:32 · 2948 阅读 · 0 评论 -
变分法
变分法弦平衡方程的导出,建立起横向位移u,张力T,外力f之间的关系:方一、根据受力平衡导出推导时用的技巧或假设:1.泰勒展开近似 同理2.3.小变形假设,张力均匀,即4.方程推导中忽略二阶小量,即只要某一项带就忽略掉他。方二、根据总势能最小导出u是x的函数,J是u的函数,求J的min,即求极值,求泛函的极值即变分问题v为试探函数,此时将问题从对泛函求极值转化为对函数求极值,即J(u(x))求极值转化为J(λ),这种方法为变分法推导时用到的技巧:1、2、泰勒展...原创 2021-03-08 00:51:09 · 3309 阅读 · 0 评论 -
《结构分析的有限元法与MATLAB程序设计》笔记
学习资源目前求解PDE的数值方法:FEM(有限单元法)、DEM(离散元)、FVM(有限体积法)、X-FEM(扩展有限元)、FDM(有限差分法)、LBM(玻尔兹曼格子法)、SPH(光滑粒子流体动力学)朱伯芳《有限单元法原理与应用》王勖成《有限单元法》Logan《A first course in finite element method》Jacob Fish和Ted Belytschko《A first course in finite elements》Bathe《Finite Eleme原创 2021-02-21 10:56:57 · 6961 阅读 · 2 评论 -
材料力学——弯曲内力、弯曲应力
弯曲内力弯曲内力:受弯杆件横截面上的内力由给定的梁与外力,列出梁的剪力方程、弯矩方程,绘制剪力图、弯矩图剪力方程:Fs = Fs(x)弯矩方程:M = M(x)载荷集度q(x)弯曲应力纯弯曲及弯曲正应力剪切弯曲:梁上既有弯矩也有剪力,此时横截面有正应力、切应力纯弯曲:梁上只有弯矩,此时横截面只有正应力讨论几何方程、本构方程、平衡方程。几何方程:纯弯曲后,中性层的曲率为ρ,与中性层距离y的纵向纤维的应变为ε本构方程:平衡方程:内力矩与外力矩M平衡,内力矩由正应力积分原创 2021-02-21 10:29:47 · 5216 阅读 · 0 评论 -
特征值、特征根、本征值
线性代数狭义特征值问题 Ax = λx 广义特征值问题Ax = λBxλ为特征值,x为λ对应的特征向量在求解特征值时|A-λE|=0 转化为求解特征多项式的特征根高数特征根的概念在求解二阶常微分方程中也有二阶微分方程x'' + px' + qx = 0 对应的特征方程 r^2 + p*r + q = 0控制方程在控制方程中也有特征根二阶微分方程x'' + px' + qx = 0 经过拉氏变换 得到特征方程 s^2 + p*s+ q = 0特征方程就是传递函数...原创 2020-09-20 08:34:39 · 14400 阅读 · 0 评论 -
数学物理方程 第五章 勒让德多项式
Legendre多项式分离变量法中使用的本征函数是三角函数系,本章使用的本征函数是勒让德多项式构成的本征函数系。1. 对比Sturn-Liouville本征值问题分离变量法一章讨论的是 X''+λX=0 根据不同的边值条件得到不同的本征值与本征函数使得解X不为0的λ称为本征值求解关于u=u(x,t)的齐次微分方程utt = a^2*uxx时,使用了分离变量法u=X(x)*T(t)勒让德方程使得解y不为0的λ称为本征值求解微分方程时,设,即直接设解的形式为幂级数可看出以上两个..原创 2020-08-24 22:45:33 · 18935 阅读 · 0 评论 -
数学物理方程—积分复习
三维直角坐标系 dxdydz柱坐标系 drdθdz球坐标系 drdθdφ令x=rcosθ y=rsinθdxdydz =令x=rsinφcosθ y=rsinφsinθ z=rcosφdxdydz =原创 2020-08-24 22:44:33 · 5140 阅读 · 0 评论 -
数学陈述
公理:没有经过证明,但被当作不证自明的命题公设:本学科内的,没有经过证明,但被当作不证自明的命题公理是许多科学分支所共有的,而各个科学分支中的公设则是不同的。定理:在公理的基础上推导出来,并被证明为真的命题命题:可以判断真假的陈述句,有真命题和假命题,也有悖论推论:一个从定理随之而即时出现的命题引理:某个定理的证明的一部分的命题假说:根据已知的科学事实和科学原理,对所研究的自然现象及其规律性提出的推测和说明定律:定义:定则:法则:原理:有时间会补充...原创 2020-08-12 23:45:04 · 348 阅读 · 0 评论 -
《现代控制理论》绪论
西北工业大学课程1.控制理论的发展经典控制理论/自动控制理论现代控制理论鲁棒控制理论2.经典控制理论的发展指南车:开环自动调节系统水运仪象台:闭环自动调节系统 开环、闭环的区别就是有没有反馈十八世纪,自动控制理论用到工业当中十九世纪—二十世纪,发展判断自动控制系统是否稳定的理论马克斯维尔提出稳定性判据劳斯与赫尔维兹(胡尔维兹)奈奎斯特提出频域的频率响应法伊万斯提出复数域的根轨迹法维纳 1948《控制论—关于动物和机器中控制与通讯的科学》建立控制理论一般方法...原创 2020-08-04 09:25:50 · 1590 阅读 · 0 评论 -
控制工程——传递函数与状态方程
1.基本概念状态方程:信号理论、控制理论与 热学都有状态方程,我现在指的是前者的。状态:控制系统中 可决定系统状态的最小数目的一组变量组成的集合状态向量:状态变量在某一时刻的值。在t=0时刻的的状态向量也称初始状态向量。如该系统的状态变量为{变量1 变量2 变量3}那么 t=t0时刻的 状态向量为{值1 值2 值3}状态空间:状态={变量1,变量2,...} 状态空间={状态1,状态2,...} 状态空间就是一个集合的集合现代控制理论以状态空间理论为基础,经典控制理论以传递函数理论为基...原创 2020-08-02 19:11:27 · 14771 阅读 · 0 评论 -
多自由度系统矩阵形式的动力学方程
1. 二自由度圆盘自由扭振原创 2020-07-17 22:31:31 · 2222 阅读 · 0 评论 -
线性振动近似计算方法(三)
用传递矩阵法求解梁的横向振动的固有频率1. 对支座、质量点、梁段进行编号支座也算在质量点中,共4个质量点共3个梁段2.边界条件即第0个质量点(支座) 右端的挠度与弯矩为0即第3个质量点(支座) 左端的挠度与弯矩为03.无量纲形式X上面加bar,这里表示这是一个无量纲状态变量X = (y θ M Fs)_X = (_y θ _M _Fs)_y = y/L_M = ML/(EI)_Fs = Fs*L**2/(EI)将边界条件也转化为无量纲的将点传递...原创 2020-07-17 19:24:51 · 302 阅读 · 0 评论 -
线性振动近似计算方法(二)
圆盘扭转传递矩阵法三圆盘系统的自由扭振1.边界条件边界条件:即第1个圆盘左侧的扭矩与第3个圆盘右侧的扭矩为02.每个圆盘的状态第i个圆盘左端与右端的传递关系为:点传递矩阵第i个轴左端与右端的传递关系为:场传递矩阵第i个圆盘右端与第i+1个圆盘右端的传递关系为:单元传递矩阵3.模态频率由边界条件得第3个圆盘右端的扭矩为0得从而求出前三阶固有频率w1 w2 w34.模态振型将w1依次带入 第1个圆盘右侧的转角,第2个圆盘右侧的转角,第3...原创 2020-07-16 22:38:40 · 454 阅读 · 0 评论 -
线性振动近似计算方法(一)
1.dunkerley method计算系统的基频λ1 ≈ _λ1+_λ2 = f11m1 + f22m2f11:给m1施加单位力时,m1所产生的位移原创 2020-07-16 14:50:38 · 1156 阅读 · 0 评论 -
目前掌握的来自Gabriel Lamé的一些知识
加布里埃尔·拉梅(Gabriel Lamé,1795年7月22日-1870年5月1日),法国数学家。1. Lamé公式1883年Lamé首次对厚壁圆筒进行应力分析,提出的应力计算公式叫做Lamé公式。这是在压力容器设计中提到的。2. Lamé常数σ_ij = λ*e_aa*δ_ij+ 2μ*e_ijλ μ就是Lamé常数,第二个常数通常写作G,称之为剪切模量。δ_ij是克罗内克函数e_aa = e_xx +e_yy+e_zz 称为单位体积的体积改变量;体应变;wiki上面称 体积变...原创 2020-07-10 21:03:50 · 1475 阅读 · 1 评论 -
克罗内克δ与利威尔-奇维塔符号
1. 克罗内克δi = j时 δij = 1i != j时 δij = 02. 利威尔-奇维塔符号也称置换符号二维当ijk不是12的顺序时,ε ij = -1当ijk重复时,ε ij = 0当ijk是12的顺序时,ε ij = 1三维当ijk不是123的顺序时,ε ijk = -1当ijk重复时,ε ijk = 0当ijk是123的顺序时,ε ijk = 1则 ε 123 = 1,ε 312 = 1, ε 213 = -1, ε 112 = 0疑问:为什么前者称为原创 2020-07-07 00:33:40 · 5217 阅读 · 0 评论 -
惯性力——举例说明
1. 旋转的圆盘上,放置一个物体,物体在圆盘上静止以地面为参考系,在地面建立坐标系此时物体受重力,摩擦力,支撑力。这三种力的合力使得物体运动。以圆盘为参考系,在圆盘上建立坐标系此时物体受重力,摩擦力,支撑力,惯性力。这四种力使得物体静止。这里指的惯性力包括什么?如果圆盘非匀速转动,那么惯性力包括离心力与欧拉力如果圆盘匀速转动,那么惯性力就是离心力2. 你静止站在行驶的公交车上以地面为参考系此时你受摩擦力,重力,支撑力。这三种力使你运动。以汽车为参考系此时你受摩擦力,重力,支原创 2020-07-05 13:53:48 · 6805 阅读 · 0 评论 -
材料力学与弹性力学中讲到的 扭转
材料力学圆轴扭转、非圆截面扭转、开口薄壁杆件扭转、闭口薄壁杆件扭转圆轴的塑性扭转1.非圆截面扭转矩形梁进行扭转,边缘处有最大切应力。在边缘上,最大切应力位于长边的中点处。短边处的最大切应力τ1与长边处的最大切应力τ2都位于他们的中点处,但后者大于前者。当h/b > 10为狭长矩形,τ max = T/(h*δ**2/3),φ = TL/(G*h*δ**3/3),实际上可视作长边上的切应力相等,都等于τ1,短边上的切应力相等,都等于τ2,只有在四个点附近才迅速减小为02.开口薄壁原创 2020-07-02 14:16:03 · 7703 阅读 · 0 评论 -
C0型单元刚度矩阵的性质
1. 刚度矩阵的性质注:Ke表示单元刚度矩阵,K表示整体刚度矩阵(1) Ke中的元素Ke_ij表示:单元的第j个节点产生单位位移,单元的其余节点位移为0,则需要在第i个节点施加的力为Ke_ij其中,对角线上的元素Ke_ii > 0(2) Ke是对称的,即Ke_ij = Ke_ji(3) Ke是半正定的。若单元无刚体位移,Ke正定。若单元有刚体位移,det(Ke) = 0 ,即Ke是奇异的(4) 刚度矩阵的任一行或列表示一个平衡力系,当节点位移全部为线位移时(C0单元),任一行或列的代原创 2020-06-30 14:15:50 · 9957 阅读 · 1 评论 -
加百利号角/托里拆利小号
根据宗教传说,天使长加百利吹号角以宣布审判日的到来。下图就是加百利号角这个号角可看作曲线y=1/x,x>=1的部分绕x轴旋转而成,求这个号角的内表面积以及体积?使用微积分求解该问题这是托里拆利发现的,那会微积分还没发明,所以他使用了祖暅(gèng)原理来得出这一结论:小号的体积为有限,面积却无限。这就是一个佯谬(悖论 paradox),虽然结果如此,但实际中这种小号是不存在的,我们无法制作出这样的小号,因为它长为无穷,所以这样的小号永远都处于被制作中,也就是现实中永远都不存在.原创 2020-06-27 22:44:09 · 3229 阅读 · 0 评论 -
理论力学中的 动量定理、动量矩定理、动能定理
1. 质点的动量定理(mv)' = F质点动量对时间的变化率等于质点所受的力(mv)2 - (mv)1 = I质点动量的增量等于质点所受的力的冲量2.1 质点系的动量定理每个质点的动量的变化率 等于 所受的内力与外力之和 (mivi)' = Fi+Fi* 式1 注:Fi*表示第i个质点受到的内力质点系 就是对式1求和,质点系的内力为0 则Σ(mivi)' = ΣFi 式2质点系的动量定理:质点系所受外力的矢量和 等于 质点系动量的变量率将式2积分,得 [Σ(miv...原创 2020-06-20 20:33:53 · 8178 阅读 · 0 评论 -
二阶偏微分方程
一般形式:注:A B C D E F G均为x y的函数二阶偏微分方程对应的特征方程:对二阶PDE进行分类:B**2-AC>0 如果是在x0,y0处满足,则该点处 方程为双曲型。若任意点都满足,则该方程为双曲型。B**2-AC=0 同理,抛物型B**2-AC<0 同理,椭圆型===========================================双曲型:波动方程抛物型:热传导方程椭圆型:位势方程一维波动:utt=a**2*uxx+f(x,t)...原创 2020-06-17 09:42:03 · 6593 阅读 · 0 评论 -
物理学史上的三次扭秤实验
Torsion balance experiment1.库伦扭秤实验证明带同种电荷的小球之间的排斥力与两个小球中心之间的距离的平方成反比早在1785年,约翰·罗比逊,卡文迪许就已经做实验发现了这个现象,只不过前者未在意,后者未公布。2.卡文迪许扭秤实验1797年8月5日至1798年5月完成,测出了地球引力常数与地球质量3.厄特沃什扭秤实验证明物体的惯性质量与引力质量等价。...原创 2020-06-10 20:56:46 · 3863 阅读 · 0 评论 -
数学物理方程课程感想
数理学院的课程真的离谱,我本以为数学物理方程会侧重讲物理意义,没想到只是推导微分方程。早知道我就不报了,还不如报个概率论也好过这个。联想到上学期上的矩阵论。如果面对工科生更应该讲物理含义,而不是如何解题,照着书一步一步来解。...原创 2020-05-16 12:04:06 · 1367 阅读 · 0 评论 -
数学物理方程 第三章 行波法
3.1 一维波动3.1.1 无界弦的自由振动----------------------定解问题:范定方程+条件先利用公式求解范定方程,得到u的通解再带入条件,得到u的定解----------------------------定解问题求出的定解该定解称为无限长弦的自由振动的达朗贝尔公式3.1.2 齐次化原理/冲量原理/外力化初速度原理 用于求解非齐次式3.1.3 无界弦的受迫振动 利用杜哈梅原理求解将该定解问题分解为...原创 2020-05-12 23:45:37 · 10327 阅读 · 0 评论 -
数学物理方程 第二章 线性偏微分方程的通解
n阶线性偏微分方程的求解对于我目前的问题来说没什么用,我只是解二阶线性偏微分方程。下面我只讲二阶方程如何解2.1 一维波动方程的解2.2常系数线性齐次偏微分方程的通解 Lu=0 (当然2.2.1简单的讨论了非齐次的情形)2.2.1 L(Dx,Dy)为Dx与Dy的齐次式 只考二阶代数特征方程的解都不相同 α=a1,a2 a1 != a2那么 u=Φ1*(y+a1*x) + Φ2*(y+a2*x)代数特征方程的解都相同 α=a1,a2 a1=a2=a...原创 2020-05-12 23:37:36 · 3900 阅读 · 0 评论 -
数学物理方程 第一章 数学物理方程及其定解问题
微分方程包括:常微分方程ode与偏微分方程pde偏微分方程阶数:方程中存在的最高阶数线性偏微分方程:方程中未知函数u及其导数都是线性的齐次偏微分方程:齐次+线性偏微分方程,齐次是指f(x,t)=0拟线性偏微分方程/半线性偏微分方程:最高阶项是线性的非线性pde完全非线性偏微分方程:最高阶是非线性的范定方程:描写一类物理现象的,不含定解条件的pde定解条件包括初值条件,边值条件,衔接条件范定方程+定解条件来构成定解问题1.1 波动方程及其定解问题1.1...原创 2020-05-12 23:33:06 · 3917 阅读 · 0 评论