经管之家App
让优质教育人人可得
立即打开
全部版块
我的主页
›
论坛
›
数据科学与人工智能
›
数据分析与数据科学
›
R语言论坛
求助!!为什么R语言计算状态空间模型的波动矩阵算出来不对称~
楼主
少年的乌托邦
911
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语言论坛
经管文库(原现金交易版)
经管高考
金融学(理论版)
教师之家与经管教育
市场行情分析
热门文章
北美PHD,焦虑的一批~
通用指标与场景指标:CDA数据分析师的核心分 ...
2024年合集 ESG评级数据大全(彭博 华证 Wi ...
技术趋势2026
人工智能赋能应用实践指南
CAS EXAM6 2023版TIA教材
中国企业高质量发展展望穿越周期与聚力创新
中国企业高质量发展展望穿越周期与聚力创新
人工智能赋能应用实践指南
多模态大语言模型技术发展报告
推荐文章
2026JG学术冬训营:从Stata初高到Python机器 ...
【必看】【本版版规,欢迎发悬赏贴求助】
【新课】26年3月|Gemini辅助论文写作与数据 ...
关于如何利用文献的若干建议
关于学术研究和论文发表的一些建议
关于科研中如何学习基础知识的一些建议 (一 ...
一个自编的经济学建模小案例 --写给授课本科 ...
AI智能体赋能教学改革: 全国AI教育教学应用 ...
2025中国AIoT产业全景图谱报告-406页
关于文献求助的一些建议
说点什么
分享
微信
QQ空间
QQ
微博
扫码加好友,拉您进群
各岗位、行业、专业交流群