应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 用matlab实现一个偏最小二乘法运算
211/3返回列表123下一页
查看: 2707  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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矩阵, 重复上述手续计算潜变量矩阵的第二列, 直列所选定的主因子
数为止。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-10-20 13:42:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

Rachellhom

铜虫 (小有名气)
引用回帖:
8楼: Originally posted by yijincao0123 at 2012-10-25 12:08:48
具体什么方法？...

你下载一个SPSS的软件。这个软件可以用于很多数据分析。然后你再在网上下载一个PLS的插件。这样的话，你就可以进行偏最小二乘的运算了。
赞一下(1人) 回复此楼
9楼2012-10-26 10:40:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

Rachellhom

铜虫 (小有名气)
假如我表述的那个原理看不的太清楚的话，麻烦大家查一下这篇文献：甘子琼等 偏最小二乘—分光光度法同时测定重油中的铁、镍、钒。
这个里面有详细的原理。。。。。。我实在不懂这个。拜托大家啦。。。。。。。感激涕零。。。。。。。
赞一下 回复此楼
2楼2012-10-20 13:44:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazewei

木虫 (小有名气)
感谢参与，应助指数 +1
xiegangmai: 应助指数-1, 非应助请勿选择“应助回帖” 2012-10-24 19:48:57
我也不会，太复杂了
赞一下 回复此楼
3楼2012-10-22 15:02:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Rachellhom

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by mazewei at 2012-10-22 15:02:05
我也不会，太复杂了

55。亲，我都抑郁了。。你能帮我找个会的人么？帮帮我。。。大侠。
赞一下 回复此楼
4楼2012-10-22 17:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yijincao0123

新虫 (初入文坛)
同求
回复此楼
没人逼你努力，但不能抱怨不如人
5楼2012-10-24 13:42:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Rachellhom

铜虫 (小有名气)
引用回帖:
5楼: Originally posted by yijincao0123 at 2012-10-24 13:42:07
同求

可能因为我这个太复杂了。不见有人回复。不过我最近找到另一种SPSS来完成这个PLS分析。你也可以试试，不需要编程什么的。
赞一下 回复此楼
6楼2012-10-25 10:35:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazewei

木虫 (小有名气)
7楼2012-10-25 10:40:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yijincao0123

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by Rachellhom at 2012-10-25 10:35:17
可能因为我这个太复杂了。不见有人回复。不过我最近找到另一种SPSS来完成这个PLS分析。你也可以试试，不需要编程什么的。...

具体什么方法？
回复此楼
没人逼你努力，但不能抱怨不如人
8楼2012-10-25 12:08:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yijincao0123

新虫 (初入文坛)
引用回帖:
9楼: Originally posted by Rachellhom at 2012-10-26 10:40:08
你下载一个SPSS的软件。这个软件可以用于很多数据分析。然后你再在网上下载一个PLS的插件。这样的话，你就可以进行偏最小二乘的运算了。...

PLS插件太难找了，你要有的话发我一份，感激不尽
赞一下 回复此楼
没人逼你努力，但不能抱怨不如人
10楼2012-11-06 13:42:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Rachellhom 的主题更新
211/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 317求调剂 +8 申子申申 2026-03-19 14/700 2026-03-21 16:29 by 我的船我的海
[考研] 求调剂 +4 十三加油 2026-03-21 4/200 2026-03-21 16:20 by barlinike
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 15/750 2026-03-21 10:16 by 脱颖而出
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 311求调剂 +5 冬十三 2026-03-18 5/250 2026-03-21 00:16 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 328求调剂，英语六级551，有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 有没有道铁/土木的想调剂南林，给自己招师弟中～ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
信息提示
关闭
请填处理意见
关闭