svd小结（包含程序，测试数据，数据说明，检验程序和说明）

Legendary

大神，我有两个小问题。
1.在做蒙特卡洛检验的时候，你引用的图片 （魏凤英老师的书上给出的公式）上方差贡献Ck=入/（入的累加），但是此书前一页上公式（8.81）却说，每对奇异向量的方差贡献率=入平方/(入平方的累加)，给出的程序貌似也是算的带平方的。这是公式错了吗？
2.如果选用5个主模态，蒙特卡罗检验会在  tezhengzhi_mengtekaluo.txt  输出500个U,是不是每100个排一次大小顺序，然后选取第95或99个U，跟对应模态的方差贡献率进行比较？

770809

看来很不错，下载下来学习

泰达巴神

下来试试，最近也在做类似的分析，希望多多交流

更大华

l楼主辛苦

心语

楼主 弱弱的问句 一般实矩阵进行SVD分解  矩阵值的数量级波动较大 在10的-6到10的-13次方之间   这种矩阵能进行SVD分解么？

家肥猫

小蝌蚪 发表于 2013-10-1 17:02
楼主，国庆快乐啊~今天在研究SVD显著性检验的问题，有一些问题不太明白，想请教一下。你的分享里用svd_test ...

小蝌蚪，你之前提出的问题解决了吗？我也很想知道怎么办，迫切想知道，这个检验怎么用呀{:5_207:}{:5_207:}

家肥猫

斑主，我特别想知道关于那个检验程序要修改什么呀？比如Source1.f90中需要修改什么呢？还有svd_test1.f中是不是只修改一下文件路径和名称就可以了呢？还有最后生成的tezhengzhi_mengtekaluo.txt 里的结果要怎么用，才能知道我做的svd是否通过检验了呢？{:5_275:}{:5_275:}拜托拜托解答一下吧，真心是菜鸟，不要嫌弃啊

小蝌蚪

家肥猫 发表于 2014-9-16 17:40
小蝌蚪，你之前提出的问题解决了吗？我也很想知道怎么办，迫切想知道，这个检验怎么用呀{:5_207 ...

我后来就没再做检验了。现在看来，你编个fortran小程序把tezhengzhi_mengtekaluo.txt把这100个试验中同一模态的方差贡献从大到小排列，实测资料做SVD得到的奇异值的方差贡献就用svd_info.txt里面 The singular values:下面的几个奇异值粗略代表所有奇异值之和，再根据魏凤英的计算方法计算就好了。

家肥猫

小蝌蚪 发表于 2014-9-17 09:39
我后来就没再做检验了。现在看来，你编个fortran小程序把tezhengzhi_mengtekaluo.txt把这100个试验中同一 ...

魏凤英的哪种计算方法哟？咋比较的？

家肥猫

Legendary 发表于 2014-9-7 17:58
大神，我有两个小问题。
1.在做蒙特卡洛检验的时候，你引用的图片 （魏凤英老师的书上给出的公式）上方差 ...

“2.如果选用5个主模态，蒙特卡罗检验会在  tezhengzhi_mengtekaluo.txt  输出500个U,是不是每100个排一次大小顺序，然后选取第95或99个U，跟对应模态的方差贡献率进行比较？”
是这样比较吗？？？？