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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

dingyy1987

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1088
紅花: 2
帖子: 56
在線: 35.7小時
蟲號: 2203357
注冊: 2012-12-24
性別: GG
專業(yè): 爆炸與沖擊動力學(xué)

[求助] 哪位大神有istopt，我網(wǎng)上下了一個不能運(yùn)行sharedmodel，求幫我運(yùn)行一下已有2人參與

Title "123";
Parameter c1,c2,c3,c4,c5;
Variable x,z,f;
SharedModel;
Function
z=(x^c2/c3/c4)^(1/c4);
f=c1*(x^c2/c3/c4)^(1/c4)*exp(-1/c4)+c5*(1-exp(-1/c4));
data;
14857.77758 0.0414 8.94345E6
18553.17237 0.05291 8.99335E6
20174.23991 0.06071 8.98066E6
24127.20834 0.07239 8.98738E6
24273.43798 0.07442 9.01031E6
26951.77543 0.08156 9.05148E6
20564.9827 0.07419 8.59135E6
19910.72847 0.07655 8.38573E6
18344.07188 0.07431 8.19985E6
16706.67637 0.07078 8.1676E6
16478.50121 0.0693 8.33608E6
22360.23753 0.08006 8.7983E6
21293.07227 0.07569 8.8397E6
20308.9913 0.07284 8.78015E6
23946.8144 0.08034 8.8947E6
20146.16809 0.07223 8.61393E6
17522.67615 0.06795 8.41823E6
17196.59303 0.06577 8.45751E6
17575.12863 0.06553 8.43412E6
15275.12691 0.0583 8.35688E6
18133.81317 0.06291 8.55274E6
18676.61758 0.06373 8.62984E6
16488.34249 0.05812 8.47193E6
16964.95475 0.06073 8.4801E6
17575.16675 0.06188 8.50543E6
17177.05496 0.05957 8.52406E6
15914.92064 0.05696 8.523E6
14213.56581 0.05257 8.42825E6
15066.24199 0.05465 8.42138E6
15927.43786 0.05661 8.44597E6
13486.53921 0.0505 8.27752E6
13013.96774 0.04854 8.32841E6
13728.11816 0.05051 8.39552E6
14144.65298 0.05275 8.41224E6
14954.59782 0.05585 8.40791E6
13582.45324 0.05299 8.31854E6
12648.72128 0.0508 8.27208E6
11698.11831 0.04785 8.21884E6
11582.30157 0.04853 8.22806E6
10380.45227 0.04524 8.15268E6
11064.31308 0.04744 8.19212E6
11121.39599 0.04783 8.15213E6
11078.61898 0.04852 8.08062E6
9802.8059 0.0446 8.00574E6
7929.85965 0.04011 7.80904E6
9503.21012 0.047 7.85323E6
8308.91314 0.04326 7.79567E6
8082.84269 0.04295 7.78816E6
7749.99714 0.04217 7.72199E6
6839.21332 0.03864 7.56499E6
6910.27878 0.03997 7.51553E6
6940.63115 0.04122 7.44532E6
6401.43753 0.04093 7.32896E6
5752.16969 0.0386 7.21224E6
4958.64579 0.03575 7.08191E6
4917.3898 0.03696 7.00645E6
4279.60599 0.03513 6.83156E6
3229.51003 0.02891 6.70232E6
2500.59496 0.02539 6.50041E6
2086.75039 0.02345 6.29762E6
2273.74314 0.02434 6.3008E6

回復(fù)此樓

» 猜你喜歡

2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有6人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有88人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有65人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)

1樓 2016-10-25 13:09:45

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

dingd

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

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

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

感謝參與，應(yīng)助指數(shù) +1

不能運(yùn)行是版本問題吧。
z和f相差的太大了，很難保證兩個擬合的都好。

贊一下

回復(fù)此樓

2樓2016-10-25 14:39:41

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

yz457694

鐵桿木蟲 (正式寫手)

計算強(qiáng)帖: 1
應(yīng)助: 28 (小學(xué)生)
金幣: 7192
散金: 1000
紅花: 16
帖子: 721
在線: 168.7小時
蟲號: 2820058
注冊: 2013-11-22
性別: GG
專業(yè): 傳遞過程

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
dingyy1987: 金幣+50, ★★★★★最佳答案, 很有幫助，而且也很熱心 2016-10-25 20:48:03

