D**s 发帖数: 6361 | 1 2016年1月20日12:01,印度IRNSS-1E由PSLV-XL(C31)自萨迪什亚万航天中心发射成功
印度今年第一射,PSLV还是比较可靠的,顺利完成任务,把一颗印度区域导航系统卫星
送入轨道 |
|
w**********5 发帖数: 1741 | 2 中国四大逼王: 目不识美刘强东,从不碰钱马云爹,1E小目标王健林,房子大点马化腾
SOD新年脑洞大开,推出全新作品《声を出さない女が必死で快感を伝える笔谈SEX》,
整部片几乎听不到任何声音,被封为“史上最安静A片”。 SOD找来爱华美玲饰演“哑
巴素人妻”,只靠纸笔交流。做爱过程中女优全程静音,却不时振笔疾书“再来!再来
!好爽”。 |
|
C********g 发帖数: 9656 | 3 http://my.cnd.org/modules/wfsection/article.php?articleid=31137
■土立方~earththree
俗话说,“最安全的地方往往最危险”,这句用在256号上再恰当不过。
上篇《贰、解开256号航线之迷》分析了,256号强行在山海关机场起飞时并没有准
备好如何飞到外蒙。这个航线只能在起飞后的10分钟时间里,即“问号”阶段根据地标
用航位推测法临时计算出来。所以飞出了270度、290度、310度和340度、325度几个不
同的航向转弯。
256号最终飞出的航线,成功的避开了中国空军的拦截可能,也避开了苏蒙的防空
阵地。但是,此航线也是条从没人走过的“绝版”航线。沿途除了承德市可作为航线检
查点外,直到温都尔汗的全程1011公里(从325度航向转弯点起)都几乎没有再可以参
考的检查点。同时,避开了苏蒙军的空防威胁也就意味着回绝了苏蒙可能的“接待”。
这,对于256号是致命的!
叁、揭开林彪256号坠毁真相
在之前博文《壹、谁有能力击落256号?》中,有读者回复“中国最好的飞行员驾
驶着当代最高档,最安全的飞机,自己掉下来,鬼才信?”。... 阅读全帖 |
|
n**********u 发帖数: 77 | 4 在附上一句,刚才做完了第三次实验,无奈仍然是一样的结果:
0.1 5.787840843
0.1 5.888357639
0.1 5.543690681
0.01 6.571349144
0.01 6.605658531
0.01 6.556361198
0.001 6.80411005
0.001 6.599585533
0.001 6.589509487
1E-04 6.875332832
1E-04 6.537699699
1E-04 6.553117275
1E-05 25.38688278
1E-05 6.44115591
1E-05 6.519612789
1E-06 28.60762405
1E-06 28.89318085
1E-06 28.75268936
1E-07 30.38112831
1E-07 30.44132042
1E-07 30.48270035
样品:
sample 1 6.308719635
sample 1 ... 阅读全帖 |
|
|
n**********u 发帖数: 77 | 6 这个是我的qPCR标准曲线, 怎么是这样的?不是线性。 共7个梯度,一个梯度3个点:
1 4.496082306
1 4.484809875
1 4.525729656
0.1 5.561143398
0.1 5.41778326
0.1 5.438686371
0.01 5.596251965
0.01 5.594720364
0.01 5.552800655
0.001 5.570918083
0.001 5.653218746
0.001 5.576761723
1E-04 5.676965237
1E-04 5.579387188
1E-04 21.28894424
1E-05 24.55447006
1E-05 24.60056686
1E-05 24.69136047
1E-06 28.44901848
1E-06 27.91756248
1E-06 28.0588932
重复两次了,一样的。 有什么可能的原因呢? |
|
g******s 发帖数: 733 | 7 多谢!用[n,x]=hist(),再对n进行处理好象更方便些.
直接用hist画图,再用你说的命令应该也可以,不过自己都不知道是怎么对的:) 用hist(a)出来的纵坐标从1e2到1e6,想让纵坐标的值除以1e6,用
set(gca, 'Yscale', 'log', 'YTick', [1e-4 1e-3 1e-2 1e-1 1], 'YTickLabel', [
1e-4 1e-3 1e-2 1e-1 1]);之后纵坐标没了,再用
set(gca, 'Yscale', 'log', 'YTick', [1e2 1e3 1e4 1e5 1e6]);纵坐标就除以1e6了.
是怎么回事,自己都不知道.
还有,用对数之后曲线成虚线了.
谢了! |
|
S*********4 发帖数: 5125 | 8 DF-21D具体是怎麽打航母的
来源: 文一城 于 09-10-03 19:21:28 [档案] [博客] [旧帖] [转至博客] [给我悄悄
话]
http://blog.huanqiu.com/?uid-6885-action-viewspace-itemid-2009
********************************************************************
作为一支近海防御力量,中国海军建立了以空潜快炮为主的海防体系,然而共和国的
决策者们在海湾战争后猛然发现,苦心孤诣40年的心血已经无法继续信任下去,世纪末
战争中美国海军以射程1000海里的巡航导弹告诉全世界,你们的海岸——不设防。如果
说海湾战争对中国军队还只是震撼,96年的台海冲突则给了海军切肤之痛。没有远洋打
击力量,更不具备远海防御能力,堂堂三百万人民解放军竟然无法对抗美国海军的两个
航母战斗群。近海防御海军的窘迫就在于无法脱离陆地独立行动,不具备远洋攻防作战
体系,不能构筑海上钢铁长城御敌于国门之外,只能作为陆军的附属保卫区区十二海里
领海线。
知... 阅读全帖 |
|
d******e 发帖数: 7844 | 9 X=0 人没病
X=1 人有病
Y=0 检测没病
Y=1 检测有病
P(X=1,Y=0)=1e-10。
P(X=1, Y=1)=1e-6(1-1e-4)
P(X=0, Y=1)=(1-1e-6)1e-4
P(X=0, Y=0)=(1-1e-6)(1-1e-4)
P(X=1|Y=1) = P(X=1,Y=1)/P(Y=1) = P(X=1,Y=1)/[P(X=1,Y=1)+P(X=0,Y=1)] = 0.
00990002
说明基本没谱 |
|
发帖数: 1 | 10 今年以来,中国航天事业屡屡传出好消息,先是在今年4月将我国首颗微重力科学实验
卫星实践十号送上预定轨道,又在今年8月发射了全球首颗量子卫星墨子号。不久前,
神舟11号与天宫二号完成了对接和分离。在11月1日,中国航天科技集团在珠海航展现
场召开新闻发布会上宣布我国首个北斗全球“厘米级”定位系统——“夔龙系统”建设
工作全面启动,夔龙系统通过计算从全球多达300个以上的多系统卫星导航参考站所获
取的观测数据,将卫星导航终端定位精度提高到“厘米级”。
在这些人造卫星和神舟飞船中,有一样电子元件发挥这至关重要的作用,那就是宇航级
CPU,这些装载在人造卫星和神舟飞船上的CPU的作用相当于人类的大脑。那么,宇航级
CPU相对于日常商用的CPU有何不同?宇航级CPU是怎么做到在太空中正常工作的?中国
在宇航级CPU上和美国还有多大差距呢?
CPU在太空中要面临恶劣的工作环境
宇航级CPU构成了人造卫星的大脑,为了能在星际空间这样的恶劣条件下工作,不仅要
应对极端苛刻的高温和低温,还要能应对无处不在的宇宙辐射。
在太空环境中,物体的温度取决于太阳的光照,由于不存在空气散热,受光面和被光面
温差非常... 阅读全帖 |
|
c****3 发帖数: 10787 | 11 抗辐射芯片上中国和美国有多大差距?
在很长一段时间,中国的星载计算机处理器大多依赖进口,国产的人造卫星也大多使用
进口CPU。而美国为首的西方国家严格限制高性能的宇航级和军品级CPU的出口,使我国
星载计算机的研制受到很大的限制。在十多年前,就有美籍华裔学者因购买了几十片军
品级Intel 486 CPU并出售给中国,被美国司法机关和国防部犯罪调查局调查指控的案
例。
在为数不多的军品级、宇航级国产CPU中,其实大多也是对国外产品的逆向工程,比如
应用于航天的386EX,以及被应用于军用电子设备的486DX、ARM7,还有由国内某些单位
仿制过PowerPC603e和SM1750、SM1753、SM1754,国产版的P1750还曾被用于遥感一号、
风云三号、试验四号等卫星。某所还有基于SPARC开源代码修改设计的BM3802RH和
BM3803MGRH。但这些仿制或基于开源代码修改设计的国产CPU,大多存在性能偏低的问
题,P1750、386EX 、486DX、ARM7性能都低于50MIPS,BM3802RH、PowerPC603e虽然超
过了100MIPS,已经强于RAD6000,与... 阅读全帖 |
|
o****o 发帖数: 8077 | 12 要做一个很简单的数值积分问题,其实就是求如下一个简单幂函数的期望
F(U)=U^(a-1), 实数 a > 0。
U的分布是一个Normal(\mu, \sigma) truncated at 0.
在对数据集中每个点进行这个计算的时候,出现了很奇怪的问题。我用的是quad,高斯
quadrature法
当abs(\mu) 远离0点, \sigma又较小的时候,即使\mu1 \mu2很接近计算出来的值也相差
很大。例如
>> quad(@numF1_e,1e-10,1-1e-10,0.00001,0, 3.0692,1.2014,0.5740,9.389589)
ans =
139.9154
>> quad(@numF1_e,1e-10,1-1e-10,0.00001,0, 3.0692,1.2014,0.5740,9.789589)
ans =
2.5851e-005
>> quad(@numF1_e,1e-10,1-1e-10,0.00001,0, 3.0692,1.2014,0.5740,9.689589)
ans =
149.5704
@numF1_e 是那个函数 |
|
c*********2 发帖数: 3402 | 13 21、256号三叉戟飞到北戴河上空了吗?
9月13日凌晨,北戴河林彪别墅的工作人员听到了飞机响声。这可能吗?北戴河在山海
关机场的东北50公里,上空没有航线,34师师长时念堂说:三叉戟起飞后,再也没有返
回山海关机场。飞机时速500公里,早超过了北戴河的范围。但是如果飞机建立起落航
线,北戴河就能听见。9月12日晚上林立果乘三叉戟到北戴河,在外面乘凉的林办工作
人员就听到了飞机落地的声音。256三叉戟从山海关机场强行起飞,先向西南,右转弯
,对着北京方向,到迁安(这是飞行员潘景寅的老家)改向西北。
时任北戴河北京军区疗养院外科护士的杨桂兰回忆:1971年9月13日凌晨她起床给儿子
换尿布,忽然听见房顶轰隆隆响,大得好像房顶上开过去拖拉机。杨桂兰是飞行员家属
,她对飞机起飞并不陌生,但她说在北戴河半夜听见飞机声是从来没有过的事情,所以
她根本没有当成是飞机声。北戴河北京军区疗养院和林彪别墅是邻居,住在林彪别墅附
近的人听到256三叉戟起飞也是真实的。
林彪警卫秘书李文普也清楚地听到了256三叉戟起飞的声音。李文普自伤后,被八三四
一部队的卫生员简单包扎,经外科医生出身的张清霖检查,之后很... 阅读全帖 |
|
u*******r 发帖数: 864 | 14 3月31日,中国北斗卫星导航系统多了一位新成员。作为第17颗北斗卫星,这位新成员
接过了“自主可控”的接力棒。据该卫星总设计师、来自上海微小卫星工程中心的林宝
军介绍,这颗由中科院负责研制的新一代导航卫星大量使用国产化器部件,并首次使用
中国制造的“龙芯”中央处理器。
“此次搭载的处理器是‘龙芯’的两款芯片——龙芯1E和龙芯1F抗辐照处理器。”中科
院国家空间科学中心研究员、中科院复杂航天系统电子信息技术重点实验室副主任安军
社在接受《中国科学报》记者采访时说。
安军社表示,新一代北斗卫星上有3个被称为“单机”的黑盒子,每个约有4本400页的
32开图书摞起来那么大。其中两个黑盒子里,每个装了2片龙芯1E芯片和4片龙芯1F芯片。
“这两个黑盒子的主要任务是用于完成星间链路的数据处理,即支持北斗各颗卫星之间
的通讯数据处理任务。龙芯1E负责进行常规运算,龙芯1F完成数据采集、开关控制、通
讯等处理功能。”安军社说。
“目前,与我们在国际上能买来的芯片相比,龙芯1E和龙芯1F的性能还是要更好一些,
处理能力和可靠性也与国际上能买到的持平。”安军社告诉记者。 |
|
f*******3 发帖数: 206 | 15 3是典型的bayes formula, 求conditional probability。
X:0 坏人, 1 好人
Y:0 报坏人, 1 报好人
已知 P(X=0)= 100/1M = 1e-4, P(X=1) = 1-1e-4 \approx 1, P(Y=1|X=0)=P(Y=0|X=1
)=0.01
求 P(X=0|Y=0) = P(X=0,Y=0)/P(Y=0)
= P(Y=0|X=0)P(X=0)/
P(Y=0|X=0)P(X=0)+P(Y=0|X=1)P(X=1)
= 0.99*1e-4/
0.99*1e-4 + 0.01*1
= 0.01
说明坏人绝对数量少,报错有概率的话,报警器可信型很低。 |
|
l*****e 发帖数: 276 | 16 我在dual core Opetron 64位 Ubuntu 机器上用g++ -O3 ./test.cpp -o test 编译的
两个程序。time显示第一个只用0.195s 但第二个却用了0.668s。有人知道为什么么?
难道说如果用64位机器把所有用float的地方换成double反而会节省时间? 谢谢。
int main()
{
for(double b=0;b<1;b+=1e-4) for(double a=0;a<1;a+=1e-4);
}
int main()
{
for(float b=0;b<1;b+=1e-4) for(float a=0;a<1;a+=1e-4);
} |
|
n******7 发帖数: 12463 | 17 是差10倍
一个是1e-3,一个是1e-4
但是这对结果有多大影响?
比如目标是call出AF>=0.01的variant site
这个就是1e-2,跟1e-3的错误率都差了10倍了
况且Q30算比较保守的了,实际大部分base的质量要高很多
另外,quality score也不过是infer出来的,又不是ground truth
没必要太当回事
我自己也搞过计算quality score的模型
很多时候不过给人一些安慰罢了
即使对WGS,我觉得把quality score离散化,分成几个类
比如不靠谱,有点靠谱,靠谱,非常靠谱
就够了,即减少数据体积,也不会过度解读这个score
要不你给我个简单的例子
说说Q30和Q40怎么区别大了
我给你发包子? |
|
f**l 发帖数: 2041 | 18 I would choose c.
reason: f(x) is almost like x^8 since x is around 3.
then let's assume the error behave also like in the function x^8.
the problem then change to estimate (x+1e-7)^8 - x^8.
it's estimated as: 8*x^7*1e-7 just like u take the linear approximation.
ok, now 8*3^7*1e-7 is at the magnitude of 1e-3. so choose c. |
|
b*******r 发帖数: 1130 | 19 例如TSMC.25um的spice model
.MODEL CMOSN NMOS ( LEVEL = 49
+VERSION = 3.1 TNOM = 27 TOX = 5.7E-9
+XJ = 1E-7 NCH = 2.3549E17 VTH0 = 0.3951677
+K1 = 0.4881677 K2 = -1.465714E-6 K3 = 1E-3
+K3B = 1.2222493 W0 = 1E-7 NLX = 1.927361E-7
+DVT0W = 0 DVT1W = 0 DVT2W = 0
+DVT0 = 0.4751286 DVT1 = 0.5660833 DVT2 = -0.3026132... 阅读全帖 |
|
a**r 发帖数: 224 | 20 这是他们不同方法的结果
Micromeritics给出的第一个数据点,对应的相对压力可能是 ~1E-7
而Quantachrome的仪器,你要是从1E-7开始的话,可能要take for ever
所以我一般都从1E-6 or 1E-5开始,对应的吸附量假设是~100 cc/g
这个吸附量要是用Micromeritics的话,就是第20个数据点(if 5 cc/dose)
至于哪种方法更好,各有各的说法
反正micropore pore size distribution and micropore volume strongly depend on
the models you are going to use.
Quantachrome给的程序包里有很多不同的models
不知到Micromeritics现在有多少models
我以前用的ASAP2010好像就一个DFT
BTW:包子收到,笑纳了 |
|
a***r 发帖数: 420 | 21 处理一堆1E-300左右的数,统统log之
为了需要相加的时候,写了个小function:
logsum <- function(obj) {
vect=as.vector(obj)
logsumm = -Inf
for (i in 1:length(vect)) {
temp=vect[i]
if (temp>-Inf) {
logsumm=temp+log(1+exp(logsumm-temp))
}
}
return (logsumm)
}
#ie. a=c(1E-10,1E-20,1E-200); b=log(a); then: logsum(b)=log(sum(a))
似乎运行得很好
可是今天发现一个bug(?),以下e是在实际计算中产生的array,需对里面每个element求
和,不辞
辛劳地抄下来,give you a flavor of how puzzled I am...:
>e=array(c(-384.1895,-299.9747,-512.4600,-591.3290,-299.9747,-278.6... 阅读全帖 |
|
v*******e 发帖数: 11604 | 22
这太容易了。Pvalue can be as small as you want.
N=c(100,1000,10000,100000)
PValue=rep(as.numeric(NA),length(N))
RSquared=rep(as.numeric(NA),length(N))
for (i in 1:length(N)){
sigmma=rnorm(N[i],0,1)
x=rnorm(N[i],0,1)
y=0.1*x+sigmma
RSquared[i]=summary(lm(y~0+x))$"r.squared"
PValue[i]=summary(lm(y~0+x))$"coefficients"[1,4]
}
data.frame(N,RSqaured,PValue)
N RSquared PValue
1 1e+02 0.007324693 3.947863e-01
2 1e+03 0.014429159 1.392340e-04
3 1e+04 0.010822982 1.792036e-25
4 1e... 阅读全帖 |
|
s*******d 发帖数: 370 | 23 我看节目时,挺佩服周炜心算出结果的准确度,但是同时对他只写出整数或最多小数点
后一位感到有些小失望。因为我期待的天才应该是,如果结果为无限小数,主持人没喊
停的话,他应该会固执的一直写下去….
看到肘子和板上一些人的质疑,本人尝试看看能否通过背诵准确推测出精确到小数点后
一位的结果。我以节目中的15位数为例,开14次方到16次方,发现只要背诵 29个数值
(精确到小数点后两位)就能达到目的,即准确推测出任何一个15位数,开14次方到16
次方,精确到小数点后一位的结果。
下面是excel输出的需要记住的30个结果
14√ 15√ 16√
1E+14 10 8.57 7.49
2E+14 10.50 8.98 7.83
3E+14 10.81 9.22 8.03
4E+14 11.04 9.40 8.17
5E+14 11.21 9.54 8.29
6E+14 1... 阅读全帖 |
|
s*****r 发帖数: 43070 | 24 D68成都东-武汉,1E成都局
D328成都东-上海虹桥,1E成都局
D368成都-武汉,1E成都局,这个和D68套跑
D2202成都东-郑州,2A郑州局
D2208成都东-上海虹桥,2A上海局
D2224成都东-杭州,1B上海局
D2238成都东-南昌,2A南昌局
D2244成都东-福州,1A南昌局
D2256成都东-南京,1A成都局
D2260成都东-汉口,2A武汉局
D2218重庆北-上海虹桥,2A上海局
D2228重庆北-福州,2A成都局
D2234重庆北-厦门,2A成都局
D2252重庆北-汉口,2A武汉局
D2266重庆北-宜昌,2A成都局
D2274重庆北-武汉,2A成都局
D2278重庆北-武汉,2A成都局,这个和2274套跑 |
|
h**********g 发帖数: 3962 | 25 有25匹马,任意两匹马的速度都不同。一次比赛可有不多于5匹马参加。
问:最少需要几次比赛参能选出最快的五匹马?
下面是一个最多使用8次比赛的方法。
第一步:
把25匹马分成5组,每一组有5匹马。每一组赛一次。
不失一般性,假设结果如下。
1A < 1B < 1C < 1D < 1E
2A < 2B < 2C < 2D < 2E
3A < 3B < 3C < 3D < 3E
4A < 4B < 4C < 4D < 4E
5A < 5B < 5C < 5D < 5E
第二步:各组的第二名在一起比赛。不失一般性,假设结果如下
1B < 2B < 3B < 4B < 5B
现在我们可以知道
(1) 1A一定入选
(2) 2D, 2E; 3B, 3C, 3D, 3E; 4B, 4C, 4D, 4E; 5B, 5C, 5D, 5E被淘汰
第三步:1D, 2B, 3A, 4A, 5A比赛。
我们分几种情况考虑。
(a) 1D < 2B < {...}
答案是1A, 1B, 1C, 1D加上min{1E, 2A}
(b) 1D < iA < {...}, i=3, 4, 5
答案是1A... 阅读全帖 |
|
Q****a 发帖数: 296 | 26 我用牛顿迭代法和二分法都试过, precision 最小能设成 1e-15 (或者9e-16),再
小(比如1e-16) 就死循环了 = =
可是我用precision = 1e-15做,返回结果比java里面Math.sqrt() 精度少一位 (比如
sqrt(2.0) 我return 1.414213562373095, 而java math return 1.
4142135623730951)。
请问大牛Math。sqrt怎么能精度多一位的? |
|
l*******s 发帖数: 26303 | 27 再帮你把你的想法具体一下,其实第二题靠简单记忆完全能应付,不需要速算。
发信人: scorfield (飞鱼), 信区: Military
标 题: Re: 周玮的最强大脑就是个娱乐噱头
发信站: BBS 未名空间站 (Tue Jan 21 22:53:46 2014, 美东)
我看节目时,挺佩服周炜心算出结果的准确度,但是同时对他只写出整数或最多小数点
后一位感到有些小失望。因为我期待的天才应该是,如果结果为无限小数,主持人没喊
停的话,他应该会固执的一直写下去….
看到肘子和板上一些人的质疑,本人尝试看看能否通过背诵准确推测出精确到小数点后
一位的结果。我以节目中的15位数为例,开14次方到16次方,发现只要背诵 29个数值
(精确到小数点后两位)就能达到目的,即准确推测出任何一个15位数,开14次方到16
次方,精确到小数点后一位的结果。
下面是excel输出的需要记住的30个结果
14√ 15√ 16√
1E+14 10 8.57 7.49
2E+14 10.50 8.98 ... 阅读全帖 |
|
t********t 发帖数: 5415 | 28 再这么下去最先倒霉的应该是theo。offseason叫嚣什么"run prevention",结果搞来
的scutaro现在3E,beltre 1E,cameron 1E, bill hall 2E, 今天hermida又添1E...目
前只有lackey达到大家的期望... |
|
s*******d 发帖数: 370 | 29 我看节目时,挺佩服周炜心算出结果的准确度,但是同时对他只写出整数或最多小数点
后一位感到有些小失望。因为我期待的天才应该是,如果结果为无限小数,主持人没喊
停的话,他应该会固执的一直写下去….
看到肘子和板上一些人的质疑,本人尝试看看能否通过背诵准确推测出精确到小数点后
一位的结果。我以节目中的15位数为例,开14次方到16次方,发现只要背诵 29个数值
(精确到小数点后两位)就能达到目的,即准确推测出任何一个15位数,开14次方到16
次方,精确到小数点后一位的结果。
下面是excel输出的需要记住的30个结果
14√ 15√ 16√
1E+14 10 8.57 7.49
2E+14 10.50 8.98 7.83
3E+14 10.81 9.22 8.03
4E+14 11.04 9.40 8.17
5E+14 11.21 9.54 8.29
6E+14 1... 阅读全帖 |
|
l*********y 发帖数: 142 | 30 下面这个方程计算 num! 的结果在base进制下的位数
比如说 5!的结果在10进制下是3位数,在16进制下是2位数。
题目限定 1<= num <= 2^20, and 1 <= base <= 800.
int GetNumDigit(int num, int base)
{
double result = 0;
for (int i = 1; i <= num; i++) {
result += log(i);
}
int digit = ceil(result / log(base) + 1e-10);
return digit;
}
我想请教1e-10在这里的使用? 我读过很多float and double精度的文章,所以知道result/log(base)的结果会有精度loss,但是我对精度补偿值的用法不是很懂。
包括何时用 还有 如何定值。以前还碰到一个例子是atof 时加了1e-10,
现在结合这个例子就完全没方向了。
谢谢了! |
|
s*******4 发帖数: 164 | 31 server的log如下:
Jun 26 12:12:46 os dhcpd: Listening on LPF/eth0/00:22:6b:bb:bb:df/192.168.1/
24
Jun 26 12:12:46 os dhcpd: Sending on LPF/eth0/00:22:6b:bb:bb:df/192.168.1/
24
Jun 26 12:12:46 os dhcpd: Sending on Socket/fallback/fallback-net
Jun 26 12:13:03 os dhcpd: DHCPDISCOVER from 00:1e:8c:33:db:ff via eth0
Jun 26 12:13:04 os dhcpd: DHCPOFFER on 192.168.1.254 to 00:1e:8c:33:db:ff
via eth0
Jun 26 12:13:07 os dhcpd: DHCPREQUEST for 192.168.1.254 (192.168.1.1) from
00:1e:8c:33:db:ff via eth0
Ju |
|
g*********s 发帖数: 1782 | 32 #include
#include
#include
#include
int main(){
float f = 0.0001;
printf("f = %f\n", f);
printf("f = %g\n", f);
printf("f = %e\n", f);
std::cout <<"f = "<
std::cout.setf(std::ios_base::fmtflags(0), std::ios_base::floatfield);
std::cout <<"f = "<
}
RH3:
f = 0.000100
f = 1e-04
f = 1.000000e-04
f = 1e-04
f = 1e-04
RH4:
f = 0.000100
f = 0.0001
f = 1.000000e-04
f = 0.0001
f = 0.0001
有啥办法统一格式么?机器是RH3/RH4都有。改程序也 |
|
b*****y 发帖数: 163 | 33 C has long double, but may not be implemented by particular compiler,
since there is little numerical advantage to do long double rather
than just double precision. In double precision, eps = 1e-16 which
is enough to handle your 1e-10. For single precision, eps = 1e-07.
Most numerical computation is carried under either single or double
precision, which are IEEE standard, using 32 or 64 bits coding and very
easy to handle with modern computer architecture. One can
achieve arbitrary eps such as 1 |
|
h***o 发帖数: 539 | 34 来自主题: Computation版 - 菜鸟问题 打个比方
1.0 + 1e-27 = 1.0
1e-27 + 1e-27 = 2e-27
计算机是这么作算术的 |
|
h***z 发帖数: 233 | 35 With a condition number of 4.6328e+020, your matrix might still technically
be full rank, but for all practical purposes, it is essentially rank-
deficient. For example, consider the matrix
A = [ 1e-20 0; 0 1 ]
A is technically full rank with a condition number of 1e+20, but in
practice, it really is very similar to the rank deficient matrix [ 0 0;
0 1]. For example, consider the vectors
u = [ 0 , 1 ]' and v = [ 1e-10 , 1 ]'
u and v are very similar and should |
|
g******s 发帖数: 733 | 36 多谢了!按照你说的,写了几行程序,但是算出的电压都有无穷大和复数了。
由方程 V(x)^2 dy/dx=a x,得V(x)=sqrt(a x dx/dy) = sqrt(alpha.*2.*omega_r.^2.*I./N./dc_dalpha)。请问程序错误在什么地方?所有参数都取自原论文。
N=120;
omega_r=240*2*pi; %omega_r=240;
density=2330;
l=1e-3;
w=1e-3;
h=50e-6;
I=density*h*l*w^3/12;
alpha=[0 1 2 3 4 6 7].*pi./180; % 取自图6b
dc_dalpha=[-0.25 0.495 0.75 0.74 0.495 -0.505 -1].*1e-10; % 取自图6b
V=sqrt(alpha.*2.*omega_r.^2.*I./N./dc_dalpha); % 公式 (4)
谢了! |
|
g******s 发帖数: 733 | 37 多谢了!我用公式V(x)=sqrt(a x dx/dy) = sqrt(alpha.*2.*omega_r.^2.*I./N./dc_dalpha)算出的电压和转角的关系,原文计算出的电压和转角的关系,都附在附件上。看横坐标,电压相差一个数量级。
N=120;
omega_r=240*2*pi; %omega_r=240;
density=2330;
l=1e-3;
w=1e-3;
h=50e-6;
I=density*h*l*w^3/12;
alpha_rad=[1 2 3 4].*pi./180; % 取自图6b
alpha_degree=[1 2 3 4]; % 取自图6b
dc_dalpha=[0.495 0.75 0.74 0.495].*1e-10; % 取自图6b
V=sqrt(alpha_rad.*2.*omega_r.^2.*I./N./dc_dalpha); % 公式 (4)
plot(V,alpha_degree)
还有什么可能的错误?原文作者不理我。
多谢了! |
|
l*****e 发帖数: 276 | 38 算2PI和Asin(x)在[0,2PI]的曲线长度的比,当A=0时Asin(x)=0,是一条直线,比值应
该是1。
在Maple里用plot画出的结果在A=0时是1.2..,如果手动用pointplot画那么结果是对的。
d := 2*3.1415926; T := proc (A) options operator, arrow; d/(int(sqrt(1+A^2*
cos(x)^2), x = 0 .. d)) end proc;
with(plots); pointplot({seq([0+0.1e-1*xx, T(0+0.1e-1*xx, 6)], xx = 0 ..
1000)});
T(0.1e-3); plot(T(C), C = 0 .. 10);
同样的东西在Mathematica里就可以给出正确的结果
d = 2*3.1415926
B[A_] = d/\!\(
\*SubsuperscriptBox[\(\[Integral]\), \(0\), \(d\)]\(
SqrtBox[\(1 +
\*SuperscriptBox[\(A\), \(2\)]*Cos[ |
|
v***y 发帖数: 25 | 39 The material is lumina. The input file is like:
** MATERIALS
**
*Material, name=Material-1
*Density
780.,
*Elastic, type=Lamina
1.2e+10, 1.56e+10, 0.01, 1e+08, 1e+08, 1e+08
*Fail Stress
100.e6,-100.e6,100.e6,-100.e6,10.e6
但是即使有远远超过100.e6Pa的应力,也观察不到failure。同时abaqus也不删除相应
的单元。
有哪位大侠知道是怎么回事吗? 着急要用。 多谢多谢。 |
|
e**********n 发帖数: 359 | 40 扯到哪里去了?这里没有什么弱相互作用。高激发态的电子云扩散得很开,其外围与其它原子有相互作用。这个原理费米早就研究
过了,无非是个微扰展开式的计算,得到一个原子在另一个原子高激发态电子云中运动的等效势能。电子云密度大的地方势能是吸
引的,所以就有了束缚态。现在实验里观察到了这些束缚态,写文章的就提出了第五种化学键的说法来吸引眼球。
实验里用冷原子气,先用激光激发一部分原子,再测量激发态能谱,在单原子激发态峰一侧有一个小峰,正好对应与双原子的束缚
态。
这些人是打量子计算的旗号做冷原子气的研究,BEC已经得过炸药奖了,所以就在与BEC接近得系统中找一些好推销的东东,否则
花了纳税人的钱也不好交代。要问有什么应用,看看实验参数就知道了,1e-6K 量级的温度,每立方厘米1e+10 到1e+11个原
子的密度,与原子高激发态差不多的寿命。要用这个东西,现置办一个冷原子气实验室,办成后当然会有助于量子计算机的物理实
现,不过还需要继续投很多银子下去。 |
|
l**n 发帖数: 908 | 41 一般来说Nanovoltmeter (象Keithley 2182) 达不到1E-9 V的精度,实际上1E-8 V都
达不到。不过1E-7 V是可以的。这些都是DC信号。一个刚calibrated的表,我看过它在
4K的drifting,基本上是在-8的位置,跳动范围大概是二十到三十个nV。
AC信号,个人经验是低于uV就开始不行了。digital的lock-in普遍比不上以前普林斯顿
那款analog的。
提高精度最好的办法还是改善屏蔽并用preamplifier来放大信号(当然还要保证信号不
失真)。
有兴趣可以看看普林斯顿N.P. Ong他们组做的Nernst Effect的工作,好像是把噪声控
制在-8以下了。 |
|
j******n 发帖数: 1568 | 42 简单做个半定量分析。
摄像头分辨本领的极限是瑞利判据。角分辨率约等于1.22 * 波长/通光孔径。取可见光
波长0.6微米, 摄像头通光孔径6毫米, 角分辨率的极限是 1.22 * 1E(-4)弧度。考虑
到ccd的分辨率以及各种像差,光照的影响, 实际的角分辨率一般要下降一个数量级,
也就是1.22 * 1E(-3)弧度。 摄像头与被摄物体的距离大多在10米这个量级, 所以摄
像头能分辨得两个点的最小间隔就是1.22厘米。这和戒指的大小类似。也就是说摄像头
能看到的戒指就是一个点。我不知道法官长了对什么样的眼睛能够从录像里看出那是个
戒指。 |
|