由买买提看人间百态

topics

全部话题 - 话题: 曲线拟合
1 2 3 4 5 下页 末页 (共5页)
a***g
发帖数: 2761
1
来自主题: Statistics版 - 请教一个曲线拟合的问题 (转载)
【 以下文字转载自 Computation 讨论区 】
发信人: akoug (xuantai), 信区: Computation
标 题: 请教一个曲线拟合的问题
发信站: BBS 未名空间站 (Fri Dec 4 12:27:58 2009, 美东)
这个问题的表述是我自己想的:
假设有n个数据:x1,.....,xn;
现在把1,.....,n视作自变量,把x1,.....,xn视作因变量
可以假设xi=f(i)
如果x1,....,xn真的服从这样的曲线f
这个时候尚不知道f的表达式
只能先按数据拟合一条光滑的曲线g
那么用什么方法拟合g同时计算其相关微分可以比较准确的估计f在{1,....,n}这n个点上的导数呢
谢谢
i*****r
发帖数: 65
2
来自主题: Mathematics版 - 请教一个曲线拟合的问题
发到这个版块来,不知道哪位仁兄又做曲线拟合的经历? 多谢了!
附件说明: 要求用那个方程来拟合x,y组成的曲线,求出a 和b. x和y在excel里面。
email: w******[email protected]
K*****2
发帖数: 9308
3
曲线拟合都是扯蛋的,你想拟合出什么都行
不过破10万是极有可能的
a***g
发帖数: 2761
4
来自主题: Mathematics版 - 请教一个曲线拟合的问题 (转载)
【 以下文字转载自 Computation 讨论区 】
发信人: akoug (xuantai), 信区: Computation
标 题: 请教一个曲线拟合的问题
发信站: BBS 未名空间站 (Fri Dec 4 12:27:58 2009, 美东)
这个问题的表述是我自己想的:
假设有n个数据:x1,.....,xn;
现在把1,.....,n视作自变量,把x1,.....,xn视作因变量
可以假设xi=f(i)
如果x1,....,xn真的服从这样的曲线f
那么有什么方法可以比较准确的估计f在{1,....,n}这n个点上的导数呢
谢谢
a***g
发帖数: 2761
5
来自主题: Computation版 - 请教一个曲线拟合的问题
这个问题的表述是我自己想的:
假设有n个数据:x1,.....,xn;
现在把1,.....,n视作自变量,把x1,.....,xn视作因变量
可以假设xi=f(i)
如果x1,....,xn真的服从这样的曲线f
这个时候尚不知道f的表达式
只能先按数据拟合一条光滑的曲线g
那么用什么方法拟合g同时计算其相关微分可以比较准确的估计f在{1,....,n}这n个点上的导数呢
谢谢
H*******5
发帖数: 83
6
来自主题: Physics版 - 请问关于Matlab曲线拟合
这里牛人多,我就发这问了。
请问诸位大侠,
Matlab的polyfit和lsqcurvefit对数据拟合只给参量的拟合结果,不给误差。请问如何
才能得到相应的拟合误差(或置信区间)呢?(指拟合参量的误差,不是总体数据的\chi^
2)
或者Matlab是否有其他函数能实现呢?
感激涕零!!!
x*****8
发帖数: 10683
7
来自主题: Military版 - 逻辑斯蒂曲线
哪位有graphing calculator的同学帮我把美国这几个月数据在不知道carrying
capacity的情况下用logistic regression拟合一下吧。
本虾在1月份用逻辑斯蒂曲线拟合中国数据是很成功的,预测拐点只差了2天。
t**********a
发帖数: 7
8
来自主题: Computation版 - 如何比较两条曲线的相似程度?
谢谢大家回帖.
是啊,我前面帖子里已经作了 cross correlation,但是不能反映曲线的大小相差程度,即
index = Sum( (xi-average(X))*(yi-average(Y)) )/ sqrt( sum((xi-average(X))^2 )
* sum ((xi-average(Y))^2) )
但是cross correlation不能反应两列值的大小相差程度,例如将X乘2,结果一样.
另外一个指标,"曲线拟合度"(coefficient of effciency) 能够较好反映出而列数的值的
大小关系
index = 1- sum((xi-yi)^2)/sum( (yi-average(Y))^2 )
但是对于shape反映不好.我要做的曲线里除掉杂音后,一般有1到3个峰. 类似的曲线应该
峰的数量一致.

)
是Y,sorry
f***o
发帖数: 171
9
谢谢回答
你的意思是,没法象*PLASTIC那样直接用表单功能加入新的piece-wise应力应变曲线,
只能写自己的UMAT,用拉/压实验曲线拟合所选本构的参数?
我看了ABAQUS的casting iron和修正DP material model,都无法描述sofening,我用
的材料sofening不可忽略,有没有什么更好点的方法
谢谢!
h***s
发帖数: 226
10
来自主题: Mathematics版 - 请教一个曲线拟合的问题
看来这里的“数学家”是不会帮你解决这么具体的问题的。
方法可能有几个
1.弄一个数值计算工具,比如matlab等,里面有回归的部分,选择需要拟合的函数,和
xy数值,得到ab。但是否能有你需要的函数种类,不得而知;
2.弄一个已经编写好的非线性回归程序控件,然后写你的输入输出函数调用即可:http://download.csdn.net/source/142460
3.自己写程序,比如利用非线性最小二乘方法编写程序,不过这个程序还不是那么容易
编写的,如果谁白手起家不看任何资料就能写出来,是比较牛的。
g****t
发帖数: 31659
11
来自主题: Mathematics版 - 请教一个曲线拟合的问题
matlab有专门的非线性最小二乘法的函数.
查一下help lsqnonlin
或者用fminsearch