用matlab給你編了個程序
主程序:

CODE:

%共享參數(shù)擬合

%原理：非線性最小二乘擬合

clear,clc

format long

data=[

14857.77758  0.0414  8.94345E6

18553.17237  0.05291  8.99335E6

20174.23991  0.06071  8.98066E6

24127.20834  0.07239  8.98738E6

24273.43798  0.07442  9.01031E6

26951.77543  0.08156  9.05148E6

20564.9827  0.07419  8.59135E6

19910.72847  0.07655  8.38573E6

18344.07188  0.07431  8.19985E6

16706.67637  0.07078  8.1676E6

16478.50121  0.0693  8.33608E6

22360.23753  0.08006  8.7983E6

21293.07227  0.07569  8.8397E6

20308.9913  0.07284  8.78015E6

23946.8144  0.08034  8.8947E6

20146.16809  0.07223  8.61393E6

17522.67615  0.06795  8.41823E6

17196.59303  0.06577  8.45751E6

17575.12863  0.06553  8.43412E6

15275.12691  0.0583  8.35688E6

18133.81317  0.06291  8.55274E6

18676.61758  0.06373  8.62984E6

16488.34249  0.05812  8.47193E6

16964.95475  0.06073  8.4801E6

17575.16675  0.06188  8.50543E6

17177.05496  0.05957  8.52406E6

15914.92064  0.05696  8.523E6

14213.56581  0.05257  8.42825E6

15066.24199  0.05465  8.42138E6

15927.43786  0.05661  8.44597E6

13486.53921  0.0505  8.27752E6

13013.96774  0.04854  8.32841E6

13728.11816  0.05051  8.39552E6

14144.65298  0.05275  8.41224E6

14954.59782  0.05585  8.40791E6

13582.45324  0.05299  8.31854E6

12648.72128  0.0508  8.27208E6

11698.11831  0.04785  8.21884E6

11582.30157  0.04853  8.22806E6

10380.45227  0.04524  8.15268E6

11064.31308  0.04744  8.19212E6

11121.39599  0.04783  8.15213E6

11078.61898  0.04852  8.08062E6

9802.8059  0.0446  8.00574E6

7929.85965  0.04011  7.80904E6

9503.21012  0.047  7.85323E6

8308.91314  0.04326  7.79567E6

8082.84269  0.04295  7.78816E6

7749.99714  0.04217  7.72199E6

6839.21332  0.03864  7.56499E6

6910.27878  0.03997  7.51553E6

6940.63115  0.04122  7.44532E6

6401.43753  0.04093  7.32896E6

5752.16969  0.0386  7.21224E6

4958.64579  0.03575  7.08191E6

4917.3898  0.03696  7.00645E6

4279.60599  0.03513  6.83156E6

3229.51003  0.02891  6.70232E6

2500.59496  0.02539  6.50041E6

2086.75039  0.02345  6.29762E6

2273.74314  0.02434  6.3008E6];

x=data(:,1);

ze=data(:,2);

fe=data(:,3);

fun=@(c)myfun(x,c,ze,fe);

c0=[7.55629E8,0.18679,46.54595,0.35525,6.06252E6]';%初值

options = optimset('Display','on','MaxFunEvals',1000,'MaxIter',1000);

c=lsqnonlin(fun,c0,[],[],options)

zc=(x.^c(2)/c(3)/c(4)).^(1/c(4));

fc=c(1)*(x.^c(2)/c(3)/c(4)).^(1/c(4)).*exp(-1/c(4))+c(5)*(1-exp(-1/c(4)));

%相關(guān)系數(shù)

ye=[ze;fe];

yc=[zc;fc];

R=corrcoef(ye,yc)

%繪圖

subplot(2,1,1);

scatter(x,ze,'o');

hold on

scatter(x,zc,'+');

subplot(2,1,2);

scatter(x,fe,'o');

hold on

scatter(x,fc,'+')

函數(shù)

CODE:

