版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

wangoec

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 519.8
帖子: 59
在線: 52.6小時(shí)
蟲號(hào): 1925630
注冊(cè): 2012-08-05
性別: GG
專業(yè): 電磁環(huán)境與電磁兼容

[求助] 含分布源的雙導(dǎo)體傳輸線在電磁脈沖激勵(lì)下響應(yīng)的計(jì)算FDTD程序請(qǐng)教

含分布源的雙導(dǎo)體傳輸線在電磁脈沖激勵(lì)下響應(yīng)的計(jì)算FDTD程序
程序如下：
clear;clc;
tic
Z1=50;Z2=50;%端接負(fù)載
L=2;h=1e-3;%線長(zhǎng)度和兩線間距
T=4*1e-7;%仿真時(shí)間
e0=8.85e-12;u0=12.57e-7;
v=1/sqrt(u0*e0);%模式傳播速度
C0=14.845*1e-12;L0=0.7485*1e-6;%分布參數(shù)
dez=0.01;det=1/3*1e-10;%空間和時(shí)間步長(zhǎng)
a=90*pi/180;thi=90*pi/180;fi=90*pi/180;  %入射平面波角度
max_space=fix(L/dez);%空間點(diǎn)數(shù)
max_time=fix(T/det);%時(shí)間點(diǎn)數(shù)
V=zeros(max_time+1,max_space+1);%初始化電壓
I=zeros(max_time+1,max_space+1);%初始化電流
E0=zeros(max_time+201,1);%激勵(lì)源初始化
for n=1:max_time
ex=sin(a)*sin(fi);ez=-sin(a)*cos(fi)*sin(thi)+cos(a)*cos(thi);
vx=-v/cos(fi);vz=-v/(sin(fi)*sin(thi));
E0(n)=38*(exp(-1.2*n*det*1e7)-exp(-8.8*n*det*1e7));%激勵(lì)源
      V(n+1,1)=1/((dez/det)*Z1*C0+1)*(((dez/det)*Z1*C0-1)*V(n,1)-2*Z1*I(n,1)-1*h*ex*Z1*C0*dez/det*(E0(n+1)-E0(n)));%始端電壓
for k=2:max_space
      if fix(n+1-(k-1)*dez/det/vz)<=1;
         V(n+1,k)=V(n,k)-det/dez/C0*(I(n,k)-I(n,k-1));%中間點(diǎn)電壓
      else
         V(n+1,k)=V(n,k)-det/dez/C0*(I(n,k)-I(n,k-1))-1*h*ex*(E0(fix(n+1-(k-1)*dez/det/vz))-E0(fix(n-(k-1)*dez/det/vz)));%中間點(diǎn)電壓
      end
end

for s=1:max_space
      if fix(n+1-(s-1)*dez/det/vz)<=1;
         I(n+1,s)=I(n,s)-det/dez/L0*(V(n+1,s+1)-V(n+1,s));%中間點(diǎn)電流
      else
         I(n+1,s)=I(n,s)-det/dez/L0*(V(n+1,s+1)-V(n+1,s))+1*h*(ex/vz-ez/vx)/L0*(E0(fix(n+1-(s-1)*dez/det/vz))-E0(fix(n-(s-1)*dez/det/vz)));%中間點(diǎn)電流
      end
      if fix(n+1-(max_space*dez/det/vz))<=1;
         V(n+1,max_space+1)=1/((dez/det)*Z2*C0+1)*(((dez/det)*Z2*C0-1)*V(n,max_space+1)+2*Z2*I(n,max_space));%末端電壓
      else
         V(n+1,max_space+1)=1/((dez/det)*Z2*C0+1)*(((dez/det)*Z2*C0-1)*V(n,max_space+1)+2*Z2*I(n,max_space)-(1*h*ex*(dez/det)*Z2*C0)*(E0(fix(n+1-(max_space*dez/det/vz)))-E0(fix(n-(max_space*dez/det/vz)))));%末端電壓
      end
end
end
%%%%%%%%%%%%%%%%plots progression of V and I%%%%%%%%%%%%%%%%%%%%%
toc
p=0:max_time;
plot(p*det*1e6,2*V(p+1,1),'r',p*det*1e6,2*V(p+1,201),'b');

1.JPG

2.JPG

3.JPG

4.JPG

5.JPG

回復(fù)此樓

» 本帖@通知

@您好

» 猜你喜歡

328求調(diào)劑已經(jīng)有7人回復(fù)
340求調(diào)劑已經(jīng)有5人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑推薦材料 304 已經(jīng)有15人回復(fù)
291求調(diào)劑已經(jīng)有19人回復(fù)
331環(huán)境科學(xué)與工程求調(diào)劑已經(jīng)有3人回復(fù)
289求調(diào)劑已經(jīng)有8人回復(fù)
315分求調(diào)劑已經(jīng)有7人回復(fù)
295求調(diào)劑已經(jīng)有6人回復(fù)
291求調(diào)劑已經(jīng)有7人回復(fù)

愛(ài)拼才會(huì)贏

1樓 2012-09-26 16:03:58

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

dongge67896

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 51.2
帖子: 29
在線: 2.9小時(shí)
蟲號(hào): 1891660
注冊(cè): 2012-07-13
性別: MM
專業(yè): 計(jì)算機(jī)網(wǎng)絡(luò)

支持樓主帖子太高了沒(méi)爬到頭

贊一下

回復(fù)此樓

http://www.ltbaopi.com http://www.hbyjbt.com/

2樓2012-09-27 09:30:36

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 wangoec 的主題更新