24小时热门版块排行榜    

查看: 2553  |  回复: 18
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

tjuxh

金虫 (小有名气)

[求助] POD正交分解求助 已有2人参与

诸位虫友,我最近在想用POD正交分解来分解一个随时间变化的压力场,想按能量分解然后看起压力模态,但是我是机械专业的,木有学过概率论,但是最近看文献发现大多数通过Karhunen-Loeve 展开来进行压力场的分析和简化的,但是不知其具体的原理,与具体的操作方法。所以向各位虫虫求助,推荐相关的学习资料,最好是那种深入浅出的,学一段时间就能上手做的,谢谢啦~~~
回复此楼
一切皆缘分
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
tjuxh: 金币+5, ★★★很有帮助 2016-01-11 09:41:47
引用回帖:
3楼: Originally posted by tjuxh at 2015-12-02 23:17:47
你好,非常感谢你的回答!我看了一些相关的例子,确实是按照您所说的这样进行分解。目前的情况是,我可以提取出随时间变化的一系列随机向量,然后据此可以得出一个协方差矩阵,但是我不知道如何通过此矩阵求出一组 ...

最近又有了一些新收获,跟楼主交流一下。
KL、POD、PCA思想都是一样的,最根本的应该是KL,离散化之后就是PCA(主成分分析)。和对称方阵的特征值分解以及一般矩阵的SVD(奇异值分解)关系非常密切。
1. 首先以对称方阵的特征值分解为例来说明。(特征值按从大到小排列。由于对称阵的不同特征值的特征向量正交,故可将特征向量矩阵S标准正交化。)



a) 基变换角度。
AX的变换效果从几何上可理解为将X旋转(S'作用)——拉伸(作用)——旋转回来(S作用)。因此,S中的每一个列向量可以理解为原始变量空间(实际上为A的行空间)的基通过旋转(因为S为正交阵)所得的一组新基。第一个基s1上包含的信息最多(特征值最大),之后的依次减少。所以任给一个d维列向量X,可通过在矩阵S代表的基上投影,依次提取其第一大分量、第二大分量、……直到第d大分量。根据实际需要,选择前若干个就可以充分表达X所包含的信息了,如果选择全部d个分量,那就是完全无损表示了。数据压缩、图像处理、数据降维等都是这样的原理,提取前几个较大的分量就可以以较小的损失表示原来的信息。
具体来说,可以定义新的列向量:


即Y的每一个分量都是X分量的线性组合,或者理解为将X沿着特征向量展开。选择前r个较大特征值对应的Y分量,就可以将原始的d维向量X降为r维向量Y。(注:A的特征向量又称作主成分,主成分分析PCA大概就是这个意思,提取最主要的那几个成分。)

b) 也可以从对称矩阵的谱分解角度理解

即对称矩阵可表达为d个秩1矩阵之和,每个秩一矩阵的权重为相应特征值。数据压缩就是选择前少数几个较大特征值对应的秩一矩阵来表示原矩阵A,即


2. SVD角度。
假设数据有d个特征或属性,例如d维变量X。假设观测N次,得N个样本,用N×d矩阵A表示,将A进行SVD分解

V中所含的A的右奇异向量就是原始数据的主成分,根据奇异值大小排列,也称作第一主成分、第二主成分、……第d主成分。之所以这么叫,是因为A的右奇异向量实际上是A’A的特征向量

而A’A又可看作d维变量X的协方差矩阵cov(X, X)。如果X均值不为0,可通过平移得到新的零均值变量。令 R = cov(X, X) ,则有

这里的矩阵V相当于第1节中的矩阵S。前r个特征向量表示对X方差贡献最大的r个基,即前r个主成分。类似的,可以定义新向量Y = V'X,提取主成分,压缩数据等等。

3. 小结
1) 将观测矩阵A进行处理,比如平移使得其变为零均值、标准化处理等;
2) 对A'A进行特征值对角分解;
3) 定义新向量:Y=V'X。根据处理能力及实际需要,选择前若干个较显著的Y分量。

4. 参考资料
浅谈对主成分分析(PCA)算法的理解 http://blog.csdn.net/passball/article/details/24037593
机器学习中的数学(4)-线性判别分析(LDA), 主成分分析(PCA) http://www.cnblogs.com/LeftNotEa ... chine-learning.html
机器学习中的数学(5)-强大的矩阵奇异值分解(SVD)及其应用 http://www.cnblogs.com/LeftNotEa ... d-applications.html
知其然,知其所以然。
11楼2016-01-09 14:29:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
tjuxh: 金币+5, 有帮助 2015-12-06 15:08:26
KL展开就是把一个随机场分解为若干随机系数和确定性函数的线性组合,随机系数就是一些随机变量,确定性函数需要求随机场的协方差函数的特征值和特征函数,大概思路是这样。
资料我也没有合适的,网上找找吧,挺经典的方法了。
知其然,知其所以然。
2楼2015-12-01 12:55:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjuxh

金虫 (小有名气)

引用回帖:
2楼: Originally posted by mygt_hit at 2015-12-01 12:55:10
KL展开就是把一个随机场分解为若干随机系数和确定性函数的线性组合,随机系数就是一些随机变量,确定性函数需要求随机场的协方差函数的特征值和特征函数,大概思路是这样。
资料我也没有合适的,网上找找吧,挺经典 ...

你好,非常感谢你的回答!我看了一些相关的例子,确实是按照您所说的这样进行分解。目前的情况是,我可以提取出随时间变化的一系列随机向量,然后据此可以得出一个协方差矩阵,但是我不知道如何通过此矩阵求出一组正交的特征函数,并由这组正交特征函数得出相对应的随机系数,不知可否赐教~或者指明一个方法,谢谢啦~~~
一切皆缘分
3楼2015-12-02 23:17:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by tjuxh at 2015-12-02 23:17:47
你好,非常感谢你的回答!我看了一些相关的例子,确实是按照您所说的这样进行分解。目前的情况是,我可以提取出随时间变化的一系列随机向量,然后据此可以得出一个协方差矩阵,但是我不知道如何通过此矩阵求出一组 ...

具体我也不清楚了,理论上要解一个积分方程,求出特征函数和特征值,好像有一些数值方法。

发自小木虫Android客户端
知其然,知其所以然。
4楼2015-12-02 23:58:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 23/1150 2026-08-13 00:35 by 要乐观耀哥
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 关于Filecode分析方法 +8 majunge000 2026-08-10 9/450 2026-08-12 22:17 by 虫友是什么虫
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +10 Tide man 2026-08-11 10/500 2026-08-12 20:59 by 陈晨晨陈啊
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
请填处理意见