goodisok 的博客

本文是《复变函数完全入门》系列的第五篇（完结篇），覆盖复势的概念、用复变函数描述基本势流解、圆柱绕流的复势形式和儒科夫斯基翼型的完全解。本文串联起前三篇的所有数学工具——将它们应用到空气动力学势流理论的实际问题中。 第五篇：复变函数在空气动力学中的应用——复势与翼型绕流一、复势——统一描述二维势流1.1 概念引入在第二篇中我们学过：任意解析函数 f(z)=ϕ+iψf(z) = \phi + i\psif(z)=ϕ+iψ 的实部和虚部都是调和函数——即满足拉普拉斯方程。 在空气动力学中： ϕ\phiϕ = 速度势 ψ\psiψ = 流函数 f(z)=ϕ+iψf(z) = \phi + i\psif(z)=ϕ+iψ 称为**复势**（Complex Potential） 所以：每一个解析函数对应一种可能的二维不可压缩无旋流动。 1.2 从复势得到流场信息速度场从复势的导数获得——复速度： dfdz=∂ϕ∂x+i∂ψ∂x=u−iv \frac{df}{dz} = \frac{\partial\phi}{\partial x} + i\frac{\partial\psi}{\p...

本文是《复变函数完全入门》系列的第四篇，覆盖保角变换的定义、性质、儒科夫斯基变换及其在翼型空气动力学中的应用。阅读本文前需掌握解析函数的定义和CR方程。 第四篇：保角变换——从圆到翼型的数学桥梁一、概念引入——为什么需要保角变换1.1 问题势流理论中，圆柱绕流已经有了精确解析解。但工程师关心的是翼型绕流，不是圆柱绕流。 有没有一种方法，能把圆柱绕流”变形”成翼型绕流？ 有——保角变换。 核心思想： 在 ζ\zetaζ 平面（计算平面）上，我们有一个圆——绕流解已知 通过一个解析函数 z=f(ζ)z = f(\zeta)z=f(ζ)，把圆映射为翼型形状 圆柱绕流解同步映射——就得到了翼型绕流解 1.2 历史的五分钟 时间 人物 贡献 1851 黎曼 黎曼映射定理：任何单连通区域都可以保角映射到单位圆 1910 儒科夫斯基（Zhukovsky） 提出儒科夫斯基变换——从圆到翼型的精确映射 1922 库塔（Kutta） 应用儒科夫斯基变换计算翼型升力 1933 西奥多森（Theodorsen） 提出一般翼型的数值保角变换方法 二、保角变换的定义2....

本文是《复变函数完全入门》系列的第三篇，覆盖复积分、柯西积分定理、柯西积分公式和留数定理。阅读本文前需了解解析函数的定义和CR方程。 第三篇：复积分与留数定理——环路积分的威力一、概念引入——从实积分到复积分1.1 定义复变函数的积分定义为： ∫Cf(z) dz=∫C[u(x,y)+iv(x,y)][dx+i dy]=∫Cu dx−v dy+i∫Cv dx+u dy \int_C f(z)\,dz = \int_C [u(x,y) + iv(x,y)][dx + i\,dy] = \int_C u\,dx - v\,dy + i\int_C v\,dx + u\,dy ∫C​f(z)dz=∫C​[u(x,y)+iv(x,y)][dx+idy]=∫C​udx−vdy+i∫C​vdx+udy 其中 CCC 是复平面中的一段曲线。 若曲线参数化为 z(t)=x(t)+iy(t)z(t) = x(t) + iy(t)z(t)=x(t)+iy(t)，a≤t≤ba \le t \le ba≤t≤b： ∫Cf(z) dz=∫abf(z(t)) z′(t) dt \int_C f(z)\...

本文是《复变函数完全入门》系列的第二篇，覆盖解析函数的定义、柯西-黎曼方程、初等解析函数和调和函数。阅读本文前需熟悉复数的极坐标表示和欧拉公式。 第二篇：解析函数与柯西-黎曼方程——复可微的苛刻条件一、概念引入——复变函数是什么1.1 定义一个复变函数 w=f(z)w = f(z)w=f(z) 将一个复数 z=x+iyz = x + iyz=x+iy 映射到另一个复数 w=u+ivw = u + ivw=u+iv。 f(z)=u(x,y)+i v(x,y) f(z) = u(x,y) + i\,v(x,y) f(z)=u(x,y)+iv(x,y) 其中 u(x,y)u(x,y)u(x,y) 和 v(x,y)v(x,y)v(x,y) 都是实函数。从几何上看，复变函数就是一个从 C\mathbb{C}C 到 C\mathbb{C}C 的映射——把复平面上的点搬到另一个复平面上。 最简单的例子：f(z)=z2f(z) = z^2f(z)=z2。 令 z=x+iyz = x + iyz=x+iy： f(z)=(x+iy)2=x2−y2+i(2xy) f(z) = (x+iy)^2...

本文是《复变函数完全入门》系列的第一篇，覆盖复数的定义、运算、几何表示、欧拉公式和复平面。阅读本文仅需掌握高中数学中的实数运算和三角函数。本系列的目标是为理解势流理论中的保角变换打下完整的数学基础——保角变换是连接复变函数和空气动力学势流理论的数学桥梁。 第一部分：为什么要有复数一、概念引入——一个简单方程引发的”数字危机”1.1 问题解方程 x2+1=0x^2 + 1 = 0x2+1=0。 实数范围内无解——因为任何实数的平方都非负。但如果我们”强行”定义一个数 iii，使得 i2=−1i^2 = -1i2=−1，那么方程的解就是 x=±ix = \pm ix=±i。 这个看起来很简单的”强行定义”，在16世纪引起了数学界激烈的争论。当时的数学家将 iii 称为”虚数”（imaginary number）——这个名称本身就带有贬义。 1.2 历史的五分钟 时间 人物 贡献 1545 卡尔达诺（Cardano） 在《大术》中首次使用 −15\sqrt{-15}−15​ 计算三次方程的解——“虚数”的首次出场 1572 邦贝利（Bombelli） 系统地使用了 ...

本文是《微分方程完全入门》系列的第六篇（完结篇），覆盖偏微分方程的基本概念、分类、分离变量法和拉普拉斯方程。阅读本文前需掌握二阶线性ODE。本篇最终目标是理解势流理论中的拉普拉斯方程——nabla2phi=0\\nabla^2\\phi=0nabla2phi=0。 第六篇：偏微分方程引论——从热方程到拉普拉斯方程一、概念引入——从一维到多维常微分方程中的未知函数只有一个自变量（通常是时间 ttt）。但现实世界中的场量——温度 T(x,y,z,t)T(x,y,z,t)T(x,y,z,t)、压力 p(x,y,z,t)p(x,y,z,t)p(x,y,z,t)、速度势 ϕ(x,y,z,t)\phi(x,y,z,t)ϕ(x,y,z,t)——都是多个自变量的函数。 偏微分方程（PDE）描述的就是这种多自变量函数的演化或分布规律。 1.1 历史的5分钟 时间 人物 贡献 1747 达朗贝尔 提出波动方程——最早的PDE之一 1777 拉普拉斯 拉普拉斯方程 ∇2ϕ=0\nabla^2\phi=0∇2ϕ=0——势流的控制方程 1807 傅里叶 提出热传导方程和傅里叶级数解...

本文是《微分方程完全入门》系列的第五篇，覆盖ODE数值解法的基本方法：欧拉法、改进欧拉法、四阶龙格-库塔法（RK4），以及自适应步长和刚性方程的概念。阅读本文前需理解一阶ODE的基本形式。本文重点在方法的几何直觉和误差分析，而非纯数值分析视角。 第五篇：数值解法——从欧拉法到 RK4一、概念引入——为什么要数值解法在之前四篇中，我们学习了多种解析解法——分离变量、积分因子、特征根法、拉普拉斯变换。 但这些方法都有一个致命的限制：大多数现实中的ODE没有解析解。 考虑一个稍微修改的摆方程——加入非线性项和非正弦驱动力： d2θdt2+gLsin⁡θ=F0cos⁡(ωt) \frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = F_0\cos(\omega t) dt2d2θ​+Lg​sinθ=F0​cos(ωt) 这个方程——受迫阻尼摆——无法用初等函数表示其解。但它描述的就是真实世界的物理。 数值解法就是回答这个问题的：用计算机逐”小步”逼近解。 二、欧拉法——最直观的起点2.1 思想一个一阶 ODE：dydt=f(t,y)\f...

本文是《微分方程完全入门》系列的第四篇，覆盖一阶常微分方程组的相平面分析、平衡点分类和李雅普诺夫稳定性。阅读本文前需掌握特征根的概念——我们将把特征根从”单个方程”推广到”方程组”。 第四篇：微分方程组与相平面——从一维到多维一、概念引入——为什么需要方程组一个四旋翼无人机的动力学需要 12 个一阶 ODE 来描述——位置(3)+速度(3)+姿态(3)+角速度(3)。这不是单个方程能处理的。 一阶 ODE 方程组的一般形式： {dx1dt=f1(x1,x2,…,xn)dx2dt=f2(x1,x2,…,xn)⋮dxndt=fn(x1,x2,…,xn) \begin{cases} \frac{dx_1}{dt} = f_1(x_1, x_2, \dots, x_n) \\ \frac{dx_2}{dt} = f_2(x_1, x_2, \dots, x_n) \\ \vdots \\ \frac{dx_n}{dt} = f_n(x_1, x_2, \dots, x_n) \end{cases} ⎩⎨⎧​dtdx1​​=f1​(x1​,x2​,…,xn​)dtdx2​​=f2​(...

本文是《微分方程完全入门》系列的第三篇，覆盖拉普拉斯变换的定义、性质、逆变换及其在解ODE中的应用。阅读本文前需掌握一阶和二阶ODE的基础知识。掌握拉普拉斯变换后，你会发现：解微分方程变成了做代数运算——这是控制理论中传递函数概念的数学起点。 第三篇：拉普拉斯变换——从时域到频域的桥梁一、概念引入——为什么需要一个”变换”1.1 问题解一个二阶ODE需要： 求特征方程的根 确定齐次解形式 猜测特解（待定系数法，常常还要处理”共振”的例外） 代入初始条件确定任意常数 如果是高阶系统（四阶、六阶），这个流程变得极其繁琐。更糟的是，待定系数法的每一步都需要针对 f(t)f(t)f(t) 的形式特殊处理。 有没有一种方法，可以让求解 ODE 变得像解代数方程一样简单？ 有：拉普拉斯变换。它的核心思想是： 把时域中的微分方程，通过一个积分变换，变成频域中的一个代数方程。解完代数方程后，再变换回时域。 1.2 历史 时间 人物 贡献 1744 欧拉（Euler） 首次使用了类似拉普拉斯变换的积分形式求解ODE 1785 拉普拉斯（Pierre-Simon Lapl...

本文是《微分方程完全入门》系列的第二篇，覆盖二阶线性常微分方程的完整理论。阅读本文前需熟悉一阶ODE的解法。掌握本篇后你将理解四旋翼飞行动力学的”阻尼比”和”自然频率”这两个核心参数。 第二篇：二阶线性方程——弹簧、阻尼与共振一、二阶 ODE 从哪来——三个物理问题1.1 概念引入我们身边的很多物理系统都可以用一个统一的方程来描述： 弹簧-质量-阻尼系统： md2xdt2+cdxdt+kx=F(t) m\frac{d^2x}{dt^2} + c\frac{dx}{dt} + kx = F(t) mdt2d2x​+cdtdx​+kx=F(t) RLC电路： Ld2qdt2+Rdqdt+1Cq=E(t) L\frac{d^2q}{dt^2} + R\frac{dq}{dt} + \frac{1}{C}q = E(t) Ldt2d2q​+Rdtdq​+C1​q=E(t) 四旋翼俯仰角动力学（简化）： Id2θdt2+cθdθdt+kθθ=τ(t) I\frac{d^2\theta}{dt^2} + c_\theta\frac{d\theta}{dt} + k_\theta\t...