 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 912 | 回復(fù): 0 | ||
xxz903金蟲 (小有名氣)
|
[求助]
matlab數(shù)據(jù)擬合求助！求出來的是復(fù)數(shù) 應(yīng)該怎么修改？
|
|
function F=unique(x,xdata) %x=[] t12=x(1); t21=x(2); q1=x(3); q2=x(4); r1=x(5); r2=x(6); h=x(7); m=xdata; z=10; x1=m./(m+h); x2=1-x1; o1=x1*r1./(x1*r1+x2*r2); o2=x2*r2./(x1*r1+x2*r2); b1=x1*q1./(x1*q1+x2*q2); b2=x2*q2./(x1*q1+x2*q2); l1=0.5*z*(r1-q1)-(r1-1); l2=0.5*z*(r2-q2)-(r2-1); y1=log(o1./x1)+0.5*z*q1*log(b1./o1)+o2*(l1-r1*l2/r2)-q1*log(b1+b2*t21)... +b2*q1.*(t21./(b1+b2*t21)-t12./(b2+b1*t12))-log(r1/r2)... -05*z*q1*log(q1*r2/(q2*r1))-(l1-r1*l2/r2)+q1*log(t21)-q1*(1-t12); y2=log(o2./x2)+0.5*z*q2*log(b2./o2)+o1*(l2-r2*l1/r1)-q2*log(b2+b1*t12)... +b1*q2.*(t12./(b2+b1*t12)-t21./(b1+b2*t21))-log(r2/r1)... -05*z*q2*log(q2*r1/(q1*r2))-(l2-r2*l1/r1)+q2*log(t12)-q2*(1-t21); F=[y1,y2]; end function xxz2 data=[ 0.1 -0.673344553 -0.003214304 0.2 -0.783071888 -0.005447496 0.3 -0.841647189 -0.008713882 0.4 -0.869884359 -0.010990339 0.5 -0.884307686 -0.014309599 0.6 -0.889162064 -0.018687569 0.7 -0.889162064 -0.02207116 0.8 -0.88673193 -0.025484701 0.9 -0.879476759 -0.029975034 1 -0.869884359 -0.038731848 1.2 -0.850971266 -0.042636996 1.4 -0.830113036 -0.051995295 1.6 -0.807436327 -0.061533872 1.8 -0.783071888 -0.071257081 2 -0.759286983 -0.081169452 2.5 -0.701179352 -0.107972916 3 -0.650087691 -0.136153773 3.5 -0.599656837 -0.165810242 4 -0.553385238 -0.197052419 4.5 -0.512493681 -0.227286842 5 -0.47320876 -0.260553522]; %xdata=data(:,1); %Im %ydata=data(:,2:3); % %x0=ones(1,7); %t12,t21,q1,q2,r1,r2,h %options= optimset('MaxFunEvals',5000,'Maxiter',1000); %lb=[10,20,0,0]; %ub=[10000,10000,100,1000]; %lb=[]; %ub=[]; %[x,resnorm]=lsqcurvefit(@unique,x0,xdata,ydata,lb,ub,options) %y=unique(x,xdata); %plot(xdata,ydata,'g',xdata,y,'r') %end xdata=data(:,1); %Im ydata=data(:,2:3); % options= optimset('MaxFunEvals',5000,'Maxiter',1000); lb=-ones(1,7); ub=ones(1,7); t=1; for i=1:1:101 x0=rand(1,7); %t12,t21,q1,q2,r1,r2,h [x(t,  ,resnorm(t)]=lsqcurvefit(@myfun,x0,xdata,ydata,lb,ub,options);t=t+1; end m=find(resnorm==min(resnorm)); xm=x(m, y=myfun(xm,xdata); plot(xdata,ydata,'g',xdata,y,'r') |
回復(fù)此樓» 猜你喜歡
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有81人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有137人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
跪求1stOpt2.5以上版本，或者請(qǐng)高人幫忙擬合包含復(fù)數(shù)的函數(shù)的參數(shù)
已經(jīng)有8人回復(fù)
- matlab擬合求值
已經(jīng)有17人回復(fù)
- 需要用matlab擬合數(shù)據(jù)的外行人急需幫助
已經(jīng)有10人回復(fù)
- 求助大俠幫忙，有關(guān)matlab模擬數(shù)據(jù)
已經(jīng)有29人回復(fù)
- 求各位前輩幫忙！matlab解復(fù)數(shù)方程，出現(xiàn)圖中的警告，怎來修改呢？謝謝了啊
已經(jīng)有5人回復(fù)
- 【求助】使用MATLAB擬合出按照公式的系數(shù)
已經(jīng)有19人回復(fù)
- 【求助】matlab曲面擬合表達(dá)式
已經(jīng)有10人回復(fù)
- 【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過一個(gè)點(diǎn)
已經(jīng)有8人回復(fù)
- 【求助】使用MATLAB怎么實(shí)現(xiàn)擬合力場參數(shù)的程序？【已解決】
已經(jīng)有15人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲開始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xxz903 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +11 | Faune 2026-03-21 | 11/550 |
|
|
|
[考研] 07化學(xué)280分求調(diào)劑 +7 | 722865 2026-03-23 | 7/350 |
|
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
|
[考研] 291求調(diào)劑 +3 | HanBeiNingZC 2026-03-24 | 3/150 |
|
|
|
[考研] 321求調(diào)劑 +4 | Ymlll 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +5 | kotoko_ik 2026-03-23 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +7 | vv迷 2026-03-22 | 7/350 |
|
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|