24小时热门版块排行榜    

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

Rachellhom

铜虫 (小有名气)

引用回帖:
8楼: Originally posted by yijincao0123 at 2012-10-25 12:08:48
具体什么方法?...

你下载一个SPSS的软件。这个软件可以用于很多数据分析。然后你再在网上下载一个PLS的插件。这样的话,你就可以进行偏最小二乘的运算了。
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求助调剂,跨调 +6 X十甫寸Y 2026-04-11 6/300 2026-04-11 10:43 by qingpingzhu
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 284求调剂 +19 梵@@ 2026-04-06 21/1050 2026-04-10 21:12 by zhouxiaoyu
[考研] 085404 298分求调剂 +10 呼啦呼啦呼呼呼 2026-04-10 11/550 2026-04-10 16:44 by wangy0907
[考研] 071000生物学调剂求助 +17 zzzzwww 2026-04-09 20/1000 2026-04-10 15:55 by 求调剂zz
[考研] 一志愿中南大学物理学,英一66,求调剂 +4 长烟旖旎 2026-04-08 5/250 2026-04-10 10:31 by 颖果儿
[考研] 调剂 +24 不逢春 2026-04-07 25/1250 2026-04-10 09:39 by 690616278
[考研] 085400电子信息类(川大控制工程)求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 一志愿武理车辆 281 求调剂 +5 上岸研究生. 2026-04-07 5/250 2026-04-09 15:56 by only周
[考研] 求调剂 +8 吃口冰激凌 2026-04-07 8/400 2026-04-09 08:03 by 5268321
[考研] 0703化学调剂 348分 +14 唉我超真没招了 2026-04-06 15/750 2026-04-08 19:16 by 我减肥1
[考研] 304求调剂 +16 c297914 2026-04-05 17/850 2026-04-08 13:00 by grayjzr
[考研] 312求调剂 +4 Say Never 2026-04-04 4/200 2026-04-08 08:41 by barlinike
[考研] 一志愿北交大材料工程总分358求调剂 +10 cs0106 2026-04-05 12/600 2026-04-06 19:41 by 无际的草原
[考研] (调剂)一志愿报考哈尔滨工业大学0857资源与环境专业378分考生 +7 狠狠加油 2026-04-05 8/400 2026-04-06 16:52 by momo皓
[考研] 0817化学工程与技术求调剂,一志愿中海洋319 +14 lv945 2026-04-04 14/700 2026-04-06 10:20 by 蓝云思雨
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
信息提示
请填处理意见