 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1280 | 回復(fù): 11 | |||
autumnfallen銅蟲 (初入文壇)
|
[求助]
麻煩大神幫忙用1stopt擬合得到未知參數(shù) 已有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 以上是輸入和代碼，但是我的不是正版的，能夠麻煩運(yùn)行一下，不能運(yùn)行的話能否在幫忙修改一下，在此謝過了 公式 |
回復(fù)此樓» 猜你喜歡
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有32人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有225人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有76人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有1人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab 擬合求解參數(shù)值
已經(jīng)有4人回復(fù)
- 這樣的動(dòng)力學(xué)方程能夠有1stopt擬合出參數(shù)來么
已經(jīng)有4人回復(fù)
- 看了下1stOpt的介紹，感覺略顯浮夸
已經(jīng)有14人回復(fù)
- 求助 1stopt低版本 方程的擬合（含若干未知參數(shù)）
已經(jīng)有10人回復(fù)
- 1stopt自變量超多時(shí)是否有簡(jiǎn)單語(yǔ)法？AND 自變量多 數(shù)據(jù)量大的話1stopt都會(huì)算的特慢么
已經(jīng)有7人回復(fù)
- 求大神跑下1stOpt，一個(gè)sharemodel，謝謝！
已經(jīng)有7人回復(fù)
- 求大神給看看CFTOOL工具箱非線性參數(shù)擬合的
已經(jīng)有5人回復(fù)
- 求助 EC50 怎么算�。。。� 求大神幫忙
已經(jīng)有3人回復(fù)
- 求大神幫忙擬合一個(gè)非線性方程，求出模型參數(shù)
已經(jīng)有15人回復(fù)
- 用1stOpt進(jìn)行復(fù)數(shù)曲線擬合時(shí)，結(jié)果出錯(cuò)。急求教！��！
已經(jīng)有5人回復(fù)
- 一個(gè)1stOpt的初級(jí)問題：復(fù)數(shù)方程組計(jì)算問題（有代碼）
已經(jīng)有5人回復(fù)
- 1stopt能進(jìn)行隱函數(shù)非線性擬合嗎？誰(shuí)能幫擬合一下新型液體的粘度混合計(jì)算公式
已經(jīng)有5人回復(fù)
- 1stopt做最小二乘法回歸參數(shù)，請(qǐng)求大家?guī)兔χ笇?dǎo)下
已經(jīng)有3人回復(fù)
- 用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大
已經(jīng)有5人回復(fù)
- 大家在使用1stOpt進(jìn)行數(shù)據(jù)擬合的時(shí)候有沒有遇到過這種情況？
已經(jīng)有8人回復(fù)
- 用origin或1stopt擬合一條隱函數(shù)方程曲線
已經(jīng)有4人回復(fù)
autumnfallen
銅蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 155.6
- 帖子: 36
- 在線: 6.9小時(shí)
- 蟲號(hào): 3611302
- 注冊(cè): 2014-12-24
- 性別: GG
- 專業(yè): 控制理論與方法
2樓2015-03-17 11:36:46
feixiaolin
榮譽(yù)版主 (文壇精英)
-
專家經(jīng)驗(yàn): +518 - 應(yīng)助: 942 (博后)
- 貴賓: 1.275
- 金幣: 3430
- 散金: 58785
- 紅花: 532
- 沙發(fā): 11
- 帖子: 24215
- 在線: 2601.8小時(shí)
- 蟲號(hào): 2139575
- 注冊(cè): 2012-11-21
- 專業(yè): 光學(xué)信息獲取與處理
- 管轄: 數(shù)學(xué)
3樓2015-03-17 21:57:05
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +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
autumnfallen
銅蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 155.6
- 帖子: 36
- 在線: 6.9小時(shí)
- 蟲號(hào): 3611302
- 注冊(cè): 2014-12-24
- 性別: GG
- 專業(yè): 控制理論與方法
5樓2015-03-18 18:01:35
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
autumnfallen: 金幣+5, ★有幫助 2015-03-18 21:58:28
autumnfallen: 金幣+5, ★有幫助 2015-03-18 21:58:28
|
均方差(RMSE):0.102319169558223 殘差平方和(SSE):0.0942229121317601 參數(shù) 最佳估算 -------------------- ------------- a1 6.35449569400801E-5 b1 0.167518107515883 c1 52.8945994593321 d1 -0.304188153936499 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No 目標(biāo) y 計(jì)算 y 1 1 1.0746489425274 2 1 0.912041978907172 3 1 0.977075337840864 4 1 0.847834098373018 5 1 1.083406434405 6 1 0.920870246098658 7 1 1.03781299546324 8 1 0.882465450911845 9 1 1.16962160291593 |
6樓2015-03-18 21:30:46
autumnfallen
銅蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 155.6
- 帖子: 36
- 在線: 6.9小時(shí)
- 蟲號(hào): 3611302
- 注冊(cè): 2014-12-24
- 性別: GG
- 專業(yè): 控制理論與方法
|
能不能麻煩您再幫我運(yùn)行一下，上次的參數(shù)不知道怎么編輯錯(cuò)了，再次感謝�。。�！ //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
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
autumnfallen8樓2015-03-19 11:35:47
autumnfallen
銅蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 155.6
- 帖子: 36
- 在線: 6.9小時(shí)
- 蟲號(hào): 3611302
- 注冊(cè): 2014-12-24
- 性別: GG
- 專業(yè): 控制理論與方法
9樓2015-03-19 14:27:06
autumnfallen
銅蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 155.6
- 帖子: 36
- 在線: 6.9小時(shí)
- 蟲號(hào): 3611302
- 注冊(cè): 2014-12-24
- 性別: GG
- 專業(yè): 控制理論與方法
10樓2015-03-20 14:26:53
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 autumnfallen 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
