




下載本文檔
版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、#include H-.88 #include j53Td. #include wwHV;i PKM8MY vo #define NUM_of_ZAXIS 400 nE3m) 2U;ImC1g float gauss_pulse(float T,float t0,float spread); IZkQmA= L: ;L/ int file_save(float* data,char* filename); m:7$oq| ,!s;o6|*y sSMcF2I WU -_Y void main() _8CE|Cn (hn+hc float exNUM_of_ZAXIS,hyNUM_of_ZAXIS
2、; O!D0 hW4 float obj_parametersNUM_of_ZAXIS4;/*模型參數(shù)設(shè)置 */ ng,qEI;Z ?X1vU0 c float dt,ddz; %8)GuxG* float source,T; 4|ryt4B int pos_driv_source; P| o_/BS int i,n,Nsteps; |_/i/H int e_low,e_high; ZqHh$QBD 9 HG2iy Vu DSjh 4+t9SD FILE* fp; #:Bd8PS 4#$gTc NGD?. (G /*初始化各個(gè)變量 */ 83dwr;- Pi=3.14159; / c1=OJ
3、Epsilon=8.85e-12; =U sl Mu=4*Pi*1e-7; NPws uFDJRQJ var_cb=dt/Epsilon; lP*_dt9 var_cp=0.5*h_sigma*dt/Mu; okfGd= & var_cq=dt/Mu; :_;9&H9ha hsIC5s3 ca=(1-var_ca)/(1+var_ca); 4.|-?qG cb=var_cb/(1+var_ca); sWw cp=(1-var_cp)/(1+var_cp); #9aB3C cq=var_cq/(1+var_cp); 07WZ w1(; 4.)hCb var_border=(C*dt-ddz)/(
4、C*dt+ddz); /i;/$ o ex_high_m1=0; :7 .a._WZF ex_low_s1=0; =ML6jr ex_low_s2=0; 4naL2 Y! ex_high_s1=0; z(14250 ex_high_s2=0; 6.? Ke8iC 8z&7wO / 上 ,下 邊界位置 S#+Dfa8X e_low=20; p%) 1(R8qM e_high=380; RRh0G* k.nq, pos_driv_source=200; FHpS?htRy ihCa $ +92/0 for(i=0;iNUM_of_ZAXIS;i+) sE87Lz pyA;%vJn ex=0; tr
5、8a_CV hy=0; pv$ o obj_parameters0=1; ETk4I obj_parameters1=1; 0) u#,#z_ /2 V printf(Nsteps-); nPlg5&E scanf(%d,&Nsteps); fc,LFvZ +5AWX,9,- for(n=0;nNsteps;n+) jL3/0o .W;,.l T+=1; d v8q&_ P*=M?:Jb, / 計(jì)算 Ex Fa |Y for(i=e_low+1;ie_high;i+) iXtL:,i x:wv#Wh:l7 / 采用葛德彪書中的一維FDTD 公式 TR7j? ex=CA*ex-CB*(hy-hy
6、i-1)/ddz; 5t=7- 7XM:4whw $ kAal26z if(T200) b8FSVV 7 obzdH:S expos_driv_source+=gauss_pulse(T,40,12); K;L6a A# $K12Wcp- Xtk3 b)(si/ TxA%0 /Mur 邊界吸收 OCN: 52P0T 7X$+S ex_low_s2=ex_low_s1; Q$)z_jai ex_low_s1=exe_low+1; 0Qr|!B:+9) 36QZV*P ex_high_m2=ex_high_m1; 1$n! Lj=5 ex_high_m1=exe_high; n2N:rP &BOG
7、&ot ex_high_s2=ex_high_s1; ,O|#?G ex_high_s1=exe_high-1; 0RF:9x2 fq(3uEnC 2 #_ i_j exe_low=ex_low_s2+var_border*(exe_low+1-ex_low_m2); HV*=Qi exe_high=ex_high_s2+var_border*(exe_high-1-ex_high_m2); -TgA*_5 MTkrv(G /行波時(shí)延法 mNQ*YCq. /exe_high=0.5*(ex_high_m2+ex_high_s2); Ps7cl /exe_low=0.5*(ex_low_m2+ex_low_s2); (iB+c +3K)9H / 計(jì)算 Hy t+z r3 ;|D8D6 for(i=e_low;ie_high;i+) 2 X5YN; e$4$G / 存儲(chǔ)數(shù)據(jù) 95-ptO1 (6 e$j fp=fopen(Ex,w); V38v2LI for(i=0;iUl0&Mf? fclose(fp); vw /X PtT$#hx fp=fopen(Hy,w); kz?m 1 for(i=0;i= Mg8 3kS M*#UV2 float gauss_pulse(float T,float t0,float spread)
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 運(yùn)動(dòng)與健康促進(jìn)工作計(jì)劃
- 管理者的影響力與說(shuō)服力計(jì)劃
- 2025年雙氰胺合作協(xié)議書
- 評(píng)優(yōu)評(píng)先活動(dòng)與實(shí)施方案計(jì)劃
- 包裝行業(yè)設(shè)計(jì)創(chuàng)新包裝材料計(jì)劃
- 2025年微循環(huán)測(cè)試系統(tǒng)項(xiàng)目建議書
- 2025年激光轉(zhuǎn)速測(cè)量?jī)x合作協(xié)議書
- 2025年中國(guó)消毒柜行業(yè)發(fā)展環(huán)境、供需態(tài)勢(shì)及投資前景分析報(bào)告(智研咨詢發(fā)布)
- 七年級(jí)下冊(cè)《平行線》課件與練習(xí)
- 電子政務(wù)系統(tǒng)權(quán)限設(shè)置指南
- 國(guó)庫(kù)集中支付培訓(xùn)班資料-國(guó)庫(kù)集中支付制度及業(yè)務(wù)操作教學(xué)課件
- 詩(shī)詞寫作入門課件
- 2023年上海青浦區(qū)區(qū)管企業(yè)統(tǒng)一招考聘用筆試題庫(kù)含答案解析
- 植物之歌觀后感
- 理發(fā)店個(gè)人門面轉(zhuǎn)讓合同
- 空氣能熱泵安裝示意圖
- 建筑工程施工質(zhì)量驗(yàn)收規(guī)范檢驗(yàn)批填寫全套表格示范填寫與說(shuō)明
- 2020年中秋國(guó)慶假日文化旅游市場(chǎng)安全生產(chǎn)檢查表
- 03J111-1 輕鋼龍骨內(nèi)隔墻
- 資產(chǎn)負(fù)債表模板范本
- 人教版高中數(shù)學(xué)選擇性必修二導(dǎo)學(xué)案
評(píng)論
0/150
提交評(píng)論