看来这里的“数学家”是不会帮你解决这么具体的问题的。
方法可能有几个
1.弄一个数值计算工具,比如matlab等,里面有回归的部分,选择需要拟合的函数,和
xy数值,得到ab。但是否能有你需要的函数种类,不得而知;
2.弄一个已经编写好的非线性回归程序控件,然后写你的输入输出函数调用即可:http://download.csdn.net/source/142460
3.自己写程序,比如利用非线性最小二乘方法编写程序,不过这个程序还不是那么容易
编写的,如果谁白手起家不看任何资料就能写出来,是比较牛的。
n******g
发帖数: 17225
12
请问哪个是这个二次曲线的quadratic term?
k**o
发帖数: 3006
13
☆─────────────────────────────────────☆
apolitical (无期徒刑) 于 (Thu Apr 11 12:52:03 2013, 美东) 提到:
趁板神在黑屋里关着,我来行为艺术啦。用希大师的eos 3拍了一场比赛。说明书上写
连拍是每秒4.7张,但实际操作达不到,不过我算了一下,一卷24张,洗扫花8块多,基
本按动快门每秒就扔出去一块多钱,连拍体育有肾上腺激素的涌动,这次跟着流动的还
有现金。
重点难题在对焦上,虽然能在胶片机上找到佳能现在数码对焦操作的影子,但我还没找
到感觉。对上焦很快很准,但追焦跟踪的难度很大。都说佳能技术储备强,我时光倒流
追抚奸商的过去,还需要更多时间摸索。我应该还会再谈谈对焦问题。
fuji superia 400
kodak ektar 100
顺便对比一下用7d拍的,光照条件相似。
对比胶片和数码,除了明显的不同色调外,胶片在明暗过度的处理上会不会好一些?条
件是大太阳直射,无云,7d给了明显的数码味,对比度很强,显得光很硬,胶片则显得
在细节上不是那么对比强烈,能看出是晴天,但不是那种让人睁不开眼... 阅读全帖
m*l
发帖数: 507
14
来自主题: Environmental版 - GIS在环境工程中的应用(9)
1.基于地理的研究方法
2.与环境相关的地理研究
3.GIS在市政工程中的应用
4.GIS在水文水利中的应用, 2
~~~~~~~~~~~~~~~
5.GIS在其它环境工程中的应用
首先说一下地形是怎么测出来的,就是一个点一个点测,要想准确就测的密一些。所以
,地形的最原始的数据就是点。有了这个点以后,就可以画等高线了。
在gis中的做法有两种,第一种是根据已知的点通过插值的方法来得到一个连续的平面。
三维的插值不太好理解,可以考虑一下二维的插值,其实就是曲线拟合。
大家都知道只要是曲线拟合就会有误差,而且不同情况要用不同的曲线来拟合。对三维
地形也是同样的,简单的可以是线性的,复杂就要用到统计方法来决定样本和其权值了。
另一种是triangulated irregular network (TIN),与上一种方法不同的是,当已知的
点已经给定之后,算法会自动生许多小的三角形,从而形成一个矢量的曲面。
这两种方法互有利弊,所以在gis中大部分情况都是结合在一起使用。
在实际工作中并不需要了解gis具体实现这些功能的细节,在esri arcgis中,通过使用
3d analyst+spa
q*********n
发帖数: 103
15
借楼主的帖子问下三次样条拟合。
Matlab里用的spline函数是not a knot的边界条件,不知道实际效果和自然边界(两头二阶
导数为零)有多大的区别?
还有一种 Kochanek-Bartels Hermite Spline interpolation,能通过调节三个参数
得到不同的拟合效果。这种spline是怎么算的?常用么?
(1) tension: 调节峰的尖锐程度
(2) continuity: 节点两边的平滑度
(3) bias: 好像是关于节点两侧对称程度
http://www.shaneaherne.com/research/splines.html
请问怎么为不同类型的曲线选spline的算法? 比如,是不是光滑曲线要用自然三次样
条, 有很多峰值的曲线则应该用其他的方法?
多谢多谢~
s***e
发帖数: 911
16
来自主题: Science版 - Re: how to 数值拟合?

