24小时热门版块排行榜    

查看: 2737  |  回复: 1

qlpeter

银虫 (初入文坛)

[求助] matlab多元非线性拟合 已有1人参与

希望虫友能给与帮助
采用多元非线性回归:
data.txt
168.7        202.9        145.9        12.1
142.0        206.4        152.9        15.8
181.9        205.8        144.2        10.7
119.7        206.8        156.2        20.9
243.9        203.9        130.6        5.0
170.1        205.8        141.9        17.2
167.3        200.1        138.7        16.3
269.2        203.2        132.0        9.7
254.1        205.4        128.9        10.6
325.2        208.8        123.4        6.7
83.9        210.0        153.4        17.1
26.7        206.4        176.6        92.9
-1.3        207.1        190.8        149.8
79.3        209.2        157.2        21.4
68.2        210.6        167.9        44.7
41.1        205.2        162.9        61.3
8.4        198.8        176.1        115.7
50.1        209.8        167.3        64.9
189.7        216.8        160.2        20.9
50.7        200.4        158.2        46.4
35.9        207.1        167.9        64.0
function main
clc
clear
data=load('data.txt');
y=data(:,1);
x=data(:,2:4);
beta=;            % 回归系数的初值,任意取
yy=fun(beta,x);
=nlinfit(x,y,@fun,beta);   
beta_out
mse
betaci=nlparci(beta_out,r,'Jacobian',J);   
betaa=        
=nlpredci(@fun,x,beta_out,r,'Jacobian',J);
alpha=0.05;
nlintool(x,y,'fun',beta,alpha)  
r=corrcoef(y,yy)       % 求相关系数r
plotregression(y,yy)   % 并且画图r
regstats(y,yy)
function yy=fun(beta,x)
b1 = beta(1);
b2 = beta(2);
b3 = beta(3);
b4 = beta(4);
yy=b1*(x(:,1).^b2).*(x(:,2).^b3).*(x(:,3).^b4);


运行时显示??? Error: File: D:\新建文件夹\main.m Line: 9 Column: 1
The expression to the left of the equals sign is not a valid target for an assignment.
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
均方差(RMSE):29.1845118412956
残差平方和(SSE):17886.4503597092
相关系数(R): 0.950674501500641
相关系数之平方(R^2): 0.903782007803491
确定系数(DC): 0.899149439045897
卡方系数(Chi-Square): 116.489341933952
F统计(F-Statistic): 56.1887471230383

参数        最佳估算
----------        -------------
b1        0.0690396368659353
b2        4.68095500291962
b3        -3.29925868731259
b4        -0.310074339479526

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

No        实测值y        计算值y
1        168.7        145.834326238406
2        142        124.609811365089
3        181.9        168.291349599024
4        119.7        107.453774499923
5        243.9        282.878567758331
6        170.1        153.172181912177
7        167.3        147.231225521178
8        269.2        218.821981736768
9        254.1        242.150731537058
10        325.2        348.111406472481
11        83.9        130.430380974856
12        26.7        44.7214275898647
13        -1.3        30.3560270325071
14        79.3        110.242995026258
15        68.2        72.8394588615087
16        41.1        64.6176535703849
17        8.4        35.3797172907226
18        50.1        64.4980548296721
19        189.7        123.308050070583
20        50.7        69.4516535737008
21        35.9        60.2511608251563
2楼2014-02-27 11:46:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qlpeter 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 00:23 by YXM1FDCSE8Bn
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 00:23 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 00:22 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-10-31 23:23 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-10-31 23:18 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-10-31 23:17 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-10-31 23:10 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-10-31 23:10 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 24/1200 2026-10-31 22:39 by YXM1FDCSE8Bn
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-10-31 21:59 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-10-31 21:54 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-10-31 21:51 by YXM1FDCSE8Bn
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-10-31 21:51 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 9/450 2026-10-31 21:34 by YXM1FDCSE8Bn
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[基金申请] 论文的挂名怎么评判贡献?基金申请中,排名第一的挂名文章 +3 biily 2026-10-25 3/150 2026-10-30 18:33 by 旁观2020
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见