function y=myfun(x,c,ze,fe)

    zc=(x.^c(2)/c(3)/c(4)).^(1/c(4));

    fc=c(1)*(x.^c(2)/c(3)/c(4)).^(1/c(4)).*exp(-1/c(4))+c(5)*(1-exp(-1/c(4)));

    n=length(zc);

    i=1:n;

    y=(zc-ze).^2+(fc-fe).^2;

結(jié)果：
c =
1.0e+008 *
7.556291199986442
0.000000001867029
0.000000465455226
0.000000003562897
0.060625164104285
R=0.999425581659652
見圖1

哪位大神有istopt，我網(wǎng)上下了一個不能運(yùn)行sharedmodel，求幫我運(yùn)行一下

圖1.jpg

贊一下(4人)

回復(fù)此樓

3樓2016-10-25 16:17:27

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

dingyy1987

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1088
紅花: 2
帖子: 56
在線: 35.7小時
蟲號: 2203357
注冊: 2012-12-24
性別: GG
專業(yè): 爆炸與沖擊動力學(xué)

引用回帖:

3樓: Originally posted by yz457694 at 2016-10-25 16:17:27
用matlab給你編了個程序
主程序:

%共享參數(shù)擬合
%原理：非線性最小二乘擬合
clear,clc
format long
data=;
x=data(:,1);
ze=data(:,2);
fe=data(:,3);
fun=@(c)myfun(x,c,ze,fe);
c0=';%初值
options ...

很厲害，請問一下，你初始參數(shù)怎么定的？
每次我改初始參數(shù)后，擬合效果就很差，就是說本程序嚴(yán)重依賴初始參數(shù)的預(yù)定。
不知道是誤差規(guī)定的，還是迭代不夠？?

贊一下

回復(fù)此樓

4樓2016-10-25 20:00:36

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

yz457694

鐵桿木蟲 (正式寫手)

計算強(qiáng)帖: 1
應(yīng)助: 28 (小學(xué)生)
金幣: 7192
散金: 1000
紅花: 16
帖子: 721
在線: 168.7小時
蟲號: 2820058
注冊: 2013-11-22
性別: GG
專業(yè): 傳遞過程

你在這個值附近修改結(jié)果變差是因為達(dá)到了一些邊界條件，比如迭代過程MaxFunEvals，MaxIter等參數(shù)超過了預(yù)設(shè)的值，運(yùn)行結(jié)束的時候命令窗口會提示擬合結(jié)束的原因，你可以相應(yīng)的在options那一行代碼里面增強(qiáng)一下迭代參數(shù)，結(jié)果就還是會收斂在這個值(你的數(shù)據(jù)z,f之間量級差太大了，所以迭代參數(shù)可以調(diào)高一些)。這個初值開始也很頭大，后來觀察發(fā)現(xiàn)z中只含有c2,c3,c4，f中全都包含，就想試了試先通過z-x擬合確定c2,c3,c4，然后固定c2,c3,c4，擬合f-x,結(jié)果收斂還不錯(通過origin實(shí)現(xiàn))，然后作為初值用matlab來整體擬合。

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

贊一下(1人)

回復(fù)此樓

5樓2016-10-25 20:31:44

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

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

返回列表

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

[考研] 317求調(diào)劑 +4	申子申申 2026-03-19	8/400	2026-03-20 11:20 by 申子申申
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +5	Ncdx123456 2026-03-19	5/250	2026-03-20 11:15 by wangy0907
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +3	lemonzzn 2026-03-17	4/200	2026-03-20 11:04 by lemonzzn
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +7	Doleres 2026-03-19	7/350	2026-03-20 10:39 by EBSD
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	14/700	2026-03-20 09:30 by kkcoco25
[考研] 304求調(diào)劑 +5	曼殊2266 2026-03-18	5/250	2026-03-20 09:00 by ZHANG0tao
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 307求調(diào)劑 +9	冷笙123 2026-03-17	9/450	2026-03-19 22:44 by 學(xué)員8dgXkO
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實(shí)用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

dingyy1987

[求助] 哪位大神有istopt，我網(wǎng)上下了一個不能運(yùn)行sharedmodel，求幫我運(yùn)行一下 已有2人參與

» 猜你喜歡

dingd

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

yz457694

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

dingyy1987

yz457694

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 哪位大神有istopt，我網(wǎng)上下了一個不能運(yùn)行sharedmodel，求幫我運(yùn)行一下已有2人參與