| 3 | 1/1 | 返回列表 |
| 查看: 2213 | 回復: 2 | |||
| 【懸賞金幣】回答本帖問題,作者束文生將贈送您 25 個金幣 | |||
[求助]
求助介電常數(shù)(復數(shù))怎么用一階德拜方程擬合 已有1人參與
|
|||
|
介電常數(shù)在每個頻點的數(shù)據(jù)為復數(shù),現(xiàn)在知道部分頻點的介電常數(shù),需要用一階德拜方程來擬合,如何得到德拜方程的系數(shù)呢?用matlab還是origin? 一階德拜方程: y=(a+(b-a)/(1+(2*pi*f*c)^2))-j*(d/(2*pi*f*(1/(36*pi)*10^-9))+((b-a)*2*pi*f*c)/(1+(2*pi*f*c)^2));其中f為頻率,a b c d 為需要擬合的系數(shù); 頻率: 3000000 4000000 5000000 6000000 7000000 8000000 9000000 10000000 15000000 20000000 25000000 30000000 40000000 50000000 60000000 70000000 80000000 90000000 100000000 介電常數(shù)實部: 13.5 10.61 8.8 7.58 6.7 6.04 5.53 5.12 3.92 3.35 3 2.78 2.52 2.37 2.26 2.18 2.11 2.048 2.01 介電常數(shù)虛部: 17.85 15.28 13.35 11.86 10.67 9.72 8.93 8.27 6.07 4.83 4.02 3.46 2.72 2.25 1.92 1.67 1.47 1.3 1.17 |
|
用OpenLu(可從www.forcal.net下載)求解。 !!!using["luopt","math","win"]; //使用命名空間 g(a,b,c,d,f)= (a+(b-a)/(1+(2*pi*f*c)^2))-1.0i*(d/(2*pi*f*(1/(36*pi)*1e-9))+((b-a)*2*pi*f*c)/(1+(2*pi*f*c)^2)); 目標函數(shù)(a,b,c,d : i,s, e11,e22 : max, f, e1, e2)= { i=-1, s=0, while{++i<max, toreal[g(a,b,c,d,f(i,0)), &e11,&e22], s=s+[e1(i,0)-e11]^2+[e2(i,0)-e22]^2 }, s }; main(: tArray : max, f, e1, e2)= { tArray=matrix{ //存放實驗數(shù)據(jù) //f e1[realPart] e2[imagPart] "3000000 13.5 17.85 4000000 10.61 15.28 5000000 8.8 13.35 6000000 7.58 11.86 7000000 6.7 10.67 8000000 6.04 9.72 9000000 5.53 8.93 10000000 5.12 8.27 15000000 3.92 6.07 20000000 3.35 4.83 25000000 3 4.02 30000000 2.78 3.46 40000000 2.52 2.72 50000000 2.37 2.25 60000000 2.26 1.92 70000000 2.18 1.67 80000000 2.11 1.47 90000000 2.048 1.3 100000000 2.01 1.17 " }, len[tArray,0,&max], f=tArray(all:0), e1=tArray(all:1),e2=tArray(all:2), //用len函數(shù)取矩陣的行數(shù),f等取矩陣的列 Opt1[@目標函數(shù)] //Opt1函數(shù)全局優(yōu)化 }; 結(jié)果(a,b,c,d,最小值): 2.393920870707986 18.24067405097755 -3.758382749525807e-008 -1.610640545623026e-003 3.654009362095743 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +5 | 手機用戶 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 344求調(diào)劑 +5 | knight344 2026-03-16 | 5/250 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 281求調(diào)劑(0805) +4 | 煙汐憶海 2026-03-16 | 11/550 |
|
|
[考研] 0703化學求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
[考研] 求調(diào)劑,總分315,考的生物醫(yī)藥,一志愿湖南師范大學。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進步 2026-03-11 | 5/250 |
|
|
[考研] 材料專碩274一志愿陜西師范大學求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 0703化學調(diào)劑,求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|