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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助：多变量曲线拟合，没有经验公式
101/1返回列表
查看: 2336  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

shidakfy

新虫 (初入文坛)

[交流] 求助：多变量曲线拟合，没有经验公式 已有5人参与

各路大神：
           小弟遇到一个多元曲线拟合的问题，请大神们多多指教！一个因变量，三个自变量，在未知其经验模型（公式）的情况下，希望通过拟合找出因变量和三个自变量之间的关系方程。本人数学小白！！各路软件实在是搞不明白，跪求大神帮忙！！
这是我的实验数据：
X       Y      Z        m
10        80        297        139
10        40        297        154
10        20        297        156.5
20        80        297        169
20        40        297        182
20        20        297        187
30        80        297        200
10        80        340        211
50        80        297        211
30        40        297        212
30        20        297        215
70        80        297        215
50        40        297        218
10        40        340        223
50        20        297        224
70        40        297        226
10        20        340        231
70        20        297        231
20        80        340        234
20        40        340        242
20        20        340        248
30        80        340        253
30        40        340        256
10        80        384        258
50        80        340        258
10        40        384        264
50        40        340        264
70        80        340        269
30        20        340        271
10        20        384        274
10        80        393        274
70        40        340        276
50        20        340        278
20        80        384        279
10        40        393        281
20        40        384        283
20        20        384        289
70        20        340        290
10        20        393        294
20        80        393        294
30        80        384        295
20        40        393        297
30        40        384        298
20        20        393        305
30        20        384        306
50        80        384        308
30        80        393        310
30        40        393        314
50        40        384        314
50        20        384        319
30        20        393        321
50        80        393        321
70        80        384        325
70        40        384        328
50        40        393        329
70        20        384        331
50        20        393        334
70        80        393        338
70        40        393        344
70        20        393        349
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-03-19 19:14:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
m=a*x+b*y+c*z+d+ e*x^2+f*y^2+g*z^2+ h*x*y+k*x*z+l*y*z

相关系数(R): 0.993488624488174
相关系数之平方(R^2): 0.987019646987404
决定系数(DC): 0.987019646987405
卡方系数(Chi-Square): 4.52272468670753
F统计(F-Statistic): 4410.291420404

参数        最佳估算
----------        -------------
a                  2.865001615092
b                  -0.934855325951024
c                  2.61500660142316
d                  -436.463275749684
e                  -0.0181087864234095
f                  0.0031874999923929
g                  -0.00205952018605737
h                  0.000825123153197607
k                  -0.00127126387431178
l                  0.000969309457521069
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2015-03-19 19:23:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖仅楼主可见
3楼2015-03-22 16:47:42
已阅   申请计算强帖   回复此楼   编辑   查看我的主页

月只蓝

主管区长 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by dyhseu at 2015-03-22 16:47:42
您好！我是一名大二的学生，用1stopt 1.5有一个问题想请教您，因为刚刚注册不能发帖和密语，所以找到您，麻烦您帮忙解答下谢谢！~代码如下：
Constant     W=25000,b=0.3,D=1,R=0.5;                             ...

MATLAB代码：
CODE:
function solve_nonleqs

clear all;clc
format long
global i

ispan=0.1:0.05:0.9;
m=length(ispan);
for j=1:m
    i=ispan(j);
x0=0.5;
options=optimset('TolFun',1e-6);
[o(j) fval(j)] =fsolve(@eqs,x0,options);


W=25000;b=0.3;D=1;R=0.5;                                                      
c=5000;phi=(30/180*pi);n=0.8;kc=10000;kphi=1000000;k=0.015;  %土壤参数

                                                            %滑转率
                                                    %接触角，挂钩牵引力

sigma=(kc/b+kphi)*(0.5*R*o(j)^2)^n*(1-n/3);                    %驱动轮正应力
tau=c+sigma*tan(phi);                                                     %驱动轮剪应力
eta=sigma/tau;
zeta=i*R/k;


