matlab codant pour adi FDTD

Firma ART powiększyła swoją ofertę o akcesoria dla osób prowadzących aktywny tryb życia. Bezprzewodowe, douszne słuchawki AP-B21 łączą się z urządzeniami mobilnymi za pomocą modułu Bluetooth. Słuchawki zostały przygotowane z myślą o osobach, które lubią aktywnie spędzać czas i równocześnie cieszyć się ulubioną muzyką. Niewielkie rozmiary i elastycz...

Read more...
 
ces mai être utile, ce
n'est 2d

----------------------
fonction FDTD2_2 = fdtd2_2 (ER1, ER2, MR1, MR2, o1, o2, ET1, ET2, ET3, Et4, Et5, HT1, HT2)
% FDTD2_2 Implements FDTD algorithme pour un domaine de Hy pour (k) = 0 à la
%
right boundary-Part2-Question2
% Attire les graphiques de l'ex en temps de ET1, ET2, ET3, Et4, Hy Et5 et en temps de
% HT1, HT2 où le temps est un entier
% Si vous ne voulez pas d'un graphique à un moment, juste attribuer zéro pour que
%
du temps (Et5 doit être attribué et il doit être le plus grand moment)
% ER1: constante diélectrique du premier moyen
% ER2: constante diélectrique du deuxième moyen
MR1%: permittivité relative du premier moyen
Mr2%: permittivité relative de la deuxième moyen
% o1: Conductivité du premier moyen
% o2: Conductivité de la deuxième moyen

% Initialisations
e1 = ER1 * 10 ^ -9/36/pi;
e2 = ER2 * 10 ^ -9/36/pi;
M1 = MR1 * 4 * pi * 10 ^ -7;
M2 = mr2 * 4 * pi * 10 ^ -7;
v1 = 1 / (sqrt (M1 * e1));
v2 = 1 / (sqrt (M2 * e2));
K = 500;
K_D = 250;
z = (0: K);
Ex = zeros (1, K 1);
Hy = zeros (1, K 1);
dz = 0,015;
dt = 2.5 * 10 ^ -11;
T = 1 * 10 ^ -9;
N = 2500;
h1 = figure;
h2 = figure;

Computations%
pour i = 1: Et5,

Ex1 = Ex (1);
Ex2 = Ex (2);
ExK = Ex (K);
Ex ExK_plus1 = (K 1);
j = 1: K,
si j <= K_D
Hy (j) = Hy (j) - (dt/M1/dz) * (Ex (j 1)-Ex (j));
elseif j == K
Hy (j) = 0;
autre
Hy (j) = Hy (j) - (dt/M2/dz) * (Ex (j 1)-Ex (j));
fin

fin
j = 1: K 1,

if ((i-1) * dt) <= (6 * T)
Ex (1) = exp (-((( i-1) * dt * T-3) ^ 2) / (T ^ 2));
autre
Ex (1) = 0;

fin
si j ~ = 1 & j <= K_D
Ex (j) = (1-dt o1/2/e1 *) / (1 dt * o1/2/e1) * Ex (j) - (dt/dz/e1) / (1 dt * o1 / 2 / e1) * (Hy (j)-Hy (j-1));
elseif j == K_D
1Ex (j) = (1-dt * o1 / 2 / ((E1 E2) / 2)) / (1 dt * o1 / 2 / ((E1 E2) / 2)) * Ex (j) -- (dt / dz / ((E1 E2) / 2)) / (1 dt * o1 / 2 / ((E1 E2) / 2)) * (Hy (j)-Hy (j-1));
elseif j> K_D
1 & j <= K
Ex (j) = (1-dt o2/2/e2 *) / (1 dt * o2/2/e2) * Ex (j) - (dt/dz/e2) / (1 dt * o2 / 2 / e2) * (Hy (j)-Hy (j-1));
fin
fin

% Attire les graphiques
si i == ET1 | i == ET2 | i == ET3 | i == Et4
figure (h1);
fin
si i == ET1
plot (z, Ex, 'k');
tenir;
fin
si i == ET2
plot (z, Ex, 'b');
tenir;
fin
si i == ET3
plot (z, Ex, 'g');
tenir;
fin
si i == Et4
plot (z, Ex, 'r');
tenir;
fin
si i == HT1 | i == HT2
figure (h2);
fin
si i == HT1
plot (z, Hy, 'c');
tenir;
fin
si i == HT2
plot (z, Hy, 'm');
xlabel ( 'z (0.015m)');
ylabel ( 'Hy (A / m)');
title ( 'FDTD2_2 Implements FDTD algorithme pour un domaine de Hy pour (k) = 0 à la frontière droite-Part2-Question2');
legend ( 'Hy à t1', 'Hy à t2');
grille;
tenir;
fin
fin
figure (h1);
plot (z, Ex, 'y');
xlabel ( 'z (0.015m)');
ylabel ( 'Ex (V / m)');
title ( 'FDTD2_2 Implements FDTD algorithme pour un domaine de Hy pour (k) = 0 à la frontière droite-Part2-Question2');
legend ( 'Ex à t1', 'Ex à t2 "," Ex à t3', 't4 à Ex »,« Ex à t5');
grille;
retour

 
c'est matlab code pour FDTD dans 1d.il
n'est pas adi FDTD.
l'un peut
s'il vous plaît télécharger matlab codant pour adi FDTD.Je trouve difficile dans la mise en œuvre.il prend plus de temps pour la simulation.Je trouve aussi bien difficile à mettre en œuvre pml condition.
du moins
s'il vous plaît
quelqu'un peut me guider.

 

Welcome to EDABoard.com

Sponsor

Back
Top