| 10 | 1/1 | 返回列表 |
| 查看: 2476 | 回復(fù): 9 | ||
一直在糾結(jié)鐵桿木蟲 (正式寫手)
|
[求助]
用MATLAB回歸方程的參數(shù),是個(gè)隱函數(shù),請(qǐng)附上程序 已有3人參與
|
|
數(shù)據(jù)和方程如下, x= 53.428 54.076 55.046 57.571 78.753 79.154 79.485 80.022 80.488 80.938 81.162 81.387 81.613 82.068 82.296 82.524 83.209 83.436 84.116 84.344 84.809 85.048 85.544 85.989 86.577 87.222 87.987 88.951 89.272 89.832 91.809 101.409 138.539 y= 1.29E+01 1.30E+01 1.32E+01 1.37E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.91E+01 1.91E+01 1.91E+01 1.91E+01 1.91E+01 1.89E+01 1.85E+01 1.84E+01 1.83E+01 1.83E+01 1.83E+01 1.83E+01 1.84E+01 1.88E+01 2.14E+01 3.36E+01 方程如下: logy=a1+a2/(a3+x)+a4*x+a5*x^2+a6*logx+a7*y/x^2 其中a1~a7是方程要回歸的7個(gè)參數(shù) 能幫忙回歸出來(lái)就可以,不限于用MATLAB 但是請(qǐng)把解決的程序附上,多謝!! |
金蟲 (正式寫手)
鐵桿木蟲 (正式寫手)
金蟲 (小有名氣)

木蟲 (正式寫手)
| 參考 http://www.gaoyang168.com/bbs/viewthread.php?tid=3240154 一樓有完整代碼,只要修改一下函數(shù)和數(shù)據(jù)就可以 |
鐵桿木蟲 (正式寫手)
|
呃 仿照他的輸入了一下有個(gè)錯(cuò)誤 但是不知道是哪里不行 能幫忙看一下嗎 p=[1.29E+01 1.30E+01 1.32E+01 1.37E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.91E+01 1.91E+01 1.91E+01 1.91E+01 1.91E+01 1.89E+01 1.85E+01 1.84E+01 1.83E+01 1.83E+01 1.83E+01 1.83E+01 1.84E+01 1.88E+01 2.14E+01 3.36E+01]; c=[53.428 54.076 55.046 57.571 78.753 79.154 79.485 80.022 80.488 80.938 81.162 81.387 81.613 82.068 82.296 82.524 83.209 83.436 84.116 84.344 84.809 85.048 85.544 85.989 86.577 87.222 87.987 88.951 89.272 89.832 91.809 101.409 138.539]; cc=c'; pp=p'; F=@(q,x)(q(1)+q(2)/(q(3)+x(:,1))+q(4)*x(:,1)+q(5)*(x(:,1)^2)+q(6)*log(x(:,1))+q(7)*x(:,2)/x(:,1)^2)-log(x(:,2)); x=[cc pp]; q0=[1 1 1 1 1 1 1]; warning off q=nlinfit(x,zeros(size(x,1),1),F,q0); disp(num2str(q)); 顯示錯(cuò)誤,應(yīng)該是方程輸入錯(cuò)誤,但是不知道錯(cuò)誤出在了哪里 |
|
c=[53.4280000000000,54.0760000000000,55.0460000000000,57.5710000000000,78.7530000000000,79.1540000000000,79.4850000000000,80.0220000000000,80.4880000000000,80.9380000000000,81.1620000000000,81.3870000000000,81.6130000000000,82.0680000000000,82.2960000000000,82.5240000000000,83.2090000000000,83.4360000000000,84.1160000000000,84.3440000000000,84.8090000000000,85.0480000000000,85.5440000000000,85.9890000000000,86.5770000000000,87.2220000000000,87.9870000000000,88.9510000000000,89.2720000000000,89.8320000000000,91.8090000000000,101.409000000000,138.539000000000;]; p=[12.9000000000000,13,13.2000000000000,13.7000000000000,19.3000000000000,19.3000000000000,19.3000000000000,19.3000000000000,19.3000000000000,19.3000000000000,19.3000000000000,19.2000000000000,19.2000000000000,19.2000000000000,19.2000000000000,19.2000000000000,19.2000000000000,19.1000000000000,19.1000000000000,19.1000000000000,19.1000000000000,19.1000000000000,18.9000000000000,18.5000000000000,18.4000000000000,18.3000000000000,18.3000000000000,18.3000000000000,18.3000000000000,18.4000000000000,18.8000000000000,21.4000000000000,33.6000000000000;]; cc=c'; pp=p'; F=@(q,x)(q(1)+q(2)./(q(3)+x(:,1))+q(4)*x(:,1)+q(5)*x(:,1).^2+q(6)*log(x(:,1))+q(7)*x(:,2)./x(:,1).^2-x(:,2)); x=[cc pp]; q0=[1,1,1,1,1,1,1]; warning off q=nlinfit(x,zeros(size(x,1),1),F,q0); disp('H A B分別為') disp(num2str(q)); plot(x(:,1),x(:,2),'r^');hold on; ezplot(@(x,y)F(q,[x,y])) |
金蟲 (著名寫手)

金蟲 (著名寫手)

金蟲 (著名寫手)
|
p=[1.29E+01 1.30E+01 1.32E+01 1.37E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.93E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.92E+01 1.91E+01 1.91E+01 1.91E+01 1.91E+01 1.91E+01 1.89E+01 1.85E+01 1.84E+01 1.83E+01 1.83E+01 1.83E+01 1.83E+01 1.84E+01 1.88E+01 2.14E+01 3.36E+01]; c=[53.428 54.076 55.046 57.571 78.753 79.154 79.485 80.022 80.488 80.938 81.162 81.387 81.613 82.068 82.296 82.524 83.209 83.436 84.116 84.344 84.809 85.048 85.544 85.989 86.577 87.222 87.987 88.951 89.272 89.832 91.809 101.409 138.539]; cc=c'; pp=p'; x=[cc pp]; F=@(q,x)(q(1)+q(2)./(q(3)+x(:,1))+q(4).*x(:,1)+q(5).*(x(:,1).^2)+q(6).*log(x(:,1))+q(7).*x(:,2)./x(:,1).^2)-log(x(:,2)); q0=[1 1 1 1 1 1 1]; q=nlinfit(x,zeros(size(x,1),1),F,q0); disp(num2str(q)); |

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-30 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +25 | 阿斯蒂芬2004 2026-03-28 | 25/1250 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 12/600 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +17 | Li李魚 2026-03-26 | 17/850 |
|
|
[考研] 286求調(diào)劑 +6 | Faune 2026-03-30 | 6/300 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +8 | 哇呼哼呼哼 2026-03-27 | 10/500 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-30 | 6/300 |
|
|
[論文投稿]
chinese chemical letters英文版投稿求助
130+3
|
Yishengeryi 2026-03-30 | 3/150 |
|
|
[考研] 070300化學(xué)279求調(diào)劑 +7 | 哈哈哈^_^ 2026-03-31 | 9/450 |
|
|
[考研] 合肥區(qū)域性重點(diǎn)一本招收調(diào)劑 +4 | 6266jl 2026-03-30 | 4/200 |
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
[考研] 310求調(diào)劑 +10 | 爭(zhēng)取九點(diǎn)睡 2026-03-30 | 10/500 |
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|