9 Aug

三次方程求根器(VB程序+源码,“低手”拙作)

三次方程求根器-界面

三次方程求根器-界面

点击阅读全文...

29 Aug

计算夏至的精确时刻2——提高精确度

之前曾经得到过一条计算夏至精确时间的公式,现在检验一下(之前推导是根据了2009年的数据)

公元Y年的夏至日期为该年的6月
$$21.9938+0.2422Y-\lfloor Y/4 \rfloor-\lfloor Y/400 \rfloor+\lfloor Y/100 \rfloor$$
其中$\lfloor x \rfloor$表示整数部分。

点击阅读全文...

10 Sep

级数求和——近似的无穷级数

级数是数学的一门很具有实用性的分支,而级数求和则是级数研究中的核心内容之一。很多问题都可以表示成一个级数的和或积,也就是$\sum_{i=1}^n f(i)$或者是$\prod_{i=1}^n f(i)$类型的运算。其中,$\ln(\prod_{i=1}^n f(i))=\sum_{i=1}^n \ln(f(i))=k$,因此$\prod_{i=1}^n f(i)=e^k$,也就是说,级数求积也可以变为级数求和来计算,换言之我们可以把精力放到级数求和上去。

为了解决一般的级数求和问题,我们考虑以下方程的解:
$$f(x+\epsilon)-f(x)=g(x)\tag{1}$$其中g(x)是已知的以x为变量的函数式,$\epsilon $是常数,初始条件是$f(k)=b$,要求f(x)的表达式。

点击阅读全文...

30 Oct

太阳帆技术的粗浅分析(补充)

上星期,BoJone凭借简陋的物理知识,发表了《太阳帆技术的粗浅分析》一文,并转到了牧夫天文论坛上,希冀能够抛砖引玉。很幸运得到了牧夫上的高手的指正。他们指出了我的文章中$a=a_{ray}-a_G > 0$这一条件过于苛刻。因为,除了太阳光压外,还有另外一种力量能够战胜太阳引力——惯性离心力

重新把上篇文章的一个结果列出来:
$$a=a_{ray}-a_G=(\frac{L}{2\pi c (\rho h+{m'}/S)}-GM_{sun})\frac{1}{r^2} $$

点击阅读全文...

28 Nov

《自然极值》系列——3.平衡态公理

黄果树大瀑布

黄果树大瀑布

光学定律无疑是一个美妙的原理,而自然界中还存在另外一个我们随处可见的“公理”。平时的生活中,我们总能看见“水往低处流”的现象,这是因为水处于地球重力场的结果(也正因为如此,某些轻生者的自杀活动才得以顺利进行;当然,我们并不需要为了验证这一点而亲自试验。)。由此我们可以联想到一个名词:重力势能。“水往低处流”意味着什么呢?高度变低了。高度更低意味着什么呢?重力势能降低了!换句话说,自然界中物体有趋于势能最低的倾向。我们可以从这个角度来解释:体系总有趋于稳定的倾向,而拥有的能量(势能)越高,则越不稳定。

点击阅读全文...

10 Dec

《自然极值》系列——6.最速降线的解答

通过上一小节的小故事,我们已经能够基本了解最速降线的内容了,它就是要我们求出满足某一极值条件的一个未知函数,由于函数是未知的,因此这类问题被称为“泛分析”。其中还谈到,伯努利利用费马原理巧妙地得出了答案,那么我们现在就再次回顾历史,追寻伯努利的答案,并且寻找进一步的应用。

最速降线-1

最速降线-1

为了计算方便,我们把最速降线倒过来,把初始点设置在原点。在下落过程中,重力势能转化为动能,因此,在点(x,y)处有$\frac{1}{2} mv^2=mgy\Rightarrow v=\sqrt{2gy}$,由于纯粹为了探讨曲线形状,所以我们使g=0.5,即$v=\sqrt{y}$。在点(x,y)处所走的路程为$ds=\sqrt{dy^2+dx^2}=\sqrt{\dot{y}^2+1}dx$,所以时间为$dt=\frac{ds}{v}=\frac{\sqrt{\dot{y}^2+1}dx}{\sqrt{y}}$,于是最速降线问题就是求使$t=\int_0^{x_2} \frac{\sqrt{\dot{y}^2+1}dx}{\sqrt{y}}$最小的函数。

点击阅读全文...

26 Feb

线圈感抗和电容容抗的计算

形形色色的电容

形形色色的电容

学到人教版高二物理选修3-2的同学们,眼前会出现许多新的名词,如楞次定律、自感(电感)、感抗、容抗等等。其中对于电感,在中文维基百科给予的解释为:当电流改变时,因电磁感应而产生抵抗电流改变的电动势(EMF,electromotive force)。电路中的任何电流,会产生磁场,磁场的磁通量又作用于电路上。依据楞次定律,此磁通会借由感应出的电压(反电动势)而倾向于抵抗电流的改变。磁通改变量对电流改变量的比值称为自感,自感通常也就直接称作是这个电路的电感

自感的计算公式为:$U=-L\frac{dI}{dt}$,U是自感电动势,I是电流,负号表示自感电动势反抗原来的电流。L是比例系数,就称为电感,对于同一个线圈来说,L是常数,单位是$V\cdot t//A=\Omega \cdot t$,同时也简记为$H$(亨利)。

点击阅读全文...

30 Apr

蘑菇的最优形状模型

淡白口蘑

淡白口蘑

达尔文的进化学说告诉我们,自然界总是在众多的生物中挑出最能够适应环境的物种,赋予它们更高的生存几率,久而久之,这些物种经过亿万年的“优胜劣汰”,进化成了今天的千奇百怪的生物。无疑,经过长期的选择,优良的形状会被累积下来,换句话讲,这些物种在某些环境适应能力方面已经达到最优或近乎最优的状态(又是一个极值问题了)。好,现在我们来考虑蘑菇。

蘑菇是一种真菌生物,一般生长在阴暗潮湿的环境中。喜欢湿润的它自然也不希望散失掉过多的水分,因此,它努力地调整自身的形状,使它的“失水”尽可能地少。假设单位面积的蘑菇的失水速度是一致的,那么问题就变成了使一个给定体积的立体表面积尽可能少的问题了。并且考虑到水平各向同性生长的问题,理想的蘑菇形状应该就是一个平面图形的旋转体。那么这个旋转体是什么呢?聪明的你是否想到了是一个球体(的一部分)呢?

点击阅读全文...