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

登录注册

应《网络安全法》要求，自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
求助我自己仿照其他帖子写了个分段函数非线性拟合的程序一直出错求高手！新手刚刚接触非常感谢

回复此楼

» 猜你喜欢

求助中医院临床中医中标的青年标书学习，谢谢已经有6人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有259人回复
物理学I论文润色/翻译怎么收费? 已经有282人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复
25届硕士申请，材料计算方向已经有2人回复
课题组招1-2名硕士已经有1人回复
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上已经有0人回复
斯坦福大学与Elsevier世界前2%学者新排名查询已经有1人回复
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有19人回复
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计已经有1人回复
apl投稿状态已经有2人回复

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

关于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.')

赞一下