☆皓月當(dāng)空☆

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 73.7
散金: 10
紅花: 1
帖子: 38
在線: 6.4小時(shí)
蟲號: 2173071
注冊: 2012-12-04
專業(yè): 激光

[求助] 求一個(gè)用FDTD（時(shí)域有限差分法）模擬電磁波三維RCS（散射）的源程序已有2人參與

求一個(gè)用FDTD（時(shí)域有限差分法）模擬電磁波三維RCS（散射）的源程序。本人現(xiàn)在急用一個(gè)三維散射的程序，用于仿真模擬室內(nèi)常用材料對電磁波的輻射和散射還有透射過程，希望有哪位好心人給個(gè)源程序，matlab或者C都行（本人不會c++，fortran等語言）。一定要是三維的，而且是分析散射的。

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有52人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

求教關(guān)于球形納米粒子局域表面等離子體共振的散射光強(qiáng)度分布問題已經(jīng)有7人回復(fù)

1樓 2013-11-02 16:49:58

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

feixiaolin

榮譽(yù)版主 (文壇精英)

專家經(jīng)驗(yàn): +518
應(yīng)助: 942 (博后)
貴賓: 1.275
金幣: 3430
散金: 58785
紅花: 532
沙發(fā): 11
帖子: 24215
在線: 2601.8小時(shí)
蟲號: 2139575
注冊: 2012-11-21
專業(yè): 光學(xué)信息獲取與處理
管轄: 數(shù)學(xué)

百度: FDTD 時(shí)域有限差分法源程序

贊一下

回復(fù)此樓

2樓2013-11-02 18:04:26

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

☆皓月當(dāng)空☆

新蟲 (初入文壇)

引用回帖:

2樓: Originally posted by feixiaolin at 2013-11-02 18:04:26
百度: FDTD 時(shí)域有限差分法源程序

3樓2013-11-03 10:58:58

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

sumin84

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5
帖子: 1
在線: 3小時(shí)
蟲號: 3398582
注冊: 2014-09-04

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

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

贊一下

回復(fù)此樓

4樓2014-10-07 16:11:09

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

雪影舞動

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 229.5
帖子: 6
在線: 1.4小時(shí)
蟲號: 4139955
注冊: 2015-10-13
性別: MM
專業(yè): 光學(xué)

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

請問你也是做RCS么
我看到網(wǎng)上有一個(gè)C++的源程序就是RCS的
http://www.codeforge.cn/article/214605
我也是做RCS的需要仿真完全不理解可不可以交流一下？扣扣：397099837

贊一下

回復(fù)此樓

你要不斷努力，以保持原地奔跑

5樓2015-10-21 15:27:14

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

1099797132

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 16.5
帖子: 10
在線: 3.6小時(shí)
蟲號: 3490643
注冊: 2014-10-22

matlab,FDTD的散射交流，QQ1099797132

贊一下

回復(fù)此樓

6樓2015-11-02 21:33:42

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

雪影舞動

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 229.5
帖子: 6
在線: 1.4小時(shí)
蟲號: 4139955
注冊: 2015-10-13
性別: MM
專業(yè): 光學(xué)

引用回帖:

6樓: Originally posted by 1099797132 at 2015-11-02 21:33:42
matlab,FDTD的散射交流，QQ1099797132

找不到這個(gè)群呢

回復(fù)此樓

你要不斷努力，以保持原地奔跑

7樓2015-11-03 09:36:37

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

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

[基金申請] 山東省面上項(xiàng)目限額評審 +3	石瑞0426 2026-03-19	3/150	2026-03-21 18:59 by RRRKKK
[考研] 材料學(xué)碩301分求調(diào)劑 +6	Liyouyumairs 2026-03-21	6/300	2026-03-21 17:42 by JourneyLucky
[考研] 286求調(diào)劑 +8	Faune 2026-03-21	8/400	2026-03-21 17:16 by 學(xué)員8dgXkO
[考研] 0856材料專碩353求調(diào)劑 +3	NIFFFfff 2026-03-20	3/150	2026-03-21 10:23 by luoyongfeng
[考研] 08工科 320總分求調(diào)劑 +6	梨花珞晚風(fēng) 2026-03-17	6/300	2026-03-21 03:40 by JourneyLucky
[考研] 一志愿中國石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-18	6/300	2026-03-21 00:32 by JourneyLucky
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	5/250	2026-03-20 22:11 by 云游重陽
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky
[考研] 環(huán)境工程調(diào)劑 +9	大可digkids 2026-03-16	9/450	2026-03-20 17:38 by 醉在風(fēng)里
[考研] 08工學(xué)調(diào)劑 +5	用戶573181 2026-03-20	5/250	2026-03-20 15:47 by xia_2003
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky