dlionel

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3424.6
散金: 109
紅花: 4
帖子: 842
在線: 156.5小時
蟲號: 761805
注冊: 2009-05-02
專業(yè): 固體力學(xué)

[求助] 1stopt 非線性約束擬合已有1人參與

請教大家一個關(guān)于多非線性約束下擬合的問題，

自變量： theta
應(yīng)變量：V
參數(shù)： c11, c12, c13, c33, c44
常數(shù)： ro

1stopt (v1.5) 的代碼如下：

Parameter c11=53.872,c12=10,c13=15,c33=47.424,c44=19.968;
Constant ro=2080;
Function y=sqrt((c44+c11*sin(x)^2+c33*cos(x)^2+sqrt(((c11-c44)^2*sin(x)^2+(c44-c33)^2*cos(x)^2)^2+(c13+c44)^2*sin(2*x)^2))/(2*ro));
      c11>c12;
      c12>0;
      c44>0;
      (c11+c12)*c33>2*c13^2;
Data;
0                            22.8
0.382227106186758    22.2
0.911061869541040 29.2
1.22696646415201 25.0
1.54461638801498 22.8
1.57079632679490 25.9

擬合結(jié)果：
迭代數(shù): 16
計算用時(時:分:秒:微秒): 00:00:00:63
優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法
計算結(jié)束原因: 達到收斂判斷標(biāo)準(zhǔn)
均方差(RMSE): 23.7738301499779
殘差平方和(SSE): 3391.17
相關(guān)系數(shù)(R): 0
相關(guān)系數(shù)之平方(R^2): 0
決定系數(shù)(DC): -95.2443593018306
卡方系數(shù)(Chi-Square): 130.843543616615
F統(tǒng)計(F-Statistic): 0

參數(shù) 最佳估算
---------- -------------
c11         53.872
c12         10
c13         15
c33         47.424

====== 結(jié)果輸出 =====

No 實測值y 計算值y
1 22.8 1.0000000
2 22.2 1.0000000
3 29.2 1.0000000
4 25 1.0000000
5 22.8 1.0000000
6 25.9 1.0000000

參數(shù)擬合結(jié)果明顯是錯誤的，而且沒有參數(shù)c44的估算，計算值y都是1.00。
請問是多個約束條件命令不對嗎？

回復(fù)此樓

» 猜你喜歡

計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有108人回復(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)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名已經(jīng)有0人回復(fù)

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

matlab 非線性擬合已經(jīng)有16人回復(fù)
請教一個非線性方程非線性項轉(zhuǎn)化為多項式的方法已經(jīng)有31人回復(fù)
關(guān)于有約束非線性函數(shù)的優(yōu)化問題已經(jīng)有8人回復(fù)
怎么做matlab多元非線性擬合已經(jīng)有8人回復(fù)
不可線性化的非線性曲線方程擬合已經(jīng)有5人回復(fù)
求大神幫忙擬合一個非線性方程，求出模型參數(shù) 已經(jīng)有15人回復(fù)
求助1stopt含約束非線性擬合已經(jīng)有5人回復(fù)
]緊急求助復(fù)介電常數(shù)的模型擬合-Matlab or 1stOpt 已經(jīng)有9人回復(fù)
在進行多元非線性擬合的時候，用matlab擬合的值殘差都大于0 已經(jīng)有14人回復(fù)
非線性擬合得到的模型可靠性怎么驗證？已經(jīng)有8人回復(fù)
用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)
大家在使用1stOpt進行數(shù)據(jù)擬合的時候有沒有遇到過這種情況？已經(jīng)有8人回復(fù)
用1stopt或者origin非線性擬合磁導(dǎo)率已經(jīng)有5人回復(fù)
非線性數(shù)據(jù)擬合的數(shù)學(xué)原理�。。。。�？？？？已經(jīng)有30人回復(fù)
求指點 langmuir線性擬合與非線性擬合不一致已經(jīng)有12人回復(fù)
自定義非線性擬合中兩個自變量，如何擬合？已經(jīng)有4人回復(fù)

Good times, hard times, but never bad times.

1樓 2015-01-06 10:31:12

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

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

★ ★ ★ ★ ★
dlionel: 金幣+5, ★★★很有幫助 2015-01-06 13:52:09

