应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab非线性拟合模型参数，大神你在哪！
742/8返回列表上一页123456下一页末页
查看: 25189  |  回复: 73
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

stwangkaiyan
好

[ 发自手机版 http://muchong.com/3g ]
11楼2014-01-15 04:29:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2012200838

新虫 (正式写手)
引用回帖:
2楼: Originally posted by jerkwin at 2013-11-05 11:23:25
不劳matlab, origin, 1stopt足矣

2013-11-04_212307.png

大神，请问Matlab有没有直接求决定系数的命令。
赞一下 回复此楼
生命就是一场旅行，不要错过了路边的风景！
12楼2014-10-18 16:01:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)
引用回帖:
7楼: Originally posted by 月只蓝 at 2013-11-05 14:09:49
（1） MATLAB代码及结果

function fit_nonl
clear all;clc
format long

data=;
xspan=data(:,2);    %x的数据，在此输入
Texp=data(:,1);    %T的数据，在此输入

k0=;
lb=-*1e9;
ub=*1e9;


%--- ...

能不能请您也帮我用最小二乘法拟合下我的三元二次多项式回归方程？公式请见附件照片当然授人以鱼不如授人以渔，如果能知道编程命令，更是求之不得。可是作为MATLAB盲，命令可能一时半会儿也学不会，数据先发您好吗？谢谢！
matlab非线性拟合模型参数，大神你在哪！
QQ图片20150120163713.png


matlab非线性拟合模型参数，大神你在哪！-1
QQ截图20150120161540.png
赞一下 回复此楼
13楼2015-01-20 16:40:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
13楼: Originally posted by 梦落在草地上 at 2015-01-20 16:40:18
能不能请您也帮我用最小二乘法拟合下我的三元二次多项式回归方程？公式请见附件照片当然授人以鱼不如授人以渔，如果能知道编程命令，更是求之不得。可是作为MATLAB盲，命令可能一时半会儿也学不会，数据先发您好吗 ...

请把公式和数据以文本的形式给出，数据点很多，一个一个地敲比较麻烦。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
14楼2015-01-20 17:10:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)
引用回帖:
14楼: Originally posted by 月只蓝 at 2015-01-20 17:10:07
请把公式和数据以文本的形式给出，数据点很多，一个一个地敲比较麻烦。...

好的：
自变量X1=[45 45 40 40 35 35 30 30 25 25 20 20 15 15]
         X2=[4 5 6 7 8 9 10 4 5 6 7 8 9 10]
         X3=[125 75 25 150 100 50 0 150 100 50 0 125 75 50]
因变量Y=[20.8 3.8 41.2 4.51 6.33 45.9 36.2 53.3 62.6 8.74 34.2 63 32.5 30.8]
公式Y=bo+b1*X1+b2*X2+b3*X3+b4*X1*X2+b5*X1*X3+b6*X2*X3+b7*X1^2+X2^2+X3^2
请您教我命令代码 同时给出回归方程结果成吗？多谢！
赞一下 回复此楼
15楼2015-01-20 18:12:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
15楼: Originally posted by 梦落在草地上 at 2015-01-20 18:12:06
好的：
自变量X1=
         X2=
         X3=
因变量Y=
公式Y=bo+b1*X1+b2*X2+b3*X3+b4*X1*X2+b5*X1*X3+b6*X2*X3+b7*X1^2+X2^2+X3^2
请您教我命令代码 同时给出回归方程结果成吗？多谢！...

用1stopt软件做比较方便，该软件参见：http://muchong.com/bbs/viewthread.php?tid=4586018&fpage=1&target=blank

求解该问题的代码：
CODE:
Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9;
Variable X1,X2,X3,Y;
Function Y=b0+b1*X1+b2*X2+b3*X3+b4*X1*X2+b5*X1*X3+b6*X2*X3+b7*X1^2+b8*X2^2+b9*X3^2;
Data;
   45.0000    4.0000  125.0000   20.8000
   45.0000    5.0000   75.0000    3.8000
   40.0000    6.0000   25.0000   41.2000
   40.0000    7.0000  150.0000    4.5100
   35.0000    8.0000  100.0000    6.3300
   35.0000    9.0000   50.0000   45.9000
   30.0000   10.0000         0   36.2000
   30.0000    4.0000  150.0000   53.3000
   25.0000    5.0000  100.0000   62.6000
   25.0000    6.0000   50.0000    8.7400
   20.0000    7.0000         0   34.2000
   20.0000    8.0000  125.0000   63.0000
   15.0000    9.0000   75.0000   32.5000
   15.0000   10.0000   50.0000   30.8000