我这个办法肯定不很好,不过你这个例子不对, hehe...针对你这个例子我的办法一定
给出a=1. 这个问题里面你有三个实验点. 其中任取两个都可以定出a,b. 你这个例子
无论取那两个都是a=1, b=0.于是最后平均还是a=1, b=0.
我猜是我没说清楚,你误解了我的算法. 我这个算法的思想是:
平面上有N个点, 我用一个曲线模型R(t,P)来拟合.加入n P, 则我们根据实验点共可以产生C_N^n条可能的曲线. 最后对参数平均, 其实就是对这个
C_N^n条曲线作了平均.
假如我们的曲线根据的model是正确的, 假设真正的那个曲线是R(t,P0). 那么实验点
确定出来了曲线系宗就看成围绕R(t,P0)的随机波动. 假设这个分布是Gaussian的,
那么你知道平均正好对应最可几.
以上是偶的办法根据的一点儿道理乐:P
n******g
发帖数: 17225
17
来自主题: WaterWorld版 - 半年连续12跳,富士康中邪了?
数学和统计建模从来都没有证明模型正确描述了"真理"这个步骤啊.
原文首先从这次事件是一种情绪diffusion的假设开始,然后用于拟合diffusion最常用
的工具就是logistic曲线,所以得到了后续的拟合结果。
如果从阴谋论或者化学毒素造成幻觉等其他基本假设出发,应该不会去用logistic曲线
拟合,甚至不会有适合的模型去拟合。
R*j
发帖数: 248
18
来自主题: Pharmaceutical版 - 有人做LC/MS/MS吗?问个标准曲线问题
经常看到有些LC/MS/MS的标准曲线不是线性的。然后作者用一些奇怪的方程去拟合。(
我理解加权是必要的)
谁能讲一下如果标准曲线不是线性的,那末用什末方程去拟合啊?
z*****4
发帖数: 45
19
我试图将数据拟合到一个模型中: Signal = a*exp(-R1*t)+b*exp(-R2*t).要拟合的参
数是R1,R2,a和b。 R1和R2更重要。有12个数据点: 12个t,对应12个Signal。我想用
贝叶斯概率理论来求R1,R2,a和b。 我在看贝叶斯分析的逻辑回归模型: https://www
.mathworks.com/help/stats/examples/bayesian-analysis-for-a-logistic-
regression-model.html。 但是,我还没想出如何写the prior distribution, the
posterior distribution, or maximum likelihood for my question。 这个网页也可
能是相关的: https://www.mathworks.com/help/stats/examples/curve-fitting-and-
distribution-fitting.html。
有人可以教我贝叶斯方法的实现数据拟合么,或是分享Matlab code实现贝叶斯方法... 阅读全帖
f**********r
发帖数: 91
20
【 以下文字转载自 Chemistry 讨论区 】
发信人: flyingwisher (flyingwisher), 信区: Chemistry
标 题: 紧急求助:关于origin 的线性拟合,谢谢!
发信站: BBS 未名空间站 (Thu Jul 26 12:26:21 2007)
现在有一条曲线,由几个点组成,类似抛物线,现在想把它分别进行线性拟合,就是作
成两条相交的直线,然后确定交点的值,但是我发现,分别拟合出来的红线,显示数值
范围太小,没有相交,请问如何延长这个红线,谢谢。
r******n
发帖数: 88
21
来自主题: Mathematics版 - 急问非线性拟合
祝大家新年好!
我得到的扩散偶曲线在界面附近斜率很陡(因为扩散系数依赖于浓度),所以常规的最
小方差拟合(y的最小方差)得到的结果明显不是很好。不知道对非线性拟合还有很么
别的判断标准?总有什么方法定量地描述视觉上的最佳拟合吧?
多谢!
I*M
发帖数: 937
22
【 以下文字转载自 Mathematics 讨论区 】
发信人: ILM (工业光魔), 信区: Mathematics
标 题: 请问Origin怎样对直方图进行双峰高斯拟合?
发信站: BBS 未名空间站 (Wed Jan 6 15:30:55 2010, 美东)
如附件所示。直方图显示了在某个长度区间内所具有的分子数。在40-50和80-90附近有
两个峰。但是用origin右击data选plot details-data-curve里选normal,只能拟合出
一条高斯曲线。怎样同时对这两个峰进行拟合呢?非常感谢!!
z*****4
发帖数: 45
23
借贵宝地请教:
我试图将数据拟合到一个模型中: Signal = a*exp(-R1*t)+b*exp(-R2*t).要拟合的参
数是R1,R2,a和b。 R1和R2更重要。有12个数据点: 12个t,对应12个Signal。我想用
贝叶斯概率理论来求R1,R2,a和b。 我在看贝叶斯分析的逻辑回归模型: https://www
.mathworks.com/help/stats/examples/bayesian-analysis-for-a-logistic-
regression-model.html。 但是,我还没想出如何写the prior distribution, the
posterior distribution, or maximum likelihood for my question。 这个网页也可
能是相关的: https://www.mathworks.com/help/stats/examples/curve-fitting-and-
distribution-fitting.html。
有人可以教我贝叶斯方法的实现数据拟合么,或是分享Matlab cod... 阅读全帖
z*****4
发帖数: 45
24
signal不是正态分布。如果把signal vs t 画出来,是exponential decay。这个问题
其实是把这12个数据点,fit成bi-exponential curve。
来做一下simulation: 假设a=0.5, b=0.5, R1=0.002, R2=0.1,根据已知的12个t,可以
算出12个Signal.然后Signal再加上随机的noise,在matlab可以用randn实现。这样我
们有了12个t,12个Signal(加了noise的),来fit a, b, R1和R2。在matlab里,可用fit
这个function来做。但每运行一次fit,得到的a,b,R1和R2的值,都是不一样的,经常
会和原始值a=0.5, b=0.5, R1=0.002, R2=0.1相差很远。
一篇文章上写着(是一幅图):用Bayesian方法,算出来的R1和R2是一个分布,像是正态
分布,R1的值是横坐标,次数是纵坐标。然后取 出现次数最多的R1最为最后的拟合值
。可惜这篇文章根本不告诉你这个过程是怎么实现的,开头是Bayesian概率论的基本公
式,后面直接上结果。所... 阅读全帖
g****h
发帖数: 481
25
来自主题: Stock版 - 波浪理论的高手有么?
重点是曲线拟合,基本上所有曲线拟合的理论在股市里面都是没用的。
b*******n
发帖数: 1267
26
来自主题: WaterWorld版 - 数学与战争 zt
数学与战争 [ 晨枫 ] 于:2010-09-22 10:11:42 主题帖
军人讲究的是孔武之道,写写算算是文人的事情。在连队里,和数学最接近的大概就是
司务长了,要把粮饷之事计算清楚。但数学和战争其实相距不远。
第一次世界大战前夕,多才多艺的英国人兰切斯特用数学开创了半经验的作战模拟方法
,建立了经典的兰切斯特方程。兰切斯特用平方律定量地解释了特拉法尔加海战中纳尔
逊各个击破的成功诀窍(人称Nelson Touch),恩格尔在1954年用线性律精确地复现了
硫磺岛中美军伤亡情况。经典兰切斯特方程对士气、地形、机动、增援和撤退等没有考
虑,但对战斗的一般规律仍有指导意义。
兰切斯特把战斗简化为两种基本情况:远距离交火杀伤和近距离集中火力杀伤。远距离
交火时,一方损失率既和对方兵力成正比,也和己方兵力成正比。换句话说,敌人越多
,自己损失越大;另一方面,自己人越多,目标越大,损失也越大。这个情况以微分方
程表示即为
dy/dt=-axy
dx/dt=-bxy
其中x和y分别为红军和蓝军的战斗单位数量,a和b分别为红军和蓝军的平均单位战斗力
,因此双方实力相等的条件为
ax=by
即任一... 阅读全帖
M********t
发帖数: 5032
27
来自主题: Joke版 - 引力波 带人类倾听星辰大海之声 新浪科技 ugmbbc 8小时9分钟前 在一波又一波传言后,终于!北京时间2016年2月11日23:40左右,激光干涉引力波天文台(LIGO)负责人、加州理工学院教授David Reitze宣布,人类首次发现了引力波。发布会上究竟揭开了哪些谜团?与之前的传言有何不同吗?中国科学院国家天文台黑洞来客团队(团队主要成员苟利军 研究员为国家天文台恒星级黑洞研究创新小组负责人,而这次被探测到的引力波正是来自双黑洞系统),全面解读这场激动人心的发布会。 正如发布会所言,在被预言将近百年、苦苦追寻几十年之后,首个位于地球之外13亿光年的引力波源GW150914被人类直接探测到,这是一个值得纪念的伟大时刻,一个新时代的序幕正在拉开——地球人,欢迎你来到引力波时代! 激动人心的发布会结束了,其重点内容可以被简单归纳为三点: (1)引力波终于被探测到了。 (2)引力波产生于两个恒星量级黑洞的合并(merger)。 (3)引力波是美国的激光干涉引力波天文台(英文简称LIGO)发现的。 接下来,让我们逐个分析和解释一下以上三点,从而对这项具有划时代意义的科学发现做一个稍微深入的了解: http://static.cnbetacdn.com/article/2016/0212/41080c0cddd85b8.jpg 双黑洞系统想象图。(来源于LIGO网站)) 天外飘来引力波 对于“波”,我们并不陌生,生活中时常会听到无线电波、电磁波、声波、光波等等,引力波也是波的一种。 既然称之为引力波,它必然与引力有关。所以,在更深一步了解引力波之前,我们需要了解一下人类对于引力的认识过程。17世纪末的物理学家牛顿看到了下落的苹果,意识物体之间普遍存在的一种力,称之为“引力”,并且将其数学化,这就是我们熟知的万有引力。万有引力认识的精髓是物体质量的存在导致了引力,这在牛顿之后的两百多年里被认为是宇宙间的绝对真理。直到1905年狭义相对论发表,再到1915年广义相对论的发表,爱因斯坦提出了一种完全不同的对于引力的看法,引力是因为质量对于时空造成了变形所导致,而非质量之间的吸引。这就意味着,我们时空可被当做一种可以变形的介质来认识。所以引力波,简单来说,就是时空自身的波动。相比较我们熟知的无线电波(或者电磁波),它仅仅是在时空之中传播的,时空是它的媒介。 人们常说“星辰大海”,如果将时空视作海洋,那么天体就如同海洋生物一般。可以想象,如果大海中的某个生物摇了摇尾巴、或是晃了晃头,海水由此所产生的波动就会向外传播。与此类似,宇宙中某个天体的剧烈活动,会对所在的时空产生扰动,时空自身的波动也会向远处传播,如果足够强,就能够为地球上的我们所感知。 在引力的世界中,我们的宇宙通常是平静的。可是在北京时间2015年9月14日17点50分45秒,地球上的LIGO探测器却探测到了来自于宇宙深处距离地球13亿光年之外的一场引力风暴,来自于一个双黑洞系统的合并,以它的探测日期命名为 GW150914。 http://static.cnbetacdn.com/article/2016/0212/12c0afa9d0fc4fb.jpg LIGO的两个观测站探测到了同一个引力波事件。上面为观测得到的曲线,下面是和理论相比较之后的拟合结果。(来源于LIGO所发文章) 两个黑洞的“火并” 此次发布会的另外一个亮点就是双黑洞。这也是人们首次直接发现双黑洞,这两个黑洞的质量分别为26和39太阳质量,属于恒星级黑洞。或许你已经听腻了黑洞,生活中时不时的会听到某某黑洞爆发了,某某黑洞吞吃恒星了等等。但是这此发现却有些不同,两个天体都是黑洞,互相绕转,最后合并。这听起来像是一场黑吃黑的火拼,甚至有点儿像港片里熟悉的火爆场面。黑洞合并产生了非常强烈的时空振荡,所以我们遥远的地球才观测到了。 黑洞通常认为是大质量(超过25个太阳质量,请注意这是前身星的质量)的恒星在其演化的最终阶段,恒星中心会形成我们了解的恒星级黑洞。它们的质量通常预计在3个太阳质量到100个太阳质量之间。因为黑洞本身没有任何的辐射(不考虑量子效应下的霍金辐射,它的电磁辐射也是异常微弱),所以我们不能直接看到黑洞。我们的银河系或者其他类似的星系当中当中,每个星系都预计存在着上千万个恒星级黑洞,但是绝大多数的黑洞都是孤独存在的天体,如同幽灵般,没有任何辐射或者辐射及其微弱,所以很难被看到。 所幸,有的黑洞处于双星系统当中,就像发布会中提及的系统,而且另外一个天体是正常的恒星(也称之为伴星)。在这种情形之下,黑洞会从正常恒星上吸积气体,在其周围产生一个吸积盘,以至于某些时候吸积气体的量过多,不能被黑洞直接吞掉,这时还会沿着黑洞的两个转轴将多余的气体抛射出去,从而产生非常壮观的喷流。正是因为吸积盘和喷流的存在,他们都能够产生我们非常熟知的电磁辐射(也就是有我们熟知的光子产生),从而我们利用传统的探测方式,比如地面或者太空的望远镜,就可以间接地探测黑洞的存在。 大约在50年前,人类就是利用此方法发现了第一个黑洞候选体,天鹅座X-1 (Cygnus X-1)。在1974年,地球上两颗聪明的大脑、理论物理学家霍金和好朋友基普·索恩就这个候选体是不是黑洞而打了个赌,他们以一年的成人杂志《阁楼》作为赌注。后来的观测是利用天鹅座X-1中的伴星运动测得了黑洞质量,大约为15个太阳质量,从而霍金认输并且在两人的赌书上签名按上了自己的手印。基普·索恩从那时起就赢了。 对于双黑洞系统,他们几乎不会产生能够为传统方式所观测到的光子。所以,即使它们存在,仅凭借传统的观测方式,我们也无法发现他们。况且,很多的人都怀疑它的存在。但是,在双黑洞绕转,尤其是合并之时,会产生很强的引力波。只要引力波探测器足够灵敏,我们就可以发现它们的踪影。面对大家的怀疑,LIGO的发现用事实证明了它的存在。而此次发现引力波的天文台的创建人之一就是基普·索恩——索恩教授又一次赢了! 引力波为黑洞做名片 在观测到了完整的引力波形之后,利用一种叫做匹配滤波(waveform matching)的方法,理论上就可以推断出系统的性质信息,包括合并之前和之后。比如,对于双黑洞系统,可以推断出合并之前的黑洞质量,自旋和轨道,以及合并之后的质量和自旋。此次新闻发布会中所提及的引力波事件,就得到黑洞的质量在合并之前是26个和39个太阳质量,合并之后是62个太阳质量(合并之前的两个黑洞自旋参数值限制的并不是很好),合并后黑洞是一个克尔黑洞,其自旋参数值为0.67。 你或许会有疑问,合并之后怎么少了3个太阳质量,它跑到什么地方去了?引力波也是携带能量的,在黑洞合并之时,它的形状非常不对称,不是我们看到的单个黑洞的球形,所以在振荡恢复的过程当中,一部分质量就通过引力波的方式辐射出去,从而为我们所接收。合并的时间非常之短,仅仅在大约0.05秒的时间,就将3个太阳质量(大约6.0E30公斤)的能量就通过引力波的方式释放出去,也就是说在一秒钟可以释放出大约10^32公斤的能量。相比较之下,我们的整个宇宙包含了大约1.0E22个太阳,而每个太阳每秒钟向外辐射大约4.0E9公斤的物质能量,所以黑洞合并的最后释放出比整个宇宙每秒钟辐射出的电磁总能量还要高出3倍。 http://static.cnbetacdn.com/article/2016/0212/3aa5437d0cd140c.jpg 双黑洞系统在不同阶段所产生波形随时间演化图。双黑洞系统的演化包括三个阶段:旋近(inspiral), 合并(merger)和铃振(ring down)阶段。(来源于LIGO所发文章) 我们注意到,合并前黑洞的质量为26个和39个太阳质量——这两个质量都比目前已知的恒星级黑洞要重许多。我们现在已经确认了大约20多个恒星级黑洞,最重的恒星级黑洞位于M33 X-7系统中,为16个太阳质量。而其它所有确定的黑洞质量均比这个小,大多数是七八个太阳质量左右。尽管理论之前预言了更大质量的黑洞的存在,但是之前从来还没有发现过,所以此次发现更大质量黑洞,对于黑洞的形成渠道也有着重要影响。 我们所看到的所有恒星,都是处于旋转当中。所以对于那些由恒星塌缩而形成的黑洞而言,它们也是处在旋转之中。对于黑洞的旋转,天文学家用一个叫做“自旋参数”的量来表示,它在0和1之间变化——0代表了没有转动的黑洞(也叫“史瓦西黑洞”),而1代表了理论上黑洞所具有的转动最大值(也叫“极端克尔黑洞”)。对于此次观测中合并以后的黑洞,它的自旋参数为0.67。如果我们以转速来描述,它在以每秒钟100转的速率转动。相较而言,我们刚才提到的人类第一个黑洞“天鹅座X-1”,它是一个极端克尔黑洞,差不多在以超过1000转每秒的速率转动。 如果你觉得转速依旧很难理解,那我们可以想象一下《星际穿越》电影的一个情节,主人公在卡刚都亚黑洞的行星上只呆了3个小时,但是后来却发现飞船上已经过了23年,时间相差了6万倍。要有这样巨大的时间差,其中条件就行星极度靠近黑洞的同时,黑洞也以最大速度转动。按照相对论理论所言,这样行星上的时间就会被极大的拉长。天鹅座黑洞的转速就具有类似于电影中卡刚都亚黑洞那样的转速。对于此次引力波探测到黑洞而言,即使有某个行星在其周围最靠近的稳定存在,那么对于它的时间也会流逝的很慢,不过不会有6万倍那么大的差别,最多也就2倍左右。 对于黑洞而言,有着非常著名的无毛定理,也就是说黑洞只需要简单的几个量就可以描述。对于宇宙当中的黑洞,只需要我们上面所说的质量和自旋参数就可以完整的描述。当我们知道了黑洞的质量和自旋参数一些性质时,我就可以很容易的对黑洞本身的全貌做出一个描绘,就如同给出了一个人的完整自画像。而引力波的方法可以快速给出黑洞的完整信息,这相比较传统的观测方式,更为有效,尽管在观测方面有些困难。 氢原子的百亿分之一 从预言到探测,物理学家等待引力波的到来已有一百年之久,为什么引力波这么难得一见?主要原因是,相比较其他的几种力(强力,弱力,电磁力),引力是最弱的,相应的引力波效应也就很弱。想当初爱因斯坦在发表自己新的理论之后,就估算了引力波的强度。引力波的强度通常利用相对变形大小来表示,但是结果往往是小的可怜,几乎无法探测到。引力波是时空的自身变形,在一个方向上被拉伸,在其垂直的另外一个方向上就会被压缩。如果我们有一天,我们被同样的双黑洞系统在合并时所产生的引力波(变化强度为1.0E-21)所击中的话,理论上来说,我们同样会经历一个稍微变高变瘦,然后变胖变矮一些的过程。实际上,对我们身高不超过2米的人类来说,导致的变化大约为2E-21,为一个氢原子的五百亿分之一(一个氢原子的大小大约为1.0E -10米)。 引力波的效应是如此之小,所以一方面需要增加探测的长度,来增强变化的效应,另外一方面通过巧妙的方法来探测到微小的变化。这也是此次新闻发布会中提到的激光干涉引力波天文台(LIGO)在建造之初所考虑的。它有两个观测点,分别建在美国华盛顿州的列文斯顿和路易斯安那州的汉福德。每个观测点都有两个互相垂直,每条长达4公里的臂。长臂中间是高度真空的管子,而在长臂两端,悬挂着大约直径34厘米、重达40kg的反射镜。LIGO的探测器利用激光干涉技术,不间断地测量每对反射镜之间的距离。每当引力波通过探测器时,人们会探测到两对反射镜之间的距离呈现此消彼长的周期性变化。 即使对于LIGO天文台4公里的长臂,引力波所造成的变化也是极其微小的。对于此次新闻报道中的双黑洞合并,其可能产生的尺度相对变化最高可为1.0E-21,意味着4公里的长度也仅仅只变化了一个氢原子直径的2.5千万分之一。为了达到这个精度,LIGO的科学家做了许多精密的设计,保证探测系统的稳定,保证LIGO反射镜的位置随机涨落小于一个氢原子大小的百亿分之一,从而保证可以相对比较容易的探测到可能的引力波源。 LIGO在1999年建成并且开始运行。但是在进行升级之前(也就是2010年),没有探测到任何确定的引力波事件。从2010年开始,LIGO对探测器进行了第二阶段的升级,2015年6月进行测试运行,2015年9月开始正式运行,第一阶段的科学运行一直持续到2016年1月,升级后的版本被称为Advanced LIGO (简称aLIGO)。而此次新闻发布会的结果其实就是在刚刚升级完之后,由还在进行测试中的的aLIGO所探测到的。相比较之前,aLIGO的探测灵敏度提高了10倍。而且此次的双黑洞所产生的引力波强度就在仅仅比最初LIGO的灵敏度低一些,所以当LIGO的升级刚刚完成,在试运行的阶段就发现了所报道的双黑洞系统。探测到引力波似乎就在本来的预料当中。就像发布会中所言,这或许是大自然给我们苦盼许久的一份礼物。2012年,LIGO天文台创建人基普·索恩在《科学》杂志的一篇评论文章中说,预计在2017年,第一例黑洞合并的事例将会被发现。当时预计LIGO的升级会在2016年底完成,结果是升级进度超前,让我们提前听到了宇宙时空的声音。 http://static.cnbetacdn.com/thumb/article/2016/0212/a4e034c670c7a62.jpg_600x600.jpg LIGO和aLIGO灵敏度比较(左);LIGO和aLIGO探测范围比较(右)。(图片来源于LIGO网站) 谁将撼动时空 那么在我们的星辰大海中,什么样的天体才能够撼动我们的这个宇宙时空,让位于遥远地球上的LIGO探测到呢?现在通常认为有如下四种: (1)旋进(In-spiral)或者合并的致密星双星系统。比如中子星或者黑洞的双星系统。非常类似于发布会当中的系统。 (2)快速旋转的致密天体。这类天体会通过周期性的引力波辐射损失掉角动量,它的信号的强度会随着非对称的程度增加而增加。可能的候选体包括非对称的中子星之类的。在《星际穿越》电影当中,教授说它发现了引力波,而它的其中一个产生机制很可能就是由一个快速转动的中子星,其表面大约2厘米的凸起所产生的(具体分析可以参阅由基普·索恩撰写的星际穿越一书)。 (3)随机的引力波背景。非常类似于我们通常熟知的宇宙背景辐射,这一类背景引力波,也通常叫做原初引力波,它是早期宇宙暴涨的遗迹。2014年由加州理工、哈佛大学等几个大学的研究人员所组成的BICEP2团队曾宣称利用南极望远镜找到了原初引力波,但是后来证实为银河系尘埃影响的结果。原初引力波的探测将是对暴胀宇宙模型的直接验证,对于它的探测依旧在努力寻找之中。 (4)超新星或者伽马射线暴爆发。恒星爆发时非对称性动力学性质也会产生引力波。而直接探测到来自于这些天体的引力波,将是提供对这些天体最直接而且最内部的信息。 除过LIGO之外,还有意大利的VIRGO引力波探测器,日本的KAGRA探测器以及英国德国的联合GEO-HF探测器。相信在不远的未来,引力波的探测事例会如同春笋般爆发,越来越多。 http://static.cnbetacdn.com/article/2016/0212/8dad832e823b89e.jpg 全球引力波天文台分布(来源于LIGO网站) 为一窥宇宙初生 毫无疑问引力波是对广义相对论的一个最直接的验证。另外它在弱场中已经得到验证,但是对已强引力之下的验证,之前却从来没有验证过。所以此次的观测,是对广义相对论的一个非常好的检验。 引力波以光速传播,它与物质的相互作用非常非常的弱,所以引力波可以给我们提供我们宇宙几乎无阻挡的图景,而这个几乎是无法利用我们熟知的电磁波来达到的。比如,利用引力波,我们可以看到宇宙的最早期,宇宙大爆炸之后的1.0E-36秒开始的宇宙形成过程,而对于电磁波而言,它最早只能看到大爆炸后的大约300,000之后的宇宙历史,在此之前,电磁波是不能给我们提供的。所以引力波是我们了解我们宇宙形成的最好工具。 如果还记得,在《星际穿越》电影中的结尾之时,主人公库珀身处一个5维时空的超体方体中,为了将从黑洞中心所提取出来的信息传递给身处4维时空的女儿墨菲,人为的制造引力波效应,成功将信息传递,从而人类得以解救。引力波从目前物理学家的认识来看,是唯一一种可以在不同维度传播的波。不同宇宙之间的碰撞,会产生引力波。说不定在不远的将来,我们也可以依靠引力波来判断多重宇宙的存在与否。 就如同一个天生的聋哑人,一直在听别人说声音的存在,突然有一天听力恢复了。我想我们此刻的心情也是差不多如此。引力波给我们打开了一扇全新的窗口。引力波是一种方式,是一种看待世界的方式。历史的发现轨迹告诉我们,每一扇新的窗口被打开,都会有令人称奇的发现。虽然LIGO的探测能力还是有限,一旦这个引力波的世界被撬开了一道小的裂缝,让我们看到了春天的种子,相信硕果累累的引力波丰收季节也不会太远。
引力波
带人类倾听星辰大海之声
新浪科技
ugmbbc
8小时9分钟前
在一波又一波传言后,终于!北京时间2016年2月11日23:40左右,激光干涉引力波天
文台(LIGO)负责人、加州理工学院教授David Reitze宣布,人类首次发现了引力波。
发布会上究竟揭开了哪些谜团?与之前的传言有何不同吗?中国科学院国家天文台黑洞
来客团队(团队主要成员苟利军 研究员为国家天文台恒星级黑洞研究创新小组负责人
,而这次被探测到的引力波正是来自双黑洞系统),全面解读这场激动人心的发布会。
正如发布会所言,在被预言将近百年、苦苦追寻几十年之后,首个位于地球之外13亿光
年的引力波源GW150914被人类直接探测到,这是一个值得纪念的伟大时刻,一个新时代
的序幕正在拉开——地球人,欢迎你来到引力波时代!
激动人心的发布会结束了,其重点内容可以被简单归纳为三点:
(1)引力波终于被探测到了。
(2)引力波产生于两个恒星量级黑洞的合并(merger)。
(3)引力波是美国的激光干涉引力波天文台(英文简称LIGO)发现的。
接下来,让我们逐个分析和解释一下以上三点,从而对这项具有划时代意义的科学发现
... 阅读全帖
b****o
发帖数: 360
28
来自主题: Stock版 - 靠 TA 炒股
曾经看过一篇文章,有个街上大牛说,每天真正看fa只不过15分钟,大部分时间是做ta
,做模型。
不过话说回来,我们都做过曲线拟合这种东西,历史数据是可以拟合的,未来却是无法
拟合出来的。所以fa 和ta都是无奈强做之举。
m*1
发帖数: 5
29
来自主题: Mathematics版 - 问一个关于函数拟合的问题 (转载)
【 以下文字转载自 Engineering 讨论区 】
发信人: m81 (好好的干活), 信区: Engineering
标 题: 问一个关于函数拟合的问题
发信站: BBS 未名空间站 (Sun Nov 18 13:25:13 2007)
===问题是===
先考虑二维情况,已知一系列点(x1,y1),(x2,y2),...(xn,yn),但是不知道他们背后的
函数是什么。如何在某种误差限度内,找到一种带宽比较小的函数较好地拟合这些点。
所谓带宽较小,举例就是,(1,4) (3,-2) (4,10) 可以找到穿过这三个点的
曲线很多,希望找到穿过y=0坐标轴次数较少的解答。
=========
比较急,在线等,感谢先 :)
I*M
发帖数: 937
30
如附件所示。直方图显示了在某个长度区间内所具有的分子数。在40-50和80-90附近有
两个峰。但是用origin右击data选plot details-data-curve里选normal,只能拟合出
一条高斯曲线。怎样同时对这两个峰进行拟合呢?非常感谢!!
B******u
发帖数: 23763
31
【 以下文字转载自 PKU 讨论区 】
发信人: Benbenxu (悠悠), 信区: PKU
标 题: 学术一下,怎么求三次样条的拟合函数阿?
发信站: BBS 未名空间站 (Wed Oct 19 21:02:09 2011, 美东)
知道一组数值,需要寻找其拟合样条曲线函数。
通过软件最好拉,
比如matlab,或者其他啥。
包子谢哈。
btw,
千万别发我信箱,打不开的。
发我邮箱c*******[email protected]可以
或者版面指导哈。
多谢多谢。
B******u
发帖数: 23763
32
【 以下文字转载自 PKU 讨论区 】
发信人: Benbenxu (悠悠), 信区: PKU
标 题: 学术一下,怎么求三次样条的拟合函数阿?
发信站: BBS 未名空间站 (Wed Oct 19 21:02:09 2011, 美东)
知道一组数值,需要寻找其拟合样条曲线函数。
通过软件最好拉,
比如matlab,或者其他啥。
包子谢哈。
btw,
千万别发我信箱,打不开的。
发我邮箱c*******[email protected]可以
或者版面指导哈。
多谢多谢。
r*******n
发帖数: 131
33
来自主题: Joke版 - 小崔的这些话我赞成
傻逼,老爱提相对论,包括拿诺奖的光电效应,都是实验上出现了当时的理论解释不了
的现象,才催生出来的
理论本质上都不是对的,只是在可观测范围内的一个模型而已,说白了和曲线拟合一样
。你有两个点的时候只能猜是线性的,点多了可能就会发现这是一条二次曲线,多测了
几个点发现好像是个高斯。但只有两个点的时候你在那嚷嚷这是二次曲线是高斯曲线,
可能会被你猜对,这跟猴子在打字机上打出莎士比亚的著作是一回事
最后,不是说要反对文科生,也不是说文科生没资格说话,但是没有科学精神不事实求
是的人的确是没有资格的,这是态度问题不是知识问题
b*****e
发帖数: 762
34
来自主题: Mathematics版 - 这种曲线最初是怎么搞出来的?
(-1 + x^2 + y^2)^3 == x^2 y^3
除了天才能想出来。现在人如果也想做是不是要先画一个桃,然后拟合拟出来的吗?但
是用机器 fit的话,是不是会fit出形
式非常复杂的曲线?还是限制了一个形式然后用机器去fit? 有没有简单的手工的,拟
合出这种曲线,曲面的方法?
c**********9
发帖数: 135
35
来自主题: _Auto_Fans版 - 发动机的共况图是怎么得来的

