1128阅读
5回复

[求助]端接非线性负载和外场激励传输线传输线FDTD编程

UID ：69145

注册：2010-11-11
登录：2015-09-01
发帖：146
等级：仿真二级

0楼发表于: 2012-09-19 15:57:39

各位前辈，很不好意思，还劳您在百忙之中帮我看看这两个问题。我编写了两个程序，一个用于计算多导体端接非线性负载情况下的电磁暂态响应，另一个用于计算平面波激励下双导体传输线的响应。不知道是什么原因，计算结果和算例就是不一致，反复检查就是查不出原因。我把算例和程序传上来，请大家 ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

描述:请输入描述端接非线性负载和外场激励传输线传输线FDTD编程
附件:

端接非线性负载和外场激励传输线传输线FDTD编程.rar (474 K) 下载次数:44

离线wangever

UID ：69145

注册：2010-11-11
登录：2015-09-01
发帖：146
等级：仿真二级

1楼发表于: 2012-09-22 22:16:05

请懂行的前辈帮一下吧，寸心焦虑，诚不知所止

离线wangever

UID ：69145

注册：2010-11-11
登录：2015-09-01
发帖：146
等级：仿真二级

2楼发表于: 2012-09-23 17:24:40

各位前辈，请在百忙之中驻足一下，帮小弟看看这个程序问题出在哪里，这是一个用于计算平面波激励下双导体传输线响应的FDTD程序。不知道是什么原因，计算结果和算例就是不一致，反复检查就是查不出原因。我把算例和程序传上来，请大家施以援手。不胜感激！
含分布源的双导体传输线在电磁脉冲激励下响应的计算FDTD程序

程序如下：
clear;clc;
tic
Z1=50;Z2=50;%端接负载
L=2;h=1e-3;%线长度和两线间距
T=4*1e-7;%仿真时间
e0=8.85e-12;u0=12.57e-7;
v=1/sqrt(u0*e0);%模式传播速度
C0=14.845*1e-12;L0=0.7485*1e-6;%分布参数
dez=0.01;det=1/3*1e-10;%空间和时间步长
a=90*pi/180;thi=90*pi/180;fi=90*pi/180;  %入射平面波角度
max_space=fix(L/dez);%空间点数
max_time=fix(T/det);%时间点数
V=zeros(max_time+1,max_space+1);%初始化电压
I=zeros(max_time+1,max_space+1);%初始化电流
E0=zeros(max_time+201,1);%激励源初始化
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));%激励源
        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));%中间点电压
        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)));%中间点电压
        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));%中间点电流
        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)));%中间点电流
        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');