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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 用matlab实现一个偏最小二乘法运算
211/3返回列表123下一页
查看: 2760  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +17 努力奋斗112 2026-04-06 20/1000 2026-04-11 00:31 by wangjihu
[考研] 284求调剂 +19 梵@@ 2026-04-06 21/1050 2026-04-10 21:12 by zhouxiaoyu
[考研] 材料工程日语考生求调剂 +5 0856?调剂 2026-04-10 5/250 2026-04-10 20:45 by yukihao
[考研] 材料工程085601，270求调剂 +28 @ASDF1234 2026-04-08 30/1500 2026-04-10 19:28 by dick_runner
[考研] 本9一志愿2 0854低分专硕286求调剂 +10 芒种111 2026-04-04 10/500 2026-04-10 12:31 by luosha500
[考研] 调剂 +24 不逢春 2026-04-07 25/1250 2026-04-10 09:39 by 690616278
[考研] 化工求调剂！ +34 RichLi_ 2026-04-06 34/1700 2026-04-09 20:39 by zhouxiaoyu
[考研] 312求调剂 +3 李鸿飞飞 2026-04-06 3/150 2026-04-09 17:32 by wp06
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
[考研] 327求调剂 +10 Xxjc1107. 2026-04-06 11/550 2026-04-09 01:21 by lature00
[考研] 材料科学与工程320求调剂，080500 +12 黄瓜味薯片 2026-04-06 12/600 2026-04-08 16:26 by luoyongfeng
[考研] 274求调剂求调剂 +10 Jachenbingoo 2026-04-06 13/650 2026-04-08 14:25 by zhq0425
[考研] 288环境专硕,求调材料方向 +35 lllllos 2026-04-04 39/1950 2026-04-07 23:24 by 一只好果子?
[考研] 22408 一志愿双一流人工智能300分 四六级，数据分析国奖 +4 zzfeng123 2026-04-06 6/300 2026-04-07 21:02 by zzfeng123
[考研] 316求调剂 +7 yyx想调剂 2026-04-05 7/350 2026-04-07 14:31 by shdgaomin
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
[考研] 307求调剂 +3 所念及所望 2026-04-06 3/150 2026-04-06 17:30 by 土木硕士招生
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-05 5/250 2026-04-06 15:40 by lin-da
[考研] 一志愿沪9，求生物学调剂，326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
信息提示
关闭
请填处理意见
关闭