 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2742 | 回復: 20 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
Rachellhom銅蟲 (小有名氣)
|
[求助]
用matlab實現(xiàn)一個偏最小二乘法運算
|
|
|
各位大俠，我想用matlab寫一個最小二乘法的運算代碼。 需要實現(xiàn)以下運算： 原理是： 設配制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之間可用下式關聯(lián): U= BT U, T中的列, P, Q中的行都分別相互正交, 維數a 在統(tǒng)計意義上要使T、P盡可能接近X. 迭代過程中導入一權重矢量W, 迭代循環(huán)可按下列公式分步進行: 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新+小于指定的收斂標準時, 迭代收斂. 退出迭代循環(huán)后進行下列運算: Paxp（轉置）= t（轉置）X / t（轉置）t Enxp= Xnxp—t P（轉置） b=U（轉置）t / t（轉置）t Fnxq=Ynxq — btq（轉置）. 括號里面的轉置出現(xiàn)在哪個字母后面就是代表那個矩陣的轉置, 大寫字母表示矩陣, 小寫字母表示列矢量. 將計算所得的各潛變量矩陣第一列貯存. 用新算出的E與F矩陣當作修正的X與Y矩陣, 重復上述手續(xù)計算潛變量矩陣的第二列, 直列所選定的主因子 數為止。 |
回復此樓» 猜你喜歡
環(huán)境工程 085701，267求調劑
已經有9人回復
-
085601一志愿西北工業(yè)大學初試346
已經有4人回復
-
材料與化工（0856）304求B區(qū)調劑
已經有10人回復
-
本2一志愿C9-333分，材料科學與工程，求調劑
已經有4人回復
-
一志愿中海洋320化學工程與技術學碩求調劑
已經有6人回復
-
070300化學279求調劑
已經有7人回復
-
085601一志愿中山大學深圳材料工程330求調劑
已經有5人回復
-
293分求調劑，外語為俄語
已經有4人回復
-
0856調劑
已經有3人回復
-
一志愿大連理工大學材料求調劑
已經有6人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 遞推最小二乘法在什么書本中有詳細的講解和matlab仿真介紹？
已經有5人回復
- matlab用多元線性最小二乘法求解多元線性方程的問題
已經有11人回復
Rachellhom
銅蟲 (小有名氣)
- 應助: 6 (幼兒園)
- 金幣: 1095.4
- 散金: 65
- 紅花: 4
- 帖子: 221
- 在線: 91.1小時
- 蟲號: 2054481
- 注冊: 2012-10-11
- 性別: MM
- 專業(yè): 環(huán)境工程
4樓2012-10-22 17:51:26
Rachellhom
銅蟲 (小有名氣)
- 應助: 6 (幼兒園)
- 金幣: 1095.4
- 散金: 65
- 紅花: 4
- 帖子: 221
- 在線: 91.1小時
- 蟲號: 2054481
- 注冊: 2012-10-11
- 性別: MM
- 專業(yè): 環(huán)境工程
2樓2012-10-20 13:44:32
mazewei
木蟲 (小有名氣)
- 應助: 4 (幼兒園)
- 金幣: 1754.4
- 帖子: 251
- 在線: 89.6小時
- 蟲號: 1783394
- 注冊: 2012-04-27
- 性別: GG
- 專業(yè): 機械測試理論與技術
3樓2012-10-22 15:02:05
Rachellhom
銅蟲 (小有名氣)
- 應助: 6 (幼兒園)
- 金幣: 1095.4
- 散金: 65
- 紅花: 4
- 帖子: 221
- 在線: 91.1小時
- 蟲號: 2054481
- 注冊: 2012-10-11
- 性別: MM
- 專業(yè): 環(huán)境工程
6樓2012-10-25 10:35:17
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調劑 +3 | 圖鑒212 2026-03-30 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 合肥區(qū)域性重點一本招收調劑 +4 | 6266jl 2026-03-30 | 4/200 |
|
|
[考研] 285求調劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
|
[考研] 11408總分309，一志愿東南大學求調劑，不挑專業(yè) +5 | 天賦帶到THU 2026-03-29 | 6/300 |
|
|
|
[考研] 281求調劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
|
[考研] 286分調劑 +10 | Faune 2026-03-30 | 11/550 |
|
|
|
[考研] 一志愿華中師范化學332分求調劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
|
[考研] 303求調劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 289求調劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
|
[考研] 一志愿北京理工大學本科211材料工程294求調劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研] 349求調劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
|
[考研] 0856材料化工調劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
|
[考研] 藥學105500求調劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 08開頭275求調劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 316求調劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 一志愿吉大071010，316分求調劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[碩博家園] 北京林業(yè)大學碩導招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 調劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 281求調劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
|
[考研] 303求調劑 +6 | 藍山月 2026-03-25 | 6/300 |
|