28 Sep

开始学习数学软件Scilab

其实很早之前我就想学习一款数学软件的使用,以前很感兴趣的是mathematica,也玩弄过一阵子,但毕竟在高中没有多大需要,也就没有坚持下来。更重要的是,这些软件都是要收费的。上了大学后,听了师兄姐对数学建模的讲述,发现他们基本上也是用mathematica或者matlab的,但这两个软件都是要收费的,我不大想用破解版本。既然我都已经用上了ubuntu了,那么我就该好好利用它。据说命令跟matlab很相似的软件是scilab,还有octave,不同的是这些都是开源免费的。

出于熟悉代码操作和数学软件编程的目的,我选择了学习scilab。虽然网上说octave与matlab的相似程度更高,但是我感觉scilab比octave用的更广一些,所以就用它。所谓“一理通百理明”,先专心学好一个。

下面是我编写的第一个scialb程序,利用威尔逊方法来进行素性测试。这个代码的主要目的是练习条件语句和循环语句,以及一些输出输入的技巧而已。程序本身比较丑陋。

//我的第一个scilab程序
//完成于2012.09.27

label1=['p:';];  //定义标签
B=x_mdialog(['本程序使用威尔逊方法判断进行素数测试。';'请输入要判断的数'],label1,['127';]);  //输入框
p=evstr(B(1));  //提取输入框里边的数字进行赋值
i=1;
j=1;
q=p-1;
while i<q
    j=j*i;
    j=modulo(j,p);//这个是模函数。
    i=i+1;
end
if j==1
    messagebox(['这是一个素数';],['测试结果']);  //输出,其中后边的“测试结果”是输入框的标题
else
    messagebox(['这是一个合数';],['测试结果']);
end

点击阅读全文...

30 Sep

中秋节快乐!

中秋2012

中秋2012

明月几时有,把酒问青天,不知相聚时刻,它夕是何年?我欲乘风归去,又带月饼佳酿,唯恐不尽兴。起舞弄月影,胜似在天宫!

转楼台,低月映,照无眠。不应有愁,何事月在秋心圆。人有萍聚萍散,月有半圆半缺,此事古难求。但愿人长久,千里共婵娟。

——BoJone修改版《水调歌头》

我的好朋友,祝你中秋快乐,快乐一生,事业爱情如意!让科学空间与你我情谊地久天长!^_^

11 Oct

2012诺贝尔奖...

又是一年诺奖公布时......每年的这个时候,诺贝尔奖又会被热门地提及到,现在三个自然科学方面的奖项都已经公开了。简略收集如下:

诺贝尔生理学或医学奖
京都大学物质-细胞统合系统据点iPS细胞研究中心主任长山中伸弥(Shinya Yamanaka)、英国发育生物学家约翰-戈登因(John B. Gurdon)

2012诺贝尔生理学或医学奖

2012诺贝尔生理学或医学奖

原因:在细胞核重新编程研究领域的杰出贡献而获奖。所谓细胞核重编程即将成年体细胞重新诱导回早期干细胞状态,以用于形成各种类型的细胞,应用于临床医学。细胞核重编程指细胞内的基因表达由一种类型变成另一种类型。通过这一技术,可在同一个体上将较容易获得的细胞(如皮肤细胞)类型转变成另一种较难获得的细胞类型(如脑细胞)。这一技术的实现将能避免异体移植产生的排异反应。

点击阅读全文...

18 Oct

证明光速不变的一个理想实验??

在狭义相对论发表之前和之后,都有不少实验从不同角度论证了它的正确性。这些实验大多数是实际测量得出结果的,当然也存在着一些“理想实验”,这些实验只需要一定的逻辑推理,而实际上是无法完成的。下面就是我很久之前在某本书(很抱歉,我真的忘记书名了)看到的一个用来推翻光速可叠加的伽利略变换的理想实验。它只用寥寥几句,就好像已经证明了“c+c=c”(c是真空中的光速)的事实。可是“c+c=c”在狭义相对论上是作为原理出现的,是不可能通过逻辑推理来证明的。事实究竟如何?我们先来看这个实验。

