當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

anarky

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 150
帖子: 16
在線: 8.2小時(shí)
蟲號(hào): 2655893
注冊: 2013-09-15
專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)

[求助] 請教matlab非線性擬合結(jié)果不唯一的問題已有2人參與

我用matlab擬合一組實(shí)驗(yàn)數(shù)據(jù)，函數(shù)是y=a*（b+x)^c+d*(1-f*exp(-g*x)),現(xiàn)在的問題是發(fā)現(xiàn)只要改一次初值就會(huì)得到不同的擬合結(jié)果，或者改一下上下限，也會(huì)得到不同的擬合結(jié)果，這么一來，我覺得就不可能找到最優(yōu)解，雖然上下限可以取得很大，認(rèn)為是無界，但是初值怎么取？因?yàn)槲也豢赡車L試所有的參數(shù)初值，因?yàn)槲乙獢M合多組實(shí)驗(yàn)數(shù)據(jù)，獲得公式中參數(shù)a(1)~a(6)的變化規(guī)律，如果參數(shù)結(jié)果不唯一，我怎么確定哪個(gè)參數(shù)才是最優(yōu)的參數(shù)呢？事實(shí)上，如果從物理意義來限定參數(shù)取值也可以，但是那樣得到的參數(shù)并不是最優(yōu)解，而最優(yōu)解的參數(shù)又不一定符合物理意義。我該怎么處理呢？哪位能幫忙看一下，非常感謝！
實(shí)驗(yàn)數(shù)據(jù)和代碼如下：
主程序
clc
xdata=load('78028-ta.txt');  %read section1
size1=size(xdata);             %section1 行、列數(shù)
sizen(1)=size1(1);             %section1數(shù)據(jù)量
for i=1:size1;
xd(i)=xdata(i,1);
yd(i)=xdata(i,2);
end

set(gca,'Fontsize',15);
plot(xd,yd,'ro');
xlabel('True Strain');
ylabel('True Stress/MPa');

hold on
lb=[-100000, 0, 0, -10000, -10000, -10000];
ub=[ 100000, 0.02, 1,  10000,  10000,  10000];
a1=1898;a2=0.00002;a3=0.136;a4=542;a5=0.784;a6=62.4;
a0=[a1,a2,a3,a4,a5,a6];
options = optimset('TolFun',1e-6,'LargeScale','on','TolCon',1.0e-6,'MaxFunEvals',2.e10,'MaxIter',2000000,'TolX',1e-6);
exitflag=0;
while (exitflag<=0)                               %反復(fù)迭代
[a,resnorm,residual,exitflag]=lsqnonlin(@SV,a0,lb,ub,options);    %nonlinear fitting, solve coefficient a
a0=a;
end
R2=1-resnorm/sum((yd-mean(yd)).^2);
standard=sqrt(resnorm/size(residual,2));
stp=1;
out(stp,1)=R2;
out(stp,2)=standard;
out(stp,3)=a(1);
out(stp,4)=a(2);
out(stp,5)=a(3);
out(stp,6)=a(4);
out(stp,7)=a(5);
out(stp,8)=a(6);

hold on
fid=fopen('SV-fit-78028-ta.txt','wt');
for i=1:size1;
   G(i)=a(1)*(a(2)+xd(i))^a(3)+a(4)*(1-a(5)*exp(-a(6)*xd(i)));
fprintf(fid,'%8.6f %10.6f\n',xd(i),G(i));
end

fprintf(fid,'    R2       standard    s1          s2          s3       V1       V2       V3\n');
fprintf(fid,'%10.6f %10.6f %10.6f %10.6f %10.6f %10.6f %10.6f %10.6f\n',out(i,1),out(i,2),out(i,3),out(i,4),out(i,5),out(i,6),out(i,7),out(i,8));
fclose(fid);
plot(xd,G,'b','LineWidth',2)
hold on
legend('EXP','SV');
hold off
clear
函數(shù)
function F=SV(a)
xdata=load('78028-ta.txt');
size1=size(xdata);
sizen(1)=size1(1);

for i=1:sizen(1);
xd(i)=xdata(i,1);
yd(i)=xdata(i,2);
end

for i=1:sizen(1);
F(i)=a(1)*(a(2)+xd(i))^a(3)+a(4)*(1-a(5)*exp(-a(6)*xd(i)));
F(i)=F(i)-yd(i);
end

