经管之家App
让优质教育人人可得
立即打开
全部版块
我的主页
›
论坛
›
数据科学与人工智能
›
数据分析与数据科学
›
R语言论坛
求助!!为什么R语言计算状态空间模型的波动矩阵算出来不对称~
楼主
少年的乌托邦
827
2
收藏
2021-09-03
如题,根据卡尔曼滤波计算出来的矩阵第二列的数据小数点位数比其他列的少一位,做乔氏分解一直显示非正定,有哪位大神可以指点一下么~感激不尽!!
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
全部回复
沙发
llb_321
2021-9-4 08:32:05
是要计算马氏距离吗?
第二列的问题,确实很奇怪,不知道你用什么函数做的。
library(Rmpfr)
prec=8*log2(10) #浮点精度8位
#比较一下两种计算的结果
exp(2)
#[1] 7.389056
exp(mpfr(2, prec))
#1 'mpfr' number of precision 26 bits
#[1] 7.38905609
你参考上面的例子,浮点数指定精度后计算,看看出来的结果是不是可以进行Cholesky分解。
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
藤椅
少年的乌托邦
2021-9-4 09:33:12
llb_321 发表于 2021-9-4 08:32
是要计算马氏距离吗?
第二列的问题,确实很奇怪,不知道你用什么函数做的。
不是,是卡尔曼滤波,状态空间模型中需要计算波动V_t_T这个矩阵,但是R算出来的结果不对称,小数点位数第二列总是少一位,没法做乔氏分解。而且有的矩阵计算的没问题,有的矩阵却不行,太玄乎了。<br>
这位大神,问下这是什么情况
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
栏目导航
R语言论坛
外文文献专区
文献求助专区
经管高考
哲学与心理学版
区域经济学
热门文章
【全美经典】离散数学
understanding climate change perceptions ...
【同程商旅】中国企业出海差旅研究报告
“十四五”能源发展成就报告
智算无界AIDC的超越和重构2025
2000离散数学习题精解
当社科基础理论重大理论发现的时候
【24重磅,自用整理!】2000-2024上市公司投资 ...
2025年我国医药航空冷链发展现状与趋势展望 ...
中国数字经济规模数据、报告(2005-2023年) ...
推荐文章
AI狂潮席卷学术圈,不会编程也能打造专属智 ...
10月重磅来袭|《打造Coze/Dify专属学术智能 ...
最快1年拿证,学费不足5W!热门美国人工智能 ...
关于如何利用文献的若干建议
关于学术研究和论文发表的一些建议
关于科研中如何学习基础知识的一些建议 (一 ...
一个自编的经济学建模小案例 --写给授课本科 ...
AI智能体赋能教学改革: 全国AI教育教学应用 ...
2025中国AIoT产业全景图谱报告-406页
关于文献求助的一些建议
说点什么
分享
微信
QQ空间
QQ
微博
扫码加好友,拉您进群
各岗位、行业、专业交流群