| 2 | 1/1 | 返回列表 |
| 查看: 277 | 回復: 1 | |||
[交流]
【求助】lsqcurvefit非線性擬合的問題 已有1人參與
|
|
function F=myfun(p,data) F=p(1)*(xda(1).^p(2))*(xda(2).^p(3))*(xda(3).^p(4))+p(5)*(xda(4).^p(6)) y=[41.894 115.712 58.319 57.324 67.9495 130.542 71.109 74.447 93.078 128.794 128.573 124.831 107.041 141.092 174.634 217.774 183.355 177.148]; xda(1)=[0.1 0.1 0.1 0.1 0.1 0.1 0.2 0.2 0.2 0.2 0.2 0.2 0.3 0.3 0.3 0.3 0.3 0.3]; xda(2)=[0.1 0.1 0.15 0.15 0.2 0.2 0.1 0.1 0.15 0.15 0.2 0.2 0.1 0.1 0.15 0.15 0.2 0.2 ]; xda(3)=[200 250 200 250 200 250 200 300 200 250 200 300 200 300 250 200 250 300]; xda(4)=[0.095 0.425 0.228 0.124 0.127 0.427 0.110 0.148 0.133 0.302 0.125 0.149 0.134 0.309 0.312 0.315 0.128 0.235]; data=[xda(1);xda(2);xda(3);xda(4)]; p0=[1,2,3,4,5,6] [p,resnorm]=lsqcurvefit(@myfun,p0,data,y) 請教各位大俠為什么我的程序運行不了,提示 “??? Input argument "x" is undefined.” |
銀蟲 (初入文壇)
|
我看完你的程序,提出以下幾點看法: 第一,程序的前兩行應該是單獨的一個m文件; 第二,你有xda(1),...,xda(4)分別表示矩陣,這里就有問題; 第三,data和y矩陣的形式應該是一致的,在你的程序里data是12*6,而y是3*6,這肯定不行。 你用xda(1),...,xda(4)是否要做4次計算?如果是這樣的話,可以先給出一個12*6矩陣,然后每三行取值做計算,這樣循環(huán)4次就可以了。 下面是我自己試著做的,僅供參考: function F=myfun(x,xdata) F=x(1)*(xdata.^x(2)).*(xdata.^x(3)).*(xdata.^x(4))+x(5)*(xdata.^x(6)); clear;clc ydata=[41.894 115.712 58.319 57.324 67.9495 130.542 71.109 74.447 93.078 128.794 128.573 124.831 107.041 141.092 174.634 217.774 183.355 177.148]; data=[0.1 0.1 0.1 0.1 0.1 0.1 0.2 0.2 0.2 0.2 0.2 0.2 0.3 0.3 0.3 0.3 0.3 0.3 0.1 0.1 0.15 0.15 0.2 0.2 0.1 0.1 0.15 0.15 0.2 0.2 0.1 0.1 0.15 0.15 0.2 0.2 200 250 200 250 200 250 200 300 200 250 200 300 200 300 250 200 250 300 0.095 0.425 0.228 0.124 0.127 0.427 0.110 0.148 0.133 0.302 0.125 0.149 0.134 0.309 0.312 0.315 0.128 0.235]; x0=[1,2,3,4,5,6]; for i=1:4 fprintf('Case %d:',i) xdata=[data(3*i-2,: );data(3*i-1,: );data(3*i,: )]; [x,resnorm]=lsqcurvefit(@myfun,x0,xdata,ydata) end [ Last edited by anyuezhiji on 2010-7-5 at 17:30 ] |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請求調劑愿意跨考! +3 | lemonzzn 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 0856調劑,是學校就去 +6 | sllhht 2026-03-19 | 7/350 |
|
|
[考研] 【考研調劑】化學專業(yè) 281分,一志愿四川大學,誠心求調劑 +5 | 吃吃吃才有意義 2026-03-19 | 5/250 |
|
|
[考研]
|
.6lL 2026-03-18 | 6/300 |
|
|
[考研] 求調劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
[考研] 281求調劑(0805) +9 | 煙汐憶海 2026-03-16 | 19/950 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 一志愿985,本科211,0817化學工程與技術319求調劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
[考研] 332求調劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[考研] 0703化學 305求調劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 26調劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 085601專碩,總分342求調劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 考研化學學碩調劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 一志愿南京大學,080500材料科學與工程,調劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 一志愿,福州大學材料專碩339分求調劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 283求調劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 274求調劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
[考研] 333求調劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 297求調劑 +4 | 學海漂泊 2026-03-13 | 4/200 |
|