m*******s
发帖数: 3142 | 1 【 以下文字转载自 Mathematics 讨论区 】
发信人: manifolds (流形), 信区: Mathematics
标 题: 近似计算matrix exponential的代码？
发信站: BBS 未名空间站 (Wed Aug 10 23:44:38 2011, 美东)
给定一个non-Hermitian matrix A, exp(A)v可以用Arnoldi algorithm很好的近似，不
知道大家可否推荐一些比较可靠的fortran code。
我google了一下，有几个现成的matlab script file ，不过似乎没有找到什么fortran
的。
谢谢！ |
|
发帖数: 1 | 2 物理都是近似计算
不会有人真以为这个世界是理想化的吧 |
|
g****t
发帖数: 31659 | 3 我想到一个办法.你看行不行.
x*exp(-x^2*a)的Fourier变换不需要广义函数就可以搞定.
那么如果把a取的很小,例如a=0.00001,这就近似于x的Fourier变换.
这个过程中,Dirac函数被1/a^2类似的东西近似了.
在你的问题中,可以把B(x)展开成正交多项式,然后求B(x)*exp(-x^2*0.00000001)
的Fourier变换. 计算中出现的x^n*exp(-a*x^2)之类的东西可以查表搞定.
(
如果不乘exp(-x^2*a),无论按什么多项式展开,总有Dirac函数出现.
hermite,Laguerre的第一项都是1,第二项都是a+bx形式的,这样,如果用Dirac函数,
最后很可能因为数值误差出现消不掉的Dirac函数---虽然原理上来讲到最后这些
剩下的Dirac函数的系数应该非常小.
)
B(x)如果是指数函数的线性组合,那有现成的方法做.
function |
|
m*******s
发帖数: 3142 | 4 我现在遇到一个数值计算问题，大概如下
一个形如x-A(x)矩阵的逆矩阵，左乘右乘一个不同的常数矩阵，得到一个矩阵B，所有
这都是数值计算
得到的。
当然矩阵B也可以视为以x为自变量的算符。
然后我需要求这个矩阵B(x)的Fourier变换，
\int_{-\infty}^ {\infty}dx \boldsymbol{B} \left (x\right )e^{-ixt}
对每个t，都要重复做相同的数值积分程序，比较耗时。
我想问的是，有没有可能用一种形如
f_1(x)B_1+f_2(x)B_2+f_3(x)B_3+\cdots
的式子来拟合这个矩阵B(x)，
其中f_1,f_2,f_3.....都是x的常见函数,B_1,B_2,B_3....是常数矩阵。
使得可以使用复变函数的residue theorem来方便地求得Fourier变换?
难就难在如何确定这些f_1,f_2,f_3,B_1,B_2,B_3....
似乎常规的least square方法派不上用场。
我感觉可以用B(x)的trace来找f_1,f_2,f_3.....但是这样也定不了B_1,B_2,B_3..
有没有更... 阅读全帖 |
|
g****t
发帖数: 31659 | 5 B(x)每个元素都展开成chebyshev 多项式,
然后 Chebyshev多项式的Fourier变换应该有表可以查?
我现在遇到一个数值计算问题，大概如下
一个形如x-A(x)矩阵的逆矩阵，左乘右乘一个不同的常数矩阵，得到一个矩阵B，所有
这都是数值计算
得到的。
当然矩阵B也可以视为以x为自变量的算符。
然后我需要求这个矩阵B(x)的Fourier变换，
\int_{-\infty}^ {\infty}dx \boldsymbol{B} \left (x\right )e^{-ixt}
对每个t，都要重复做相同的数值积分程序，比较耗时。
我想问的是，有没有可能用一种形如
f_1(x)B_1+f_2(x)B_2+f_3(x)B_3+\cdots
的式子来拟合这个矩阵B(x)，
其中f_1,f_2,f_3.....都是x的常见函数,B_1,B_2,B_3....是常数矩阵。
使得可以使用复变函数的residue theorem来方便地求得Fourier变换?
难就难在如何确定这些f_1,f_2,f_3,B_1,B_2,B_3....
似乎常规的least square方法派不上用场... 阅读全帖 |
|
B*Q
发帖数: 25729 | 6 物理搞不明白时
近似一下
忽悠一下
一门学科就出来了 |
|
发帖数: 1 | 7 那为啥还要统计专业
：物理搞不明白时
：近似一下 |
|
n****4
发帖数: 12553 | 8 所有的物理，包括牛逼的爱斯坦因物理，都是描述宇宙的某种近似 |
|
m*******s
发帖数: 3142 | 9 给定一个non-Hermitian matrix A, exp(A)v可以用Arnoldi algorithm很好的近似，不
知道大家可否推荐一些比较可靠的fortran code。
我google了一下，有几个现成的matlab script file ，不过似乎没有找到什么fortran
的。
谢谢！ |
|
q********g
发帖数: 10694 | 10 Sobereva
Department of Chemistry, University of Science and Technology Beijing,
Beijing 100083, China
前言：本文主要介绍过渡态、反应路径的计算方法，并讨论相关问题。由于这类算法极
多，可以互相组合，限于精力不可能面面俱到展开，所以只介绍常用，或者实用价值有
限但有启发性的方法。文中图片来自相关文献，做了一定修改。由于本文作为帖子发布
，文中无法插入复杂公式，故文中尽量将公式转化为文字描述并加以解释，这样必然不
如公式形式严谨，而且过于复杂的公式只能略过，但我想这样做的好处是更易把握方法
的梗概，有兴趣可以进一步阅读原文了解细节。对于Gaussian中可以实现的方法，文中
对其在Gaussian中的使用进行了一些讨论，希望能纠正一些网上流传的误区。虽然绝大
多数人不专门研究计算方法，其中很多方法也不会用到，但多了解一下对开阔思路是很
有好处的。
文中指的“反应”包括构象变化、异构化、单分子反应等任何涉及到过渡态的变化过程
。“反应物”与“产物”泛指这些过程的初态和末态。“优化”若未注明，... 阅读全帖 |
|
l***a
发帖数: 198 | 11 原文：
http://blog.sciencenet.cn/blog-268546-1154610.html
量子计算的理论基础
量子计算的理论基础是量子力学，相比传统计算，它声称的巨大性能提升来自于量子态
的叠加性。
我发现很难跟量子计算（量子信息）的研究者解释清楚，他们说的量子力学实际上是非
相对论量子力学，而非相对论量子力学只是一套近似理论，并不代表真实的物理。在量
子信息研究者的眼中，只有一套神奇的理论叫做量子力学，谁也说不清楚怎么回事，但
是它就是对的。谁反对都没有用，爱因斯坦也不行。无论爱因斯坦有多么伟大，在与量
子力学的交锋中全都败下阵来。
其中最重要的一个争论，爱因斯坦认为世界是局域实在的，局域意味着满足相对论，不
能超光速，实在意味着客观世界独立于人的意识而存在。而量子力学认为，世界（量子
）是非局域的，量子现象可以超光速；世界（量子）在测量之前讨论它是不是客观实在
是没有意义的；量子的非局域性已经得到了众多实验的严格证明。（这里说的量子力学
观点存在重大争议，量子非局域性的实验证明，也没有得到广泛的认可。）
在非相对论量子力学中，作用势是唯象的，作用量传播速度无穷大... 阅读全帖 |
|
b****a
发帖数: 4465 | 12 为了落实第二条原则，分析总结了力学所遇到阻碍的技术原因。由于涡轮轴形状奇特，
与经常见到的形状匀称规则的各种轴不同，涡轮轴直径沿轴线不断增大，成为喇叭形状
，而且喇叭口附近还有用于气体密封的几圈尺寸严格的环形台阶，断裂往往就从台阶根
底发生。
李敏华教授试图寻找应力分布的解析解，对于这样奇形怪状的轴，形状曲线都不连续，
求解当然很难。当时在计算结构力学领域刚刚诞生了有限元方法。不去寻求解析解，而
是通过网格迭代近似计算的方法进行数值计算。有限元法将连续的求解域离散为一组单
元的组合体，用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数
，近似函数通常由未知场函数及其导数在单元各节点的数值插值函数来表达。
从而使一个连续的无限自由度问题变成离散的有限自由度问题。举一个简单的比喻，解
析计算是要通过解微分方程求出函数精确解，有限元计算则是弹性力学的一种近似解法
，它将连续体变换为离散化结构，然后再利用分片插值技术与虚功原理或变分方法寻求
近似数值解。网格足够多、足够细小，数值解的精度也就越高。网格的划分设定，能够
达到工程需要的计算精度就行了。那个时候，有限元领域只有个别院校... 阅读全帖 |
|
s********k
发帖数: 6180 | 13 【 以下文字转载自 Stock 讨论区 】
发信人: silverhawk (silverhawk), 信区: Stock
标 题: 请教有没有可以实时计算时间学列correlation的近似算法？ (转载)
发信站: BBS 未名空间站 (Fri May 21 18:39:21 2010, 美东)
发信人: silverhawk (silverhawk), 信区: EE
标 题: 请教有没有可以实时计算时间学列correlation的近似算法？
发信站: BBS 未名空间站 (Wed May 19 20:42:59 2010, 美东)
假设我有两个时间序列x(i),y(i),i=1:n. 一般的correlation计算要等到全部n个采样
完成之后再计算。我现在想用一种online的方法来实时计算，从i=2开始就开始计算两
个序列的相关性，每次时间序列有一个新的采样correlation更新一次，当然由于采样
不完全，所以允许correlation的误差存在，但是希望大体走势应该是越来越趋近最后
的准确值，不知道有没有这样的实时算法存在。谢谢 |
|
发帖数: 1 | 14 日前，美国总统科学技术咨询委员会（President’s Council of Advisors on
Science and Technology，简称PCAST）发表了名为《确保美国半导体的领导地位》的
报告。在报告中提到，中国的半导体的崛起，对美国已经构成了“威胁”，委员会建议
政府对中国产业加以限制。
在详细介绍半导体产业的全球现状之前，委员会在报告里面写了一封给美国总统奥
巴马的信，简要翻译如下：
敬爱的奥巴马总统，
这个报告是由工作组内的相关行业领袖、杰出研究员和曾经的政策制定者拟定
的，PSAST已经做了审核并通过。报告主要是谈及半导体产业在创新、竞争和安全方面
面临的挑战和机会。
半导体是现代生活的重要组成部分，在半导体领域取得的进展已经将机基于其
打造的设备和服务提升到一个新的阶段。与此同时，还开拓了很多新的业务模式和产业
，为美国相关从业人员与消费者带来了巨大的收益，对促进全球经济的发展也产生了重
大的影响。我们也要明白到，尖端的半导体技术对于美国的国防系统和军队实力来说，
也是重要的保证。无处不在的半导体使得我们还同时面临了网络安全... 阅读全帖 |
|
|
f********t
发帖数: 6999 | 16 【 以下文字转载自 JobHunting 讨论区 】
发信人: mudhoof (正在长牙的羊), 信区: JobHunting
标 题: 这么热闹, 我也报Google offer
发信站: BBS 未名空间站 (Tue Feb 23 12:32:47 2010, 美东)
今天刚刚通知的, 特别感谢一起讨论的krone, geniusxsy, hnm, 特别是blaze教了我很
多, 还要特别感谢mitbbs59的总结帖
一起报offer, 好事成三, 大吉大利, 包子分光为止
贴下我的复习材料
题目大全:
http://www.spellscroll.com/viewquestions/?tag=algorithm
http://www.thecareerplus.com/?page=resources&cat=10
http://interviewcyclopedia.blogspot.com/
http://www.doctorinterview.com/A.html
http://toptechnotes.blogspot.com/search/label/algorith... 阅读全帖 |
|
x******i
发帖数: 3022 | 17
经典和量子电动力学的共同点在于，作用量是j*A，j是电流，A是电磁势。
在经典电动力学，这个项导致了电子的加速运动放出电磁波，以及电磁场
对电子的洛伦兹力。电磁波可以用Lineard-Wiechert势来计算，而这个势
可以简单的用库仑场加洛伦兹变幻得到。
在量子场论，这个项相应的过程是电子动量的改变放出光子，或者入射光子
被电子吸收改变电子的动量。可以用有一个结点的费曼图表示。虽然这个概
念貌似是从辐射场出发的，这只是因为数学处理的时候在动量空间比较简单。
非相对论电子受到原子核的库伦场，也可以用电子和原子核之间交换光子来
得到。
如果这个费曼图表达了量子电动力学的本质，那“加速电子放出光子”就是个
很本质的过程。
经典和量子的区别在于电子的波动性，还在于不确定原理。电子的波动性导
致了原子能级的分立。不确定原理导致了电磁场的零点振动。
如果忽略电磁场的量子涨落，则电子有多个分立的能量本征态，其波函数是
一些驻波。在数学上，这相当于哈密顿量里面的j*A变成j*Ac，把算符矢量A
变成复数矢量Ac。这个哈密顿量是《量子力学》的出发点。由于Ac是经典电
磁场，它已经没法被电子的量子化的... 阅读全帖 |
|
发帖数: 1 | 18 GPS需要相对论是谣言，这种谣言一到具体，很容易不攻自破
GPS 不是广义相对论的实验证据
王令隽
有人宣称GPS全球定位系统必须进行相对论修正，否则就不干活。许多人也信了。
最近和孙南生教授通信，他推荐了一篇业内专家R Li 的不同意见：
作者：R Li
链接：http://www.zhihu.com/question/24796597/answer/45911638
来源：知乎
摘要如下：
GPS不是相对论的应用
相对论改正只是卫星钟误差中的一个影响很小的系统误差改正项。
GPS系统有两个用途
主要用途， to drop 5 bombs in the same hole – 定位用途（这也可以看出设计这
个系统时军方目标精度是米级）相对论影响对目标精度的影响可忽略不计。
–主要用途中衍生出来的高精度测绘级定位用差分观测法消掉卫星钟误差，不需考虑相
对论影响（直接被绕过了）。
GPS定位是从未知点观测观测已知点从而确定未知点的坐标方位（专业说法叫resection
后方交会），观测量是未知点，你手持的GPS到已知坐 标的卫星的距离（卫星发射基于
卫星钟时钟信号，其实是有长周... 阅读全帖 |
|
j***h
发帖数: 4412 | 19 见附图GB50011-2001《建筑抗震设计规范》第255页的条文说明之第6.2.4、6.2.5、6.2
.8条对89年规范“强剪弱弯”近似计算公式的简化说明，89年规范的近似计算公式中已
经考虑了轴压比。
GB50011-2001对89年规范的近似计算公式的简化是通过数百个算例的统计分析而得出的
。因此1989年的建筑抗震设计规范生效后所建设的楼房，如果是正确设计及施工的，均
可满足89年规范已经提出的“强剪弱弯”要求，框架的抗侧力（水平地震力）构件的抗
剪能力应可大于其抗弯能力， 在大震（罕遇地震）下应发生延性的（有一定延时的）
弯曲破坏，而不是发生脆性的（突然的）剪切破坏。
箍筋不足从而造成剪切破坏是可能的。还有可能是轴向力的作用，或者轴向力和弯矩的
共同作用。
1990年的国标我不懂，但是无论是剪切还是轴向力的问题，大多都是Northridge以后研
究出来的。Northridge中的房屋和桥梁很多都是由于箍筋不足而破坏，所以后来的规范
有很大修改。以前的规范只要求箍筋间距12feet，后来发现远远不够。后来的巨大项目
retrofit都是为了弥补箍筋的不足。 |
|
s********k
发帖数: 6180 | 20 【 以下文字转载自 EE 讨论区 】
发信人: silverhawk (silverhawk), 信区: EE
标 题: 请教有没有可以实时计算时间学列correlation的近似算法？
发信站: BBS 未名空间站 (Wed May 19 20:42:59 2010, 美东)
假设我有两个时间序列x(i),y(i),i=1:n. 一般的correlation计算要等到全部n个采样
完成之后再计算。我现在想用一种online的方法来实时计算，从i=2开始就开始计算两
个序列的相关性，每次时间序列有一个新的采样correlation更新一次，当然由于采样
不完全，所以允许correlation的误差存在，但是希望大体走势应该是越来越趋近最后
的准确值，不知道有没有这样的实时算法存在。谢谢 |
|
r******r
发帖数: 700 | 21 海量数据处理：十道面试题与十个海量数据处理方法总结
作者：July、youwang、yanxionglu。
时间：二零一一年三月二十六日
说明：本文分为俩部分，第一部分为10道海量数据处理的面试题，第二部分为10个海量
数据处理的方法总结。
本文之总结：教你如何迅速秒杀掉：99%的海量数据处理面试题。有任何问题，欢迎随
时交流、指正。
出处：http://blog.csdn.net/v_JULY_v。
------------------------------------------
第一部分、十道海量数据处理面试题
1、海量日志数据，提取出某日访问百度次数最多的那个IP。
首先是这一天，并且是访问百度的日志中的IP取出来，逐个写入到一个大文件中
。注意到IP是32位的，最多有个2^32个IP。同样可以采用映射的方法，比如模1000，把
整个大文件映射为1000个小文件，再找出每个小文中出现频率最大的IP（可以采用hash
_map进行频率统计，然后再找出频率最大的几个）及相应的频率。然后再在这1000个最
大的IP中，找出那个频率最大的IP，即为所求。
或者如下阐述（雪... 阅读全帖 |
|
r******r
发帖数: 700 | 22 海量数据处理：十道面试题与十个海量数据处理方法总结
作者：July、youwang、yanxionglu。
时间：二零一一年三月二十六日
说明：本文分为俩部分，第一部分为10道海量数据处理的面试题，第二部分为10个海量
数据处理的方法总结。
本文之总结：教你如何迅速秒杀掉：99%的海量数据处理面试题。有任何问题，欢迎随
时交流、指正。
出处：http://blog.csdn.net/v_JULY_v。
------------------------------------------
第一部分、十道海量数据处理面试题
1、海量日志数据，提取出某日访问百度次数最多的那个IP。
首先是这一天，并且是访问百度的日志中的IP取出来，逐个写入到一个大文件中
。注意到IP是32位的，最多有个2^32个IP。同样可以采用映射的方法，比如模1000，把
整个大文件映射为1000个小文件，再找出每个小文中出现频率最大的IP（可以采用hash
_map进行频率统计，然后再找出频率最大的几个）及相应的频率。然后再在这1000个最
大的IP中，找出那个频率最大的IP，即为所求。
或者如下阐述（雪... 阅读全帖 |
|
s**********o
发帖数: 14359 | 23 【 以下文字转载自 JobHunting 讨论区 】
发信人: rongxuer (蓉儿), 信区: JobHunting
标 题: 如何秒杀99%的海量数据处理面试题
发信站: BBS 未名空间站 (Thu Apr 5 02:08:57 2012, 美东)
海量数据处理：十道面试题与十个海量数据处理方法总结
作者：July、youwang、yanxionglu。
时间：二零一一年三月二十六日
说明：本文分为俩部分，第一部分为10道海量数据处理的面试题，第二部分为10个海量
数据处理的方法总结。
本文之总结：教你如何迅速秒杀掉：99%的海量数据处理面试题。有任何问题，欢迎随
时交流、指正。
出处：http://blog.csdn.net/v_JULY_v。
------------------------------------------
第一部分、十道海量数据处理面试题
1、海量日志数据，提取出某日访问百度次数最多的那个IP。
首先是这一天，并且是访问百度的日志中的IP取出来，逐个写入到一个大文件中
。注意到IP是32位的，最多有个2^32个IP。同样可以采用映射的... 阅读全帖 |
|
l*****n
发帖数: 608 | 24 如果你只有两个序列的话, 假设他们是对齐的 (No time delay between them), 你可
以用一个moving window 去计算他们每一段的hamming distance, 那么你会有足够多的
样本去判断. 但是我觉得您是否需要先判断你关心数据的哪一个方面? 这两个数据之间
的关系是什么?
比如, 要比较序列的趋势,同时数据的绝对大小并不重要的话, 你可以计算他们的
Kendall tau distance,
http://en.wikipedia.org/wiki/Kendall_tau_distance
如果是模拟出来的数据, 要比较是否完全重复了实验的数据,也可以使用Granger
Causality 比较他们两者之间的causality (实际也是一种delay embadding做
regression的方法) 计算他们的distance. 或者用time embadding的方法, 将时间序列
投影到N维空间内, 计算他们对应点的距离.
但是, 如果数据是独立的,那么你可能只需要判断他们的统计参数, 比如傅立叶变换的
频谱, 关联长度等等, 来决定... 阅读全帖 |
|
n********g
发帖数: 6504 | 25 Original AI前线小组 AI前线 今天
作者｜Mikhail Dyakonov
译者｜李志
编辑｜Debra
AI 前线导读：爱好科幻小说的你，一定对量子计算并不陌生。而如今，备受瞩目的量
子计算到底进展如何呢？本文将带你走进当今的量子计算世界，听听理论物理学界大师
Mikhail Dyakonov 为什么不相信量子计算。
更多干货内容请关注微信公众号“AI 前线”（ID：ai-front）
插画作者：Christian Gralingen
当今，量子计算备受瞩目，风靡一时。每天，各种新闻媒体都争相报道这项高大上的技
术能为人类社会带来的神奇改变。然而，大多数评论家都忘记了——或者说有意掩盖了
这样一个事实：人们几十年来一直殚精竭虑地进行着量子计算研究，却到现在仍没有什
么实质性的应用展示出来。
有人说（相关文章：https://www.research.ibm.com/ibm-q/learn/what-is-quantum-
computing/）， 量子计算机可以“在许多学科中提供突破，包括材料和药物发现，复
杂人造系统的优化和人工智能”。还有人说（相关文章：https... 阅读全帖 |
|
h*j
发帖数: 393 | 26 不透明的“透明计算”——张尧学是如何骗取国家自然科学一等奖的
·方舟子·
教育部原司长、中南大学现校长、中国工程院院士张尧学领衔的“透明计算”
课题获得国家自然科学一等奖，舆论大哗。多名IT专业人士在网上一面倒地揭露
这个课题低劣、造假、骗人的实质，这些揭露文章发表在科学网、知乎、新浪博
客等处，随后即被删除得干干净净。财新网在1月17日发了一篇报道《“透明计
算”获大奖引发诸多争议》，随后也被删除。对这些批评、揭露，张尧学没有给
出任何回应，而这些试图压制批评的删帖行为，不知是张尧学公关的结果，还是
政府部门为遮羞采取的行动。不管怎样，对学术争议只能以删帖来回应，其实质
如何，已昭然若揭了。
我不是IT专业人士，不想从学术的角度来评价张尧学的成果。我想从别的角
度，让公众了解一个被IT专业人士一致评为低劣、造假的课题何以能够获得国家
自然科学一等奖。
我们先来看看张尧学的简历：
1976年10月至1978年10月，在湖南澧县余家台水电站工作，测绘员；
1978年10月至1982年07月，在陕西西安电子科技大学（原陕西西北电讯工... 阅读全帖 |
|
m********3
发帖数: 3280 | 27 教育部原司长、中南大学现校长、中国工程院院士张尧学领衔的“透明计算”课题获得
国家自然科学一等奖，舆论大哗。多名IT专业人士在网上一面倒地揭露这个课题低劣、
造假、骗人的实质，这些揭露文章发表在科学网、知乎、新浪博客等处，随后即被删除
得干干净净。财新网在1月17日发了一篇报道《“透明计算”获大奖引发诸多争议》，
随后也被删除。对这些批评、揭露，张尧学没有给出任何回应，而这些试图压制批评的
删帖行为，不知是张尧学公关的结果，还是政府部门为遮羞采取的行动。不管怎样，对
学术争议只能以删帖来回应，其实质如何，已昭然若揭了。
我不是IT专业人士，不想从学术的角度来评价张尧学的成果。我想从别的角度，让
公众了解一个被IT专业人士一致评为低劣、造假的课题何以能够获得国家自然科学一等
奖。
我们先来看看张尧学的简历：
1976年10月至1978年10月，在湖南澧县余家台水电站工作，测绘员；
1978年10月至1982年07月，在陕西西安电子科技大学（原陕西西北电讯工程学院）
电子工程系学习，获学士学位；
1982年07月至1983年02月，在西安电子科技大学出国研究生预备班学习；
1983年02月... 阅读全帖 |
|
S******8
发帖数: 24594 | 28 不透明的“透明计算”——张尧学是如何骗取国家自然科学一等奖的
教育部原司长、中南大学现校长、中国工程院院士张尧学领衔的“透明计算”课题
获得国家自然科学一等奖，舆论大哗。多名IT专业人士在网上一面倒地揭露这个课题低
劣、造假、骗人的实质，这些揭露文章发表在科学网、知乎、新浪博客等处，随后即被
删除得干干净净。财新网在1月17日发了一篇报道《“透明计算”获大奖引发诸多争议
》，随后也被删除。对这些批评、揭露，张尧学没有给出任何回应，而这些试图压制批
评的删帖行为，不知是张尧学公关的结果，还是政府部门为遮羞采取的行动。不管怎样
，对学术争议只能以删帖来回应，其实质如何，已昭然若揭了。
我不是IT专业人士，不想从学术的角度来评价张尧学的成果。我想从别的角度，让
公众了解一个被IT专业人士一致评为低劣、造假的课题何以能够获得国家自然科学一等
奖。
我们先来看看张尧学的简历：
1976年10月至1978年10月，在湖南澧县余家台水电站工作，测绘员；
1978年10月至1982年07月，在陕西西安电子科技大学（原陕西西北电讯工程学院）
电子工程系学习，获学士学位；
1982年07月至1983年02... 阅读全帖 |
|
发帖数: 1 | 29 关于水星近日点的移动问题
王令隽 2014年8月
水星近日点的移动问题是一个具有历史重要性的问题。其重要性不是因为
这一现象本身对天体物理有多么重要，而是因为相对论学者们把这问题当成支持相对论
反对经典力学的重要实验证据，被认为是“广义相对论的三大实验检验”之一。
所谓水星的近日点，就是水星轨道离太阳最近的那一点。所有的行星轨道都是椭圆，只
不过有的椭圆偏心率很小，几乎是个圆。地球轨道就非常接近一个圆。椭圆有两个焦点
，太阳就落在其中一个焦点上，所以椭圆轨道的长轴的两端一个离太阳最近，叫近日点
，一个离太阳最远，叫远日点。近日点和远日点的移动就是椭圆轨道长轴（或短轴）的
移动。如果把太阳和行星都看成为质点，并且在计算行星轨道时，忽略所有其他行星的
存在，那末行星轨道的长短轴是不会移动的。也就是说，近日点是不会移动的。可是这
种理想的条件不符合实际情况。因为第一，太阳本身不是一个质点，而是一个扁球体（
oblate sphere）；第二，太阳系里有九个行星和其他物质。它们的相互作用会影响到
任何一个行星的运动。所以行星的近日点都会移动，水星近日点的移... 阅读全帖 |
|
f*******2
发帖数: 717 | 30 不透明的“透明计算”——张尧学是如何骗取国家自然科学一等奖的
·方舟子·
教育部原司长、中南大学现校长、中国工程院院士张尧学领衔的“透明计算”
课题获得国家自然科学一等奖，舆论大哗。多名IT专业人士在网上一面倒地揭露
这个课题低劣、造假、骗人的实质，这些揭露文章发表在科学网、知乎、新浪博
客等处，随后即被删除得干干净净。财新网在1月17日发了一篇报道《“透明计
算”获大奖引发诸多争议》，随后也被删除。对这些批评、揭露，张尧学没有给
出任何回应，而这些试图压制批评的删帖行为，不知是张尧学公关的结果，还是
政府部门为遮羞采取的行动。不管怎样，对学术争议只能以删帖来回应，其实质
如何，已昭然若揭了。
我不是IT专业人士，不想从学术的角度来评价张尧学的成果。我想从别的角
度，让公众了解一个被IT专业人士一致评为低劣、造假的课题何以能够获得国家
自然科学一等奖。
我们先来看看张尧学的简历：
1976年10月至1978年10月，在湖南澧县余家台水电站工作，测绘员；
1978年10月至1982年07月，在陕西西安电子科技大学（原陕西西北电讯工程
学院... 阅读全帖 |
|
w*********a
发帖数: 9279 | 31 【 以下文字转载自 Military 讨论区 】
发信人: hpj (江湖), 信区: Military
标 题: 方舟子： 不透明的“透明计算?
发信站: BBS 未名空间站 (Tue Jan 20 11:37:16 2015, 美东)
不透明的“透明计算”——张尧学是如何骗取国家自然科学一等奖的
·方舟子·
教育部原司长、中南大学现校长、中国工程院院士张尧学领衔的“透明计算”
课题获得国家自然科学一等奖，舆论大哗。多名IT专业人士在网上一面倒地揭露
这个课题低劣、造假、骗人的实质，这些揭露文章发表在科学网、知乎、新浪博
客等处，随后即被删除得干干净净。财新网在1月17日发了一篇报道《“透明计
算”获大奖引发诸多争议》，随后也被删除。对这些批评、揭露，张尧学没有给
出任何回应，而这些试图压制批评的删帖行为，不知是张尧学公关的结果，还是
政府部门为遮羞采取的行动。不管怎样，对学术争议只能以删帖来回应，其实质
如何，已昭然若揭了。
我不是IT专业人士，不想从学术的角度来评价张尧学的成果。我想从别的角
度，让公众了解一个被IT专业... 阅读全帖 |
|
c*******n
发帖数: 1648 | 32 教育部原司长、中南大学现校长、中国工程院院士张尧学领衔的“透明计算”课题获得
国家自然科学一等奖，舆论大哗。多名IT专业人士在网上一面倒地揭露这个课题低劣、
造假、骗人的实质，这些揭露文章发表在科学网、知乎、新浪博客等处，随后即被删除
得干干净净。财新网在1月17日发了一篇报道《“透明计算”获大奖引发诸多争议》，
随后也被删除。对这些批评、揭露，张尧学没有给出任何回应，而这些试图压制批评的
删帖行为，不知是张尧学公关的结果，还是政府部门为遮羞采取的行动。不管怎样，对
学术争议只能以删帖来回应，其实质如何，已昭然若揭了。
我不是IT专业人士，不想从学术的角度来评价张尧学的成果。我想从别的角度，让
公众了解一个被IT专业人士一致评为低劣、造假的课题何以能够获得国家自然科学一等
奖。
我们先来看看张尧学的简历：
1976年10月至1978年10月，在湖南澧县余家台水电站工作，测绘员；
1978年10月至1982年07月，在陕西西安电子科技大学（原陕西西北电讯工程学院）
电子工程系学习，获学士学位；
1982年07月至1983年02月，在西安电子科技大学出国研究生预备班学习；
1983年02月... 阅读全帖 |
|
z*****3
发帖数: 1793 | 33 本人CS PHD，专攻方向reinforcement learning （RL），恰好AlphaGo用到很多用到很
多RL。看不下去一群文科生+下室索南YY，乱猜测。
第一个问题，AlphaGO涉及到机器学习中的supervised learning（SL）和RL。这是机器
学习下的两个分支。deep learning是目前兴起的另外一个分支。Deep Mind正是运用
deep learning的deep neural network作为基础架构来解决 prediction （涉及SL）和
action selection 问题（涉及RL）。
第二个问题，AlphaGO 基础构造
1）利用SL学习了一个policy network。这里解释一下policy，policy是RL，control
thoery中的一个重要概念，简单来讲就是当遇到各种环境反馈时，应该如何选择下一步
的action。这里主要是利用SL学习一个policy neural network。
2）利用RL中的policy gradient算法，优化上面学到的policy network
3）利用deep neura... 阅读全帖 |
|
z*****3
发帖数: 1793 | 34 【 以下文字转载自 Military 讨论区 】
发信人: zyszys3 (zyszys3), 信区: Military
标 题: 专业科普AlphaGO
发信站: BBS 未名空间站 (Fri Mar 11 13:54:23 2016, 美东)
本人CS PHD，专攻方向reinforcement learning （RL），恰好AlphaGo用到很多用到很
多RL。看不下去一群文科生+下室索南YY，乱猜测。
第一个问题，AlphaGO涉及到机器学习中的supervised learning（SL）和RL。这是机器
学习下的两个分支。deep learning是目前兴起的另外一个分支。Deep Mind正是运用
deep learning的deep neural network作为基础架构来解决 prediction （涉及SL）和
action selection 问题（涉及RL）。
第二个问题，AlphaGO 基础构造
1）利用SL学习了一个policy network。这里解释一下policy，policy是RL，control
thoery中的一个重要概念，简单来讲就是当... 阅读全帖 |
|
e***d
发帖数: 8248 | 35 呵呵，高幂次的乘方开方，比低幂次的乘方开方，容易太多。
比如说，用泰勒级数的方法做高幂次，取到一阶近似，就能有很好的结果，
如果是一阶近似，就是解个一元一次方程的问题。
取到二阶近似，基本就能得到“超人”的答案了，归结为解个一元二次方程的问题。
不适用于低幂次的乘方开方，因为需要用到高阶近似，计算量太大。
同样的数字，给他做个平方或三次方，三倍的时间内能算出来是真牛。
后面那个2的乘方，我们平常人都能脱口而出。
当然心算能力还是不错的了，相对于我们平常人。
就是个电视秀，电视秀要的是收视率，理解这一点就可以了。 |
|
O****s
发帖数: 463 | 36 很赞同。科学计算中很多是用Fortran, C, 和C++的混合编程。
Fortran虽然不够灵活，但是数值计算的性能很好：
http://annwm.lbl.gov/bench/talk.ppt;
http://stackoverflow.com/questions/13078736/fortran-vs-c-does-f
因为对于物理中的大体系，是几乎不可能用量子力学去精确求解的，所以实际处理只能
用微扰近似去计算
（即便在漂亮的string theory里，想得到实际有用的信息，当前也多是用微扰去处理
，还没有较好的non-perturbative的理论）
以前做过几年的第一性原理分子动力学模拟计算,因为需要根据薛定谔方程去求解波函
数，用平面波展开去拟合做自洽计算，体系每个周期性结构原子数几十个左右还可以在
32节点每节点16GB内存的并行机上在几天内跑出结果，体系稍大点，上百个原子，就要
个把月了，，，而且很容易出问题，特别是IO是个很大的瓶颈，虽然总体数据量不大，
但计算量非常大，现在国内各大HPC中心不少任务就是跑这个的
类似的，在粒子物理中，Lattice ... 阅读全帖 |
|
O****s
发帖数: 463 | 37 很赞同。科学计算中很多是用Fortran, C, 和C++的混合编程。
Fortran虽然不够灵活，但是数值计算的性能很好：
http://annwm.lbl.gov/bench/talk.ppt;
http://stackoverflow.com/questions/13078736/fortran-vs-c-does-f
因为对于物理中的大体系，是几乎不可能用量子力学去精确求解的，所以实际处理只能
用微扰近似去计算
（即便在漂亮的string theory里，想得到实际有用的信息，当前也多是用微扰去处理
，还没有较好的non-perturbative的理论）
以前做过几年的第一性原理分子动力学模拟计算,因为需要根据薛定谔方程去求解波函
数，用平面波展开去拟合做自洽计算，体系每个周期性结构原子数几十个左右还可以在
32节点每节点16GB内存的并行机上在几天内跑出结果，体系稍大点，上百个原子，就要
个把月了，，，而且很容易出问题，特别是IO是个很大的瓶颈，虽然总体数据量不大，
但计算量非常大，现在国内各大HPC中心不少任务就是跑这个的
类似的，在粒子物理中，Lattice ... 阅读全帖 |
|
发帖数: 1 | 38
爆发了又如何？苏联诺奖也没少拿。 so fucking what?
没有宪政、不尊重私人财产的国家，最后下场好不了
苏联时期研究成果学术成果获诺贝尔奖者如下，一共14人
1.尼;谢苗诺夫(1896-1986)著名化学物理学家，开辟了有关燃烧、爆炸、火焰传播的独
立研究领域，1934年创建了链反应的数量通论，研究了混合气体的热爆炸理论，著有《
化学反应速度与链锁反应》与《化学反应论》，1956年与美国科学家C;欣谢尔伍德共获
诺贝尔化学奖。
2.巴;切连科夫(1904-1958)物理学家，1934年在苏联科学院物理研究所作研究生时
发现了“切连科夫效应”。所谓切连科夫效应是指当带电粒子在某些透明介质中以大于
光在介质中的速度传播时，这种带电粒子就会发出一种特殊的波。切连科夫由于发现和
解释了切连科夫效应，于1958年与苏联物理学家塔姆、弗兰克分享了诺贝尔物理学奖。
3.伊;弗兰克(1908-1990)物理学家，1937年与塔姆一起，对切连科夫效应提
出了理论解释，三人因此同获1958年度诺贝尔物理学奖。
4.伊;塔姆(1895-1971)理论物理学家，主要创立了快速电子的作用和各种物
... 阅读全帖 |
|
发帖数: 1 | 39 这个标签是摘不掉了。
貌似是卢瑟福说的「所有科学要么是物理，要么是集邮」
还有狄拉克也说所有化学理论问题都解决了，剩下的就是怎么近似+计算
但是近似+计算的事情，貌似物理系的人不怎么感兴趣。只有搞化学的人来搞了。 |
|
v***r
发帖数: 1046 | 40 用finite difference 近似解PDE时，采用extrapolation方法来简化计算复杂度，发现
快很多，不知道有什么缺点吗？请指教，多谢 |
|
h****8
发帖数: 49 | 41 在超凡越圣的伟人眼中无穷大n总≈0
————符合实际的全新数学必取代几千年井底蛙数学
黄小宁
通讯：广州市华南师大南区9-303第二信箱 邮编510631
（此文公开发表于：科技信息，2008年第2期）
[摘要]近似计算常识n + 1000n ≈ 0 + 1000n揭示式中由1→∞的n相比下总≈定量0，
其所取各数n相比下全都是可忽略的极小正数，更谈不上n→∞能任意变大取一切非0自
然数n。由此可见书上各取正数的无穷大y均相比下≈定量0，相应的1/y→0也有相比下
总距0极远的另一面。本文揭示N内暗含有>“任给定正数”M的无穷大自然数n>M。
关键词 重大中学数学错误； 语文、物理与近似计算常识；数学以外的更“无理”数
；搞错变量的变域；推翻数学公理和无穷集论；数学革命
科学史上那些千载难逢的重大革命发现造福全人类，但发现的方法是“渔”，远比发现
本身更有价值。思想方法上的革命能使人的科学洞察力一下子提高无穷大倍，从而获超
凡越圣的革命发现。
本文揭示数学中，用而不知的“骨干”数远远多于已知数。例如，如无>任何标准正数
的非标准数及其倒数就绝无非标准微积分一样，若无>“任给定正数” |
|
h****8
发帖数: 49 | 42 在超凡越圣的伟人眼中无穷大n总≈0
————符合实际的全新数学必取代几千年井底蛙数学
黄小宁
通讯：广州市华南师大南区9-303第二信箱 邮编510631
（此文公开发表于：科技信息，2008年第2期）
[摘要]近似计算常识n + 1000n ≈ 0 + 1000n揭示式中由1→∞的n相比下总≈定量0，
其所取各数n相比下全都是可忽略的极小正数，更谈不上n→∞能任意变大取一切非0自
然数n。由此可见书上各取正数的无穷大y均相比下≈定量0，相应的1/y→0也有相比下
总距0极远的另一面。本文揭示N内暗含有>“任给定正数”M的无穷大自然数n>M。
关键词 重大中学数学错误； 语文、物理与近似计算常识；数学以外的更“无理”数
；搞错变量的变域；推翻数学公理和无穷集论；数学革命
科学史上那些千载难逢的重大革命发现造福全人类，但发现的方法是“渔”，远比发现
本身更有价值。思想方法上的革命能使人的科学洞察力一下子提高无穷大倍，从而获超
凡越圣的革命发现。
本文揭示数学中，用而不知的“骨干”数远远多于已知数。例如，如无>任何标准正数
的非标准数及其倒数就绝无非标准微积分一样，若无>“任给定正数” |
|
N***l
发帖数: 52 | 43 有完没完了，数学界的frjj啊。
在超凡越圣的伟人眼中无穷大n总≈0
————符合实际的全新数学必取代几千年井底蛙数学
黄小宁
通讯：广州市华南师大南区9-303第二信箱 邮编510631
（此文公开发表于：科技信息，2008年第2期）
[摘要]近似计算常识n + 1000n ≈ 0 + 1000n揭示式中由1→∞的n相比下总≈定量0，
其所取各数n相比下全都是可忽略的极小正数，更谈不上n→∞能任意变大取一切非0自
然数n。由此可见书上各取正数的无穷大y均相比下≈定量0，相应的1/y→0也有相比下
总距0极远的另一面。本文揭示N内暗含有>“任给定正数”M的无穷大自然数n>M。
关键词 重大中学数学错误； 语文、物理与近似计算常识；数学以外的更“无理”数
；搞错变量的变域；推翻数学公理和无穷集论；数学革命
科学史上那些千载难逢的重大革命发现造福全人类，但发现的方法是“渔”，远比发现
本身更有价值。思想方法上的革命能使人的科学洞察力一下子提高无穷大倍，从而获超
凡越圣的革命发现。
本文揭示数学中，用而不知的“骨干”数远远多于已知数。例如，如无>任何标准正数
的非标准数及其倒数就绝无非 |
|
h****8
发帖数: 49 | 44 在超凡越圣的伟人眼中无穷大n总≈0
————符合实际的全新数学必取代几千年井底蛙数学
黄小宁
通讯：广州市华南师大南区9-303第二信箱 邮编510631
（此文公开发表于：科技信息，2008年第2期）
[摘要]近似计算常识n + 1000n ≈ 0 + 1000n揭示式中由1→∞的n相比下总≈定量0，
其所取各数n相比下全都是可忽略的极小正数，更谈不上n→∞能任意变大取一切非0自
然数n。由此可见书上各取正数的无穷大y均相比下≈定量0，相应的1/y→0也有相比下
总距0极远的另一面。本文揭示N内暗含有>“任给定正数”M的无穷大自然数n>M。
关键词 重大中学数学错误； 语文、物理与近似计算常识；数学以外的更“无理”数
；搞错变量的变域；推翻数学公理和无穷集论；数学革命
科学史上那些千载难逢的重大革命发现造福全人类，但发现的方法是“渔”，远比发现
本身更有价值。思想方法上的革命能使人的科学洞察力一下子提高无穷大倍，从而获超
凡越圣的革命发现。
本文揭示数学中，用而不知的“骨干”数远远多于已知数。例如，如无>任何标准正数
的非标准数及其倒数就绝无非标准微积分一样，若无>“任给定正数” |
|
k*********g
发帖数: 791 | 45 1 “计算材料科学”是大热门；
凝聚态物理实验成本低，但计算材料科学成本更低。
作坊式的、一个一个做实验的传统思想，是很落伍的；
计算材料科学 是理论、数学分析方法的 arm reach，对数值计算的不信任，也是对 数
学、理论的 不信任
2 做计算的 遵循的 方程，全是数学方程，这些方程的由来，做数值计算 都应该非常
清楚，知道他们的意义、局限性。
一个数学上很完美的东西，从数值角度看，很有可能是一堆 trash，因为数学是精确的
，不讲实现的效率的；数值是近似的，讲究效率、稳定性可靠性、flexibility的，这
个叫 ERF（efficiency、reliability、flexibility）；
最典型的例子是 lagrangian multiplier，这个强大的数学武器，在计算结构力学、计
算流体力学领域却是最臭名昭著的算法 |
|
f**********r
发帖数: 18251 | 46 【 以下文字转载自 Military 讨论区 】
发信人: ntkrnl (ntkrnl), 信区: Military
标 题: 国安集团混改涉贱卖国资 近似苏联解体瓜分国家资源
发信站: BBS 未名空间站 (Sat Aug 9 10:27:31 2014, 美东)
媒体：国安集团混改涉贱卖国资 近似苏联解体瓜分国家资源
2014年08月09日 14:36
中国证券网
18
原标题：国安“混改”：中国版世纪大拍卖？
有这么个大型国企，合并资产总额826.35亿、合并净资产近155亿，经营业务涉及信息
产业、资源开发、房地产、旅游、葡萄酒、金融业务，旗下还控股了两家市值分别为
133亿、37亿的A股上市公司。您只要花16亿现金，就可获得该国企15.8%股权。
这么赞的买卖可能吗？在这个神奇的国度，在混合所有制改革风潮下，真发生了！真是
超赞的让利于民企，超便宜的国退民进！
若都这么估值，都这么引进民企，混合所有制肯定红遍全中国，中国将进入“全民混合
时代”。
霸气侧漏的国企改制
这个诡异的国企改制案例就是中信国安[0.24% 资金 研报]集团。8月5日晚，国安集团
旗下的两家上市公... 阅读全帖 |
|
f**********r
发帖数: 18251 | 47 【 以下文字转载自 Military 讨论区 】
发信人: ntkrnl (ntkrnl), 信区: Military
标 题: 国安集团混改涉贱卖国资 近似苏联解体瓜分国家资源
发信站: BBS 未名空间站 (Sat Aug 9 10:27:31 2014, 美东)
媒体：国安集团混改涉贱卖国资 近似苏联解体瓜分国家资源
2014年08月09日 14:36
中国证券网
18
原标题：国安“混改”：中国版世纪大拍卖？
有这么个大型国企，合并资产总额826.35亿、合并净资产近155亿，经营业务涉及信息
产业、资源开发、房地产、旅游、葡萄酒、金融业务，旗下还控股了两家市值分别为
133亿、37亿的A股上市公司。您只要花16亿现金，就可获得该国企15.8%股权。
这么赞的买卖可能吗？在这个神奇的国度，在混合所有制改革风潮下，真发生了！真是
超赞的让利于民企，超便宜的国退民进！
若都这么估值，都这么引进民企，混合所有制肯定红遍全中国，中国将进入“全民混合
时代”。
霸气侧漏的国企改制
这个诡异的国企改制案例就是中信国安[0.24% 资金 研报]集团。8月5日晚，国安集团
旗下的两家上市公... 阅读全帖 |
|
l***o
发帖数: 5337 | 48 数学之美番外篇：进化论中的概率论
BY 刘未鹏 – DECEMBER 2, 2007
李笑来老师在blog上转了一篇宏文，“15 Answers to Creationist Nonsense”；然后
余晟同学（顺便推荐余晟同学译的《精通正则表达式》（第3版））把它给译了出来。
漂亮的文章加上漂亮的翻译，当然是要拜读的:-)
进化论从其诞生以来受到的非难不计其数。这里提到的这篇便收集了广为神创论者提出
以及广为大众误解的一些观点。其中有一点尤其引起了我的兴趣，如下：
8. 严格说起来，我们很难相信复杂如蛋白质的物质能偶然出现，更不用说人或是活细
胞了。
偶然性在进化中确实存在（例如，偶然性的突变可以产生新的特征），但是进化并不依
赖偶然性来产生新的器官、蛋白质或其他实体。截然相反的是，自然选择，作为进化中
已知的最主要机制，却会明确保留“需要的”（能适应的）特性，消除“不需要的”（
无法适应的）特性。只要选择的影响力存在，自然选择就能把进化向一个方向推进，在
出乎意料的短时间内产生复杂的结构。举个例子，现有由13个字母构成的序列“
TOBEORNOTTOBE”，假设有几百万只猴子，每只... 阅读全帖 |
|
s********n
发帖数: 26222 | 49 1楼
列昂尼得·康托罗维奇1912年1月19日（俄历1月6日）生于圣彼得堡。父亲维塔利·莫
伊谢耶维奇·康托罗维奇Виталий Моисеевич Канторович
是一位医生，母亲名叫帕乌琳娜（波琳娜）·格利高里耶夫娜·扎克斯Паулин
а (Полина) Григорьевна Закс。
1926年，列昂尼得考入列宁格勒大学——这一年他14岁。1930年（18岁）从数学系
毕业，随后在母校读研究生。
1932年留校任教，1934年（22岁！）成为教授，1935年他未经论文答辩就被授予物
理-数学博士学位。
1938年列昂尼得·康托罗维奇结婚。他的妻子娜塔丽娅·伊利英娜是一位医生。他
俩生了一儿一女。
1938年列昂尼得·康托罗维奇为胶合板托拉斯解决了有效使用单板镟切机床、实现
最佳剪裁方案的问题。康托罗维奇领悟到，问题可以归结为将一个有着若干变量且限定
为数量众多的线性等式或不等式形式的线性程式予以极限化的任务。他将拉格朗日 （
1736-1813，法国数学家、力学家）的乘数解法进行变形以便用来解决问题并由此领悟
到，数量巨大的经济学问题都可以简化为这一类型的任务。19... 阅读全帖 |
|