“曲线是平滑出来的”
曲线是又无数个定转速测出扭矩的点组成?
还有,什么是“曲线拟合”,好像会用到这个。
l**********g
发帖数: 6198
36
是的
我们组一个人做过一个项目,简单地用一下R,然后拟合曲线拟合得很完美啊,很多个
参数,很清晰地显示出各种因素的影响性,很好
对consultant,data analysis,data mining都挺有好处,你们cs的人学一下绝对只有好
处没有坏处的,上课吧,嗯
b*****y
发帖数: 350
37
这马说,胶片并不一定比数码动态范围大,只是在同样范围内的过渡更平滑。相当于统
计中直线拟合和曲线拟合的区别。
b*****y
发帖数: 350
38
这个偶到明白,就是不明白有人讲数码替代不了胶片的说法。通过HDR技术,完全可以
用分段直线拟合逼近曲线拟合,超出人眼的分辨范围之外,那就和胶片没区别了。
v***a
发帖数: 903
39
来自主题: Joke版 - 非笑话:求数值拟合软件
拟合了做什么用?
你知道数据之间大概是啥形式的关系不?我记得matlab 那个 fit 的tool 可以自己加
函数的。
你用样条的话也就能得到光滑的曲线s吧,还能干嘛用?
s*****e
发帖数: 115
40
就说最小二乘法,是用来找“最佳”函数拟合,曲线拟合,等等的最基本的方法。属于
大二线性代数的必讲章节。在业界混你要是讲没听说过,就好像你说自己是中国人又不
晓得北京是不是在河北省一样。
下面这句话引自wiki,关键是十八岁!
"Carl Friedrich Gauss is credited with developing the fundamentals of the
basis for least-squares analysis in 1795 at the age of eighteen"
s*****e
发帖数: 115
41
就说最小二乘法,是用来找“最佳”函数拟合,曲线拟合,等等的最基本的方法。属于
大二线性代数的必讲章节。在业界混你要是讲没听说过,就好像你说自己是中国人又不
晓得北京是不是在河北省一样。
下面这句话引自wiki,关键是十八岁!
"Carl Friedrich Gauss is credited with developing the fundamentals of the
basis for least-squares analysis in 1795 at the age of eighteen"
B******u
发帖数: 23763
42
知道一组数值,需要寻找其拟合样条曲线函数。
通过软件最好拉,
比如matlab,或者其他啥。
包子谢哈。
btw,
千万别发我信箱,打不开的。
发我邮箱c*******[email protected]可以
或者版面指导哈。
多谢多谢。
J*****n
发帖数: 4859
43

