[数学][软件] FOUR_TRAN Example 4: Stock Periodicity and Volatility Analysis [FOUR_TRAN傅里叶分析工具 案例4: 股票周期性和波动性分析]

实干、实践、积累、思考、创新。 2020 年写的题目,现在已经是2022年了…… 这两天补上笔记                       案例 ( Examples )   [01] [数学][软件] FOUR_TRAN Example 1: Filtering [FOUR_TRAN傅里叶分析工具使用案例1: 滤波]   [02] …

[几何][数学] 梁跨中下挠引起的端部水平位移有多大?(How large is the end horizontal displacement caused by the midspan deflection of the beam?)

实干、实践、积累、思考、创新。 不考虑梁本身的伸长、弯曲等自身变形,单纯考虑水平放置的梁跨中发生挠度 的情况下,水平伸缩量Δ的大小。这其实是一个纯数学的推导,已经和物理无关了(PS. 题目说梁其实不对),其实说的是一根不可伸缩的直线,跨中发生挠度 ,变为两根直线后,水平伸缩量是多大。 其中 \(L\)为直线的总长度,\(\Delta \)为水平伸缩量的大小,\(X\)为直线发生倾斜变形后,水平投影长度的一半。 由总长不变,可得水平伸缩量 \(\Delta = L – 2X\) 同时 \(L\)、\(\Delta \) 及\(X\) 之间满足以下三角函数关系 \(2X = 2\sqrt {{{\left( {\frac{L}{2}} \right)}^2} – {\omega ^2}} \) 将上式代入第一个公式,可得 …

[数学][几何] 数学之心函数曲线动图

