 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 1575 | 回復: 6 | |||
[求助]
求一個用FDTD（時域有限差分法）模擬電磁波三維RCS（散射）的源程序 已有2人參與
|
| 求一個用FDTD（時域有限差分法）模擬電磁波三維RCS（散射）的源程序。本人現(xiàn)在急用一個三維散射的程序，用于仿真模擬室內常用材料對電磁波的輻射和散射還有透射過程，希望有哪位好心人給個源程序，matlab或者C都行（本人不會c++，fortran等語言）。一定要是三維的，而且是分析散射的。 |
回復此樓» 猜你喜歡
調劑
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有63人回復
-
中德博士后2026項目5月30日截止- 國內博士后入站、德國兩年全職
已經(jīng)有0人回復
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內博后入站德國兩年全職
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復
-
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。
已經(jīng)有0人回復
-
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！
已經(jīng)有1人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調劑，10人
已經(jīng)有0人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求教關于球形納米粒子局域表面等離子體共振的散射光強度分布問題
已經(jīng)有7人回復
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
5樓2015-10-21 15:27:14
6樓2015-11-02 21:33:42
7樓2015-11-03 09:36:37
相關版塊跳轉
我要訂閱樓主 ☆皓月當空☆ 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 343求調劑 +4 | 愛羈絆 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 312，生物學求調劑 +3 | 小譯同學abc 2026-03-28 | 3/150 |
|
|
|
[考研]
 一志愿廈門大學化學學碩307求調劑
+10
|
y7czhao 2026-03-26 | 10/500 |
|
|
|
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調劑！ +4 | Zuhui0306 2026-03-25 | 4/200 |
|
|
|
[考研] 083000學碩274求調劑 +7 | Li李魚 2026-03-26 | 7/350 |
|
|
|
[考研] 291求調劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
|
[考研] 322求調劑 +6 | 宋明欣 2026-03-27 | 6/300 |
|
|
|
[考研] 315分求調劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 275求調劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 286求調劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[考研] 材料調劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
|
[考研] 081700 調劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
|
[考研] 一志愿天津大學339材料與化工求調劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
|
[考研] 考研一志愿蘇州大學初始315（英一）求調劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
|
[考研] 0854人工智能方向招收調劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
|
[考研] 285求調劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] 277分求調劑，跨調材料 +3 | 考研調劑lxh 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿國科過程所081700，274求調劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
|
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|