 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1675 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
chaiqing123金蟲 (小有名氣)
|
[求助]
MATLAB用于三元回歸參數(shù)
|
|
|
function文件為下 function rcal = fun552(x,x1,x2) for n=1:11%m為所取的數(shù)據(jù)點個數(shù) R=8.315; global x1 x2 x3 T; tao12=x(1)./R/T(n);%x(1)=g12-g22 tao21=x(2)./R/T(n); %x(2)=g21-g11 tao13=x(3)./R/T(n); %x(3)=g13-g33 tao31=x(4)./R/T(n); %x(4)=g31-g11 tao23=x(5)./R/T(n); %x(5)=g23-g33 tao32=x(6)./R/T(n); %x(6)=g32-g22 tao=[tao12,tao13,tao21,tao23,tao31,tao32]; G=exp(-0.3*tao); A(n)=(tao(3)*G(3)*x2(n)+tao(5)*G(5)*x3(n))/(x1(n)+G(3)*x2(n)+G(5)*x3(n)); B(n)=-x1(n)*(x2(n)*tao(3)*G(3)+x3(n)*tao(5)*G(5))/((x1(n)+G(3)*x2(n)+G(5)*x3(n))^2); C(n)=x2(n)*G(1)*(tao(1)-(x1(n)*tao(1)*G(1)+x3(n)*tao(6)*G(6))/(G(1)*x1(n)+x2(n)+G(6)*x3(n)))/(G(1)*x1(n)+x2(n)+G(6)*x3(n)); D(n)=x3(n)*G(2)*(tao(2)-(x1(n)*tao(2)*G(2)+x2(n)*tao(4)*G(4))/(G(2)*x1(n)+G(4)*x2(n)+x3(n)))/(G(2)*x3(n)+G(4)*x2(n)+x3(n)); r1cal(n)=exp(A(n)+B(n)+C(n)+D(n)); A1(n)=(tao(1)*G(1)*x1(n)+tao(6)*G(6)*x3(n))/(G(1)*x1(n)+G(6)*x3(n)+x2(n)); B1(n)=-x2(n)*(x1(n)*tao(1)*G(1)+x3(n)*tao(6)*G(6))/((G(1)*x1(n)+x2(n)+G(6)*x3(n))^2); C1(n)=x1(n)*G(3)/(x1(n)+G(3)*x2(n)+G(5)*x3(n))*(tao(3)-(x2(n)*tao(3)*G(3)+x3(n)*tao(5)*G(5))/(x1(n)+G(3)*x2(n)+G(5)*x3(n))); D1(n)=x3(n)*G(4)/(G(2)*x1+G(4)*x2(n)+x3(n))*(tao(4)-(x1(n)*tao(2)*G(2)+x2(n)*tao(4)*G(4))/(G(2)*x1(n)+G(4)*x2(n)+x3(n))); r2cal(n)=exp(A1(n)+B1(n)+C1(n)+D1(n)); rcal=[r1cal r2cal]; end script文件如下 xx1=[0.05 0.105 0.153 0.248 0.347 0.446 0.474 0.535 0.624 0.740 0.870];%不含離子液體的x1 xx2=1-xx1; x1=xx1/1.05;%離子液體濃度為5% x2=xx2/1.05; x3=1-x1-x2; global x1 x2 x3 T; y1=[0.397 0.483 0.521 0.601 0.635 0.676 0.685 0.713 0.745 0.825 0.897]; y2=1-y1; T=[364 360.2 358.5 356.3 355.3 354.5 354.3 354.1 353.9 353.7 353.6]; t=T-273.15; a=[7.23255 1750.286 235;7.16879 1552.601 222.419]; p1s=10.^(a(1,1)-a(1,2)./(t+a(1,3))); p2s=10.^(a(2,1)-a(2,2)./(t+a(2,3))); %A=7.23255; B=1750.286;C=235; %A1=7.16879;B1=1552.601;C1=222.419; %p1s=10.^(A-B./(t+C)) %p2s=10.^(A1-B1./(t+C1)) p3s=0; %a為aotonie方程的三個系數(shù) P=101.3%總壓 r1exp=y1.*P./(x1.*p1s); r2exp=y2.*P./(x2.*p2s); rexp=[r1exp r2exp] x0=[1,1,1,1,1,1]; xx=[xx1 xx1]; x=nlinfit(xx,rexp,'fun552',x0) 最小二乘法 這個有問題 但不知道怎么改.... |
回復(fù)此樓» 猜你喜歡
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有133人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab反求未知輸入?yún)?shù)
已經(jīng)有4人回復(fù)
- 求回歸方程參數(shù)估計的最優(yōu)化算法matlab代碼
已經(jīng)有9人回復(fù)
- matlab非線性參數(shù)擬合問題
已經(jīng)有7人回復(fù)
- 如何在matlab中繪制含有三個參數(shù)方程的三維圖形
已經(jīng)有6人回復(fù)
- 【求助】matlab 求解微分方程中的未知參數(shù)
已經(jīng)有20人回復(fù)
- 有人學(xué)過<神經(jīng)網(wǎng)絡(luò)>和《模糊數(shù)學(xué)》還有《matlab》這三門課嗎？名字都好虛幻~
已經(jīng)有20人回復(fù)
- 【求助】三元二次方程的回歸擬合
已經(jīng)有5人回復(fù)
- 三元方程用matlab畫圖，求指令
已經(jīng)有1人回復(fù)
chaiqing123
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 679.7
- 紅花: 1
- 帖子: 71
- 在線: 65.2小時
- 蟲號: 1704607
- 注冊: 2012-03-20
- 性別: GG
- 專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)
4樓2012-04-01 14:48:51
695710829
新蟲 (初入文壇)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 35.1
- 帖子: 5
- 在線: 9.2小時
- 蟲號: 1448756
- 注冊: 2011-10-18
- 專業(yè): 分離過程
2樓2012-04-01 10:52:29
chaiqing123
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 679.7
- 紅花: 1
- 帖子: 71
- 在線: 65.2小時
- 蟲號: 1704607
- 注冊: 2012-03-20
- 性別: GG
- 專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)
3樓2012-04-01 14:46:38
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +6 | Y-cap 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 材料與化工272求調(diào)劑 +12 | 阿斯蒂芬2004 2026-03-28 | 12/600 |
|
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +8 | 超級伊昂大王 2026-03-24 | 9/450 |
|
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 上海電力大學(xué)材料防護與新材料重點實驗室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
|
[有機交流]
 有機合成求助
 20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|