版塊導航: 正在加載中...

調(diào)劑小程序

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (3828)
>導師招生 (425)
>考研 (391)
>文獻求助 (126)
>蟲友互識 (93)
>考博 (70)
>博后之家 (66)
>論文投稿 (64)
>基金申請 (55)
>教師之家 (48)
>碩博家園 (39)
>招聘信息布告欄 (37)
>公派出國 (25)
>第一性原理 (16)
>找工作 (16)
>休閑灌水 (15)

返回列表

hanmei98232

木蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 2315.9
紅花: 2
帖子: 285
在線: 95.7小時
蟲號: 1357946
注冊: 2011-07-31
性別: MM
專業(yè): 催化化學

[求助] 求助我自己仿照其他帖子寫了個分段函數(shù)非線性擬合的程序一直出錯求高手！已有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
求助我自己仿照其他帖子寫了個分段函數(shù)非線性擬合的程序一直出錯求高手！新手剛剛接觸非常感謝

回復此樓

» 猜你喜歡

光學工程學碩調(diào)劑信息已經(jīng)有32人回復
歡迎加入課題組已經(jīng)有0人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有196人回復
散金幣，求好運，祝面上順利！已經(jīng)有78人回復
華東師范大學芯片設計徐瓏真實水平咋樣已經(jīng)有1人回復
調(diào)劑已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

關于matlab兩個自變量、一個因變量曲線擬合問題已經(jīng)有3人回復
求助用matlab 多元非線性擬合2 已經(jīng)有21人回復
求助用matlab 多元非線性擬合已經(jīng)有4人回復
求助該圖形用什么函數(shù)曲線擬合，已經(jīng)有3人回復
如何用已知函數(shù)去擬合實驗數(shù)據(jù)，得出參數(shù) 已經(jīng)有10人回復
matlab 擬合求助已經(jīng)有8人回復
關于matlab擬合已經(jīng)有4人回復
如何用matlab實現(xiàn)分峰擬合？已經(jīng)有18人回復
求助Matlab最小二乘法非線性擬合已經(jīng)有4人回復
Matlab非線性擬合（y=sqrt(1-(x-a).^2)）已經(jīng)有5人回復
求助origin非線性擬合的編譯（origin C）已經(jīng)有4人回復
畢不了業(yè)了，發(fā)帖求助Origin8.0多元非線性擬合的問題已經(jīng)有9人回復
用matlab里的nlinfit如何擬合分段函數(shù) 已經(jīng)有8人回復
二階非線性微分方程求解已經(jīng)有14人回復
matlab中關于非線性擬合，用nlinfit.程序總出錯。謝謝各位大俠已經(jīng)有3人回復
請教大蝦，一個分段函數(shù)的 matlab 程序?qū)崿F(xiàn) 已經(jīng)有4人回復
[求助]origin中非線性擬合中Guass和GuassAmp兩個函數(shù)的區(qū)別？已經(jīng)有4人回復
求助Origin中自定義求和函數(shù)的非線性擬合已經(jīng)有0人回復
【求助】急求關于matlab做非線性擬合的函數(shù) 已經(jīng)有8人回復
【求助】這個帶有分段函數(shù)偏微分方程怎么用for程序去解決已經(jīng)有14人回復
【求助】求助仿真中需要一個分段函數(shù)來控制電流源已經(jīng)有5人回復
【求助】關于origin進行非線性擬合的自定義函數(shù) 已經(jīng)有7人回復
【求助】關于分段函數(shù)的表示成一個式子的方法已經(jīng)有7人回復
【求助】多變量非線性擬合【已完結(jié)】已經(jīng)有8人回復
【求助】緊急求助有關MATHEMATICA的問題已經(jīng)有9人回復

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
計算強帖: 12
應助: 289 (大學生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學
管轄: 計算模擬

【答案】應助回帖

感謝參與，應助指數(shù) +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.')

贊一下

回復此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2樓2014-03-19 14:56:37

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉(zhuǎn) 我要訂閱樓主 hanmei98232 的主題更新