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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

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

回复此楼

» 猜你喜欢

调剂已经有0人回复
【新加坡】纳米电子器件项目组有“联合培养博士生”名额已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有295人回复
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职已经有0人回复
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职已经有0人回复
【新加坡】纳米电子器件项目组有“联合培养博士生”名额已经有2人回复
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。已经有0人回复
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！已经有1人回复
天津城建大学理学院化学、环境工程、材料物理等方向接收10人已经有0人回复
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人已经有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.')

赞一下