我有一些离散的点(x_i,y_i),先用某种曲线拟合的技巧(样条,多项式回归)得出其拟
合函数f(x),然后再把f(x)传到g中,供g使用。
现在的问题是f的参数个数不一定,比如用多项式回归,那么f中不但要传x_i的值,还
要传一个多项式的degree:
double f (vector& x_i, const double& degree);
如果用三次样条插值,虽然不需要degree了,但是除了x_i外还要传一个vector进去用
于存放三次样条的knots(边界点):
double f (vector& x_i, vector& knots);
作为g函数本身,并不知道用户用哪种方法去拟合,也就是不知道f所带的参数是怎么样
的。
请教该如何解决这个问题。
谢谢。
c*u
发帖数: 916
44
来自主题: Computation版 - 怎么拟合这样的曲线
我有一组实验数据,
y = y_high, if x > thr
y = y_low, if x < thr
y_high, y_low, thr 是未知的常数,thr是所需要得到的东西.好像3D bruteforce
search 寻找最低点不是很有效率.
举个通俗的例子,假设x是温度,y是水的体积,那么Thr是熔点.但是实验数据中
y的测量有很大noise, 所以必须拟合.
谢谢,
g******s
发帖数: 733
45
x=[0.1 0.2 0.4 0.6 0.8 1 2 4 6 8 10 20 40 60 80];
y=[200 100 50 35 28 22 12 6.6 5 4 3.6 2.5 1.9 1.8 1.8];
loglog(x,y,'o-')
set(gca,'FontSize',14)
请问什么函数(不能分段)可以拟合附图中的曲线?
先谢了!
c******k
发帖数: 1140
46
画的这一条曲线类似于三段直线连接在一起,但是他们有不同的斜率,就是想知道每一
段的斜率,我用的是Origin 8, 在 Analysis-->Fitting-->Fit Linear, 本想分段拟
合直线,但是每次都是把三段直线当成整体,作的直线拟合,指点一下,多谢了!!!
m*1
发帖数: 5
47
来自主题: Engineering版 - 问一个关于函数拟合的问题
===问题是===
先考虑二维情况,已知一系列点(x1,y1),(x2,y2),...(xn,yn),但是不知道他们背后的
函数是什么。如何在某种误差限度内,找到一种带宽比较小的函数较好地拟合这些点。
所谓带宽较小,举例就是,(1,4) (3,-2) (4,10) 可以找到穿过这三个点的
曲线很多,希望找到穿过y=0坐标轴次数较少的解答。
=========
比较急,在线等,感谢先 :)
c*u
发帖数: 916
48
来自主题: Mathematics版 - 这种分段曲线怎么拟合?
我有一组实验数据,理论上应当符合以下规律:
y = y_high, if x > thr
y = y_low, if x < thr
y_high, y_low, thr 是未知的常数,thr是所需要得到的东西.好像3D bruteforce
search 寻找最低点不是很有效率.
举个通俗的例子,假设x是温度,y是水的体积,那么Thr是熔点.但是实验数据中
y的测量有很大noise, 所以必须拟合.
谢谢,
D*****i
发帖数: 8922
49
中国古代的代数很发达,郭守敬就用三次、四次曲线拟合天体运行轨道,这种东西现代
只有到了大学讲样条曲线时才会讲到。一般文科生对这个根本就没有概念,个别的傻逼
文科生居然还在这个讨论中妄谈什么才是数学。
r*****g
发帖数: 9999
50
来自主题: Military版 - 哈哈 阿里造假被扒皮了
阿里的数据绝对没有造假,但也不能算100%纯天然,之所以曲线拟合的这么完美,是因
为阿里有指哪打哪的能力。阿里每年双11的销售额是有非常具体的目标的,这些目标很
有可能就是用一个曲线方程设定的,剩下的就是如何达到这个目标。阿里有无数种方法
确保在真实交易的前提下下完成目标。哥有几个朋友就是所谓的天猫上的KA大卖,每年
销售额几十个亿级别,阿里对这些大卖都下达了具体的销售目标,在双11当天更是每小
时check一次,如果与目标有差距的话,店小二直接打电话,说下一个小时如果不能达
标立刻从首页拿掉,而且未来也别想进入天猫主推位置,这帮大卖就算赔本也会保证销
售业绩的,所以阿里根本犯不着造假,不是真实交易的话,造假的数据一样要交税,阿
里不至于傻到自己赔钱赚吆喝,但下面的买家是不是赔钱赚吆喝就不知道了。
1 2 3 4 5 下页 末页 (共5页)