1611阅读
6回复

[求助]多导体传输线串扰的FDTD程序

UID ：69145

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

0楼发表于: 2012-09-16 15:15:55

各位前辈：我最近开始了FDTD的学习，想用之解决传输线耦合电磁场计算的问题，但遇到了很多困难。这几天，对照别人的算例，编写了一个计算三导体传输线串扰的FDTD程序，感觉程序也没出错误，但计算结果就是跟算例对不上，百思不得其解，故求助大伙在百忙之中耐心给我看看，指 ..

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

描述:多导体传输线串扰问题FDTD程序请教
附件:

多导体传输线串扰问题FDTD程序请教.rar (152 K) 下载次数:38

离线shangshan

无我

UID ：96634

注册：2012-07-07
登录：2015-02-11
发帖：710
等级：积极交流五级

1楼发表于: 2012-09-16 18:18:39

来看一下！

离线wbzhou

做做算法还是有点意思！

UID ：64301

注册：2010-08-03
登录：2017-10-27
发帖：525
等级：积极交流四级

2楼发表于: 2012-09-17 08:39:53

%%%%%%%%%%%2012-9-12%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%多导体传输线的FDTD程序%%%%%%
clear;clc;
tic
L0=[1.10418 0.690094;0.690094 1.38019]*1e-6;C0=[40.628 -20.314;-20.314 29.7632]*1e-12;%分布电感和电容
det=5/3*1e-11;dez=5e-3;%分别为空间和时间步长
Rs=[50 0;0 50];RL=[50 0;0 50];%首末两段端接阻抗
N=2001;K=101;M=2;%分别为仿真时间点数、空间点数、传输线数
u=zeros(K,M);i=zeros(K,M);%初始化
E=[1 0;0 1];%单位矩阵
us=[1;0];
for n=1:N-1;
    u(1,:)=(inv(dez/det*Rs*C0+E))*((dez/det*Rs*C0-E)*[u(1,1);u(1,2)]-2*Rs*[i(1,1);i(1,2)]+us);%首端电压
    for k=2:K-1;
        u(k,:)=[u(k,1);u(k,2)]-det/dez*(inv(C0))*([i(k,1);i(k,2)]-[i(k-1,1);i(k-1,2)]);%中间点电压
    end;
    u(K,:)=(inv(RL*C0*dez/det+E))*((RL*C0*dez/det-E)*[u(K,1);u(K,2)]+2*RL*[i(K-1,1);i(K-1,2)]); %末端电压
    for s=1:K-1;
        i(s,:)=[i(s,1);i(s,2)]-det/dez*(inv(L0))*([u(s+1,1);u(s+1,2)]-[u(s,1);u(s,2)]); %中间点电流
    end;

     f1(n)=u(1,1);f2(n)=u(K,1);
end
toc
p=1:N-1;
plot(p*det*1e9,2*f1,'r',p*det*1e9,2*f2,'b');
程序经过简化（只记录当前时刻的值，节约内存），红色处为你的笔误，顺序有调整，可对照自己的代码看看。