model.m

这是一个基于FDTD时域有限差分原理的Radar计算程序

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Projekt Numerische Feldberechung   %
%                                    %
%   Marco Angliker, Remo Huber       %
%                                    %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function model()

clear all;
close all;

epsilon_r1 = 2; %lossy
epsilon_r2 = 2; %Debey
epsilon_r3 = 2; %Plasma
epsilon_r4 = 2; %Lorentz
epsilon = 2;


omega1 = 1e6:1e6:1000e6;
omega2 = 200e12:1e12:4000e12; 
omega1 = (2*pi).*omega1;
omega2 = (2*pi).*omega2;

%Variablen und Konstanten
vc = 57e12;
omega_p = 2*pi*2000e12;
chii = 2;
t0 = 0.001e-6;
delta = 0.25;
omega_0 = 2*pi*100e6
epsilon_0 = 8.85419e-12;
sigma = 0.01;

% lossy Dielektrikum
epsilon_1 = (sigma./(j*omega1*epsilon_0));
epsilon_1 = epsilon_1 +epsilon_r1;

% frequency dependent media (Debye)
epsilon_2 = (sigma./(j*omega1*epsilon_0)) +(chii./(1+j*omega1*t0));
epsilon_2 = epsilon_r2 + epsilon_2;

% unmagnetized Plasma
epsilon_3 = ((omega_p).^2./(omega2.*(j*vc-omega2)*epsilon_0));
epsilon_3 = (epsilon_3 + epsilon_r3);

% Lorentz Medium
epsilon_4 = epsilon./(1+j*delta*(omega1/omega_0)-(omega1/omega_0).^2);
epsilon_4 = epsilon_r4+epsilon_4;


%Plots
plot((omega1/(2*pi)),log(abs(epsilon_1)),'r');
hold on
plot((omega1/(2*pi)),log(abs(epsilon_2)),'b');
plot((omega1/(2*pi)),log(abs(epsilon_4)),'g');
legend('Lossy Dielektrikum','Debye Medium', 'Lorentz Medium' )
xlabel('Frequenz')
title('relative Permittivit鋞 (logarithmiert)');

figure(2);
plot((omega1/(2*pi)),(abs(epsilon_1)),'r');
hold on
plot((omega1/(2*pi)),(abs(epsilon_2)),'b');
plot((omega1/(2*pi)),(abs(epsilon_4)),'g');
legend('Lossy Dielektrikum','Debye Medium', 'Lorentz Medium' )
xlabel('Frequenz')
title('relative Permittivit鋞 (linear)');

figure(3);
plot((omega2/(2*pi)),abs(epsilon_3));
title('unmagnetized Plasma');

figure(4);
plot(-atan2(imag(epsilon_1),real(epsilon_1)),'r');
hold on;
plot(-atan2(imag(epsilon_2),real(epsilon_2)),'b');
plot(-atan2(imag(epsilon_4),real(epsilon_4)),'g');
legend('Lossy Dielektrikum','Debye Medium', 'Lorentz Medium' )
xlabel('Frequenz')
title('Phase');