同样一片膜可以有几种稳定形状

先选一族能算清的曲面

给囊泡选一根对称轴，用极角 $\theta$ 从北极往南极走。径向距离 $r(\theta)$ 决定曲面。我们用两个低阶 Legendre 多项式描述它：二阶项主要区分长与扁，四阶项允许轮廓进一步调整。指数让半径始终为正。

整体尺度 $s$ 用来满足面积，系数 $a$ 和 $b$ 改变形状。这个有限变分族光滑、轴对称且相对中心呈星形：每条射线只穿过膜一次。反复穿过同一射线的凹陷、极窄颈和不对称形状超出了这个范围；受限的长形和扁形仍能提供可控的能量比较。

r(\theta)=s\exp\!\left[aP_2(\cos\theta)+bP_4(\cos\theta)\right]

这里

P_2(x)=(3x^2-1)/2

，

P_4(x)=(35x^4-30x^2+3)/8

。参数无量纲，尺度

s

以长度计。

来源：U. Seifert · W. Helfrich

把面积和体积真的固定住

先临时令整体尺度为一，沿极角积分求面积和体积。再乘一个尺度，让面积恰好等于目标值。整体缩放不会改变约化体积，因此剩下的工作是找到满足目标约化体积的形状系数。

固定一个 $b$ ，沿正、负两侧分别寻找 $a$ 。二分法每次取一个区间的中点，判断约化体积位于目标哪边，再保留夹住解的半个区间。找到两支允许形状后，才比较它们的弯曲能；不能拿体积不同的两个漂亮轮廓直接排序。

最后在四阶系数的有限范围内寻找每一支的最低能量。画出的结果是这族曲面、这个搜索范围内的受限极小值。真实囊泡的全局最低能量还要允许更多形状，并检查自接触与稳定性。

A=2\pi\int_0^\pi r\sin\theta\sqrt{r^2+r_\theta^2}\,\mathrm d\theta,\qquad V=\frac{2\pi}{3}\int_0^\pi r^3\sin\theta\,\mathrm d\theta

面积来自小圆环的周长乘以母线长度；体积来自每个立体角内的径向体积。

s=\sqrt{\frac{A_{\rm target}}{A_{s=1}}},\qquad v(a,b)=v_{\rm target}

第一式消去尺度，第二式消去一个形状自由度。剩下的参数沿受限曲线变化。

完整例题：把约束残差换成物理大小

取 $R_A=10\,\mu\mathrm m$ 、 $v=0.80$ 。目标面积为 $1.26\times10^3\,\mu\mathrm m^2$ ，目标体积约为 $3.35\times10^3\,\mu\mathrm m^3$ 。假设某个候选形状的相对体积残差为 $10^{-4}$ ，对应水量误差约 $0.335\,\mu\mathrm m^3$ 。

这一步先检查它是不是同一个问题的解。然后把角向积分节点加倍，比较能量、约化体积和高斯曲率积分。数值不断接近同一个答案，才有理由认为差别来自模型，而不是轮廓离散得太粗。

球面还给出独立标尺：零自发曲率时，平均曲率能除以 $\kappa$ 应接近 $8\pi$ ，高斯曲率积分应接近 $4\pi$ 。这两个检查分别检验局部曲率与闭合表面的整体几何。

形状照片缺少哪个尺度

假设能量只有平均曲率项，面积、体积和自发曲率都已给定。把 $\kappa$ 乘以任意正数，会把所有候选能量同时乘大，但不会移动极小值。即使轮廓测得完美，形状也无法单独确定弯曲刚度。

若再把未知张力、压力或蛋白分布一起拟合，可互相补偿的参数会更多。拟合程序返回一个数字，只表示它在当前假设下找到了一个答案。测热涨落、施加已知力或加入独立压力测量，才能补上形状里缺少的信息。

这也解释了为什么换一套材料参数有时几乎看不出轮廓差异。更可靠的反演方案先找出观测对哪些参数敏感，再选择第二种实验。下一段课程就用热运动提供一个独立能量标尺。

F_b=\kappa\,\mathcal E[\text{shape}],\qquad \arg\min F_b=\arg\min\mathcal E\quad(\kappa>0)

这个结论限于刚度作为整体正因子的模型；加入其他独立能量尺度后，需要重新分析。

同样一片膜可以有几种稳定形状

先选一族能算清的曲面

把面积和体积真的固定住

完整例题：把约束残差换成物理大小

形状照片缺少哪个尺度

把推理带进实验

比较同一组约束下的两个分支

练习与解释性答案

来源与复算范围