| 查看: 2604 | 回復(fù): 10 | |||
[求助]
需要用matlab擬合數(shù)據(jù)的外行人急需幫助
|
|
我有一函數(shù)y=f(x),里面有個(gè)未知數(shù)k,還有個(gè)x,y數(shù)值的txt文件,用這些條件求k 首先建立了一個(gè)M文件, function F=myfun(x,xdata) F=0.6518 + 0.7713./atan(xdata) - 0.0911*atan(xdata) + ((0.9121* xdata + 1.572)./(0.9121* xdata + 1.59 + 2.3026./ xdata + 0.0637*log(1 + 15.7./ xdata)) - 1)./((1 + 5.04*x* xdata).*(1 + (x^0.674)* (xdata.^0.173))); 然后就運(yùn)行下面的命令: >> load data.txt >> xdata=data(:,1); >> ydata=data(:,2); >> x=lsqcurvefit(@myfun,0.05,xdata,ydata) 然后就有錯(cuò)誤了 ??? Error using ==> qr Complex sparse QR is not yet available. Error in ==> aprecon at 57 RPCMTX = qr(TM(:,p)); Error in ==> trdog at 47 [R,permR] = feval(pcmtx,H,pcoptions,DM,DG,varargin{:}); Error in ==> snls at 346 [sx,snod,qp,posdef,pcgit,Z] = trdog(x,g,A,D,delta,dv,... Error in ==> lsqncommon at 285 [x,FVAL,LAMBDA,JACOB,EXITFLAG,OUTPUT,msg]=... Error in ==> lsqcurvefit at 186 [x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ... 這是哪里的錯(cuò)誤呢, lsqcurvefit 的運(yùn)用中有木有語(yǔ)法錯(cuò)誤,哪位大俠可以幫我解答一下,感激涕零,如果能幫我跑一下程序就最好不過(guò)了…… |

榮譽(yù)版主 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +4 |



鐵桿木蟲 (職業(yè)作家)
|
數(shù)據(jù)有很多 x=[25.9484 25.7484 25.5484 25.3484 25.1484 24.9484 24.7484 24.5484 24.3484 24.1484 23.9484 23.7484 23.5484 23.3484 23.1484 22.9484 22.7484 22.5484 22.3484 22.1484 21.9484 21.7484 21.5484 21.3484 21.1484 20.9484 20.7484 20.5484 20.3484 20.1484 19.9484 19.7484 19.5484 19.3484 19.1484 18.9484 18.7484 18.5484 18.3484 18.1484 17.9484 17.7484 17.5484 17.3484 17.1484 16.9484 16.7484 16.5484 16.3484 16.1484 15.9484 15.7484 15.5484 15.3484 15.1484 14.9484 14.7484 14.5484 14.3484 14.1484 13.9484 13.7484 13.5484 13.3484 13.1484 12.9484 12.7484 12.5484 12.3484 12.1484 11.9484 11.7484 11.5484 11.3484 11.1484 10.9484 10.7484 10.5484 10.3484 10.1484 9.9484 9.7484 9.5484 9.3484 9.1484 8.9484 8.7484 8.5484 8.3484 8.1484 7.9484 7.7484 7.5484 7.3484 7.1484 6.9484 6.7484 6.5484 6.3484 6.1484 5.9484 5.7484 5.5484 5.3484 5.1484 4.9484 4.7484 4.5484 4.3484 4.1484 3.9484 3.7484 3.5484 3.3484 3.1484 2.9484 2.7484 2.5484 2.3484 2.1498 2.003 1.8768 1.8236 1.7816 1.7138 1.632 1.5528 1.4828 1.4426 1.427 1.4128 1.4 ] y=[0.86396 0.86396 0.86384 0.86372 0.86372 0.86372 0.86348 0.86336 0.86336 0.86336 0.86324 0.86324 0.863 0.863 0.863 0.863 0.86288 0.86288 0.86288 0.86264 0.86252 0.8624 0.86228 0.86216 0.86204 0.86204 0.86192 0.8618 0.86156 0.86144 0.86144 0.86144 0.86132 0.8612 0.86108 0.86108 0.86108 0.86108 0.86084 0.86072 0.86072 0.8606 0.86036 0.86024 0.86 0.86 0.85975 0.85951 0.85951 0.85963 0.85963 0.85939 0.85939 0.85927 0.85915 0.85915 0.85915 0.85891 0.85867 0.85831 0.85843 0.85819 0.85807 0.85795 0.85783 0.85783 0.85771 0.85747 0.85747 0.85723 0.85699 0.85699 0.85687 0.85663 0.85627 0.85615 0.85579 0.85567 0.8553 0.85494 0.8547 0.85458 0.85434 0.85398 0.85374 0.8535 0.85326 0.85266 0.85242 0.85194 0.85109 0.85037 0.84965 0.84905 0.84845 0.84785 0.84664 0.8458 0.84448 0.8434 0.84243 0.84087 0.83834 0.83666 0.83462 0.83149 0.82896 0.82535 0.82295 0.82006 0.81525 0.81128 0.80611 0.79733 0.78807 0.77544 0.76004 0.73214 0.70375 0.66623 0.63375 0.58588 0.5498 0.53488 0.53007 0.5261 0.52009 0.50541 0.48677 0.4567 0.433 0.40907 ] |

金蟲 (正式寫手)


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +6 | 加一一九 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 0856材料專碩353求調(diào)劑 +5 | NIFFFfff 2026-03-20 | 5/250 |
|
|
[考研] 材料調(diào)劑 +5 | 匹克i 2026-03-23 | 5/250 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開(kāi)開(kāi)心心沒(méi)煩惱 2026-03-23 | 4/200 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +11 | hanamiko 2026-03-18 | 11/550 |
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +4 | 谷雨上岸 2026-03-23 | 5/250 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +7 | 生物工程調(diào)劑 2026-03-17 | 12/600 |
|
|
[考研]
求調(diào)劑材料學(xué)碩080500,總分289分
5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
[考研] 考研調(diào)劑 +4 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-21 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[考研] 一志愿南大,0703化學(xué),分?jǐn)?shù)336,求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|