autumnfallen

銅蟲 (初入文壇)

應助: 1 (幼兒園)
金幣: 155.6
帖子: 36
在線: 6.9小時
蟲號: 3611302
注冊: 2014-12-24
性別: GG
專業(yè): 控制理論與方法

[求助] 麻煩大神幫忙用1stopt擬合得到未知參數已有1人參與

//Parameters ;
Parameters a1,b1,c1,d1,Vc=3.33e-5,B=86e-3,l=175e-3,a=43.035e-3,nc=1.30;
//Variable ;
Variable phi,Pivc,Tivc,Vivc,N,theta_ivc,theta_soc,y;
//Function ;
Function y= int(phi^b1/(a1*exp((c1*(Pivc*(Vivc/(Vc+((pi*B^2)/4)*(l+a-a*cos(theta)-sqrt(l^2-(a^2)*sin(theta)^2))))^nc)^d1)/((Tivc*(Vivc/(Vc+((pi*B^2)/4)*(l+a-a*cos(theta)-sqrt(l^2-(a^2)*sin(theta)^2)))^(nc-1)))))*N),theta=theta_ivc,theta_soc);
Data;
0.942187500000000       101.325000000000       389.827000000000       0.000533263830555081       3000       -3.14159265358979       -0.0709441434350655       1
0.942187500000000       101.325000000000       402.817000000000       0.000487559323265132       3000       -2.44346095279206       -0.0700784601260763       1
0.753750000000000       101.325000000000       377.300000000000       0.000533263830555081       3000       -3.14159265358979       -0.0715951512460594       1
0.753750000000000       101.325000000000       394.932000000000       0.000487559323265132       3000       -2.44346095279206       -0.0702303037709998       1
0.753750000000000       111.457500000000       388.831000000000       0.000533263830555081       3000       -3.14159265358979       -0.0707556478758501       1
0.753750000000000       111.457500000000       401.750000000000       0.000487559323265132       3000       -2.44346095279206       -0.0697957167872532       1
0.753750000000000       101.325000000000       389.970000000000       0.000533263830555081       3000       -3.14159265358979       -0.0699301071396568       1
0.753750000000000       101.325000000000       403.150000000000       0.000487559323265132       3000       -2.44346095279206       -0.0685914396033772       1
0.753750000000000       101.325000000000       381.380000000000       0.000487559323265132       2000       -2.44346095279206       -0.0800233462039400       1
以上是輸入和代碼，但是我的不是正版的，能夠麻煩運行一下，不能運行的話能否在幫忙修改一下，在此謝過了

公式

回復此樓

» 猜你喜歡

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

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

matlab 擬合求解參數值已經有4人回復
這樣的動力學方程能夠有1stopt擬合出參數來么已經有4人回復
看了下1stOpt的介紹，感覺略顯浮夸已經有14人回復
求助 1stopt低版本方程的擬合（含若干未知參數）已經有10人回復
1stopt自變量超多時是否有簡單語法？AND 自變量多數據量大的話1stopt都會算的特慢么已經有7人回復
求大神跑下1stOpt，一個sharemodel，謝謝！已經有7人回復
求大神給看看CFTOOL工具箱非線性參數擬合的已經有5人回復
求助 EC50 怎么算�。。�！求大神幫忙已經有3人回復
求大神幫忙擬合一個非線性方程，求出模型參數已經有15人回復
用1stOpt進行復數曲線擬合時，結果出錯。急求教�。�！已經有5人回復
一個1stOpt的初級問題：復數方程組計算問題（有代碼）已經有5人回復
1stopt能進行隱函數非線性擬合嗎？誰能幫擬合一下新型液體的粘度混合計算公式已經有5人回復
1stopt做最小二乘法回歸參數，請求大家?guī)兔χ笇?/a> 已經有3人回復

1樓 2015-03-16 09:52:01

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

autumnfallen

銅蟲 (初入文壇)

