📄 coupleline2ndmur.m
字号:
% The main function
clear all;
Timestart=datestr(now)
% 定义仿真区，网格大小，电参数等
% [dx,dy,dz,dt,Nt,W,S,H,TT,Ls,Nx,Ny,Nz,t0,T,NH,eps,u,segema]=Boxsize(0.05e-3,0.05e-3,0.05e-3,0.05e-12,1600,0.3e-3,0.3e-3,0.25e-3,0.05e-3,2e-3,2.5e-3,0.85e-3,6e-3,4.5,1,0);
dx=0.05e-3; dy=0.05e-3; dz=0.05e-3; dt=0.05e-12; % step length
Lx=4e-3;Ly=1.15e-3; Lz=12e-3; H=0.25e-3;
Nx=Lx/dx+1; Ny=round(Ly/dy)+1; Nz=Lz/dz+1;NH=H/dy+1; Nt=1600; % 空间时间网格数
W=0.3e-3; S=0.3e-3; TT=0.05e-3; Er=4.5; Ls=4e-3; % 结构参数
t0=500*dt; T=100*dt; % parameters of Gauss inpulse
u0=pi*4e-7; ur=1; % 真空磁导率
u(1:Nx,1:Ny,1:Nz)=u0*ur;
segema0=0;
segema(1:Nx,1:Ny,1:Nz)=segema0;
eps0=8.854e-12; %真空介电常数
epsr(1:Nx,NH+1:Ny,1:Nz)=1;
epsr(1:Nx,1:NH,1:Nz)=Er;
eps=eps0*epsr;
% 差分方程跌代系数
[Exee,Exeh,Eyee,Eyeh,Ezee,Ezeh]=Coefficient_E(eps,segema,dt,Nx,Ny,Nz);
[Hxhe,Hxhh,Hyhe,Hyhh,Hzhe,Hzhh]=Coefficient_H(u,segema,dt,Nx,Ny,Nz);
% 附初值
Ex(1:Nx,1:Ny,1:Nz)=0.0; Ey(1:Nx,1:Ny,1:Nz)=0.0; Ez(1:Nx,1:Ny,1:Nz)=0.0;
Hx(1:Nx,1:Ny,1:Nz)=0.0; Hy(1:Nx,1:Ny,1:Nz)=0.0; Hz(1:Nx,1:Ny,1:Nz)=0.0;
% preEyx0=Ey(1,2:Ny-1,2:Nz-1); preEyx1=Ey(2,2:Ny-1,2:Nz-1); preEzx0=Ez(1,2:Ny-1,2:Nz-1);preEzx1=Ez(2,2:Ny-1,2:Nz-1); % 一般的表达
preEyx0=0; preEyx1=0; preEzx0=0; preEzx1=0;
preEyxL=0; preEyxL1=0; preEzxL=0; preEzxL1=0;
preExy0=0; preExy1=0; preEzy0=0; preEzy1=0;
preExyL=0; preExyL1=0; preEzyL=0; preEzyL1=0;
preEyz0=0; preEyz1=0; preExz0=0; preExz1=0;
preEyzL=0; preEyzL1=0; preExzL=0; preExzL1=0;
% 导体边界
[Exc,Ezc]=PEC(Nx,Ny,Nz,dx,dy,dz,W,S,TT,NH);
% 激励源位置
[Sxsta,Sxend,Systa,Syend,Szsta,Szend]=Source(round((Nx-(2*W+S)/dx)/2),round((Nx-S/dx)/2),1,NH,Ls/dz,Ls/dz);
% 观察电压的位置
[Oxsta1,Oxend1,Oysta1,Oyend1,Ozsta1,Ozend1]=Observe(round((Nx-S/dx-W/dx)/2),round((Nx-S/dx-W/dx)/2),1,NH-1,Szsta+40,Szsta+40);
[Oxsta2,Oxend2,Oysta2,Oyend2,Ozsta2,Ozend2]=Observe(round((Nx+S/dx+W/dx)/2),round((Nx+S/dx+W/dx)/2),1,NH-1,Szsta+40,Szsta+40);
% 观察电流的位置
[OIxsta1,OIxend1,OIysta1,OIyend1,OIzsta1,OIzend1]=Observe(round((Nx-S/dx)/2-W/dx-2),round((Nx-S/dx)/2+2),NH-2,NH+2,Szsta+40,Szsta+40);
[OIxsta2,OIxend2,OIysta2,OIyend2,OIzsta2,OIzend2]=Observe(round((Nx+S/dx)/2-2),round((Nx+S/dx)/2+W/dx+2),NH-2,NH+2,Szsta+40,Szsta+40);
for n=1:Nt
% 保存前时刻的场值，备Mur吸收边界所用
[Exy0,Exy1,Exy0addx,Exy0subx,Exy1addx,Exy1subx,Exy0addz,Exy0subz,Exy1addz,Exy1subz]=PreviousExy0(Ex,Nx,Ny,Nz);
[ExyL,ExyL1,ExyLaddx,ExyLsubx,ExyL1addx,ExyL1subx,ExyLaddz,ExyLsubz,ExyL1addz,ExyL1subz]=PreviousExyL(Ex,Nx,Ny,Nz);
[Exz0,Exz1,Exz0addx,Exz0subx,Exz1addx,Exz1subx,Exz0addy,Exz0suby,Exz1addy,Exz1suby]=PreviousExz0(Ex,Nx,Ny,Nz);
[ExzL,ExzL1,ExzLaddx,ExzLsubx,ExzL1addx,ExzL1subx,ExzLaddy,ExzLsuby,ExzL1addy,ExzL1suby]=PreviousExzL(Ex,Nx,Ny,Nz);
[Eyx0,Eyx1,Eyx0addy,Eyx0suby,Eyx1addy,Eyx1suby,Eyx0addz,Eyx0subz,Eyx1addz,Eyx1subz]=PreviousEyx0(Ey,Nx,Ny,Nz);
[EyxL,EyxL1,EyxLaddy,EyxLsuby,EyxL1addy,EyxL1suby,EyxLaddz,EyxLsubz,EyxL1addz,EyxL1subz]=PreviousEyxL(Ey,Nx,Ny,Nz);
[Eyz0,Eyz1,Eyz0addx,Eyz0subx,Eyz1addx,Eyz1subx,Eyz0addy,Eyz0suby,Eyz1addy,Eyz1suby]=PreviousEyz0(Ey,Nx,Ny,Nz);
[EyzL,EyzL1,EyzLaddx,EyzLsubx,EyzL1addx,EyzL1subx,EyzLaddy,EyzLsuby,EyzL1addy,EyzL1suby]=PreviousEyzL(Ey,Nx,Ny,Nz);
[Ezx0,Ezx1,Ezx0addy,Ezx0suby,Ezx1addy,Ezx1suby,Ezx0addz,Ezx0subz,Ezx1addz,Ezx1subz]=PreviousEzx0(Ez,Nx,Ny,Nz);
[EzxL,EzxL1,EzxLaddy,EzxLsuby,EzxL1addy,EzxL1suby,EzxLaddz,EzxLsubz,EzxL1addz,EzxL1subz]=PreviousEzxL(Ez,Nx,Ny,Nz);
[Ezy0,Ezy1,Ezy0addx,Ezy0subx,Ezy1addx,Ezy1subx,Ezy0addz,Ezy0subz,Ezy1addz,Ezy1subz]=PreviousEzy0(Ez,Nx,Ny,Nz);
[EzyL,EzyL1,EzyLaddx,EzyLsubx,EzyL1addx,EzyL1subx,EzyLaddz,EzyLsubz,EzyL1addz,EzyL1subz]=PreviousEzyL(Ez,Nx,Ny,Nz);
% 以下是电场的迭代
Ex(1:Nx,2:Ny-1,2:Nz-1)=Exee(1:Nx,2:Ny-1,2:Nz-1).* Ex(1:Nx,2:Ny-1,2:Nz-1)+Exeh(1:Nx,2:Ny-1,2:Nz-1).*((Hz(1:Nx,2:Ny-1,2:Nz-1)-Hz(1:Nx,1:Ny-2,2:Nz-1))/dy-(Hy(1:Nx,2:Ny-1,2:Nz-1)-Hy(1:Nx,2:Ny-1,1:Nz-2))/dz);
Ey(2:Nx-1,1:Ny,2:Nz-1)=Eyee(2:Nx-1,1:Ny,2:Nz-1).* Ey(2:Nx-1,1:Ny,2:Nz-1)+Eyeh(2:Nx-1,1:Ny,2:Nz-1).*((Hx(2:Nx-1,1:Ny,2:Nz-1)-Hx(2:Nx-1,1:Ny,1:Nz-2))/dz-(Hz(2:Nx-1,1:Ny,2:Nz-1)-Hz(1:Nx-2,1:Ny,2:Nz-1))/dx);
Ez(2:Nx-1,2:Ny-1,1:Nz)=Ezee(2:Nx-1,2:Ny-1,1:Nz).* Ez(2:Nx-1,2:Ny-1,1:Nz)+Ezeh(2:Nx-1,2:Ny-1,1:Nz).*((Hy(2:Nx-1,2:Ny-1,1:Nz)-Hy(1:Nx-2,2:Ny-1,1:Nz))/dx-(Hx(2:Nx-1,2:Ny-1,1:Nz)-Hx(2:Nx-1,1:Ny-2,1:Nz))/dy);
% 激励信号
Ey(Sxsta:Sxend,Systa:Syend,Szsta:Szend)=-exp(-((n-1)*dt-t0)^2/T^2)/((Syend-Systa)*dy);
% Ey(Sxsta:Sxend,Systa:Syend,Szsta:Szend)=exp(-((n-1)*dt-t0)^2/T^2);
%***** 2nd order Mur absorb boundary condition 需要改进，改进方法见西电书P47,48
c=3e8;
%% for x=0
Ey(1,2:Ny-1,2:Nz-1)=Mur2ndformula(preEyx1,Ey(2,2:Ny-1,2:Nz-1),preEyx0,Eyx0,Eyx1,Eyx0addy,Eyx1suby,Eyx1addy,Eyx1suby,Eyx0addz,Eyx0subz,Eyx1addz,Eyx1subz,dx,dy,dz,dt);
Ez(1,2:Ny-1,2:Nz-1)=Mur2ndformula(preEzx1,Ez(2,2:Ny-1,2:Nz-1),preEzx0,Ezx0,Ezx1,Ezx0addy,Ezx1suby,Ezx1addy,Ezx1suby,Ezx0addz,Ezx0subz,Ezx1addz,Ezx1subz,dx,dy,dz,dt);
%% for x=Lx
Ey(Nx,2:Ny-1,2:Nz-1)=Mur2ndformula(preEyxL1,Ey(Nx-1,2:Ny-1,2:Nz-1),preEyxL,EyxL,EyxL1,EyxLaddy,EyxL1suby,EyxL1addy,EyxL1suby,EyxLaddz,EyxLsubz,EyxL1addz,EyxL1subz,dx,dy,dz,dt);
Ez(Nx,2:Ny-1,2:Nz-1)=Mur2ndformula(preEzxL1,Ez(Nx-1,2:Ny-1,2:Nz-1),preEzxL,EzxL,EzxL1,EzxLaddy,EzxL1suby,EzxL1addy,EzxL1suby,EzxLaddz,EzxLsubz,EzxL1addz,EzxL1subz,dx,dy,dz,dt);
%% for y=0
Ex(2:Nx-1,1,2:Nz-1)=Mur2ndformula(preExy1,Ex(2:Nx-1,2,2:Nz-1),preExy0,Exy0,Exy1,Exy0addx,Exy1subx,Exy1addx,Exy1subx,Exy0addz,Exy0subz,Exy1addz,Exy1subz,dx,dy,dz,dt);
Ez(2:Nx-1,1,2:Nz-1)=Mur2ndformula(preEzy1,Ez(2:Nx-1,2,2:Nz-1),preEzy0,Ezy0,Ezy1,Ezy0addx,Ezy1subx,Ezy1addx,Ezy1subx,Ezy0addz,Ezy0subz,Ezy1addz,Ezy1subz,dx,dy,dz,dt);
%% for y=Ly
Ex(2:Nx-1,Ny,2:Nz-1)=Mur2ndformula(preExyL1,Ex(2:Nx-1,Ny-1,2:Nz-1),preExyL,ExyL,ExyL1,ExyLaddx,ExyL1subx,ExyL1addx,ExyL1subx,ExyLaddz,ExyLsubz,ExyL1addz,ExyL1subz,dx,dy,dz,dt);
Ez(2:Nx-1,Ny,2:Nz-1)=Mur2ndformula(preEzyL1,Ez(2:Nx-1,Ny-1,2:Nz-1),preEzyL,EzyL,EzyL1,EzyLaddx,EzyL1subx,EzyL1addx,EzyL1subx,EzyLaddz,EzyLsubz,EzyL1addz,EzyL1subz,dx,dy,dz,dt);
%% for z=0
Ex(2:Nx-1,2:Ny-1,1)=Mur2ndformula(preExz1,Ex(2:Nx-1,2:Ny-1,2),preExz0,Exz0,Exz1,Exz0addx,Exz1subx,Exz1addx,Exz1subx,Exz0addy,Exz0suby,Exz1addy,Exz1suby,dx,dy,dz,dt);
Ey(2:Nx-1,2:Ny-1,1)=Mur2ndformula(preEyz1,Ey(2:Nx-1,2:Ny-1,2),preEyz0,Eyz0,Eyz1,Eyz0addx,Eyz1subx,Eyz1addx,Eyz1subx,Eyz0addy,Eyz0suby,Eyz1addy,Eyz1suby,dx,dy,dz,dt);
%% for z=Lz
Ex(2:Nx-1,2:Ny-1,Nz)=Mur2ndformula(preExzL1,Ex(2:Nx-1,2:Ny-1,Nz-1),preExzL,ExzL,ExzL1,ExzLaddx,ExzL1subx,ExzL1addx,ExzL1subx,ExzLaddy,ExzLsuby,ExzL1addy,ExzL1suby,dx,dy,dz,dt);
Ey(2:Nx-1,2:Ny-1,Nz)=Mur2ndformula(preEyzL1,Ey(2:Nx-1,2:Ny-1,Nz-1),preEyzL,EyzL,EyzL1,EyzLaddx,EyzL1subx,EyzL1addx,EyzL1subx,EyzLaddy,EyzLsuby,EyzL1addy,EyzL1suby,dx,dy,dz,dt);
% (n-1)dt时刻即前前时刻的旧值
preEyx0=Eyx0; preEyx1=Eyx1; preEzx0=Ezx0;preEzx1=Ezx1;
preEyxL=EyxL; preEyxL1=EyxL1; preEzxL=EzxL;preEzxL1=EzxL1;
preExy0=Exy0; preExy1=Exy1; preEzy0=Ezy0;preEzy1=Ezy1;
preExyL=ExyL; preExyL1=ExyL1; preEzyL=EzyL;preEzyL1=EzyL1;
preEyz0=Eyz0; preEyz1=Eyz1; preExz0=Exz0;preExz1=Exz1;
preEyzL=EyzL; preEyzL1=EyzL1; preExzL=ExzL;preExzL1=ExzL1;
% 导体边界条件
Ex=Ex.*Exc;
Ez=Ez.*Ezc;
% 以下为磁场迭代
Hx(1:Nx,1:Ny-1,1:Nz-1)=Hxhh(1:Nx,1:Ny-1,1:Nz-1).*Hx(1:Nx,1:Ny-1,1:Nz-1)-Hxhe(1:Nx,1:Ny-1,1:Nz-1).*((Ez(1:Nx,2:Ny,1:Nz-1)-Ez(1:Nx,1:Ny-1,1:Nz-1))/dy-(Ey(1:Nx,1:Ny-1,2:Nz)-Ey(1:Nx,1:Ny-1,1:Nz-1))/dz);
Hy(1:Nx-1,1:Ny,1:Nz-1)=Hyhh(1:Nx-1,1:Ny,1:Nz-1).*Hy(1:Nx-1,1:Ny,1:Nz-1)-Hyhe(1:Nx-1,1:Ny,1:Nz-1).*((Ex(1:Nx-1,1:Ny,2:Nz)-Ex(1:Nx-1,1:Ny,1:Nz-1))/dz-(Ez(2:Nx,1:Ny,1:Nz-1)-Ez(1:Nx-1,1:Ny,1:Nz-1))/dx);
Hz(1:Nx-1,1:Ny-1,1:Nz)=Hzhh(1:Nx-1,1:Ny-1,1:Nz).*Hz(1:Nx-1,1:Ny-1,1:Nz)-Hzhe(1:Nx-1,1:Ny-1,1:Nz).*((Ey(2:Nx,1:Ny-1,1:Nz)-Ey(1:Nx-1,1:Ny-1,1:Nz))/dx-(Ex(1:Nx-1,2:Ny,1:Nz)-Ex(1:Nx-1,1:Ny-1,1:Nz))/dy);
% 记录电压
Vs(n)=-sum(Ey(Sxsta,Systa:Syend-1,Szsta)*dy);
V11(n)=-sum(Ey(Oxsta1,Oysta1:Oyend1,Ozsta1)*dy);
V12(n)=-sum(Ey(Oxsta1,Oysta1:Oyend1,Ozsta1+20)*dy);
V13(n)=-sum(Ey(Oxsta1,Oysta1:Oyend1,Ozsta1+40)*dy);
V14(n)=-sum(Ey(Oxsta1,Oysta1:Oyend1,Ozsta1+60)*dy);
V15(n)=-sum(Ey(Oxsta1,Oysta1:Oyend1,Ozsta1+80)*dy);
V21(n)=-sum(Ey(Oxsta2,Oysta2:Oyend2,Ozsta2)*dy);
V22(n)=-sum(Ey(Oxsta2,Oysta2:Oyend2,Ozsta2+20)*dy);
V23(n)=-sum(Ey(Oxsta2,Oysta2:Oyend2,Ozsta2+40)*dy);
V24(n)=-sum(Ey(Oxsta2,Oysta2:Oyend2,Ozsta2+60)*dy);
V25(n)=-sum(Ey(Oxsta2,Oysta2:Oyend2,Ozsta2+80)*dy);
% 记录电流
% 环路积分
I11(n)=sum(Hx(OIxsta1:OIxend1,Oysta1,Ozsta1)*dy);
I11(n)=I11(n)+sum(Hx(OIxsta1:OIxend1,Oyend1,Ozsta1)*dy);
I11(n)=I11(n)+sum(Hy(OIxsta1,Oysta1:Oyend1,Ozsta1)*dy);
I11(n)=I11(n)+sum(Hy(OIxend1,Oysta1:Oyend1,Ozsta1)*dy);
I12(n)=sum(Hx(OIxsta1:OIxend1,Oysta1,Ozsta1+20)*dy);
I12(n)=I12(n)+sum(Hx(OIxsta1:OIxend1,Oyend1,Ozsta1+20)*dy);
I12(n)=I12(n)+sum(Hy(OIxsta1,Oysta1:Oyend1,Ozsta1+20)*dy);
I12(n)=I12(n)+sum(Hy(OIxend1,Oysta1:Oyend1,Ozsta1+20)*dy);
I13(n)=sum(Hx(OIxsta1:OIxend1,Oysta1,Ozsta1+40)*dy);
I13(n)=I13(n)+sum(Hx(OIxsta1:OIxend1,Oyend1,Ozsta1+40)*dy);
I13(n)=I13(n)+sum(Hy(OIxsta1,Oysta1:Oyend1,Ozsta1+40)*dy);
I13(n)=I13(n)+sum(Hy(OIxend1,Oysta1:Oyend1,Ozsta1+40)*dy);
I14(n)=sum(Hx(OIxsta1:OIxend1,Oysta1,Ozsta1+60)*dy);
I14(n)=I14(n)+sum(Hx(OIxsta1:OIxend1,Oyend1,Ozsta1+60)*dy);
I14(n)=I14(n)+sum(Hy(OIxsta1,Oysta1:Oyend1,Ozsta1+60)*dy);
I14(n)=I14(n)+sum(Hy(OIxend1,Oysta1:Oyend1,Ozsta1+60)*dy);
I15(n)=sum(Hx(OIxsta1:OIxend1,Oysta1,Ozsta1+80)*dy);
I15(n)=I15(n)+sum(Hx(OIxsta1:OIxend1,Oyend1,Ozsta1+80)*dy);
I15(n)=I15(n)+sum(Hy(OIxsta1,Oysta1:Oyend1,Ozsta1+80)*dy);
I15(n)=I15(n)+sum(Hy(OIxend1,Oysta1:Oyend1,Ozsta1+80)*dy);
I21(n)=sum(Hx(OIxsta2:OIxend2,Oysta2,Ozsta2)*dy);
I21(n)=I21(n)+sum(Hx(OIxsta2:OIxend2,Oyend2,Ozsta2)*dy);
I21(n)=I21(n)+sum(Hy(OIxsta2,Oysta2:Oyend2,Ozsta2)*dy);
I21(n)=I21(n)+sum(Hy(OIxend2,Oysta2:Oyend2,Ozsta2)*dy);
I22(n)=sum(Hx(OIxsta2:OIxend2,Oysta2,Ozsta2+20)*dy);
I22(n)=I22(n)+sum(Hx(OIxsta2:OIxend2,Oyend2,Ozsta2+20)*dy);
I22(n)=I22(n)+sum(Hy(OIxsta2,Oysta2:Oyend2,Ozsta2+20)*dy);
I22(n)=I22(n)+sum(Hy(OIxend2,Oysta2:Oyend2,Ozsta2+20)*dy);
I23(n)=sum(Hx(OIxsta2:OIxend2,Oysta2,Ozsta2+40)*dy);
I23(n)=I23(n)+sum(Hx(OIxsta2:OIxend2,Oyend2,Ozsta2+40)*dy);
I23(n)=I23(n)+sum(Hy(OIxsta2,Oysta2:Oyend2,Ozsta2+40)*dy);
I23(n)=I23(n)+sum(Hy(OIxend2,Oysta2:Oyend2,Ozsta2+40)*dy);
I24(n)=sum(Hx(OIxsta2:OIxend2,Oysta2,Ozsta2+60)*dy);
I24(n)=I24(n)+sum(Hx(OIxsta2:OIxend2,Oyend2,Ozsta2+60)*dy);
I24(n)=I24(n)+sum(Hy(OIxsta2,Oysta2:Oyend2,Ozsta2+60)*dy);
I24(n)=I24(n)+sum(Hy(OIxend2,Oysta2:Oyend2,Ozsta2+60)*dy);
I25(n)=sum(Hx(OIxsta2:OIxend2,Oysta2,Ozsta2+80)*dy);
I25(n)=I25(n)+sum(Hx(OIxsta2:OIxend2,Oyend2,Ozsta2+80)*dy);
I25(n)=I25(n)+sum(Hy(OIxsta2,Oysta2:Oyend2,Ozsta2+80)*dy);
I25(n)=I25(n)+sum(Hy(OIxend2,Oysta2:Oyend2,Ozsta2+80)*dy);
n
Ey(Sxsta,Systa,Szsta)
Ey(Oxsta1,round((Oysta1+Oyend1)/2),Ozsta1)
end
clear n;
n=200:1:1600;
figure(1);
plot(n,Vs(200:1600));
figure(2);
plot(n,V11(200:1600));
hold on;
plot(n,V12(200:1600));
plot(n,V13(200:1600));
plot(n,V14(200:1600));
plot(n,V15(200:1600));
plot(n,V21(200:1600),'-k');
plot(n,V22(200:1600),'-k');
plot(n,V23(200:1600),'-k');
plot(n,V24(200:1600),'-k');
plot(n,V25(200:1600),'-k');
figure(3);
plot(n,I11(200:1600));
hold on;
plot(n,I12(200:1600));
plot(n,I13(200:1600));
plot(n,I14(200:1600));
plot(n,I15(200:1600));
plot(n,I21(200:1600),'-k');
plot(n,I22(200:1600),'-k');
plot(n,I23(200:1600),'-k');
plot(n,I24(200:1600),'-k');
plot(n,I25(200:1600),'-k');
% 下面进行傅立叶变换计算分布电容和电感参数
Fmax=25e9;
n=1:1*Nt;
k=1;
Nf=15;
df=Fmax/Nf;
while k<=Nf
FI11(k)=sum(I11.*exp(-i*2*pi*(k*df)*(dt*n)))*dt;
FI12(k)=sum(I12.*exp(-i*2*pi*(k*df)*(dt*n)))*dt;
FV11(k)=sum(V11.*exp(-i*2*pi*(k*df)*(dt*n)))*dt;
FV12(k)=sum(V12.*exp(-i*2*pi*(k*df)*(dt*n)))*dt;
FI21(k)=sum(I21.*exp(-i*2*pi*(k*df)*(dt*n)))*dt;
FI22(k)=sum(I22.*exp(-i*2*pi*(k*df)*(dt*n)))*dt;
FV21(k)=sum(V21.*exp(-i*2*pi*(k*df)*(dt*n)))*dt;
FV22(k)=sum(V22.*exp(-i*2*pi*(k*df)*(dt*n)))*dt;
k=k+1;
end
f=df*(1:Nf)*1e-9;
C11=abs(-j*(FI12-FI11)/(20*dz)./FV11./(2*pi*f));
C12=abs(-j*(FI12-FI11)/(20*dz)./FV21./(2*pi*f));
C21=abs(-j*(FI22-FI21)/(20*dz)./FV11./(2*pi*f));
C22=abs(-j*(FI22-FI21)/(20*dz)./FV21./(2*pi*f));
figure(4);
plot(f,C11);
hold on;
plot(f,C12);
plot(f,C21);
plot(f,C22);
Timeend=datestr(now)
Timecosted=Timeend-Timestart;
Timecosted
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -