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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者ckm0811將贈送您 10 個金幣

ckm0811

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 30.5
帖子: 29
在線: 27.6小時
蟲號: 16121443
注冊: 2019-08-22
專業(yè): 結(jié)構(gòu)工程

送紅花一朵

引用回帖:

10樓: Originally posted by wlfc at 2021-08-11 16:09:33
使用 3樓、4樓的公式進(jìn)行求解。

有兩個難點(diǎn)：
1、微分方程求解時要傳遞中間變量lamb, pa。
2、lambv為未知中間變量，需要求解。
第一個難點(diǎn)對Lu腳本來說不是問題，gsl支持的Lu擴(kuò)展數(shù)學(xué)庫中求解微分方程的函數(shù) ...

非常感謝，您博客里給出的圖形擬合效果很好，我再仔細(xì)看一下您的答復(fù)

贊一下

回復(fù)此樓

11樓2021-08-16 18:10:56

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

wlfc

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 77
帖子: 19
在線: 12小時
蟲號: 26920802
注冊: 2021-08-01

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

★ ★ ★ ★ ★
獨(dú)孤神宇: 金幣+5, 鼓勵交流 2021-08-22 08:43:45

引用回帖:

7樓: Originally posted by wlfc at 2021-08-09 08:09:41
嘗試用OpenLu（可從www.forcal.net下載）求解，與大家探討。
使用 1樓給出的第一種公式。未知中間變量lambv用解方程法求解，但不知這種方法有沒有問題？
仍類似于 4樓使用部分?jǐn)?shù)據(jù)擬合，以減少運(yùn)行時間。

Lu ...

7樓、8樓、9樓的回復(fù)，因理解有誤，故重新解答。

仍使用 1樓給出的第一種公式。仍使用部分?jǐn)?shù)據(jù)擬合。

分析：需在微分方程求解中傳遞中間變量lamb，gsl支持的Lu擴(kuò)展數(shù)學(xué)庫中求解微分方程的函數(shù)gsl_ode提供了傳遞中間變量的功能。
進(jìn)一步討論：微分方程求解傳遞中間變量時，按時間t的改變量進(jìn)行插值，可提高求解精度。例如，進(jìn)行線性插值：lamb=lamb0+(lamb1-lamb0)*[(t-t0)/(t1-t0)]。

另外，微分方程初值lambv未知，將其追加為擬合變量lambv0。

Lu代碼：

CODE:

!!!using["luopt","math"]; //使用命名空間

f(  t,lambv,dlambv,i

    : lamb

    : t_A, t_lamb, fac1,fac2,fac3,fac4,fac5,fac6,fac7) =

    lamb=t_lamb[i]+(t_lamb[i+1]-t_lamb[i])*[(t-t_A[i])/(t_A[i+1]-t_A[i])],

    dlambv=((1.0/3/fac7)*(fac1*((lamb/lambv)^fac2-(lambv/lamb)^(0.5*fac2))+fac3*((lamb/lambv)^fac4-(lambv/lamb)^(0.5*fac4))+fac5*((lamb/lambv)^fac6-(lambv/lamb)^(0.5*fac6)))),

    0;  //成功返回0

目標(biāo)函數(shù)(_fac1,_fac2,_fac3,_fac4,_fac5,_fac6,_fac7,lambv0

    : i,s,tf,pa,lamb,lambv

    : t_A, t_lamb, t_pa, t_p, max, fac1,fac2,fac3,fac4,fac5,fac6,fac7)=

{

    fac1=_fac1, fac2=_fac2, fac3=_fac3, fac4=_fac4, fac5=_fac5, fac6=_fac6, fac7=_fac7,  //傳遞優(yōu)化變量

    //最后一個參數(shù)50表示gsl_ode函數(shù)在計算時,最多循環(huán)計算50次,這樣可以提高速度

    tf=gsl_ode[@f, nil, nil, t_A, ra1(lambv0),  1e-6, 1e-6, gsl_rkf45, 1e-6,50],

    i=-1, s=0, while{++i<max,

        pa=t_pa(i), lamb=t_lamb(i), lambv=tf(i,1),

        s=s+[(pa+fac1*((lamb/lambv)^fac2/lamb-(lambv/lamb)^(fac2*0.5)/lamb)+fac3*((lamb/lambv)^fac4/lamb-(lambv/lamb)^(fac4*0.5)/lamb)+fac5*((lamb/lambv)^fac6/lamb-(lambv/lamb)^(fac6*0.5)/lamb)) - t_p(i)]^2.0

    },

    s

};

main(: tArray : t_A, t_lamb, t_pa, t_p, max)=

