版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

hanmei98232

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 2315.9
红花: 2
帖子: 285
在线: 95.7小时
虫号: 1357946
注册: 2011-07-31
性别: MM
专业: 催化化学

[求助] 求助我自己仿照其他帖子写了个分段函数非线性拟合的程序一直出错求高手！已有1人参与

clear
tic
fun1=inline('-((beta(1)*x)/(beta(2)*x+1)).*8.314*298.15*(log(x)-0.784*(1.7321*x.^0.5/(1+2.0785*x.^0.5)+1.6667*log(1+2.0785*x.^0.5))+2.6667*x.*(beta(3)+beta(4)*(1+(-1-3.464*x.^0.5+6*x).*exp(-3.464*x.^0.5))/(12*x))+2.8289*x.^2*beta(5))+0.07214','beta','x');
fun2=inline('((beta(1)*x)/(beta(2)*x+1)).*8.314*298.15*(log(x)-0.784*(1.7321*x.^0.5/(1+2.0785*x.^0.5)+1.6667*log(1+2.0785*x.^0.5))+2.6667*x.*(beta(3)+beta(4)*(1+(-1-3.464*x.^0.5+6*x).*exp(-3.464*x.^0.5))/(12*x))+2.8289*x.^2*beta(5))+0.07214','beta','x');
fun=@(beta,x)((x<1).*fun1(beta,x)+(x>=1).*fun2(beta,x));
x=[0.19727 0.39798 0.60182 0.8034 0.98454 1.78733 2.05406 2.51801 3.16362 3.69954 4.5584 ];
y=[0.07282919 0.07366338 0.07454467 0.07522262 0.07593838 0.07919294 0.08033517 0.08203459 0.08494242 0.0867815 0.09083547];
beta=nlinfit(x,y,fun,[0.1,0.2,0.1,0.1,0]);
m=beta(1);
n=beta(2);
a=beta(3);
b=beta(4);
c=beta(5);
warning off all
y1=-((beta(1)*x)/(beta(2)*x+1)).*8.314*298.15*(log(x)-0.784*(1.7321*x.^0.5/(1+2.0785*x.^0.5)+1.6667*log(1+2.0785*x.^0.5))+2.6667*x.*(beta(3)+beta(4)*(1+(-1-3.464*x.^0.5+6*x).*exp(-3.464*x.^0.5))/(12*x))+2.8289*x.^2*beta(5))+0.07214;
y2=((beta(1)*x)/(beta(2)*x+1)).*8.314*298.15*(log(x)-0.784*(1.7321*x.^0.5/(1+2.0785*x.^0.5)+1.6667*log(1+2.0785*x.^0.5))+2.6667*x.*(beta(3)+beta(4)*(1+(-1-3.464*x.^0.5+6*x).*exp(-3.464*x.^0.5))/(12*x))+2.8289*x.^2*beta(5))+0.07214;
yfit(1:5)=y1;
yfit(6:11)=y2;
plot(x,y,'*')
hold on;
plot(x,yfit,'r-')
toc
求助我自己仿照其他帖子写了个分段函数非线性拟合的程序一直出错求高手！新手刚刚接触非常感谢

回复此楼

» 猜你喜欢

如何从铁电相到顺电相。已经有1人回复
各位大佬，求一份最新Wien2k版本作为学习用途已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有86人回复
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后已经有42人回复
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有24人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有3人回复
PVA溶解已经有2人回复
基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复

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

关于matlab两个自变量、一个因变量曲线拟合问题已经有3人回复
求助用matlab 多元非线性拟合2 已经有21人回复
求助用matlab 多元非线性拟合已经有4人回复
求助该图形用什么函数曲线拟合，已经有3人回复
如何用已知函数去拟合实验数据，得出参数已经有10人回复
matlab 拟合求助已经有8人回复
关于matlab拟合已经有4人回复
如何用matlab实现分峰拟合？已经有18人回复
求助Matlab最小二乘法非线性拟合已经有4人回复
Matlab非线性拟合（y=sqrt(1-(x-a).^2)）已经有5人回复
求助origin非线性拟合的编译（origin C）已经有4人回复
毕不了业了，发帖求助Origin8.0多元非线性拟合的问题已经有9人回复
用matlab里的nlinfit如何拟合分段函数已经有8人回复
二阶非线性微分方程求解已经有14人回复
matlab中关于非线性拟合，用nlinfit.程序总出错。谢谢各位大侠已经有3人回复
请教大虾，一个分段函数的 matlab 程序实现已经有4人回复
[求助]origin中非线性拟合中Guass和GuassAmp两个函数的区别？已经有4人回复
求助Origin中自定义求和函数的非线性拟合已经有0人回复
【求助】急求关于matlab做非线性拟合的函数已经有8人回复
【求助】这个带有分段函数偏微分方程怎么用for程序去解决已经有14人回复
【求助】求助仿真中需要一个分段函数来控制电流源已经有5人回复
【求助】关于origin进行非线性拟合的自定义函数已经有7人回复
【求助】关于分段函数的表示成一个式子的方法已经有7人回复
【求助】多变量非线性拟合【已完结】已经有8人回复
【求助】紧急求助有关MATHEMATICA的问题已经有9人回复

1楼 2014-03-17 13:52:16

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

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

CODE:

fun1=inline('-((beta(1)*x)/(beta(2)*x+1)).*8.314*298.15*(log(x)-0.784*(1.7321*x.^0.5/(1+2.0785*x.^0.5)+1.6667*log(1+2.0785*x.^0.5))+2.6667*x.*(beta(3)+beta(4)*(1+(-1-3.464*x.^0.5+6*x).*exp(-3.464*x.^0.5))/(12*x))+2.8289*x.^2*beta(5))+0.07214','beta','x'); 

fun2=inline('((beta(1)*x)/(beta(2)*x+1)).*8.314*298.15*(log(x)-0.784*(1.7321*x.^0.5/(1+2.0785*x.^0.5)+1.6667*log(1+2.0785*x.^0.5))+2.6667*x.*(beta(3)+beta(4)*(1+(-1-3.464*x.^0.5+6*x).*exp(-3.464*x.^0.5))/(12*x))+2.8289*x.^2*beta(5))+0.07214','beta','x'); 

fun=@(beta,x)((x<1).*fun1(beta,x)+(x>=1).*fun2(beta,x)); 

x=[0.19727        0.39798        0.60182        0.8034        0.98454        1.78733        2.05406        2.51801        3.16362        3.69954        4.5584 ];

y=[0.07282919        0.07366338        0.07454467        0.07522262        0.07593838        0.07919294        0.08033517        0.08203459        0.08494242        0.0867815        0.09083547];

beta=nlinfit(x,y,fun,[0.1,0.2,0.1,0.1,0]);

plot(x,y,'ro',x,fun(beta,x),'k.')

赞一下