當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

☆皓月當(dāng)空☆

新蟲 (初入文壇)

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

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

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

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有199人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用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à)值貼推薦，對(duì)您同樣有幫助:

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

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

已閱回復(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的回帖

查看全部 7 個(gè)回答

feixiaolin

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

專家經(jīng)驗(yàn): +518
應(yīng)助: 942 (博后)
貴賓: 1.275
金幣: 3430
散金: 58785
紅花: 532
沙發(fā): 11
帖子: 24215
在線: 2601.8小時(shí)
蟲號(hào): 2139575
注冊(cè): 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的回帖

sumin84

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5
帖子: 1
在線: 3小時(shí)
蟲號(hào): 3398582
注冊(cè): 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的回帖

雪影舞動(dòng)

新蟲 (初入文壇)

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

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

請(qǐ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的回帖

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

[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無際的草原
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8	niko- 2026-03-19	9/450	2026-03-20 21:57 by luoyongfeng
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機(jī)用戶 2026-03-17	9/450	2026-03-20 14:24 by 無懈可擊111
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999