{

    tArray=matrix{ //存放實驗數(shù)據(jù)//t,lamb,pa,p

        "0        0.959926667        -0.218715994        -0.060966222

0.25        0.963063333        -0.200840559        -0.046836278

0.5        0.96612        -0.183553472        -0.035042256

0.75        0.969016667        -0.167289626        -0.024438988

1        0.97171        -0.152268899        -0.01706765

1.25        0.97409        -0.139075651        -0.010939712

1.5        0.976123333        -0.127862697        -0.006795087

1.75        0.977786667        -0.118729759        -0.002919842

2        0.979000001        -0.112089924        -0.001981819

2.25        0.979743333        -0.108031322        -0.001163676

2.5        0.98001        -0.106577018        -0.003491456

2.75        0.979766667        -0.107904033        -0.009723697"

    },

    len[tArray,0,&max], t_A=tArray(all:0).reshape(), t_lamb=tArray(all:1).reshape(), t_pa=tArray(all:2).reshape(), t_p=tArray(all:3).reshape(), //用len函數(shù)取矩陣的行數(shù)，t_A等取矩陣的列

    Opt1[@目標(biāo)函數(shù),  optwaysimdeep, optwayconfra] //Opt1函數(shù)全局優(yōu)化

};

結(jié)果：

7.148976005697928e-003 -14.5708961980305 -1.059296609366114 0.2342635320112431 -9.68049051324744e-005 -31.79361445870142 -2.026215042289879 1.184269268536115 3.55377979903292e-006

繪圖代碼：

CODE:

!!!using["luopt","math","win"]; //使用命名空間

f(  t,lambv,dlambv,i

    : lamb

    : t_A, t_lamb, fac1,fac2,fac3,fac4,fac5,fac6,fac7) =

    lamb=t_lamb[i]+(t_lamb[i+1]-t_lamb[i])*[(t-t_A[i])/(t_A[i+1]-t_A[i])],

    dlambv=((1.0/3/fac7)*(fac1*((lamb/lambv)^fac2-(lambv/lamb)^(0.5*fac2))+fac3*((lamb/lambv)^fac4-(lambv/lamb)^(0.5*fac4))+fac5*((lamb/lambv)^fac6-(lambv/lamb)^(0.5*fac6)))),

    0;  //成功返回0

get(t_pp,t_lambv, _fac1,_fac2,_fac3,_fac4,_fac5,_fac6,_fac7,lambv0

    : i,tf,pa,lamb,lambv

    : t_A, t_lamb, t_pa, t_p, max, fac1,fac2,fac3,fac4,fac5,fac6,fac7)=

{

    fac1=_fac1, fac2=_fac2, fac3=_fac3, fac4=_fac4, fac5=_fac5, fac6=_fac6, fac7=_fac7,  //傳遞優(yōu)化變量

    //最后一個參數(shù)50表示gsl_ode函數(shù)在計算時,最多循環(huán)計算50次,這樣可以提高速度

    tf=gsl_ode[@f, nil, nil, t_A, ra1(lambv0),  1e-6, 1e-6, gsl_rkf45, 1e-6,50],

    t_lambv.copy[tf(all,1)],

    i=-1, while{++i<max,

        pa=t_pa(i), lamb=t_lamb(i), lambv=t_lambv(i),

        t_pp[i]=(pa+fac1*((lamb/lambv)^fac2/lamb-(lambv/lamb)^(fac2*0.5)/lamb)+fac3*((lamb/lambv)^fac4/lamb-(lambv/lamb)^(fac4*0.5)/lamb)+fac5*((lamb/lambv)^fac6/lamb-(lambv/lamb)^(fac6*0.5)/lamb))

    }

};

init( main

    :  tArray, tf, i, k, t_pp,t_lambv

    : t_A, t_lamb, t_pa, t_p, max)=

