翰林苑|董翰林的个人网站

求解一元三次方程的原理与代码

By 董翰林; Published on 2023-01-29

本文记录了两种求解一元三次方程的方法：Cardano （卡丹）公式法和伴随矩阵特征值法。其中详细记录了 Cardano 公式的推导过程，以及使用三角函数对复数进行处理，避免复数运算的方法。然后简要介绍了多项式的伴随矩阵和通过伴随矩阵求特征值来获得根的方法。对各种方法分别给出了使用 C++ 以及 Python 求解单一方程和批量方程的代码。并尝试了使用 Pybind11 将 C++ 代码封装为...

使用 Python 制作倒计时视频

By 董翰林; Published on 2022-09-21

最近有一个在 Powerpoint 中设置倒计时器的需求。由于 Powerpoint 没有提供原生的功能支持，网上找到了一些插件又觉得安装起来不便，所以决定使用一个倒计时视频放置在 Powerpoint 中，即可实现倒计时的功能。由于不同用户所需要的倒计时时间不同，因此没办法找到用于倒计时的现成视频。因此想到使用 Python 的 Matplotlib 库写一个倒计时视频。代码如下。 12...

我的装机和科研软件合集

By 董翰林; Published on 2022-08-31

找软件是一个非常费时的事情，每次电脑重装，都要到处找软件，着实让人头疼。开源软件还好，如果是商业软件，就更加让人头疼。一直想着搞一个软件合集。如今我重装一台电脑，就趁机把我最常使用的软件列成表，并附上官网连接。如果读者恰好有和我一样的需求，就可以通过我的连接下载。有一些软件是有教育版的，可以使用 edu 后缀的教育邮箱注册下载。我选软件的原则比较简单：功能务实，安装纯净。也欢迎读者投稿推荐...

使用 Python 中的 os.walk 遍历文件夹

By 董翰林; Published on 2022-07-28

Python 的 os 模块提供了一系列访问系统中文件的方法。使用 getcwd() 方法可以获得当前的工作路径，使用 listdir() 方法可以列出路径中的所有文件和文件夹，至于是文件还是文件夹，需要使用 os.path 中的 isdir() 方法来判断。用这些方法遍历文件夹比较麻烦，因此提供了一个 walk 方法用于遍历文件夹。 walk 方法的设计比较巧妙。它的参数是要遍历的目录，返...

给自己理发的人

By 董翰林; Published on 2022-05-20

作为上海第一批居家健康监测的人，转眼已经在家三个月了。春节前理的发肆无忌惮地长着，爬过耳廓，就要进军到眉毛了。于是我翻出了尘封的理发推子，在自己的头上动起刀来。这个理发推子是在加拿大买的，尽管是加拿大亚马逊售卖，但收到货后可以明显地看出它是Made in China。事实上，用过几次加拿大亚马逊，买到的东西有一半以上是直接产自中国的。还有海外版的淘宝，只要不怕时间长，很多东西可以直接从国内...

土豆发芽了还能吃吗

By 董翰林; Published on 2022-04-30

我的童年在奶奶家度过。东北人有屯菜的传统，奶奶也不例外。尽管她的身体单薄，但屯起菜来，还是要以三十斤、五十斤为单位的。尤其是土豆这样容易存放的食物，必定要买整筐。尽管土豆不易变质，但是存放不当，或时间过长，还是要发芽的。于是，吃发芽土豆的技能，我也从小就练成了。通常所吃的土豆长出的芽不大，白白嫩嫩的，像一棵小豆芽。然而与豆芽是人间美味不同，土豆芽因为一字之差，变成了一种剧毒产品，真是太冤了...

编译 OpenSees 的方法

By 董翰林; Published on 2021-12-10

为了帮助大家快速学习编译 OpenSees 的方法，我制作了下面的视频。在下载源代码时，可以使用最新的代码编译，但是由于开发者把每次更新都直接推送到 master 分支，最新的代码不一定能成功编译。因此在视频中，我将代码的版本回退到了 v3.3.0 发布时的源码。编译过程中最重要的是指定 .h 头文件的路径，以及 .lib 库文件的路径。只要把这两个路径找好，就可以成功编译。每台电脑的系...

被手机奴役的人

By 董翰林; Published on 2021-09-05

人类的每一次技术革命，都是被加深奴役的过程。昨天晚上，我的华为 Mate 20 Pro 在自动更新了操作系统后，瘫痪了。手机重启后，输入锁屏密码，就陷入下一次重启。一开始我不以为意，以为手机过热，关机冷却一会儿，仍不奏效。想在电脑上联系客服，可是只提供了电话、邮件和在线聊天三种方式。电话是打不了了，邮件太慢，于是选择在线聊天。然而打开在线聊天，要同意服务协议并登录，要登录又要输入手机验证...

OpenSees 分析参数的选择方法（上）

By 董翰林; Published on 2021-08-17

在结构力学中，我们知道，结构的线性分析一般流程如下：为结构中的节点和单元编号，并建立整体坐标系；对各单元计算局部坐标系下的单元刚度矩阵 \(\bm{\bar{K}^e}\) ，并通过坐标变换矩阵 \(\bm{T}\) ，形成整体坐标系下的单元刚度矩阵 \(\bm{K^e}\) ；根据结构的约束，对自由度编号，并将单元刚度矩阵拼装为整体刚度矩阵 \(\bm{K}\) ；由节点荷载和单...

在 Bash 后台运行命令

By 董翰林; Published on 2021-07-26

在使用 OpenSees 进行批处理等耗时较长的操作时，有时希望多个命令在同一控制台同时运行。这时就涉及到在后台运行命令。在命令后面加入一个 & 符号可以实现前台与后台之间的切换。例如 1OpenSees model.tcl & 这样，OpenSees 会在后台执行。但是此时，错误信息还会打印到控制台中。如果不希望打印在控制台中，而是打印成为 log 文件，可以使用重定向...

VISITORS

次阅读

FEATURED TAGS

ABOUT ME

君子博学于文约之以礼

CONTACT

董翰林 博士
同济-多伦多大学联合培养博士
上海师范大学讲师
self@hanlindong.com