24小时热门版块排行榜

返回列表

qlpeter

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 490.8
帖子: 42
在线: 43.6小时
虫号: 1991644
注册: 2012-09-11
性别: GG
专业: 水环境与生态水利

[求助] 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.

回复此楼

» 猜你喜欢

实验室接单子已经有4人回复
全日制（定向）博士已经有4人回复
假如你的研究生提出不合理要求已经有6人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
论文终于录用啦！满足毕业条件了已经有28人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复

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

origin8.5 自定义方程的非线性拟合问题已经有10人回复
关于matlab两个自变量、一个因变量曲线拟合问题已经有3人回复
在进行多元非线性拟合的时候，用matlab拟合的值残差都大于0 已经有14人回复
求助用matlab 多元非线性拟合已经有4人回复
请帮忙拟合一组数据，求回归系数已经有3人回复
非线性拟合得到的模型可靠性怎么验证？已经有8人回复
线性拟合中的问题已经有9人回复
origin 非线性自定义拟合问题~~方程组编写及运行~~ 已经有8人回复
毕不了业了，发帖求助Origin8.0多元非线性拟合的问题已经有9人回复
用matlab里的nlinfit如何拟合分段函数已经有8人回复
关于matlab微分、及非线性拟合问题已经有5人回复
非线性数据拟合的数学原理！！！！！？？？？已经有30人回复
求助用matlab拟合三元线性回归方程及分析已经有9人回复
求助matlab解非线性方程。求代码已经有5人回复
ORINGIN的非线性拟合... 已经有10人回复
matlab多元非线性拟合已经有9人回复
matlab非线性参数拟合问题已经有7人回复
多元线性拟合，其中两组数据相关系数比较高，怎么办？已经有6人回复
【求助】多元方程组求解及拟合已经有9人回复
求助关于线性拟合的已经有1人回复
Matlab 的自定义函数拟合已经有3人回复

1楼 2014-02-27 10:44:30

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

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

感谢参与，应助指数 +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

赞一下