{

    tArray=matrix{ //存放實驗數(shù)據(jù)//t,lamb,pa,p

        "0        0.959926667        -0.218715994        -0.060966222

0.25        0.963063333        -0.200840559        -0.046836278

0.5        0.96612        -0.183553472        -0.035042256

0.75        0.969016667        -0.167289626        -0.024438988

1        0.97171        -0.152268899        -0.01706765

1.25        0.97409        -0.139075651        -0.010939712

1.5        0.976123333        -0.127862697        -0.006795087

1.75        0.977786667        -0.118729759        -0.002919842

2        0.979000001        -0.112089924        -0.001981819

2.25        0.979743333        -0.108031322        -0.001163676

2.5        0.98001        -0.106577018        -0.003491456

2.75        0.979766667        -0.107904033        -0.009723697"

    },

    len[tArray,0,&max], t_A=tArray(all:0).reshape(), t_lamb=tArray(all:1).reshape(), t_pa=tArray(all:2).reshape(), t_p=tArray(all:3).reshape(), //用len函數(shù)取矩陣的行數(shù)，t_A等取矩陣的列

    t_pp=new[real_s,max], t_lambv=new[real_s,max],

    get[t_pp, t_lambv, 7.148976005697928e-003  ,  -14.5708961980305     ,    -1.059296609366114    ,    0.2342635320112431    ,    -9.68049051324744e-005  ,  -31.79361445870142     ,   -2.026215042289879     ,   1.184269268536115  ],

    cwAttach[typeSplit], cwResizePlots(1,2,2), //左二右二分裂

    k=cwAddCurve{t_A, t_p, max, 0},    //給0子圖添加曲線

    cwSetScatter(k,0),                 //設(shè)置繪制點(diǎn)

    cwSetDataLineSize(5, k, 0),        //設(shè)置點(diǎn)的大小

    cwAddCurve{t_A, t_pp, max, 0},     //給0子圖添加曲線

    cwAddCurve{t_A, t_lamb, max, 1},   //給1子圖添加曲線

    cwAddCurve{t_A, t_pa, max, 2},     //給2子圖添加曲線

    cwAddCurve{t_A, t_lambv, max, 3}   //給3子圖添加曲線

};

ChartWnd[@init];

圖形參考：https://blog.csdn.net/wlfc/article/details/119830682

贊一下

回復(fù)此樓

12樓2021-08-20 20:33:59

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

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

返回列表

不應(yīng)助 確定回帖應(yīng)助 (注意：應(yīng)助才可能被獎勵，但不允許灌水，必須填寫15個字符以上)

普通表情龍兔虎貓

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

[考研] 一志愿314求調(diào)劑 +6	202111120625 2026-03-03	6/300	2026-03-04 14:05 by kakakapanpan
[考研] 0857求調(diào)劑 +3	CDDD1013 2026-03-04	3/150	2026-03-04 11:26 by lature00
[考研] 0856材料與化工，270求調(diào)劑 +12	YXCT 2026-03-01	14/700	2026-03-04 10:40 by winsuccess
[考博] 26申博 +5	north， 2026-02-28	5/250	2026-03-04 10:01 by StarAura
[考研] 085701環(huán)境工程求調(diào)劑 +5	xiiiia 2026-03-04	5/250	2026-03-04 09:44 by houyaoxu
[考研] 289求調(diào)劑 +8	yang婷 2026-03-02	10/500	2026-03-04 09:10 by 無際的草原
[考研] 接收調(diào)劑 +15	津萌津萌 2026-03-02	23/1150	2026-03-03 23:39 by 夢—-
[考研] 085602化學(xué)工程350，調(diào)劑，有沒有211的 +5	利好利好. 2026-03-02	9/450	2026-03-03 17:06 by 利好利好.
[考研] 283求調(diào)劑 +9	鹿沫笙 2026-03-02	9/450	2026-03-03 15:58 by tgxtgxtgx9
[考研] 材料學(xué)碩318求調(diào)劑 +11	February_Feb 2026-03-01	11/550	2026-03-03 14:24 by NUAAZXWS
[考研] 284求調(diào)劑 +6	天下熯 2026-03-02	6/300	2026-03-03 13:46 by Iveryant
[考研] 環(huán)境工程專碩307 一志愿211 四六級已過求調(diào)劑 +3	ccc! 2026-03-03	3/150	2026-03-03 12:55 by 熱情沙漠
[考研] 299求調(diào)劑 +5	kkcoco25 2026-03-02	9/450	2026-03-03 12:55 by 公瑾逍遙
[考研] 調(diào)劑材料學(xué)碩 +4	詞凝Y 2026-03-02	4/200	2026-03-03 07:53 by njzyff
[考研] 材料復(fù)試調(diào)劑 +5	學(xué)材料的點(diǎn) 2026-03-01	6/300	2026-03-02 20:01 by hypershenger
[考研] 化學(xué)，材料，環(huán)境類求調(diào)劑 +7	考研版棒棒 2026-03-02	7/350	2026-03-02 19:56 by hypershenger
[考研] 292求調(diào)劑 +7	yhk_819 2026-02-28	7/350	2026-03-02 12:43 by 無際的草原
[考研] 化工專碩342，一志愿大連理工大學(xué)，求調(diào)劑 +6	kyf化工 2026-02-28	7/350	2026-03-02 10:56 by 無際的草原
[基金申請] 本子寫完了，給DS兄弟看了，得了92分 +3	Doma 2026-03-01	7/350	2026-03-02 00:00 by jnzsy
[考研] 295復(fù)試調(diào)劑 +3	簡木ChuFront 2026-03-01	3/150	2026-03-01 14:27 by zzxw520th

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

24小時熱門版塊排行榜

ckm0811

wlfc

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

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