光速不变的理想实验

光速不变的理想实验

任意选定一个坐标原点。设想原点的正北方$c\cdot t_0$处有一架以光速$c$朝南运行的飞机1;原点的正西方$c\cdot t_0$处有一架以光速$c$朝东运行的飞机2。假设就这样匀速运动着,显然,$t_0$时间后,将会发生惨剧(飞机相撞)。

点击阅读全文...

4 Dec

[转载]复杂的机械,简单的原理

在网上收集到的一些GIF,解释了不少机器构造的原理,其中有一些甚是巧妙!比如下面的椭圆规

07椭圆规

07椭圆规

也许以后学习或者科普的时候会用到,所以我将其打包为zip,收藏在这里,并与大家分享。里面包括一些武器和发动机的构造图,还有一些巧妙的部件,一些图片我已经注明该构造的名称,没有注明的是很显然的或者是找不到名称的。

点击阅读全文...

18 Dec

黑洞融合的简单模拟

在天文爱好者眼中,黑洞是一个球体,其半径为$\frac{2GM}{c^2}$;这是广义相对论的施瓦兹黑洞的结果,也从经典力学推导推导出来,虽然用经典力学是错误的,但是对于多数的天文爱好者(包括笔者)来说,这是目前唯一的一种可行的理解方法(广义相对论那些复杂推导会让我们很崩溃的)。当然,事实上,黑洞不是一个球体,它只是一个密度很大的点。至于密度有多大,目前公认的说法是无穷大,但是严格的物理是不接受这个说法的,或者说,物理是不会接受任何无穷大的说法,所以现在积极发展量子引力理论来统一相对论和量子力学,不过这是另话了。$\frac{2GM}{c^2}$只不过是黑洞的视界,视界之内,我们就什么也不知道了。本文主要就从经典力学的角度探讨一下两个黑洞的合并过程中其视界的变化。读者将会发现,这些视界的形状相当有趣。

经典力学中的黑洞是这样定义的:天体表面的逃逸速度超过了光速,于是连光都无法逃脱,所以这个“洞”就很黑。也就是说,光子的总能量(引力势能与动能之和,经典力学意义下的)要为负,负数表示受到束缚。用数学公式来讲,就是:

$$\frac{1}{2}mc^2 - \frac{GM_1 m}{r_1}-\frac{GM_2 m}{r_2}-...-\frac{GM_n m}{r_n} \leq 0$$

点击阅读全文...

27 Dec

费曼路径积分思想的发展(三)

3、费曼图和量子电动力学的重整化

在1947年美国避难岛(Shelter Island)会议上,兰姆报导了他的重大发现,即现今所称的兰姆位移;氢原子的$2S_{\frac{1}{2}}$能级比$2P_{\frac{1}{2}}$高出约1000MHz。而按照狄拉克理论,对纯库仑相互作用的电子-质子系统,这两个能级应该是简并的。人们很快就认识到,该位移应归之于一阶近似的辐射校正[19]。贝特用一个电子的校正质量就非相对论近似得出了氢原子nS能级的位移公:

$$\frac{8}{3\pi}(\frac{e^2}{\hbar c})Ry \frac{Z^4}{n^3} Ln\frac{K}{ < E_n-E_m > _{AV}}$$

点击阅读全文...

20 Jan

评论功能修复了

一直都挺奇怪为什么这些天Blog都没有人来评论,难道我的文章质量下降了?可是统计显示访问人数比以往有所增加呀。一直郁闷但却无从找到答案。

今天收到了4anpan读者的邮件,表示评论解锁功能失效了。我自己测试才发现真的如此!!原来不一定是没有人来评论,而是评论插件出了问题。经排除,是最新的公式插件(ASCIIMathMLwFallback2)与滑动解锁的评论插件(IQapTcha)有冲突。所以只好把评论插件换回旧的验证码插件了。正在构思新一步的插件方案。

万分感谢4anpan读者的反馈,希望有更多的读者来帮忙完善科学空间。这里有你更精彩!!

期待你的声音^_^