实干、实践、积累、思考,创新。 网上看到的数学之心曲线,函数是: f(x)=x2/3 + 0.9*(3.3-x*x)1/2*sin(b*π*x) 主要参数为 x,b,其中参数x2<3.3,在整个x范围内,改变b可观察到心的形状。 可以通过EXCEL绘制几组b 值下f(x)的函数曲线。 b=1 b=3 b=5 b=8 b=10 b=20 b=40 b=60 可以发现b越大,越接近心形。 最后我们可以把曲线导入之前的编写的 GAAGM软件 [软件][工具][编程] GAAGM: Graph Animation and Animated Gif Maker [曲线图动画及Gif动图生成工具],制作漂亮的曲线动画。 相关博文 ( …

[笔记] 用Matlab进行功率谱分析为什么会出现负的功率谱密度?

实干、实践、积累、思考,创新。 最近闲暇的时候翻翻 随机振动的书,又搞到傅里叶分析,功率谱等。 在测试网上的MATLAB功率谱分析程序的时候,发现功率谱居然是负数,(又暴露了知识面的短缺, 😥 ) 后仔思考发现,主要是因为功率谱使用了 分贝 为计量单位,比如通过公式 10*log10(Pxx) 将原本的功率谱Pxx进行转换,可见当Pxx处于0~1之间的时候,通过log10(Pxx)转换出来的公式就是负数了。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[数学][几何] 费马螺线 [Fermat’s Spiral]

实干、实践、积累、思考,创新。 费马螺线 Fermat’s spiral的极坐标公式如下: 其中其中 a 为实数,θ为转角。 对于任意给定的正值θ,可以有两个r的取值,分别按以下公式: 我们可以通过以下公式将极坐标转换为笛卡尔坐标: x = r*cos(θ); y =r*sin(θ); 接下来可以用过EXCEL绘制费马螺线,如下 最后我们可以把曲线导入之前的编写的 GAAGM软件 [软件][工具][编程] GAAGM: Graph Animation and Animated Gif Maker [曲线图动画及Gif动图生成工具],制作漂亮的曲线动画。   相关博文 ( Related Topics) …

[数学][几何] 阿基米德螺旋线 (等速螺旋线) [Archimedes’ spiral]

实干、实践、积累、思考、创新。 阿基米德螺线(亦称等速螺线)(Archimedean spiral),得名于公元前三世纪希腊数学家阿基米德。阿基米德螺线是一个点匀速离开一个固定点的同时又以固定的角速度绕该固定点转动而产生的轨迹。更多的介绍可以参见:阿基米德螺旋线百度百科词条。 简单的说阿基米德螺线可以用极坐标表示为: r = a+b*θ; 其中其中 a 和 b 均为实数,θ为转角。当θ=0时,a为起点到极坐标原点的距离,b为螺旋线每增加单位角度r随之对应增加的数值。 我们可以通过以下公式将极坐标转换为笛卡尔坐标: x = (a+b*θ)*cos(θ); y = (a+b*θ)*sin(θ); 于是就可以非常简单画出图形了,废话不多说,我们用EXCEL来作图。 以下是b=0情况下, a=0及a=1000 时的螺旋线: 以下是a=0情况下, b=10及b=15时的螺旋线:   最后我们可以把曲线导入之前的编写的 GAAGM软件 [软件][工具][编程] GAAGM: Graph …

[数学][软件] FOUR_TRAN Example 3: Earthquake Ground Acceleration Frequency Spectrum Analysis [FOUR_TRAN傅里叶分析工具 案例3: 地震波频谱分析]

实干、实践、积累、思考、创新。 2020 年写的题目,现在已经是2022年了…… 这两天补上笔记                       关于软件( About the Program)   [01] [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具] 案例 ( Examples …

[数学][编程][混沌] 虫口模型的数值模拟 [Numerical simulation of insect population model]

实干、实践、积累、思考,创新。 1976 年 , 美国数学家梅 (May .R ) 在美国 《 自然 》 杂志上发表的题为 “ 具有复杂的动力学的简单模型 ” 文章中指出 , 在生态学中一些非常简单的确定性的数学模型却能产生看似随机的行为。如虫口模型:xn+ 1= μxn ( 1 -x n ),其中 xn 是第n 年的虫口数,xn+ 1 …

[数学][编程] 蒙特卡洛模拟法求圆周率π (Monte Carlo method to find PI)

实干、实践、积累、思考,创新。 蒙特· 卡罗方法是一种通过概率来得到问题近似解的方法。其实之前也测试过,用在这个投针求PI([数学][概率] Buffon’s Needle problem [蒲丰投针问题]),最近又搞起算法研究,顺带再整体各种东西都测算一次吧。 蒙特· 卡罗方法求圆周率PI 网上很多例子。直接上图吧。   微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[笔记][几何][Rhino] 抛物线与合理拱轴线 [Parabola and reasonable arch axis]

实干、实践、积累、思考、创新。 根据结构力学的知识可知,抛物线为三铰拱(Three hinged arch)在承受满跨竖向均布荷载下的合理拱轴线,即 三铰拱 按抛物线布置,拱在满跨竖向均布荷载作用下,只受轴力,不受弯矩。基本简图如下图所示。 采用Rhino结合Grasshopper进行抛物线绘制。由于Grasshopper似乎没有直接进行抛物线绘制的电池(测试了interpolate似乎不对),所以抛物线需要通过计算方程的形式进行绘制。 抛物线的基本方程为 y=ax^2+bx+c,其中,a、b、c均为待定参数。通过3个点的坐标,可获得这三个参数。此外,过3点也可以绘制一条圆弧,顺带对比圆弧 arc 和抛物线 parabola的差异。 下图为采用Grasshopper进行圆弧及抛物线绘制的脚本, 通过拖拉Slider可以控制拱的跨度及高度,获得不同矢跨比的抛物线及圆弧的拱,下图为跨度/矢高 = 2.2857时圆弧及抛物线的对比,可以看出两个曲线的明显差异。 以下对比不同跨/高比下,圆弧及抛物线的差异 由图可见,跨高比越大,圆弧及抛物线的差异越小,当跨高比在5左右时,两者相差较小。 接下来,测算拱的受力,将跨高比为2.2857的一组圆弧及抛物线导入midas gen进行计算。 以下考虑3支座条件下(支座固接,支座铰接,三角拱:支座铰接及跨中铰接)抛物线拱的受力。由图可见,在承受满跨竖向均布荷载下,只有三角拱的弯矩为0,其他支座条件下,弯矩不为零。 从这个也可知,合理拱轴线是与支座条件相关的。(简直是废话。。。。。。) 此外顺便对比在承受满跨竖向均布荷载下,抛物线拱及圆弧拱的受力差异。由图可见,圆弧拱的弯矩远大于抛物线拱。 前面例子测试中提了,合理拱轴线与边界条件有关,其实合理拱轴线是在特定的荷载和边界条件下提出来的。抛物线作为合理拱轴线,主要是适用于三铰拱在满跨竖向均布荷载的条件下。如果三铰拱荷载不是满跨竖向均布荷载,那合理拱轴线就不是抛物线了。最经典的例子是如果荷载不是满跨竖向均布荷载,而是自重作用,那么三铰拱的合理拱轴线就是倒过来的悬链线了。 相关博文 ( Related Topics) [01] [GRASSHOPPER] Grasshopper-Midas …

[笔记][几何] 圆柱与方柱侧面积与体积之比 (The ratio of the side area to the volume of the cylinder and the square column)

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 圆柱与方柱侧面积与体积之比,直接看下图。 The ratio of the side area to the volume of the cylinder and the square column.   由图可见: (1)随着直径或者边长的增大,AS/V减小。 (2)直径为D的圆柱与边长为D的方柱,其AS/V是一样的。 (3) 直径为D的圆柱的体积或侧面积 只有 边长为D的方柱的体积或侧面积的0.785倍。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[编程][算法][优化] 编程训练:函数的CONLIN线性化

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。     2018 12 03 挖的坑,现在来更新…… 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[数学][算法][编程训练] 最小二乘法曲线拟合( Least square curve fitting )

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 最小二乘法曲线拟合( Least square curve fitting ),又叫最小二乘法多项式曲线拟合,根据给定的m个点,并不要求这条曲线精确地经过这些点,而是曲线y=f(x)的近似曲线y= φ(x)。 使得近似曲线与y=f(x)的偏差最小。按偏差平方和最小的原则选取拟合曲线,并且采取二项式方程为拟合曲线的方法,称为最小二乘法。在某些情况有应用。 以下是编程测试的结果。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[参数化][笔记] 采光顶参数化建模分析案例 [Rhino/Grasshopper/Midas Gen]

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 来自小伙伴 张俊毫 的参数化建模分享。 这是之前做过的一个类似的项目,此次对其进行精简处理,使生成逻辑更加清晰有条理,以便有利于参数化入门、初级水平读者的阅读。以下根据采光顶的生成逻辑,分成几部分进行讲述: 1.设置采光顶的整体控制参数,即控制采光顶整体定位参数“采光顶中心”,采光顶的体量参数“采光顶直径D”和“采光顶高H”,采光顶外观参数“花瓣宽度W”。 2.设置通过采光顶底板圆上的左、右端点及顶点三点生成采光顶外轮廓弧线,将弧线绕中心轴旋转360°可生成采光顶外表皮(此处为结构外表皮)。 3.通过如图所示三个点画弧线,并将此弧线镜像形成一个“花瓣”,并通过环形阵列每隔15°生成一个花瓣,这样就获得了一个莲花形的图形。 4.(1)求得“花瓣弧线”各个交点,并把交点投影到采光顶外表皮上,可获得经线、纬线方向的杆件节点。 (2)删除中心最高点和重复的点,并通过数据处理,使每组数据的点按标高进行排列,通过每组点生成多段线,即得经线方向杆件。   5.将经线方向上的点进行数据翻转处理(Flip Matrix),删除重合的点,并将点按照圆周顺序排列,连接每组点即可得到纬线方向的杆件。 6.运用Tree Branch Index电池获取纬线方向上的点中标高最高的一组,再对这组数据用Cull Pattern电池每隔一个点删除一个点,得到A组点,创建与此A组点标高相同的采光顶中心O,将O点与A组点连接创建向量(相邻两个向量夹角30°),将0点沿向量移动距离L,获得B组点,连接B组点可的顶部构件的同心多边形的内多边。将A、B组点一一对应连接得到一组射线,取射线的等分点(C组点),控制两个等分点的间距不小L,并将B、C组点投影到采光顶外表皮,可得B’、C’组点。这样A、B’、C’就是顶部构件全部控制点,生成杆件的方法可参照经、纬线方向杆件的生成方法。 7.完成的整体参数化模型如下,可以导入midas计算软件进行结构计算啦! 最后看看动图演示: 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[编程训练][游戏][数学] 红酒杯与白酒杯里到底有多少红酒和白酒?

实干、实践、积累、思考、创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 一个数学题,红酒杯与白酒杯的一道数学题。从结构求真群听来的。主要是这么一个问题,有一杯白酒和一杯红酒,从白酒杯里蒯一勺白酒倒入红酒杯里,然后再从红酒杯里蒯一勺倒入白酒杯里,问:白酒杯里的红酒多,还是红酒杯里的白酒多? 于是写这个小软件测试一下,结果也十分有趣,看完结果也恍然大悟。 首先,相互蒯一勺后,白酒杯里的红酒和红酒杯里的白酒是一样的多的。 而且,只要不断相互倒酒,那么白酒杯里的红酒和红酒杯里的白酒就会一直是一样多。 其实与每次倒酒的量无关,不管倒几次,不管每次倒的量是否相等,只要最终状态还是两杯一样多,那么就是一样。最终状态控制,以白酒杯为例,出去多少白酒,就会进来多少红酒。 而且,不断相互倒酒,当倒酒次数无限多厚,白酒杯里面的红酒和白酒,红酒杯里面的红酒和白酒就一样多了。 其实也很好理解,相互不断倒酒,无限多次地倒,相当于酒杯倒均匀了。 十分有意思。感兴趣的朋友可以下载来玩玩: Download: 红酒与白酒 White_RedWine.rar 欢迎关注 “结构之旅” 微信公众号  

[数学][软件] FOUR_TRAN Example 2: Square Wave Signal Decomposition [FOUR_TRAN傅里叶分析工具使用案例2: 方波信号分解]

实干、实践、积累、思考、创新。 20200205放假期间做的,现在整理出来。整个暑假,因为疫情都待在家中。 简单的傅里叶分析例子,方波信号分解。用前面自己写的 FOUR_TRAN 傅里叶分析工具进行分析。 准备如下的方波信号数据。 打开FOUR_TRAN,导入方波数据,进行傅里叶分析。 将前2阶频率的数据累加,获得下图粉色部分曲线。曲线类似一个sin函数,与方波差异较大。 将前5阶频率的数据累加,获得下图灰色部分曲线,曲线相对有之前多,与方波接近了。 将前20阶频率的数据累加,获得下图蓝色部分曲线,发现曲线以方波为中心进行震荡,并且与方波十分接近了。 接着将前50阶频率的数据累加。 以下是,将前300阶频率的数据累加的结果,可见曲线在离散点上与方波吻合了。方波可以由这么多三角函数累加而成。 绘制成EXCEL图形如下。 关于软件( About the Program)   [01] [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具] 案例 ( Examples ) [01] …

[数学][软件] FOUR_TRAN Example 1: Filtering [FOUR_TRAN傅里叶分析工具使用案例1: 滤波]

坚持实干、实践、积累、思考,创新。 滤波,傅里叶分析的一个应用。用前面自己写的 FOUR_TRAN 傅里叶分析工具进行分析。 函数 y1: y1=10*SIN(2*PI()*0.5*x),振幅为10,频率为0.5; 函数2: y2=0.5*SIN(2*PI()*10*x),振幅为0.5,频率为10; 函数yw 为函数y1与函数y2的叠加: y=y1+y2 从图可见,由于y2振幅小,频率大,y1振幅大,频率小。叠加后y整体趋势与y1类似,相当于在y1的基础上加上了一些波动。 我么记下来通过傅里叶变化,想办法把从y中把这些波动的y2过滤掉。 打开FOUR_TRAN,导入离散后的函数y,并进行傅里叶分析。 从上图可见,傅里叶变换,可准确解析出两个主要频率成分,一个是0.5,一个时10,其中频率为0.5的傅里叶振幅谱大,而频率为10的傅里叶振幅谱小,能量成分低。与我们前面的函数一致。 接着我们把频率为0.5附近的成分累加起来(逆变换),如下图红色部分,红色曲线基本就是最初的y1,y2被成功过滤掉了。 如下,红儿曲线即是过滤掉频率10后的曲线,曲线与y1一致。 接着我们把频率为10附近的成分累加起来(逆变换),得到下图的红色部分,此时红色部分基本是函数y2,y1成功过滤掉了。 这就是简单的滤波的例子。 关于软件( About the Program)   [01] [数学][地震动][软件] FOUR_TRAN: Fourier Analysis …

[数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具]

实干、实践、积累、思考、创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 简单的离散傅里叶变换(DFT,Discrete Fourier Transform)分析工具,导入数据,进行DFT变化,并获得频率Frequency,相位谱Fourier Phase,幅值谱Fourier Amplitude,功率谱 Power Spectrum等结果。并可任意组合选定频率成分,查看组合结果,滤波功能 (filtering)。直观感性地理解傅里叶变换的魅力。当然,这个东西也可以用来分析地震波。(PS。趁着在家自我隔离的时间学习,希望新型冠状病毒肺炎疫情早日被战胜。) 程序界面 ( Program Interface) 下载 (  Download ) (  如果您需要使用这个软件,请在这给网站捐助:http://www.jdcui.com/?page_id=4813,捐助【250RMB】,捐助后在评论区留下您的评论,软件将通过网站管理员QQ:3014479529给您发送,敬请尊重劳动成果 !!!) 案例 ( Examples )   …