关于mathematica拟合拟合优度 相关系数数的问题，哪个大神能给讲解下不

来源：蜘蛛抓取(WebSpider) 时间：2017-06-11 06:01 标签： 线性拟合相关系数
小木虫 --- 600万学术达人喜爱的学术科研平台
&&查看话题
【求助】用mathematica怎么拟合这一组数据啊？
x& && && &&&1& && && && && && &2& && && && && &4& && && && &&&6& && && && & 15& && && && & 24
y& && &12.08& && && &18.6& && && & 16.84& && && & 15.24& && &&&7.6& && && &&&4.08
拟合函数为y=a*，其中a,b,c为要拟合的参数，哪位知道的话请赐教一下，谢谢！
北京学而思教育科技有限公司 地址：北京市海淀区北三环甲18号中鼎大厦A座1层102室 电话：010-1119人阅读
时间序列分析（5）
[时间序列分析][3]--自相关系数和偏自相关系数
之前在回归分析里面曾经讲过协方差和相关系数，
这里再多讲一句，协方差是会受到单位的影响的，而相关系数就是消除了量纲的影响，来看两者的相关性
这里讲的自相关系数可以说是根据最原始的定义引伸出来的。
下面分别讲一下我对自相关系数和偏自相关系数的理解：
自相关系数
其实自相关系数可以这么理解：把一列数据按照滞后数拆成两列数据，在对这两列数据做类似相关系数的操作。
看一个例子：
这组数据是求滞后数为2的自相关系数，则变成求{x1,x2,...,x8}和{x3,x4,...,x10}两者的“相关系数”，相关系数打引号是因为这个相关系数的公式和以往的有点不一样。下面看一下公式的对比：
要注意的是在计算自相管系数的时候&是使用的总体的均值, 可以看到他们除了&取得不一样,
几乎就是一样的。
所以，我们可以这么理解自相关系数, 她就是用来表达一组数据前后数据 (自己和自己) 的相关性的
在mathematica中，求自相关系数的函数为&CorrelationFunction[]
偏自相关系数
偏自相关系数在网上能查到的很少，我就详细的讲一下。
首先是定义：
但是上面这个式子不能进行计算，我们经过化简，可以得到下面的等价的式子：下面矩阵中的pi就是滞后为i的自相关系数
至于化简的过程，可以查阅一下相关的资料，用到了k阶自回归拟合，还是有点复杂的。
我们可以将上面的过程用mma实现，当然mma中是有现成的函数的，我们就全当验证一下公式是否正确。
我们来计算一下{2,3,4,3,7}的滞后系数为3的偏自相关系数
1.首先计算他的1，2，3阶滞后的自相关系数
xs = CorrelationFunction[{2, 3, 4, 3, 7}, #] & /@ {1, 2, 3}
2.接着生成如上的k*k的矩阵D和对于的Dk
x = Array[
CorrelationFunction[{2, 3, 4, 3, 7}, Abs[#1 - #2]] &, {3, 3}];
x // MatrixForm
xk[[All, 3]] =
xk // MatrixForm
3.计算Dk/D
PartialCorrelationFunction[{2, 3, 4, 3, 7}, 3]
Det[xk]/Det[x]
上面的过程其实可以帮助我们更好的理解偏自相关系数的计算，我们把上面的过程总结成一个函数
pcorr[h_, list_] := Block[{xs, x, xk, lh},
lh = Length[list];
xs = CorrelationFunction[list, #] & /@ Range[lh - 1];
x = Array[CorrelationFunction[list, Abs[#1 - #2]] &, {h, h}];
xk[[All, h]] = xs[[;; h]];
Print[&D矩阵: &, MatrixForm[x]];
Print[&Dk矩阵: &, MatrixForm[xk]];
Print[&使用自编函数: & && ToString[N@Det[xk]/Det[x]]];
Print[&使用系统函数: & && ToString[N@PartialCorrelationFunction[list, h]]];
这样在计算偏自相关系数的时候可以返回两个矩阵D和Dk，我们看一下效果
可以看到两者计算的结果是一样的，并且输出了两个矩阵。
最后欢迎关注我的微信公众号
公众号里除了会更新与博客上一样的文章外，还会增加一些我自己喜欢的电影分享之类的，欢迎大家关注。
以上，所有。
&&相关文章推荐
参考知识库
* 以上用户言论只代表其个人观点，不代表CSDN网站的观点或立场
访问：24596次
积分：1177
积分：1177
排名：千里之外
原创：91篇
评论：20条
欢迎大家关注我的公众号:
prettymath
不仅为你推荐最新的博文，还有更多惊喜和资源在等着你!
阅读：4465
(6)(1)(3)(6)(14)(9)(10)(8)(13)(7)(22)推荐这篇日记的豆列
······求助大神，mathematica调用矢量分析包的问题_mathematica吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：11,342贴子：
求助大神，mathematica调用矢量分析包的问题收藏
程序如下：Clear[&Global`*&];V = q/(4 Pi*\[Epsilon]0)*(1/r1 - 1/r2)cond = {r1 -& R - l/2*Cos[\[Theta]], r2 -& R + l/2*Cos[\[Theta]]}V = V /. cond // SimplifySeries[V, {l, 0, 6}] // NormalV1 = pCos[\[Theta]]/(4 \[Pi] (R^2) \[Epsilon]0 )Solve[V1 == c, R]到这儿都对，结果是{{R -& -(Sqrt[pCos[\[Theta]]]/(
2 Sqrt[c] Sqrt[\[Pi]] Sqrt[\[Epsilon]0]))}, {R -& Sqrt[
pCos[\[Theta]]]/(2 Sqrt[c] Sqrt[\[Pi]] Sqrt[\[Epsilon]0])}}调用适量分析包，错了，提示是Get::noopen: 无法打开 Calculus`VectorAnalysis`. &&&&Calculus`VectorAnalysis`Ee=Grad[-V1,Spherical[R,\[Theta],\[Omega]]]希望大神能够帮忙一下，万分感谢
【618有好戏】
直接用Grad就行了
登录百度帐号推荐应用
为兴趣而生，贴吧更懂你。或

我要回帖

更多关于 线性拟合相关系数 的文章

 

随机推荐