计算结果：
均方差(RMSE): 9.88659706466202
残差平方和(SSE): 1368.42722126577
相关系数(R): 0.867447328275689
相关系数之平方(R^2): 0.75246486733263
决定系数(DC): 0.752464867332631
卡方系数(Chi-Square): 41.2555999279191
F统计(F-Statistic): 36.4779670291298

参数        最佳估算
----------        -------------
b0                 1617.90831405883
b1                 -29.8254524959839
b2                 -302.054797806611
b3                 -3.20472492239783
b4                 2.09026520065958
b5                 -0.0277912442885419
b6                 0.416964792688625
b7                 0.273168393633342
b8                 14.6490169337092
b9                 0.00764543844729996

====== 结果输出 =====

No        实测值y        计算值y
1        20.8        2.3677723
2        3.8        20.4076232
3        41.2        38.0728655
4        4.51        13.0313386
5        6.33        7.3060878
6        45.9        33.0363177
7        36.2        40.4295690
8        53.3        62.4251003
9        62.6        55.2230099
10        8.74        20.8055830
11        34.2        26.7219948
12        63        52.5637974
13        32.5        35.0897160
14        30.8        36.3992237
matlab非线性拟合模型参数，大神你在哪！-2
附图1.png
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
16楼2015-01-21 09:01:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)
引用回帖:
16楼: Originally posted by 月只蓝 at 2015-01-21 09:01:30
用1stopt软件做比较方便，该软件参见：http://muchong.com/bbs/viewthread.php?tid=4586018&fpage=1&target=blank

求解该问题的代码：
Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9;
Variable X1,X2,X3,Y ...

您好，谢谢您的解答。同时请问matlab程序代码是怎样的啊？很复杂还是？
赞一下 回复此楼
17楼2015-01-21 09:13:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)
引用回帖:
16楼: Originally posted by 月只蓝 at 2015-01-21 09:01:30
用1stopt软件做比较方便，该软件参见：http://muchong.com/bbs/viewthread.php?tid=4586018&fpage=1&target=blank

求解该问题的代码：
Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9;
Variable X1,X2,X3,Y ...

你好，当我复制粘贴上您的代码后是这样的，然后什么步骤能得到计算结果、参数等等呢？
matlab非线性拟合模型参数，大神你在哪！-3
1.jpg
赞一下(1人) 回复此楼
18楼2015-01-21 09:54:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
17楼: Originally posted by 梦落在草地上 at 2015-01-21 09:13:21
您好，谢谢您的解答。同时请问matlab程序代码是怎样的啊？很复杂还是？...

MATLAB代码：
CODE:
function duoyuan_huigui
clear all;clc
format long

X1=[45 45 40 40 35 35 30 30 25 25 20 20 15 15]';
X2=[4 5 6 7 8 9 10 4 5 6 7 8 9 10]';
X3=[125 75 25 150 100 50 0 150 100 50 0 125 75 50]';
Y=[20.8 3.8 41.2 4.51 6.33 45.9 36.2 53.3 62.6 8.74 34.2 63 32.5 30.8]';

X=[ones(size(Y)) X1 X2 X3 X1.*X2 X1.*X3  X2.*X3  X1.^2 X2.^2 X3.^2];
[b,bint,r,rint,stats]=regress(Y,X)


figure(33)
subplot(1,2,1)
plot(Y,X*b,'ob',[min(Y)-5 max(Y)+5],[min(Y)-5 max(Y)+5],'r-'),axis([min(Y)-5 max(Y)+10 min(Y)-5 max(Y)+5]);

subplot(1,2,2)
plot([1:length(Y)],Y,'b-',[1:length(Y)],X*b,'r-'),axis([0 length(Y)+1 min(Y)-5 max(Y)+5]),legend('实验值','计算值','Location','Best')

计算结果：
b =

  1.0e+003 *

   1.617908294627979
  -0.029825452255521
  -0.302054793675629
  -0.003204724855954
   0.002090265179822
  -0.000027791244422
   0.000416964783243
   0.000273168391924
   0.014649016730840
   0.000007645438431

可见与1stopt软件的计算结果一致。
matlab非线性拟合模型参数，大神你在哪！-4
附图1.png
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
19楼2015-01-21 10:05:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
18楼: Originally posted by 梦落在草地上 at 2015-01-21 09:54:28
你好，当我复制粘贴上您的代码后是这样的，然后什么步骤能得到计算结果、参数等等呢？

1.jpg
...

管理员身份运行 1stopt 软件。
删除代码中所有的 “//”，其实你的代码如果完全复制我16楼给出的，也不会出错。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
20楼2015-01-21 10:07:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dshuanger 的主题更新
742/8返回列表上一页123456下一页末页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 3/150 2026-11-09 14:58 by 天下第一剑
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