| 10 | 1/1 | 返回列表 |
| 查看: 2480 | 回復(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ù) 能幫忙回歸出來就可以,不限于用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ā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)專業(yè)279調(diào)劑 +10 | 哈哈哈^_^ 2026-03-31 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +4 | 一樣YWY 2026-03-31 | 4/200 |
|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +15 | 材料學(xué)257求調(diào)劑 2026-03-28 | 16/800 |
|
|
[考研] 調(diào)劑申請(qǐng) +7 | 張張張張zy 2026-03-31 | 7/350 |
|
|
[考研] 267求調(diào)劑 +9 | uiybh 2026-03-31 | 9/450 |
|
|
[考研] 考研材料工程351分調(diào)劑 +4 | 整個(gè)好的 2026-03-31 | 4/200 |
|
|
[考研] 材料科學(xué)與工程339求調(diào)劑 +3 | hyz0119 2026-03-31 | 3/150 |
|
|
[考研] 329求調(diào)劑,一志愿西北工業(yè)大學(xué),材料工程(085601) +6 | 小小機(jī)靈蟲 2026-03-29 | 12/600 |
|
|
[考研] 311求調(diào)劑 +9 | 藍(lán)月亮亮 2026-03-30 | 9/450 |
|
|
[考博] 材料專業(yè)申博 +5 | 杜雨婷dyt 2026-03-29 | 5/250 |
|
|
[考研] 276求調(diào)劑 +3 | 趙久華 2026-03-29 | 3/150 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|