科学空间|Scientific Spaces

这两个程序都是为了天体力学服务的...

开普勒方程求根器-界面

点击阅读全文...

感谢Awank-Newton读者的来信，本文于2013.01.30作了修正，主要是弹性势能的正负号问题。之前连续犯了两个错误，导致得出了正确答案。现在已经修正。参考《平衡态公理的修正与思考》

在下面的两篇文章中，BoJone已经介绍了这个“旋转弹簧伸长”的问题，并从两个角度提供了两种解答方法。前者列出了一道积分方程，然后再转变为微分方程来解；后者直接从弹性力学的角度来列出一道二阶微分方程，两者殊途同归。
http://kexue.fm/archives/782/

http://kexue.fm/archives/826/

今天，再经过一段时间的变分法涉猎后，BoJone尝试从变分的角度（总能量最小）来给出一种新的解法。同样设r为旋转达到平衡后弹簧上一点到旋转中心的距离，该点的线密度为$\lambda =\lambda (r)$，该点到中心的弹簧质量为$m=m(r)$，旋转前的长度为$l_0$，旋转平衡后的长度为$l_1$。由于弹簧旋转后已经达到了平衡状态，由平衡态公理（参看《自然极值》系列），平衡意味着总能量“动能-势能”取极值。

点击阅读全文...

从2013.11.15开始，使用MathJax插件。主要原因是MathJax在兼容性方面比ASCIIMath Image Fallback Scripts做得好很多。而且从长远考虑，用MathJax也是应该的。
官方网站：http://www.mathjax.org/

复制数学公式：http://www.mathjax.org/demos/copy-and-paste/

-------以下内容已经过时（写于2013.01.05）--------

原来一直是使用“数学研发论坛”完善的数学公式插件来显示数学公式的，使用很简单，载入速度很快，这样一下子就用了三年了。

不过进入大学后，学习的东西越来越多，数学符号也越来越多，郭大哥的插件的不足也暴露出来了。最要命的是它居然无法显示$\hbar$，这叫我这个学习量子力学的孩子情何以堪...（不过郭大哥新版的插件已经加入了这个符号）。还有另外一个不足的地方，就是郭大哥的插件进行了大量的化简，使得数学公式的输入简单了不少，但是反而对标准的Latex代码支持不足了。久而久之，会带来一个弊端，就是迁移性不强。万一哪天这个插件无法使用了，就难以找到替代品了。考虑到这些，我写latex代码的时候总是用标准的语法而不用简化语法，后来$\hbar$的问题出来后，一下子用上了MathJax这个强大的插件（考虑过JsMath，但是发现它的行内公式显示效果不大好）。

点击阅读全文...

大概在半年前，我曾用“化圆法”解决了椭圆内定长弦中点轨迹问题，求出了轨迹方程。前几天，我收到了网名为“理想”的网友的Email，他提出了自己对这个问题的解法，并得到了形式不同的轨迹方程，因此对两者的等价性表示疑惑。经过检验，我跟他的轨迹方程基本上是等价的，不过，他求出的轨迹方程总包括了原点，这是一点不足之处。但是看起来，他的轨迹方程却感觉好看一些。这的确很让人意外，因为从他的化简过程来看，有种“化简为繁”的味道，却得出了相当简洁的答案，着实有趣。

经过网友的同意，将他的过程贴在这里与大家分享！后面附有pdf文档，欢迎下载阅读。希望在科学空间可以看到更多的读者留下的痕迹。

椭圆定长弦中点轨迹的一种解法

作者：理想

本文介绍了一种计算椭圆定长弦中点轨迹的方法。设椭圆长、短轴分别为$2a$、$2b$，弦长为$2r$，随着弦的两端在椭圆上滑动，弦的中点形成的轨迹为：
$$(\frac{x^2}{a^2} + \frac{y^2}{b^2} - 1)(\frac{x^2}{a^4} + \frac{y^2}{b^4} + \frac{r^2}{a^2b^2}) + \frac{r^2}{a^2b^2} = 0$$
它不是一个椭圆，而是一个高次曲线。

点击阅读全文...

昨天晚上一位网友与我讨论以下问题：

函数$y=\sqrt{3} x-\frac{1}{x}$的图像为双曲线，在此双曲线的两支上分别取P、Q点，求PQ的最短距离。

显然，如果双曲线是普通的$\frac{x^2}{a^2}-\frac{y^2}{b^2}=1$的形式，则这个问题是相当简单的。就是当y=0时两个点的距离，也就是2a。但是很明显这样的一条双曲线是经过旋转的。因此我们需要知道它究竟旋转了多少度$\theta$。然后列出$y=(\tan\theta) x$，联立双曲线方程就可以求出两个点了。

点击阅读全文...

2^29363731-1

很小的时候就开始对素数感兴趣了，后来是在一本《未解之谜》上看到了梅森素数、完全数、孪生素数等等东西，觉得甚是好玩。在初中买了计算机之后，就关注到了Prime 95这个梅森素数的分布式计算程序，以前也尝试过运行它，不过由于那时候计算机配置较低，一般都是运行到20%左右就没有坚持下去了。

上大学入手了一台四核的笔记本，就在去年10月份左右再次运行了这个程序，由于是四核，一次性可以同时测试四个数字。经过半年的运行，今天终于测试完了第一个数字：$2^{29363731}-1$。正如预料中的，这不是一个素数。不管怎样，它是我第一个完成的测试，也算是自己的一个独立的成果啦，呵呵，自娱自乐一番。

点击阅读全文...