这篇帖子要解决的问题通俗的讲就是，我是把更多电容放在前面效率高，还是放在后面效率高。

首先把问题变成这个形式：对于一个磁阻式电磁炮，已知弹丸初速为v0，出速为v1（v1>v0）。总加速距离为x1（x1>0）。求如何分配各级的加速度，能使效率最高。

引入一个前提条件：整个系统的电阻损耗功率Pr和加速度a的平方成正比，且与其它变量（比如弹丸正在哪一级的哪个位置）无关，即$P_r=k{a}^2$

这个前提条件对应的实际情况大概是：弹丸始终在磁饱和状态下；分级足够精细，各级紧密相邻，且各级均带关断与能量回收。

结论是：若加速度$a=(C_{3}x+C_4{)}^{\frac{1}{3}}$，则效率可能最高。其中C3和C4是和v0，v1，x1有关的常数。

求解过程如下：

以位移为零的时刻为时间起点，即$t{|}_{x=0}=0$，设位移为x1时，时间为t1（t1不是已知量）。设电阻损耗的总能量为$E_r$。则此问题可以描述为如下形式：

在满足固定边界条件$v{|}_{x=0}=v_0,v{|}_{x=x_1}=v_1$的情况下，求函数v(x)，使泛函$E_r[v(x)]=k{\int }_0^{t_1}{a}^2\mathrm{d}t$取极小值。

∵$v=\mathrm{d}x/\mathrm{d}t,a=\mathrm{d}v/\mathrm{d}t$

∴$\mathrm{d}t=\mathrm{d}x/v,a=v\cdot \mathrm{d}v/\mathrm{d}x=v\cdot v^{\prime },(v^{\prime }=\mathrm{d}v/\mathrm{d}x)$

∴$E_r[v(x)]=k{\int }_0^{x_1}v\cdot v^{\prime 2}\mathrm{d}x$

设$G(x,v,v^{\prime })=v\cdot v^{\prime 2}$

对于物理可实现的系统，$G(x,v,v^{\prime })$拥有连续的二阶偏导数

∴由欧拉—拉格朗日方程可知，Er取极值的必要条件是

$$G_v-\frac{\mathrm{d}}{\mathrm{d}x}{G}_{v^{\prime }}=0$$

即

$$v^{\prime 2}+2v\cdot v^{\prime \prime }=0$$

这是一个可降阶的二阶微分方程，可以求出

$$v^{\prime }=C_0{v}^{-\frac{1}{2}}$$

进而求出

$$v=(C_{1}x+C_2{)}^{\frac{2}{3}}$$

代入边界条件可以求出C1，C2的值，求出的结果形式上比较复杂 ，mathematica给出的结果是这样的，里面可能有些解不符实际需要舍去

$$a(x)=\frac{2}{3}{C}_1(C_{1}x+C_2{)}^{\frac{1}{3}}=(C_{3}x+C_4{)}^{\frac{1}{3}}$$

对于常见的初始速度为0的情况，易知C2=0，所以对于这种情况有

$$a=C_5{x}^{\frac{1}{3}}$$

所以初速为0时，若加速度与位置的三分之一次方成正比，则可能得到最高的效率。

至此求解完毕

然后说点闲话。记得我之前曾经发过一篇帖子说效率最高的是匀加速，当时没有发解法（当时本来是打算把另一个问题解决之后一起发出来的，后来发现另一个问题解不动……）结果前两天突然发现，我把当时的解法给忘了……然后又翻了一遍书重做了一遍，发现当时的结论不对，估计是边界条件给错了。当然这个也不保证对，毕竟欧拉方程只给出必要条件，是不是充分条件还要另求……而且关于位移，速度，加速度的那个部分我还有点迷糊……