修正亂碼：

CODE:

Constant ro=2080;

Parameters c11[10,],c12[0,],c44[0,];

Variable x,y;

Function y=sqrt((c44+c11*sin(x)^2+c33*cos(x)^2+sqrt(((c11-c44)^2*sin(x)^2+(c44-c33)^2*cos(x)^2)^2+(c13+c44)^2*sin(2*x)^2))/(2*ro));

Data;

0                              22.8

0.382227106186758     22.2

0.911061869541040        29.2

1.22696646415201        25.0

1.54461638801498        22.8

1.57079632679490        25.9

均方差(RMSE): 1.64771329728459
殘差平方和(SSE): 16.2897546602908
相關(guān)系數(shù)(R): 0.735295008053769
相關(guān)系數(shù)之平方(R^2): 0.540658748868792
決定系數(shù)(DC): 0.537682569595835
卡方系數(shù)(Chi-Square): 0.329591175354763
F統(tǒng)計(F-Statistic): 4.70812275220068

參數(shù) 最佳估算
---------- -------------
c11         1570.27237574273
c44         2.01763943470314E-9
c33         -1401.85923757196
c13         2126665.49354796

====== 結(jié)果輸出 =====

No 實測值y 計算值y
1 22.8 21.7271513
2 22.2 24.5656304
3 29.2 27.1722123
4 25 25.7695871
5 22.8 24.3644782
6 25.9 24.3537906

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

3樓2015-01-06 11:20:25

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

查看全部 7 個回答

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

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

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

CODE:

Constant ro=2080;

Parameters c11[10,],c12[0,],c44[0,];

Variable x,y;

Function y=sqrt((c44+c11*sin(x)^2+c33*cos(x)^2+sqrt(((c11-c44)^2*sin(x)^2+(c44-c33)^2*cos(x)^2)^2+(c13+c44)^2*sin(2*x)^2))/(2*ro));

Data;

0                              22.8

0.382227106186758     22.2

0.911061869541040        29.2

1.22696646415201        25.0

1.54461638801498        22.8

1.57079632679490        25.9

¾ù·½²î(RMSE): 1.64771329728459
²Ð²îÆ½·½ºÍ(SSE): 16.2897546602908
Ïà¹ØÏµÊý(R): 0.735295008053769
Ïà¹ØÏµÊýÖ®Æ½·½(R^2): 0.540658748868792
¾ö¶¨ÏµÊý(DC): 0.537682569595835
¿¨·½ÏµÊý(Chi-Square): 0.329591175354763
FÍ³¼Æ(F-Statistic): 4.70812275220068

²ÎÊý ×î¼Ñ¹ÀËã
---------- -------------
c11         1570.27237574273
c44         2.01763943470314E-9
c33         -1401.85923757196
c13         2126665.49354796

====== ½á¹ûÊä³ö =====

No Êµ²âÖµy ¼ÆËãÖµy
1 22.8 21.7271513
2 22.2 24.5656304
3 29.2 27.1722123
4 25 25.7695871
5 22.8 24.3644782
6 25.9 24.3537906

建議你看看該軟件的指導(dǎo)手冊：http://www.gaoyang168.com/bbs/viewthread.php?tid=4586018&fpage=1&target=blank

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2015-01-06 11:19:20

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

dlionel

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3424.6
散金: 109
紅花: 4
帖子: 842
在線: 156.5小時
蟲號: 761805
注冊: 2009-05-02
專業(yè): 固體力學(xué)

引用回帖:

