24小时热门版块排行榜

>论坛更新日志 (1599)
>虫友互识 (41)
>导师招生 (27)
>论文投稿 (26)
>考研 (22)
>基金申请 (21)
>硕博家园 (19)
>公派出国 (19)
>找工作 (18)
>教师之家 (17)
>考博 (17)
>休闲灌水 (16)
>博后之家 (15)
>有机交流 (14)
>论文道贺祈福 (13)
>文献求助 (9)

返回列表

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

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有270人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有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.')

赞一下