FDP(j)  =b*R*tau*(eta*(cos(o(j))-1)+sin(o(j))-1/(1+zeta^2)*(zeta*cos(o(j))+sin(o(j))-zeta*exp(-zeta*o(j)))) ;    %由o求FDP
end
o=o'  
FDP=FDP',
cancha=fval'

figure(1),plot(ispan,FDP,'MarkerEdgeColor',[1 0 0],'MarkerSize',12,'Marker','o',...
    'LineWidth',3);


% Create xlabel
xlabel('i','FontSize',24);

% Create ylabel
ylabel('FDP','FontSize',24);


function f=eqs(x)
global i

o=x;

W=25000;b=0.3;D=1;R=0.5;                                                      
c=5000;phi=(30/180*pi);n=0.8;kc=10000;kphi=1000000;k=0.015;  %土壤参数

                                                            %滑转率
                                                    %接触角，挂钩牵引力

sigma=(kc/b+kphi)*(0.5*R*o^2)^n*(1-n/3);                    %驱动轮正应力
tau=c+sigma*tan(phi);                                                     %驱动轮剪应力
eta=sigma/tau;
zeta=i*R/k;



f= 0.4*W- ( b*R*tau*((1-cos(o))+eta*sin(o)-1/(1+zeta^2)*(zeta*sin(o)-cos(o)+   exp(-zeta*o))) );       %已知W求o
            

计算结果

                   i                               o
   0.100000000000000   0.596897661759681
   0.150000000000000   0.592416842245195
   0.200000000000000   0.589392464072042
   0.250000000000000   0.587242294598273
   0.300000000000000   0.585649298230781
   0.350000000000000   0.584428784027905
   0.400000000000000   0.583467346996124
   0.450000000000000   0.582692260032518
   0.500000000000000   0.582055119345026
   0.550000000000000   0.581522653958748
   0.600000000000000   0.581071339161757
   0.650000000000000   0.580684122302388
   0.700000000000000   0.580348367320835
   0.750000000000000   0.580054526503918
   0.800000000000000   0.579795259470657
   0.850000000000000   0.579564834823698
   0.900000000000000   0.579358714965294

FDP =

  1.0e+003 *

   0.492340219836709
   1.153667339390909
   1.547889709981589
   1.799997202585284
   1.971130969270374
   2.093284993995862
   2.184196300087984
   2.254223160185328
   2.309708685008577
   2.354707936232967
   2.391915038887142
   2.423181977524017
   2.449820479830493
   2.472784217144151
   2.492782549378808
   2.510353763231505
   2.525913569602112


残差 =

  1.0e-011 *

   0.181898940354586
   0.181898940354586
   0.545696821063757
  -0.545696821063757
  -0.181898940354586
   0.181898940354586
                   0
                   0
   0.181898940354586
  -0.181898940354586
  -0.181898940354586
  -0.363797880709171
   0.545696821063757
   0.181898940354586
  -0.363797880709171
  -0.545696821063757
   0.181898940354586
求助：多变量曲线拟合，没有经验公式
附图1.png
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2015-03-22 19:04:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖仅楼主可见
赞一下
5楼2015-03-22 21:41:48
已阅   申请计算强帖   回复此楼   编辑   查看我的主页

匿名

用户注销 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖仅楼主可见
赞一下
6楼2015-03-22 21:46:13
已阅   申请计算强帖   回复此楼   编辑   查看我的主页

curly2017

禁虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽
7楼2017-04-12 10:22:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
内容已删除
回复此楼
8楼2017-04-12 10:46:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

curly2017

禁虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽
9楼2017-04-12 14:33:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
1.5版的估计都还不支持“LoopConstant”命令吧。
赞一下 回复此楼
10楼2017-04-12 16:26:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shidakfy 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯，如何抉择 +6 Isak0320 2026-10-12 9/450 2026-10-12 21:13 by Isak0320
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[教师之家] 课题组会帮你通过考核么？ +3 好运来～9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