涂涂木

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 169.3
散金: 7
紅花: 1
帖子: 322
在線: 38.6小時(shí)
蟲號(hào): 681538
注冊(cè): 2008-12-25
性別: MM
專業(yè): 催化化學(xué)

[求助] 最小二乘法擬合，發(fā)現(xiàn)，平方差和偏大，這是什么原因啊，請(qǐng)教各位可能的原因

最小二乘法擬合，發(fā)現(xiàn)，平方差和偏大，這是什么原因啊，請(qǐng)教各位可能的原因，
[~,S]=lsqnonlin(@taoget, tao0, [], [], [], C_i, t_i, beta_n, L, m);希望這個(gè)值小于10-6，但是運(yùn)算結(jié)果很大，造成這個(gè)原因是什么引起的啊，求教各位大俠了，我是菜鳥

回復(fù)此樓

» 猜你喜歡

本人考085602 化學(xué)工程專碩已經(jīng)有23人回復(fù)
求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(qǐng)(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專碩求調(diào)劑已經(jīng)有5人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【求助】最小二乘法的擬合問題已經(jīng)有4人回復(fù)

向前看

1樓 2012-09-17 16:24:33

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

baobiao007

木蟲 (職業(yè)作家)

中國特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時(shí)
蟲號(hào): 505962
注冊(cè): 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

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

★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝提示！ 2012-09-17 17:58:51

結(jié)果很大到底是多大？你確定自己的參數(shù)給的合理嗎？做擬合要求你對(duì)自己的數(shù)據(jù)十分了解，如果稀里糊涂獲得的數(shù)據(jù)十有八九做不好。

贊一下

回復(fù)此樓

我同意叔本華的觀點(diǎn)，人們投身藝術(shù)和科學(xué)領(lǐng)域的強(qiáng)烈愿望之一就是逃離痛苦、殘酷和枯燥無味的現(xiàn)實(shí)生活，逃離自己飄忽不定的七情六欲的桎梏。--愛因斯坦

2樓2012-09-17 16:46:20

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

nono2009

超級(jí)版主 (文學(xué)泰斗)

No gains, no pains.

專家經(jīng)驗(yàn): +21105
仿真EPI: 1
應(yīng)助: 28684 (院士)
貴賓: 513.911
金幣: 2555230
散金: 27828
紅花: 2148
沙發(fā): 66666
帖子: 1602255
在線: 65200.9小時(shí)
蟲號(hào): 827383
注冊(cè): 2009-08-13
性別: GG
專業(yè): 工程熱物理與能源利用
管轄: 科研家籌備委員會(huì)

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

★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝參與！ 2012-09-17 17:58:57

可能是擬合方程形式不合適吧。

贊一下

回復(fù)此樓

3樓2012-09-17 16:57:32

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

allen_baobao

至尊木蟲 (職業(yè)作家)

應(yīng)助: 102 (高中生)
金幣: 10405
散金: 3367
紅花: 31
帖子: 3801
在線: 3664.4小時(shí)
蟲號(hào): 1175925
注冊(cè): 2010-12-23
專業(yè): 控制理論與方法

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

★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝提示！ 2012-09-17 17:59:11

你的目標(biāo)函數(shù)選取得準(zhǔn)確么？

贊一下

回復(fù)此樓

Controlisahiddentechnology.——K.J.Astrom

4樓2012-09-17 17:39:59

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

allen_baobao

至尊木蟲 (職業(yè)作家)

應(yīng)助: 102 (高中生)
金幣: 10405
散金: 3367
紅花: 31
帖子: 3801
在線: 3664.4小時(shí)
蟲號(hào): 1175925
注冊(cè): 2010-12-23
專業(yè): 控制理論與方法

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

引用回帖:

4樓: Originally posted by allen_baobao at 2012-09-17 17:39:59
你的目標(biāo)函數(shù)選取得準(zhǔn)確么？

我的意思是要擬合的目標(biāo)曲線表達(dá)式。

贊一下

回復(fù)此樓

Controlisahiddentechnology.——K.J.Astrom

5樓2012-09-17 17:40:25

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

涂涂木

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 169.3
散金: 7
紅花: 1
帖子: 322
在線: 38.6小時(shí)
蟲號(hào): 681538
注冊(cè): 2008-12-25
性別: MM
專業(yè): 催化化學(xué)

一下程序是根據(jù)以上截圖編寫的，我印象中黃金分割法是求極小值或極大值的，運(yùn)行過程中，發(fā)現(xiàn)有點(diǎn)偏差，請(qǐng)教各位大俠，這里的黃金分割法求得L是否滿足圖3的運(yùn)行要求，
%% auxiliary equation %%
function y=auxiliary(x,L)
y=x*cot(x)+L-1;
%% Calculation of beta_n and F(i) %%
function S=F(C_i, t_i, L, m, n, tao0)
beta_n=ones(1,n);
   for i=1:n
      beta_n(i)=fzero(@auxiliary,[(i-1)*pi+1e-6,i*pi-1e-6], [], L);
   end
[~,S]=lsqnonlin(@taoget, tao0, [], [], [], C_i, t_i, beta_n, L, m);
%% Calculation of sum of squared error to run LSQNONLIN  %%
function f= taoget(tao, C_i, t_i, beta_n, L, m)
   f=ones(1,m);
for j=1:m
      ss=sum(2*L.*exp(-beta_n.^2.*t_i(j)./tao)./(beta_n.^2+L*(L-1)));
      f(j)=C_i(j)-ss;
end
clc,clear all;
t_i=[19.998 39.996 60 79.998 99.996 120 139.998 159.996 180 199.998 219.996 240 259.998 279.996 300 319.998 339.996 360 379.998,...
399.996 420 439.998 459.996 480 499.998 519.996 540 559.998 579.996 600 619.998 639.996 660 679.998 699.996 720 739.998 759.996 780];
C_i=[0.169560891 0.094901613 0.06222569 0.044350897 0.033182443 0.025755344 0.0205755344 0.016727315 0.01389813 0.01169039 0.009948977,...
0.008569505 0.007425149 0.006507127 0.005732515 0.005094485 0.004545233 0.0040711 0.003670136 0.003319903 0.003020399 0.002752115,...
0.002517 0.00230725 0.002129694 0.001969698 0.001819459 0.001687755 0.00157849 0.001466299 0.001364838 0.00127606 0.001191185 0.001118016,...
0.001051677 9.86313e-4 9.30704e-4 8.78023e-4 8.28268e-4];
% t(i) and C/C0(i) experimental
m=length(t_i);
% number of experimental points
n=10;
% number of beta_n(beta_1,beta_2,...,beta_n)
L1=0;
%  lower bound of L
L2=100;
%  upper bound of L
tao0=100;
%  initial guess of tao

%  goldensection method for L
c = 0.5*(3 - sqrt(5));
a = L1; b = L2; x1=L1; x2=L2;
tol=1e-6;
% tol is the error limit fof L
while (abs(x1-x2) > tol )
x1 = a+c*(b-a);
x2 = b - c*(b-a);
if F(C_i, t_i, x1, m, n, tao0)<= F(C_i, t_i, x2, m, n, tao0)
   b=x2;
else
   a=x1;
end
end

L=(x1+x2)/2;
% final result of L
S=F(C_i, t_i, L, m, n, tao0);
% final result of S(sum of squared error)
beta_n=ones(1,n);
   for i=1:n
      beta_n(i)=fzero(@auxiliary,[(i-1)*pi+1e-6,i*pi-1e-6], [], L);
   end
% final result of beta
tao=lsqnonlin(@taoget, tao0, [], [], [], C_i, t_i, beta_n, L, m);
% final result of tao

L,S,beta_n,tao

1.jpg

2.jpg

3.jpg

贊一下

回復(fù)此樓

向前看

6樓2012-09-17 20:49:45

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

涂涂木

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 169.3
散金: 7
紅花: 1
帖子: 322
在線: 38.6小時(shí)
蟲號(hào): 681538
注冊(cè): 2008-12-25
性別: MM
專業(yè): 催化化學(xué)

引用回帖:

6樓: Originally posted by 涂涂木 at 2012-09-17 20:49:45
一下程序是根據(jù)以上截圖編寫的，我印象中黃金分割法是求極小值或極大值的，運(yùn)行過程中，發(fā)現(xiàn)有點(diǎn)偏差，請(qǐng)教各位大俠，這里的黃金分割法求得L是否滿足圖3的運(yùn)行要求，
%% auxiliary equation %%
function y=auxili ...

截圖1為，擬合的函數(shù)表達(dá)式，各位大俠指點(diǎn)下我啊

贊一下

回復(fù)此樓

向前看

7樓2012-09-17 20:54:13

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

涂涂木

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 169.3
散金: 7
紅花: 1
帖子: 322
在線: 38.6小時(shí)
蟲號(hào): 681538
注冊(cè): 2008-12-25
性別: MM
專業(yè): 催化化學(xué)

我現(xiàn)在有一點(diǎn)想法，擬合的方程是一個(gè)減函數(shù)��！就是這個(gè)黃金分割法這部分程序中，縮小區(qū)間的判據(jù)，
if F(C_i, t_i, x1, m, n, tao0)<= F(C_i, t_i, x2, m, n, tao0)
改成使得abs(c_i-F(C_i,t_i,X1,m,n,tao0)<=abs(c-i-F(c_i,t_i,x2,m,n,tao0),即按照文獻(xiàn)判斷準(zhǔn)則來，各位大俠你們幫我指點(diǎn)一下啊

3.jpg

贊一下

回復(fù)此樓

向前看

8樓2012-09-18 09:46:57

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

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 286求調(diào)劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 資源與環(huán)境調(diào)劑申請(qǐng)(333分) +5	holy J 2026-03-21	5/250	2026-03-21 22:42 by Catalysis25
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 326求調(diào)劑 +5	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	8/400	2026-03-21 19:33 by ColorlessPI
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +3	石瑞0426 2026-03-19	3/150	2026-03-21 18:59 by RRRKKK
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 268求調(diào)劑 +9	簡(jiǎn)單點(diǎn)0 2026-03-17	9/450	2026-03-21 15:37 by lature00
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 316求調(diào)劑 +6	梁茜雯 2026-03-19	6/300	2026-03-21 06:32 by Ecowxq666！
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 274求調(diào)劑 +10	S.H1 2026-03-18	10/500	2026-03-20 23:51 by JourneyLucky
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 323求調(diào)劑 +3	洼小桶 2026-03-18	3/150	2026-03-20 22:54 by JourneyLucky
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	5/250	2026-03-20 22:11 by 云游重陽
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678