analy.m

时间域有限差分法（TD-FDM）源代码

%For a given omega the function computes the transform function
%V=exact(omega)
%V-->transfer function
%omega->2*pi*f
%the rest parameters must be changed inside the function

clear all;
clc;

R=0;
G=.1;
L=1;
C=1;
Rs=2;
Rl=1;
Cl=1/(4*pi);
h = 0.125;

omega = -16*pi:0.05*pi:16*pi;
number = max(size(omega));

for k=1:number
    gama=sqrt((R+i*L*omega(k))*(G+i*C*omega(k)));
    Zo=sqrt((R+i*omega(k)*L)/(G+i*omega(k)*C));
    Zl=Rl/(1+i*omega(k)*Cl*Rl);
    A=(Zl-Zo)/(Zl+Zo);
    exact_T(k)=(1+A)/((1+Rs/Zo)*exp(gama*h)+A*(1-Rs/Zo)*exp(-gama*h));
end

figure;
plot(omega,abs(exact_T),'b'); %axis([-16*pi 16*pi 0 .5]);
title('Exact Transfer function');
xlabel('\omega');