一直对欧洲粒子研究中心的印象都很好,不过仅限于计算机领域的一些管窥蠡测,比如,world wide web的创始者Tim Berners-Lee就是CERN的,再据说,CERN的科学家们喜欢用TeX排版。

去年风闻世界上最大的粒子对撞机LHC建好了,开始运行,心里蛮期待的。结果故障频繁出现,吊人胃口,直到一个多星期以前,也就是愚人节的前一天才成功进行第一次碰撞。好吧,既然成功了,我们就擦亮眼睛,看看什么结果,最近有人八卦

近日欧洲物理学会专业期刊Phys. Lett. B以最快的审稿、编辑和出版速度发表了大型强子对撞机ATLAS合作组的首个物理成果。在这篇22页的论文中,从标题到参考文献只占了10页,其他12页用于书写上千位作者和他们的单位,看起来十分滑稽。

这篇论文可以在这里找到,题目为Charged-particle multiplicities in pp interactions at sqrt(s) = 900 GeV measured with the ATLAS detector at the LHC,带电粒子的多样性,具体内容懒得仔细看,估计多半也看不太懂。大致感觉,只是打着LHC的名号,随便翻点花样,功利心居多。

花费如此巨大金钱,人力物力,是否值得,也是见仁见智了,不过感觉科研走到这种地步,丧失初衷居多。科研也许一张纸一支笔,加上兴趣就行了。当然科学家有偏理论研究的,有偏实验的,倒也无可厚非。研究是什么啊,可能也就是几个小孩对一些事物的执着兴趣,捣鼓点小小的经验规律,在大自然的面前,不都是渺小之极么?

现在感觉CERN有点像个迂腐的学术机构,急功近利,忙于出成果。不过,当一个组织变得很大的话,各方面因素影响,也许很多东西都是不可避免的吧。

Bezier曲线很常用,一般2D绘图软件里都有。比如photoshop,flash之类。
它背后的原理简单的超乎想象,体现了数学的美妙。

先从简单的开始,两个点之间进行线性插值。


很容易理解,可以得到

[latex] B(t)=P_0+t(P_1-P_0)=(1-t)P_0+tP_1,tin[0,1] [/latex]

当然这是最简单的情形,如果扩展到三个点该如何插值呢?


从上面图片上可以看到,可以分成三步,从P0到P1进行上面的一维情形,得到点Q0,再从P1到P2,得到Q1,那么就有

[latex]Q_0=(1-t)P_0+tP_1[/latex]
[latex]Q_1=(1-t)P_1+tP_2[/latex]

然后再对Q0和Q1进行线性插值,得到点B

[latex]B=(1-t)Q_0+tQ_1=(1-t)^2P_0+2t(1-t)P_1+t^2P_2[/latex]

t从0到1增加,就得到了一条曲线,如下图


同样可以推广到四个点的情形,这样的曲线中,t的最高幂是三次。三次样条曲线用的最多,因为它提供了足够的可控制性和满足大部分场合的精度,同时又保持了相对的简单。
依照上面的方法,可以得到三次的情形

[latex]B(t)=(1-t)^3P_0+3(1-t)^2tP_1+3(1-t)t^2P_2+t^3P_3,tin[0,1][/latex]

这个时候可以注意到Pi点前面的系数,是不是似曾相识?没错,就是二项式公式。可以看成是相应次幂的[latex](1-t+t)^n[/latex]的展开。这个系数叫做Bernstein多项式。

最常用的三次曲线如下图,其中中间的两个就是控制点,在一些绘图软件里用钢笔拖出来的两条调整曲线形状的直线,就是调节中间两个点的位置。


推广到任意中n的情况

[latex size=”2″]B(t)=sum_{i=0}^{n}{binom{n}{i}(1-t)^{n-i}t^iP_i}[/latex]

[latex size=”1″]=(1-t)^nP_0+binom{n}{1}(1-t)^{n-1}tP_1+..+t^nP_n,tin[0,1][/latex]

不过n大于3的时候就很少用了,除非在一些要求比较高的场合,比如飞机汽车线形的设计。

更多的资料,可以看这里
http://en.wikipedia.org/wiki/B%C3%A9zier_curve

今天又想来mathmetica来了,本来是懒得拿笔推导一个比较繁琐的公式,结果一不小心玩上瘾了,一天基本上都在搞这个软件。

呵呵,不说了,放张图片吧