 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1577 | 回復: 6 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
求一個用FDTD（時域有限差分法）模擬電磁波三維RCS（散射）的源程序 已有2人參與
|
|||
| 求一個用FDTD（時域有限差分法）模擬電磁波三維RCS（散射）的源程序。本人現(xiàn)在急用一個三維散射的程序，用于仿真模擬室內常用材料對電磁波的輻射和散射還有透射過程，希望有哪位好心人給個源程序，matlab或者C都行（本人不會c++，fortran等語言）。一定要是三維的，而且是分析散射的。 |
回復此樓» 猜你喜歡
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復
-
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。
已經(jīng)有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有51人回復
-
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！
已經(jīng)有1人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調劑，10人
已經(jīng)有0人回復
-
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調劑有論文
已經(jīng)有4人回復
-
歡迎調劑到贛南師范大學智能制造與未來能源學院
已經(jīng)有10人回復
-
0702一志愿吉大B區(qū)求調劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求教關于球形納米粒子局域表面等離子體共振的散射光強度分布問題
已經(jīng)有7人回復
6樓2015-11-02 21:33:42
feixiaolin
榮譽版主 (文壇精英)
-
專家經(jīng)驗: +518 - 應助: 942 (博后)
- 貴賓: 1.275
- 金幣: 3880
- 散金: 58785
- 紅花: 532
- 沙發(fā): 11
- 帖子: 24215
- 在線: 2601.8小時
- 蟲號: 2139575
- 注冊: 2012-11-21
- 專業(yè): 光學信息獲取與處理
- 管轄: 數(shù)學
2樓2013-11-02 18:04:26
 |
3樓2013-11-03 10:58:58
【答案】應助回帖
|
3-D FDTD code with PEC boundaries clear %*********************************************************************** % Fundamental constants %*********************************************************************** cc=2.99792458e8; %speed of light in free space muz=4.0*pi*1.0e-7; %permeability of free space epsz=1.0/(cc*cc*muz); %permittivity of free space %*********************************************************************** % Grid parameters %*********************************************************************** ie=50; %number of grid cells and Ex samples along x-direction je=24; %number of grid cells and Ey samples along y-direction ke=10; %number of grid cells and Ez samples along z-direction ib=ie+1; %grid boundary coordinates jb=je+1; kb=ke+1; is=26; %location of z-directed current source js=13; %location of z-directed current source kobs=5; dx=0.002; %space increment of cubic lattice dt=dx/(2.0*cc); %time step nmax=500; %total number of time steps %*********************************************************************** % Differentiated Gaussian pulse excitation %*********************************************************************** rtau=50.0e-12; tau=rtau/dt; ndelay=3*tau; J0=-1.0; %*********************************************************************** % Material parameters %*********************************************************************** eps=1.0; sig=0.0; %*********************************************************************** % Updating coefficients %*********************************************************************** ca=(1.0-(dt*sig)/(2.0*epsz*eps))/(1.0+(dt*sig)/(2.0*epsz*eps)); cb=(dt/epsz/eps/dx)/(1.0+(dt*sig)/(2.0*epsz*eps)); da=1.0; db=dt/muz/dx; %*********************************************************************** % Field arrays %*********************************************************************** ex=zeros(ie,jb,kb); ey=zeros(ib,je,kb); ez=zeros(ib,jb,ke); hx=zeros(ib,je,ke); hy=zeros(ie,jb,ke); hz=zeros(ie,je,kb); %figure %set(gcf,'DoubleBuffer','on') %*********************************************************************** % BEGIN TIME-STEPPING LOOP %*********************************************************************** for n=1:nmax %*********************************************************************** % Update electric fields %*********************************************************************** ex(1:ie,2:je,2:ke)=ca*ex(1:ie,2:je,2:ke)+... cb*(hz(1:ie,2:je,2:ke)-hz(1:ie,1:je-1,2:ke)+... hy(1:ie,2:je,1:ke-1)-hy(1:ie,2:je,2:ke)); ey(2:ie,1:je,2:ke)=ca*ey(2:ie,1:je,2:ke)+... cb*(hx(2:ie,1:je,2:ke)-hx(2:ie,1:je,1:ke-1)+... hz(1:ie-1,1:je,2:ke)-hz(2:ie,1:je,2:ke)); ez(2:ie,2:je,1:ke)=ca*ez(2:ie,2:je,1:ke)+... cb*(hx(2:ie,1:je-1,1:ke)-hx(2:ie,2:je,1:ke)+... hy(2:ie,2:je,1:ke)-hy(1:ie-1,2:je,1:ke)); ez(is,js,1:ke)=ez(is,js,1:ke)+... J0*(n-ndelay)*exp(-((n-ndelay)^2/tau^2)); %*********************************************************************** % Update magnetic fields %*********************************************************************** hx(2:ie,1:je,1:ke)=hx(2:ie,1:je,1:ke)+... db*(ey(2:ie,1:je,2:kb)-ey(2:ie,1:je,1:ke)+... ez(2:ie,1:je,1:ke)-ez(2:ie,2:jb,1:ke)); hy(1:ie,2:je,1:ke)=hy(1:ie,2:je,1:ke)+... db*(ex(1:ie,2:je,1:ke)-ex(1:ie,2:je,2:kb)+... ez(2:ib,2:je,1:ke)-ez(1:ie,2:je,1:ke)); hz(1:ie,1:je,2:ke)=hz(1:ie,1:je,2:ke)+... db*(ex(1:ie,2:jb,2:ke)-ex(1:ie,1:je,2:ke)+... ey(1:ie,1:je,2:ke)-ey(2:ib,1:je,2:ke)); %*********************************************************************** % Visualize fields %*********************************************************************** timestep=int2str(n); tview(:,  =ez(:,:,kobs);sview(:, =ez(:,js,;subplot('position',[0.15 0.45 0.7 0.45]),imagesc(tview'); shading flat; caxis([-1.0 1.0]); colorbar; axis image; axis xy; title(['Ez(i,j,k=5), time step = ',timestep]); xlabel('i coordinate'); ylabel('j coordinate'); subplot('position',[0.15 0.10 0.7 0.25]),imagesc(sview'); shading flat; caxis([-1.0 1.0]); colorbar; axis image; axis xy; title(['Ez(i,j=13,k), time step = ',timestep]); xlabel('i coordinate'); ylabel('k coordinate'); pause(0.05) %*********************************************************************** % END TIME-STEPPING LOOP %*********************************************************************** end |
4樓2014-10-07 16:11:09
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調劑 +10 | Micky11223 2026-03-25 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
|
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調劑！ +7 | Zuhui0306 2026-03-25 | 7/350 |
|
|
|
[考研] 一志愿南開大學0710生物學359求調劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
|
[考研] 334分 一志愿武理 材料求調劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
|
[考研] 324求調劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
|
[考研] 生物技術與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
|
[考研] 337求調劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
|
[考研] 275求調劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
|
[考研] 2026年華南師范大學歡迎化學，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學子加入！ +3 | llss0711 2026-03-28 | 6/300 |
|
|
|
[考研] 11408軟件工程求調劑 +3 | Qiu學ing 2026-03-28 | 3/150 |
|
|
|
[考研] 復試調劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研]
 一志愿廈門大學化學學碩307求調劑
+10
|
y7czhao 2026-03-26 | 10/500 |
|
|
|
[考研] 347求調劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 求調劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
|
[考研] 324求調劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
|
[考研] 081200-11408-276學碩求調劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研]
 085602化學工程求調劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 材料與化工304求B區(qū)調劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
|
[考研] 0854人工智能方向招收調劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|