| 查看: 2599 | 回复: 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矩阵, 重复上述手续计算潜变量矩阵的第二列, 直列所选定的主因子 数为止。 |
回复此楼» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
博士申请都是内定的吗？
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有12人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
15楼2013-10-08 16:33:19
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