3樓: Originally posted by 月只藍 at 2015-01-06 11:20:25
修正亂碼：
Constant ro=2080;
Parameters c11,c12,c44;
Variable x,y;
Function y=sqrt((c44+c11*sin(x)^2+c33*cos(x)^2+sqrt(((c11-c44)^2*sin(x)^2+(c44-c33)^2*cos(x)^2)^2+(c13+c44)^2*sin(2*x)^2))/(2*ro ...

多謝版主，手冊確實很有幫助。

請問約束 (c11+c12)*c33>2*c13^2 怎么沒有在你的代碼里給出啊？應(yīng)該怎么加上去呢？

根據(jù)你的回復(fù)，我又修改了一下原問題：c12是不需要擬合的，擬合函數(shù)形式不變，
自變量：theta
應(yīng)變量：V
擬合參數(shù)：c11, c13, c33, c44
已知常數(shù)：c12,ro
約束條件：c11>c12;
            c33>0;
            c44>0;
            (c11+c12)*c33>2*c13^2;

贊一下

回復(fù)此樓

Good times, hard times, but never bad times.

4樓2015-01-06 13:59:39

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

dlionel

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3424.6
散金: 109
紅花: 4
帖子: 842
在線: 156.5小時
蟲號: 761805
注冊: 2009-05-02
專業(yè): 固體力學(xué)

引用回帖:

4樓: Originally posted by dlionel at 2015-01-06 13:59:39
多謝版主，手冊確實很有幫助。

請問約束 (c11+c12)*c33>2*c13^2 怎么沒有在你的代碼里給出�。繎�(yīng)該怎么加上去呢？

根據(jù)你的回復(fù)，我又修改了一下原問題：c12是不需要擬合的，擬合函數(shù)形式不變，
自變量： ...

且已知常數(shù)
ro = 2080
c12 = 850000

參數(shù)initial guess
c11 = 1384531.2
c13 = 8000000
c33 = 1081267.2
c44 = 191692.8

贊一下

回復(fù)此樓

Good times, hard times, but never bad times.

5樓2015-01-06 14:08:04

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

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

[考研] 266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑 +4	哇呼哼呼哼 2026-04-02	4/200	2026-04-02 13:10 by yulian1987
[考研] 285求調(diào)劑 +11	AZMK 2026-04-02	11/550	2026-04-02 12:07 by 1753564080
[考研] 0856初試324分求調(diào)劑 +6	想上學(xué)求調(diào) 2026-04-01	6/300	2026-04-02 11:42 by 星空星月
[考研] 08開頭看過來！�。� +4	wwwwffffff 2026-03-31	6/300	2026-04-02 11:42 by 均值回歸
[考研] 085601 材料工程 313分求調(diào)劑 +7	Ong3 2026-03-27	7/350	2026-04-02 11:29 by 3041
[考研] 385分生物學(xué)（071000）求調(diào)劑 +6	qf626 2026-04-01	6/300	2026-04-02 10:39 by guoweigw
[考研] 337求調(diào)劑 +11	《樹》 2026-03-29	11/550	2026-04-02 10:20 by 不吃魚的貓
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +14	傾____蕭 2026-03-27	15/750	2026-04-02 09:21 by olim
[考研] 318求調(diào)劑，計算材料方向 +7	吸喵有害笙命 2026-04-01	8/400	2026-04-01 23:54 by 義文wang
[考研] 304求調(diào)劑 +12	素年祭語 2026-03-31	15/750	2026-04-01 22:41 by peike
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +10	GENJIOW 2026-03-31	11/550	2026-04-01 21:38 by GENJIOW
[考研] 284求調(diào)劑 +12	小熊～～ 2026-03-31	12/600	2026-04-01 20:23 by 花??
[考研] 材料與化工調(diào)劑一志愿大連海事085600，349 +9	吃的不少 2026-03-30	9/450	2026-04-01 11:24 by wangjy2002
[考研] 材料調(diào)劑 +10	Eujd1 2026-03-31	11/550	2026-04-01 11:23 by ivanqyq
[考研] 086000生物與醫(yī)藥298調(diào)劑求助 +4	元元青青 2026-03-31	6/300	2026-04-01 11:13 by syh9288
[考研] 254材料與化工求調(diào)劑 +3	翰冬林楠 2026-03-30	4/200	2026-03-31 17:53 by yishunmin
[考研] 086000調(diào)劑 +5	7901117076 2026-03-26	5/250	2026-03-31 17:45 by 544594351
[考研] 315求調(diào)劑 +6	akie... 2026-03-28	7/350	2026-03-31 16:48 by asdfzly
[考研] 304材料求調(diào)劑 +5	鐘llll 2026-03-26	5/250	2026-03-30 19:37 by 源_2020
[考研] 105500藥學(xué)求調(diào)劑，一志愿山東大學(xué)藥學(xué)，348分 +3	gr哈哈哈 2026-03-28	3/150	2026-03-30 18:56 by 源_2020