版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

4758188748

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 92
帖子: 1
在線: 1.6小時
蟲號: 3972860
注冊: 2015-07-15
性別: GG
專業(yè): 函數(shù)論

[求助] 求助。matlab共享參數(shù)類的自定義公式擬合修改問題已有1人參與

是共享參數(shù)類的自定義公式模擬，目的是求5個系數(shù)，y1-x1對應(yīng)q0固定系數(shù)，y2-x2對應(yīng)q1固定系數(shù)，依次類推四個系數(shù)。我這個代碼是肯定不行呀- -，想要知道這個能不能求幫忙修改，并且把這幾個曲線繪制到一張圖上，謝謝謝謝謝謝

clear
clc
data = [36.287 34.765 38.875 37.683 33.580 27.010 34.720 35.470
42.726 37.060 38.475 42.038 31.388 30.550 33.480 37.460
44.688 38.646 40.391 43.510 30.339 24.934 25.711 29.321
43.517 39.432 41.318 42.936 29.632 26.228 27.102 28.644
44.222 38.646 42.530 41.718 28.211 24.991 27.792 29.638
43.771 40.398 42.454 44.425 27.862 25.291 26.591 27.814]
q0=0;
q1=0.2;
q2=0.35;
q3=0.5;
x1 = data(:,1);
x2 = data(:,2);
x3 = data(:,3);
x4 = data(:,4);
y1 = data(:,5);
y2 = data(:,6);
y3 = data(:,7);
y4 = data(:,8);
f = @(a,b,c,d,e)[a.*x1./(b*q0^3+c*q0^2+d*q0+e);
a.*x2./(b*q1^3+c*q1^2+d*q1+e);
a.*x3./(b*q2^3+c*q2^2+d*q2+e);
a.*x4./(b*q3^3+c*q3^2+d*q3+e)]
[X,RESNORM,RESIDUAL,EXITFLAG,OUTPUT,LAMBDA] = lsqcurvefit(f, [1 1 1 1 1 ], [x1;x2;x3;x4],[y1;y2;y3;y4],[0 0 0 0 0 0],[],[])

回復(fù)此樓

» 本帖@通知

@死天 @liyaobo0509 @mycc @shgao20 @wuming524

» 猜你喜歡

311求調(diào)劑已經(jīng)有5人回復(fù)
274求調(diào)劑0856材料化工已經(jīng)有4人回復(fù)
308求調(diào)劑已經(jīng)有3人回復(fù)
271求調(diào)劑已經(jīng)有3人回復(fù)
材料工程085601調(diào)劑求老師收留已經(jīng)有10人回復(fù)
288求調(diào)劑已經(jīng)有18人回復(fù)
接受26屆調(diào)劑生已經(jīng)有22人回復(fù)
一志愿北京化工大學(xué)材料與化工296分求調(diào)劑已經(jīng)有16人回復(fù)
290求調(diào)劑已經(jīng)有8人回復(fù)
伙伴們，祝我生日快樂吧已經(jīng)有3人回復(fù)

1樓 2019-01-27 00:33:50

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

hollya

至尊木蟲 (文壇精英)

葛城?{美

應(yīng)助: 744 (博后)
金幣: 44165.3
散金: 1143
紅花: 127
沙發(fā): 780
帖子: 24910
在線: 835小時
蟲號: 458790
注冊: 2007-11-13
性別: MM
專業(yè): 計算機(jī)軟件

【答案】應(yīng)助回帖

樓主，你這個函數(shù)好像就是個線性方程，f(x)=kx, 這樣擬合結(jié)果可能不太好啊。而且就算有5個參數(shù)需要似合，5個也不是相互獨立的。
你看看下面的這個：其中 k1=b/a,k2=c/a,k3=d/a,k4=e/a

data = [36.287       34.765       38.875       37.683       33.580       27.010       34.720       35.470
      42.726       37.060       38.475       42.038       31.388       30.550       33.480       37.460
      44.688       38.646       40.391       43.510       30.339       24.934       25.711       29.321
      43.517       39.432       41.318       42.936       29.632       26.228       27.102       28.644
      44.222       38.646       42.530       41.718       28.211       24.991       27.792       29.638
      43.771       40.398       42.454       44.425       27.862       25.291       26.591       27.814]
q0=0;
q1=0.2;
q2=0.35;
q3=0.5;
q=[q0,q1,q2,q3]
Q=[q.^3;q.^2;q;ones(1,4)]
xdata=data(:,1:4)
ydata=data(:,5:8)
f = @(k,x) x./repmat([k(1),k(2),k(3),k(4)]*Q,6,1)
k0=[1 1 1 1]
k = lsqcurvefit(f,k0,xdata,ydata)
f_fit=f(k,xdata)
plot(xdata(:,1),ydata(:,1),'o',xdata(:,1),f_fit(:,1),...
xdata(:,2),ydata(:,2),'*',xdata(:,2),f_fit(:,2),...
xdata(:,3),ydata(:,3),'+',xdata(:,3),f_fit(:,3),...
xdata(:,4),ydata(:,4),'x',xdata(:,4),f_fit(:,4))

贊一下

回復(fù)此樓

でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。

2樓2019-02-02 02:45:25

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

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