| 查看: 3086 | 回复: 20 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
Rachellhom铜虫 (小有名气)
|
[求助]
用matlab实现一个偏最小二乘法运算
|
|
|
各位大侠，我想用matlab写一个最小二乘法的运算代码。 需要实现以下运算： 原理是： 设配制n组标准溶液, 其中含有q种不同量的组分. 这q种组分的吸光光谱曲线严重重叠, 无法单独测定. 在P个波长测定体系的吸光度, 构成浓度矩阵Ynxq和吸光度矩阵Xnxp. 将吸光度矩阵Xnxp与浓度矩阵Ynxq进行平均中心化和归一化, 按照偏最小二乘法(PLS), 将矩阵Xnxp与Ynxq分解成隐变量矩阵T, U与载荷矩阵P, Q. 即: X=TP+ E Y= UQ+ F 其中E, F分别为残差矩阵, Unxa与Tnxa之间可用下式关联: U= BT U, T中的列, P, Q中的行都分别相互正交, 维数a 在统计意义上要使T、P尽可能接近X. 迭代过程中导入一权重矢量W, 迭代循环可按下列公式分步进行: Waxp= U（转置）X / U（转置）U1 Waxp= W/ M1 t=XnxpWpxa/ W（转置）W1 q（转置）=t（转置）Ynxq / t（转置）t. Unxa= Yq/ q（转置）q. 其中M指欧几里德范数. 当+t新- t旧+/ +t新+小于指定的收敛标准时, 迭代收敛. 退出迭代循环后进行下列运算: Paxp（转置）= t（转置）X / t（转置）t Enxp= Xnxp—t P（转置） b=U（转置）t / t（转置）t Fnxq=Ynxq — btq（转置）. 括号里面的转置出现在哪个字母后面就是代表那个矩阵的转置, 大写字母表示矩阵, 小写字母表示列矢量. 将计算所得的各潜变量矩阵第一列贮存. 用新算出的E与F矩阵当作修正的X与Y矩阵, 重复上述手续计算潜变量矩阵的第二列, 直列所选定的主因子 数为止。 |
回复此楼» 猜你喜欢
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有109人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
2027广东省杰青
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 递推最小二乘法在什么书本中有详细的讲解和matlab仿真介绍？
已经有5人回复
- matlab用多元线性最小二乘法求解多元线性方程的问题
已经有11人回复
yijincao0123
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 189
- 帖子: 34
- 在线: 18.3小时
- 虫号: 1892872
- 注册: 2012-07-14
- 性别: MM
- 专业: 环境工程
8楼2012-10-25 12:08:48
Rachellhom
铜虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1095.4
- 散金: 65
- 红花: 4
- 帖子: 221
- 在线: 91.1小时
- 虫号: 2054481
- 注册: 2012-10-11
- 性别: MM
- 专业: 环境工程
2楼2012-10-20 13:44:32
mazewei
木虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1754.4
- 帖子: 251
- 在线: 89.6小时
- 虫号: 1783394
- 注册: 2012-04-27
- 性别: GG
- 专业: 机械测试理论与技术
3楼2012-10-22 15:02:05
Rachellhom
铜虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1095.4
- 散金: 65
- 红花: 4
- 帖子: 221
- 在线: 91.1小时
- 虫号: 2054481
- 注册: 2012-10-11
- 性别: MM
- 专业: 环境工程
4楼2012-10-22 17:51:26