實(shí)驗(yàn)數(shù)據(jù)
78028-ta.txt
0 223.634
0.0005 261.9
0.001 373.358
0.0015 425.658
0.002 462.61
0.0025 492.699
0.003 513.549
0.0035 524.142
0.004 538.681
0.0045 546.22
0.005 554.575
0.0055 563.233437
0.006 570.27162
0.0065 577.07317
0.007 583.656651
0.0075 590.037892
0.008 596.230537
0.0085 602.246461
0.009 608.09609
0.0095 613.788651
0.01 619.332371
0.0105 624.734633
0.011 630.002102
0.0115 635.14083
0.012 640.156342
0.0125 645.053702
0.013 649.837576
0.0135 654.512279
0.014 659.081814
0.0145 663.549911
0.015 667.920058
0.0155 672.19552
0.016 676.379368
0.0165 680.474496
0.017 684.483635
0.0175 688.409373
0.018 692.254165
0.0185 696.020341
0.019 699.710123
0.0195 703.325629
0.02 706.4514911
0.025 735.7217232
0.03 758.4120509
0.035 778.3482998
0.04 794.7567166
0.045 809.6868333
0.05 822.8422861
0.055 834.8128986
0.06 846.1687589
0.065 856.1567004
0.07 866.0042315
0.075 874.4658077
0.08 882.1951827
0.085 890.0445724
0.09 897.0900675
0.095 903.7413962
0.1 910.0615411
0.105 915.8285085
0.11 921.8803708
0.115 927.2834332
0.12 932.4186891
0.125 937.3558307
0.315 1073.901289

回復(fù)此樓

1樓 2015-12-07 21:45:34

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

anarky

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 150
帖子: 16
在線: 8.2小時(shí)
蟲號(hào): 2655893
注冊: 2013-09-15
專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)

引用回帖:

3樓: Originally posted by yanshuiting at 2015-12-08 12:36:50
Matlab尋優(yōu)過程得到的是局部最優(yōu)解，如果這個(gè)解能夠滿足實(shí)際需求，就可以使用，不用太過顧慮全局最優(yōu)解。當(dāng)模型復(fù)雜時(shí)，全局最優(yōu)解釋很難得到的。

感謝參與！

回復(fù)此樓

4樓2015-12-16 21:03:30

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

查看全部 5 個(gè)回答

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊: 2006-10-28

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

★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
anarky: 金幣+4, ★★★很有幫助, 非常感謝！ 2015-12-09 21:24:24

用1stOpt，不用初值，結(jié)果唯一穩(wěn)定：

CODE:

Function  y=a*(b+x)^c+d*(1-f*exp(-g*x));

Data;

0        223.634

0.0005        261.9

0.001        373.358

0.0015        425.658

0.002        462.61

0.0025        492.699

。。。。。

均方差(RMSE):6.83618029811499
殘差平方和(SSE):2944.20174730514
相關(guān)系數(shù)(R): 0.99915354759785
相關(guān)系數(shù)之平方(R^2): 0.998307811677368
確定系數(shù)(DC): 0.998307811677368
卡方系數(shù)(Chi-Square): 4.49199298531208
F統(tǒng)計(jì)(F-Statistic): 6736.83883025038

參數(shù) 最佳估算
---------- -------------
a -1838.97884865162
b 0.00549570600631707
c -0.0663243340918725
d 3047.85332211066
f 0.0810032510811905
g 817.309082573211

其它軟件很難得到上面的結(jié)果

c2.jpg

贊一下

回復(fù)此樓

2樓2015-12-08 10:18:31

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

yanshuiting

新蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1484.8
散金: 50
紅花: 7
帖子: 448
在線: 68.7小時(shí)
蟲號(hào): 2633300
注冊: 2013-09-03
專業(yè): 計(jì)算機(jī)科學(xué)

Matlab尋優(yōu)過程得到的是局部最優(yōu)解，如果這個(gè)解能夠滿足實(shí)際需求，就可以使用，不用太過顧慮全局最優(yōu)解。當(dāng)模型復(fù)雜時(shí)，全局最優(yōu)解釋很難得到的。

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

3樓2015-12-08 12:36:50

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

ybkooo

至尊木蟲 (著名寫手)

應(yīng)助: 32 (小學(xué)生)
金幣: 11360
散金: 178
紅花: 15
帖子: 2310
在線: 543.5小時(shí)
蟲號(hào): 2710519
注冊: 2013-10-09
性別: GG
專業(yè): 控制理論與方法

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

算法有局限性，可以考慮遺傳算法、粒子群算法；結(jié)果應(yīng)該穩(wěn)定

贊一下

回復(fù)此樓

abcd

5樓2015-12-16 21:54:02

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

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

[考研] 材料學(xué)碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無際的草原
[考研] 317求調(diào)劑 +9	申子申申 2026-03-19	15/750	2026-03-21 17:31 by 學(xué)員8dgXkO
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7	墨墨漠 2026-03-20	7/350	2026-03-21 16:36 by barlinike
[考研] 301求調(diào)劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