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

登錄注冊

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

24小時(shí)熱門版塊排行榜

返回列表

4758188748

新蟲 (初入文壇)

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

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

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

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

» 猜你喜歡

2026級(jí)碩士研究生招生/調(diào)劑已經(jīng)有5人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有7人回復(fù)
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生已經(jīng)有32人回復(fù)
0856材料與化工353分求調(diào)劑已經(jīng)有4人回復(fù)
282分材料專業(yè)求調(diào)劑院校已經(jīng)有13人回復(fù)
315食品工程求調(diào)劑已經(jīng)有5人回復(fù)
環(huán)境調(diào)劑已經(jīng)有3人回復(fù)
接受26屆調(diào)劑生已經(jīng)有20人回復(fù)
考研材料與化工，求調(diào)劑已經(jīng)有3人回復(fù)
歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組已經(jīng)有35人回復(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小時(shí)
蟲號(hào): 458790
注冊: 2007-11-13
性別: MM
專業(yè): 計(jì)算機(jī)軟件

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

樓主，你這個(gè)函數(shù)好像就是個(gè)線性方程，f(x)=kx, 這樣擬合結(jié)果可能不太好啊。而且就算有5個(gè)參數(shù)需要似合，5個(gè)也不是相互獨(dú)立的。
你看看下面的這個(gè)：其中 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 的主題更新