我已经给你们算好了

徐洋猛地张了张嘴，眉头微蹙，脑子飞速运转个不停。

　　他从事伺服控制领域已经有十几年了，对分数阶微积分也并不陌生，可那东西向来被业内人视作“屠龙之技”的。

　　在理论层面堪称精致完美，看似无懈可击，可落到实际工程应用里，却处处都是难以突破的壁垒，根本没法真正落地，真正把这个方法用到实际实验中的，少之又少。

　　可肖宿刚才给出的频带范围和近似阶次，精确得像是已经把仿真跑过几百遍一样。

　　“肖教授，”徐洋的声音有些发干，“分数阶项的离散化实现方案，您有更具体的……内容吗？”

　　肖宿扫了他一眼，转身在白板上写出了递归滤波器的零极点配置公式。

　　“用Oustaloup递归滤波器，五阶逼近就够了，零极点交错分布在指定频带内。”

　　“极点的频率分布可以采用几何级数，公比取10^(1/N)，N是逼近阶数，五阶逼近在你们关心的频段内幅频误差不会超过，相频误差不会超过2度。”

　　实验室里安静了几秒钟。

　　角落里一个年轻工程师小声问旁边的同事：“他刚才说的那个Ous什么东西，你听说过吗？”

　　同事摇了摇头，脸色发白。

　　徐洋顾不上管身后的动静，他脑子里正在疯狂地把肖宿给出的传递函数代入到他们现有的系统框图中去推演。

　　他越推越快，越推越心惊。

　　这个方案不是看起来可行，而是几乎不需要做任何额外的理论验证，拿过来就能往DSP里写代码了。

　　“可是肖教授，”徐洋突然想到了什么，“分数阶PID的参数整定，用相位裕度法需要用到一个参考模型，我们的参考模型又该怎么建呢？”

　　“可以用DenHartog的临界灵敏度法变形，把整数阶的Nyquist判据推广到分数阶域。”

　　肖宿把笔帽盖上，语气平淡得像在说一加一等于二，“推广的步骤很简单，构造一个Riemann曲面上的复变函数，把分数阶项映射到曲面的不同叶上，然后在主叶上做Nyquist围线积分。

　　Jacobian矩阵的结构和李群上的切空间有自然的对应关系，这一步我在另一篇论文里详细写过，回头你们可以看一下。”

　　徐洋的脑袋已经开始发胀了。

　　他搞了十几年伺服控制，自认为在国内算得上顶尖水平，可肖宿刚才说的话里，他能完全听懂的不到一半。

　　沈凌比徐洋好不到哪去。

　　他是做图像重建的，控制理论不是他的本行，但他已经敏锐地意识到了一个问题，那就是肖宿刚才解决伺服延迟的方法，和他接下来要问的成像畸变问题，在数学结构上似乎是相通的。

　　“肖教授，”沈凌深吸一口气，往前迈了一步，“成像这边的问题，您刚才也听到了。

　　我们的多视图几何重建在动态场景下会出现帧间畸变，我用梯度下降和卡尔曼滤波都试过了，效果还是不太理想。”

　　“因为你的问题不在滤波。”肖宿转身面对他，“你把你们的投影矩阵写给我看看。”

　　沈凌立刻拉过一台笔记本电脑，调出他们当前使用的三维重建算法的核心代码，把投影矩阵的定义部分展示在屏幕上。

　　肖宿只是上下扫了一遍就知道他们问题在哪儿了。