wangdeqiang

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1025.8
散金: 10
帖子: 67
在線: 25.5小時
蟲號: 1931929
注冊: 2012-08-09
專業(yè): 金屬材料的力學(xué)行為

[求助] 1stopt擬合中出現(xiàn)的問題已有1人參與

Title "DIC-K";
Parameters kI,kII,b;
constant n=3,g=8.1*10^4,u=0.26,k=2.17;
Variable x,y,z;
Function z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2)))
-sum((kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b;
Data;
376.78243 -112.14027 9.73371
374.92666 -111.4325 9.61182
373.12866 -110.7178 9.49729
371.38928 -109.99629 9.39441
369.70935 -109.26811 9.3051
368.08966 -108.5334 9.22794
366.53104 -107.79231 9.15739
365.03425 -107.04503 9.08779
363.60006 -106.29172 9.0149
本人想通過1stopt求解kI，kII，但是運(yùn)行時系統(tǒng)提示參數(shù)kI數(shù)值范圍定義有誤!
請問該如何寫程序？

回復(fù)此樓

» 猜你喜歡

計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有295人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有29人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有26人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有9人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名已經(jīng)有0人回復(fù)

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

求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有6人回復(fù)
1stopt自變量超多時是否有簡單語法？AND 自變量多數(shù)據(jù)量大的話1stopt都會算的特慢么已經(jīng)有7人回復(fù)
求大神跑下1stOpt，一個sharemodel，謝謝！已經(jīng)有7人回復(fù)
1stopt擬合出參數(shù)為0，為什么？已經(jīng)有6人回復(fù)
求助--用1stopt擬合藍(lán)寶石襯底的GaN薄膜紅外反射光譜已經(jīng)有6人回復(fù)
求助：以下非線性擬合方程可以用1stopt求解么？用什么方法已經(jīng)有16人回復(fù)
用1stOpt進(jìn)行復(fù)數(shù)曲線擬合時，結(jié)果出錯。急求教！�。� 已經(jīng)有5人回復(fù)
1stopt求最小值已經(jīng)有37人回復(fù)
在用origin擬合的過程中，程序編寫這一塊實在不會，求指教！已經(jīng)有10人回復(fù)
origin中如何自定義Birch–Murnaghan三階擬合已經(jīng)有40人回復(fù)
生產(chǎn)中出了點(diǎn)問題求解已經(jīng)有11人回復(fù)
1stopt中使用sharedmodel擬合曲線出現(xiàn)問題，求指示��！已經(jīng)有6人回復(fù)
求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有9人回復(fù)
用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)
運(yùn)算過程有復(fù)數(shù)出現(xiàn)，求高版本1stOpt 已經(jīng)有6人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
1stOpt求解常微分方程邊值問題已經(jīng)有5人回復(fù)

patience-practise-passion

1樓 2014-09-10 16:06:08

已閱回復(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
wangdeqiang: 金幣+20, ★★★很有幫助 2014-09-10 20:11:54

公式：
z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2)))
-sum((kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b;

改為：
z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2)))
-sum(i=1:n)(kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b;

結(jié)果：
均方差(RMSE):0.0334465479380339
殘差平方和(SSE):0.0100680441207408
相關(guān)系數(shù)(R): 0.989298144852019
相關(guān)系數(shù)之平方(R^2): 0.978710819407646
確定系數(shù)(DC): 0.978710819407646
卡方系數(shù)(Chi-Square): 0.000543327950684942
F統(tǒng)計(F-Statistic): 140.916649515266

參數(shù) 最佳估算
---------- -------------
ki 0.245068724509857
kii -5.08642699323206
b 9.4238457540924

贊一下

回復(fù)此樓

2樓2014-09-10 20:00:19

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

wangdeqiang

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1025.8
散金: 10
帖子: 67
在線: 25.5小時
蟲號: 1931929
注冊: 2012-08-09
專業(yè): 金屬材料的力學(xué)行為

引用回帖:

2樓: Originally posted by dingd at 2014-09-10 20:00:19
公式：
z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2)))
-sum((kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b;

改為：
z=sum(i ...

非常感謝您的答復(fù)！
其實我想求kI1,kI2,kI3,KII1,KII2,KII3.
請問能否這樣編函數(shù)：
z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2)))-sum(i=1:n)((k/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b;
參數(shù)又該如何設(shè)置？
謝謝！

贊一下

回復(fù)此樓

patience-practise-passion

3樓2014-09-10 20:16:08

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

wangdeqiang

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1025.8
散金: 10
帖子: 67
在線: 25.5小時
蟲號: 1931929
注冊: 2012-08-09
專業(yè): 金屬材料的力學(xué)行為

引用回帖:

非常感謝您的答復(fù)！
其實我想求kI1,kI2,kI3,KII1,KII2,KII3.
請問能否這樣編函數(shù)：
z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2)))-sum(i=1:n)((kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b;
參數(shù)又該如何設(shè)置？
謝謝！

贊一下

回復(fù)此樓

patience-practise-passion

4樓2014-09-10 21:14:35

已閱回復(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

公式寫成如下：
Function z=sum(i=1:n, kI)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2)))
-sum(i=1:n, KII)(kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b;

均方差(RMSE):0.000319000884196416
殘差平方和(SSE):9.15854077062857E-7
相關(guān)系數(shù)(R): 0.999999031699104
相關(guān)系數(shù)之平方(R^2): 0.999998063399146
確定系數(shù)(DC): 0.999998063399146
卡方系數(shù)(Chi-Square): 4.91879921856916E-8
F統(tǒng)計(F-Statistic): 0

參數(shù) 最佳估算
---------- -------------
b 9.37438411982016
ki1 240.991105375969
ki2 45.2897205090565
ki3 -0.00759010072910402
kii1 1285.30644313843
kii2 -40.4189608273874
kii3 -4.56724263105769

贊一下

回復(fù)此樓

5樓2014-09-10 21:53:48

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

lipenggg

專家顧問 (正式寫手)

專家經(jīng)驗: +65
應(yīng)助: 122 (高中生)
金幣: 4841.3
紅花: 43
帖子: 722
在線: 435.7小時
蟲號: 3929201
注冊: 2015-06-17
性別: GG
專業(yè): 流體力學(xué)
管轄: 計算模擬

Title "DIC-K";
Parameter kI(3),kII(3),b;
constant n=3,g=8.1*10^4,u=0.26,k=2.17;
Variable x,y,z;
Function z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2)))
-sum(i=1:n)(kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b;
Data;

贊一下

回復(fù)此樓

一團(tuán)和氣

6樓2015-06-23 15:22:05

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

lipenggg

專家顧問 (正式寫手)

專家經(jīng)驗: +65
應(yīng)助: 122 (高中生)
金幣: 4841.3
紅花: 43
帖子: 722
在線: 435.7小時
蟲號: 3929201
注冊: 2015-06-17
性別: GG
專業(yè): 流體力學(xué)
管轄: 計算模擬

贊一下

回復(fù)此樓

一團(tuán)和氣

7樓2015-06-23 15:24:33

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

lipenggg

專家顧問 (正式寫手)

專家經(jīng)驗: +65
應(yīng)助: 122 (高中生)
金幣: 4841.3
紅花: 43
帖子: 722
在線: 435.7小時
蟲號: 3929201
注冊: 2015-06-17
性別: GG
專業(yè): 流體力學(xué)
管轄: 計算模擬

贊一下

回復(fù)此樓

一團(tuán)和氣

8樓2015-06-23 15:28:21

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

lipenggg

專家顧問 (正式寫手)

專家經(jīng)驗: +65
應(yīng)助: 122 (高中生)
金幣: 4841.3
紅花: 43
帖子: 722
在線: 435.7小時
蟲號: 3929201
注冊: 2015-06-17
性別: GG
專業(yè): 流體力學(xué)
管轄: 計算模擬

將公式里的KI改為KI[I]，KII改為KII后可以用1.5版運(yùn)行。

贊一下

回復(fù)此樓

一團(tuán)和氣

9樓2015-06-23 15:31:11

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

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

[考研] 化學(xué)308分調(diào)劑 +14	你好明天你好 2026-03-30	15/750	2026-04-02 13:39 by 嘻嘻笑22
[考研] 309求調(diào)劑 +7	呆菇不是戴夫 2026-04-02	7/350	2026-04-02 12:29 by 1753564080
[考研] 292求調(diào)劑 +14	是妍子也是研子 2026-03-30	14/700	2026-04-02 10:17 by sanrepian
[考研] 372求調(diào)劑 +3	jj涌77 2026-04-02	3/150	2026-04-02 09:57 by olim
[考研] 303求調(diào)劑 +8	DLkz1314. 2026-03-30	8/400	2026-04-02 09:08 by sanrepian
[考研] 336求調(diào)劑 +5	kiyy 2026-04-01	5/250	2026-04-02 07:14 by chixmc
[考研] 342求調(diào)劑 +12	Mary Keen 2026-03-28	13/650	2026-04-01 21:02 by 流情牧豪
[考研] 310分求調(diào)劑 +4	成功上岸wang 2026-04-01	4/200	2026-04-01 20:35 by liu823948201
[考研] 070300一志愿211，312分求調(diào)劑院校 +14	小黃鴨寶 2026-03-30	14/700	2026-04-01 20:19 by 賴春艷
[考研] 379求調(diào)劑 +3	?苦瓜不苦 2026-04-01	3/150	2026-04-01 20:09 by 暮云清寒
[考研] 材料專碩322分 +9	哈哈哈吼吼吼哈 2026-04-01	9/450	2026-04-01 14:30 by 無際的草原
[碩博家園] 考研調(diào)劑 +5	駱駝男人 2026-04-01	5/250	2026-04-01 14:28 by syjjj0321
[考研] 291求調(diào)劑 +3	迷蒙木木 2026-04-01	4/200	2026-04-01 11:07 by 逆水乘風(fēng)
[考研] 070300化學(xué)專業(yè)279調(diào)劑 +10	哈哈哈^_^ 2026-03-31	10/500	2026-03-31 23:13 by liu823948201
[考研] 070300化學(xué)354求調(diào)劑 +15	101次希望 2026-03-28	15/750	2026-03-31 17:58 by jp9609
[考研] 263求調(diào)劑 +3	DDDDuu 2026-03-27	3/150	2026-03-31 16:21 by 土木碩士招生
[考研] 英一數(shù)一總分334求調(diào)劑 +4	陳陽坤 2026-03-31	4/200	2026-03-31 14:22 by 記事本2026
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6	王一一依依 2026-03-30	8/400	2026-03-30 11:52 by yjolah
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[考研] 295求調(diào)劑 +5	1428151015 2026-03-27	6/300	2026-03-28 04:04 by fmesaito