應助: 1 (幼兒園)
金幣: 155.6
帖子: 36
在線: 6.9小時
蟲號: 3611302
注冊: 2014-12-24
性別: GG
專業(yè): 控制理論與方法

引用回帖:

6樓: Originally posted by dingd at 2015-03-18 21:30:46
均方差(RMSE):0.102319169558223
殘差平方和(SSE):0.0942229121317601

參數最佳估算
-------------------- -------------
a1 6.35449569400801E-5
b1 0.167518107515883
c1 52.8945994593321
d ...

能不能麻煩您再幫我運行一下，上次的參數不知道怎么編輯錯了，再次感謝�。。�！
//Parameters ;
Parameters a1,b1,c1,d1,Vc=3.33e-5,B=86e-3,l=175e-3,a=43.035e-3,nc=1.30;
//Variable ;
Variable phi,Pivc,Tivc,Vivc,N,theta_ivc,theta_soc,y;
//Function ;
Function y= int(phi^b1/(a1*exp((c1*(Pivc*(Vivc/(Vc+((pi*B^2)/4)*(l+a-a*cos(theta)-sqrt(l^2-(a^2)*sin(theta)^2))))^nc)^d1)/((Tivc*(Vivc/(Vc+((pi*B^2)/4)*(l+a-a*cos(theta)-sqrt(l^2-(a^2)*sin(theta)^2)))^(nc-1)))))*N),theta=theta_ivc,theta_soc);
Data;
x=[0.942187500000000       101.325000000000       389.827000000000       0.000533263830555081       3000       9.42477796076938       -0.0709441434350655;
0.942187500000000       101.325000000000       402.817000000000       0.000487559323265132       3000       8.72664625997165       -0.0700784601260763;
0.753750000000000       101.325000000000       377.300000000000       0.000533263830555081       3000       9.42477796076938       -0.0715951512460594 1
0.753750000000000       101.325000000000       394.932000000000       0.000487559323265132       3000       8.72664625997165       -0.0702303037709998 1
0.753750000000000       111.457500000000       388.831000000000       0.000533263830555081       3000       9.42477796076938       -0.0707556478758501 1
0.753750000000000       111.457500000000       401.750000000000       0.000487559323265132       3000       8.72664625997165       -0.0697957167872532 1
0.753750000000000       101.325000000000       389.970000000000       0.000533263830555081       3000       9.42477796076938       -0.0699301071396568 1
0.753750000000000       101.325000000000       403.150000000000       0.000487559323265132       3000       8.72664625997165       -0.0685914396033772 1
0.753750000000000       101.325000000000       381.380000000000       0.000487559323265132       2000       8.72664625997165       -0.0800233462039400 1

贊一下

回復此樓

7樓2015-03-19 08:38:54

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

查看全部 12 個回答

feixiaolin

榮譽版主 (文壇精英)

專家經驗: +518
應助: 942 (博后)
貴賓: 1.275
金幣: 3880
散金: 58785
紅花: 532
沙發(fā): 11
帖子: 24215
在線: 2601.8小時
蟲號: 2139575
注冊: 2012-11-21
專業(yè): 光學信息獲取與處理
管轄: 數學

@dingd
dingd 請出手相助！

贊一下

回復此樓

3樓2015-03-17 21:57:05

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

dingd

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

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

【答案】應助回帖

感謝參與，應助指數 +1

因變量都為同樣的值1，視作解超越方程更合適吧。

下面一組供參考：
Parameter Best Estimate
-------------------- -------------
a1 0.00189150716596566
b1 0.183492178461146
c1 -0.0876103013145149
d1 -0.0371941222976175
vc 0.596152777032521
b 18.8251257735214
l 19.5312906804849
a 0.0460383371366824
nc -11.3790001878032

====== Output Results ======

File: Data File-1
No Target y Calculated y
1 1 1.00805909177632
2 1 0.991414168748313
3 1 1.00681046319549
4 1 0.996786505779086
5 1 0.973676198233269
6 1 0.96465215791757
7 1 1.0120386479123
8 1 1.03550477280184
9 1 1.00738120999655

贊一下

回復此樓

4樓2015-03-18 17:39:33

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

autumnfallen

銅蟲 (初入文壇)

應助: 1 (幼兒園)
金幣: 155.6
帖子: 36
在線: 6.9小時
蟲號: 3611302
注冊: 2014-12-24
性別: GG
專業(yè): 控制理論與方法

引用回帖:

4樓: Originally posted by dingd at 2015-03-18 17:39:33
因變量都為同樣的值1，視作解超越方程更合適吧。

下面一組供參考：
Parameter Best Estimate
-------------------- -------------
a1 0.00189150716596566
b1 0.183492178461146
c1 -0.087610301 ...

謝謝您的回答，能不能麻煩再幫我運行一下，因為下面這幾個參數值是固定的，但是沒有代入到方程中去，非常感謝了�。。。。�！
Vc=3.33e-5,B=86e-3,l=175e-3,a=43.035e-3,nc=1.30

贊一下

回復此樓

5樓2015-03-18 18:01:35

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

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

[考研] 321求調劑一志愿浙江工業(yè)大學生物醫(yī)藥 +4	嘿嘿HC 2026-04-01	5/250	2026-04-02 14:13 by 王亮_大連醫(yī)科大
[考研] 321求調劑 +9	y-yh 2026-04-01	10/500	2026-04-02 10:14 by 不吃魚的貓
[考研] 336求調劑 +5	kiyy 2026-04-01	5/250	2026-04-02 07:14 by chixmc
[考研] 一志愿南昌大學324求調劑 +12	hanamiko 2026-03-27	12/600	2026-04-01 22:34 by 運氣yunqi
[考研] 332求調劑 +8	Lyy930824@ 2026-03-29	8/400	2026-04-01 18:40 by 千鈞澄玉宇啊
[考研] 環(huán)境工程調劑 +9	hyzzzzzzz. 2026-04-01	9/450	2026-04-01 14:20 by salamander`
[考研] 085600，321分求調劑 +13	大饞小子 2026-03-31	13/650	2026-04-01 12:35 by chemdavid
[考研] 086000生物與醫(yī)藥298調劑求助 +4	元元青青 2026-03-31	6/300	2026-04-01 11:13 by syh9288
[考研] 08工科，295，接受跨專業(yè)調劑 +6	lmnlzy 2026-03-31	6/300	2026-04-01 11:02 by 逆水乘風
[考研] 一志愿華東師范大學有機化學專業(yè)，初試351分，復試被刷求調劑! +9	真名有冰 2026-03-29	10/500	2026-03-31 18:01 by xhai2011
[考研] 學碩274求調劑 +17	Li李魚 2026-03-26	17/850	2026-03-31 15:19 by 客爾美德
[考博] 材料專業(yè)申博 +5	杜雨婷dyt 2026-03-29	5/250	2026-03-31 11:19 by oooqiao
[考研] 274求調劑 +6	xiao愛同學 2026-03-30	6/300	2026-03-31 10:04 by cal0306
[考研] 一志愿中海洋320化學工程與技術學碩求調劑 +8	披星河 2026-03-30	8/400	2026-03-31 08:53 by lbsjt
[考研] 08工科求調劑286 +5	tgs_001 2026-03-28	5/250	2026-03-31 08:18 by 一只好果子?
[考研] 一志愿大連理工大學材料求調劑 +6	Gymno 2026-03-30	6/300	2026-03-31 07:26 by 無際的草原
[考研] 求調劑 +7	爭取九點睡 2026-03-28	8/400	2026-03-28 21:07 by 爭取九點睡
[考研] 藥學105500求調劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 一志愿吉大071010，316分求調劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw
[考研] 341求調劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid