R
ramani
Guest
моля някой да ми помогне как да извлечете S-параметър от fdtd result.if симулация някой да MATLAB кодиране за това моля да го качите.
Navigation
Install the app
How to install the app on iOS
Follow along with the video below to see how to install our site as a web app on your home screen.
Note: This feature may not be available in some browsers.
More options
You are using an out of date browser. It may not display this or other websites correctly.
You should upgrade or use an alternative browser.
You should upgrade or use an alternative browser.
C
cakalhunter
Guest
Това може да бъде полезен, това е 3D
и параметри са намерени чрез преобразувание на Фурие
------------------------
Г-н = 1, е = 1
М = г-н * 4 * пи * 10 ^ -7
E = е 8,854 * * 10 ^ -12
H = 0,075 * 10 ^ -2
DX = ч; DZ = ч;
C = 1 / (M * E) ^ 0,5
DT = 0,5 * H / C
Т = 20 * DT
а = 1,5 * 10 ^ -2
Ex = нули (370);
Ez = нули (370);
Хъ = нули (370);
C1 = 350; C2 = 350;
AB1 = (C * DT-DZ) / (C * DT DZ);
AB2 = (C * C * E * DZ * DT) / (2 * DX * (C * DT DZ));
volt_history1 = нули (1,3000);
volt_history2 = нули (1,3000);
Hy_eskilk1 = нули (20,3000);
Hy_eskilk2 = нули (20,3000);
Hy_eskison1 = нули (20,3000);
Hy_eskison2 = нули (20,3000);
син = нули (1511);
son1 = нули (1511);
fftvoltson = нули (1,4096);
fftvolthis = нули (1,4096);
bolum = нули (1,4096);
bolum1 = нули (1,4096);
bolumson = нули (1103);
за к = 1:3000
ако J <= 6 * Т / DT;
Ex (C2: C2 19,1) = Годен до (- (к * DT-3 * T) ^ 2 / T ^ 2);
% elseif J> 6 * Т / DT & & J <= 6 * Т / DT 10
% Ex (C2: C2 19,1) = 0;
% elseif J> 6 * Т / DT 10
% за B = C2: C2 19
% Хъ (B, 1) = Hy_eskilk2 (B-C2 1)-AB1 * (Hy_eskilk1 (B-C2 1)-Хъ (B, 2)) AB2 * (Ez ( 1,2 б) Ez (б 1,1)-Ez (B, 2)-Ez (B, 1));
% в края
приключвам
% ABC buraya
за х = 1:20
Hy_eskilk1 (X, J) = Хъ (C2 X-1, 1);
Hy_eskilk2 (X, J) = Хъ (C2 X-1, 2);
Hy_eskison1 (X, J) = Хъ (1, C2 X-1);
Hy_eskison2 (X, J) = Хъ (2, C2 X-1);
приключвам
за к = 1:370
за I = 1:369
ако к> 1
ако к == 350 & & I <350
Ex (I, K) = 0;
elseif к == 370
Ex (I, K) = 0;
в противен случай
Ex (I, K) = Ex (I, K) - (DT / (DZ * E)) * (Хъ (I, K)-Хъ (I, K-1));
приключвам
в противен случай
Ex (I, K) = Ex (I, K) - (DT / (DZ * E)) * Хъ (I, K);
приключвам
Ако K <370
ако к == 1 & & I> 349 & & J> 6 * Т / DT 10
Хъ (I, 1) = Hy_eskilk2 (I-C2 1, J-1)-AB1 * (Хъ (I, 1)-Хъ (I, 2)) AB2 * (Ez (I 1,2) Ez (I 1,1)-Ez (I, 2)-Ez (I, 1));
% elseif аз == 1 & & K> 349
% Хъ (1, I) = Hy_eskison2 (к 1- C1)-AB1 * (Hy_eskison1 (к 1- C1)-Хъ (2, I))-AB2 * (Ex (2, I 1) Ex (1, I 1)-Ex (2, I)-Ex (1, I));
elseif аз == 1 & & K> 349 & & J> 3
Хъ (1, к) = Hy_eskison2 (к 1- C1, J-1)-AB1 * (Хъ (1, к)-Хъ (2, к))-AB2 * (Ex (2, K 1) Ex (1, к 1)-Ex (2, K)-Ex (1, K));
elseif I ~ = 1
Хъ (I, K) = Хъ (I, K) - (DT / (DZ * m)) * (Ex (I, K 1)-Ex (I, K)) (DT / (M * DX) ) * (Ez (I 1, K)-Ez (I, K));
приключвам
ако аз == 349 & & K <350
Ez (I 1, к) = 0;
elseif аз == 369
Ez (I 1, к) = 0;
elseif I> 1
Ez (I, K) = Ez (I, K) (DT / (DX * E)) * (Хъ (I, K)-Хъ (I-1, K));
в противен случай
Ez (I, K) = Ez (I, K) (DT / (DX * E)) * Хъ (I, K);
приключвам
приключвам
приключвам
приключвам
voltage1 = 0;
за х = 1:20
voltage1 = voltage1 Ex (C2 X-1, C1-15);
приключвам
volt_history1 (к) = voltage1 * DX;
voltage2 = 0;
за х = 1:20
voltage2 = voltage2 Ez (C1-15, C2 X-1);
приключвам
volt_history2 (к) = voltage2 * DX;
% при = C1: C1 19
% Хъ (1, а) = Hy_eskison2 (а 1- C1)-AB1 * (Hy_eskison1 (а 1- C1)-Хъ (2, а))-AB2 * (Ex (2, 1) Ex (1, 1)-Ex (2, а)-Ex (1, а));
% в края
приключвам
% комплект (GCF, "Renderer", "zbuffer ');
% на окото (Ех);
% за I = 1700:2100
Volt_history1% (I) = volt_history1 (I);
Volt_history2% (I) = volt_history2 (I);
% в края
% второстепенна сюжетна линия (2,1,1);
% участък (volt_history1)
% второстепенна сюжетна линия (2,1,2);
% участък (volt_history2)
volt1% = FFT (volt_history1, 4096);
парцел% (ABS (volt1))
%------------------------------------------------- ----------------------
EX1 = нули (20,2900);
Ez1 = нули (20,2900);
Hy1 = нули (20,2900);
C1 = 350; C2 = 350;
AB1 = (C * DT-DZ) / (C * DT DZ);
AB2 = (C * C * E * DZ * DT) / (2 * DX * (C * DT DZ));
volt_history11 = нули (1,3000);
за к = 1:3000
ако J <= 6 * Т / DT;
EX1 (1:19,1) = Годен до (- (к * DT-3 * T) ^ 2 / T ^ 2);
приключвам
за к = 1:2000
за I = 1:19
ако к> 1
EX1 (I, K) = EX1 (I, K) - (DT / (DZ * E)) * (Hy1 (I, K)-Hy1 (I, K-1));
в противен случай
EX1 (I, K) = EX1 (I, K) - (DT / (DZ * E)) * Hy1 (I, K);
приключвам
Hy1 (I, K) = Hy1 (I, K) - (DT / (DZ * m)) * (EX1 (I, K 1)-EX1 (I, K)) (DT / (M * DX) ) * (Ez1 (I 1, K)-Ez1 (I, K));
ако аз == 1 | | I == 20
Ez1 (I, K) = 0;
в противен случай
Ez1 (I, K) = Ez1 (I, K) (DT / (DX * E)) * (Hy1 (I, K)-Hy1 (I-1, K));
приключвам
приключвам
приключвам
voltage11 = 0;
за х = 1:19
voltage11 = voltage11 EX1 (X, C1-15);
приключвам
volt_history11 (к) = voltage11 * DX;
приключвам
% участък (volt_history11)
% комплект (GCF, "Renderer", "zbuffer ');
% на окото (Ez);
за I = 1:3000
voltson (I) = volt_history1 (I)-volt_history11 (I);
приключвам
% участък (volt_history1);
fftvoltson = FFT (voltson, 4096);
fftvolthis = FFT (volt_history11, 4096);
fftvolthis1 = FFT (volt_history2, 4096);
за I = 1:4096
bolum (I) = fftvoltson (I) / fftvolthis (I);
приключвам
парцел% (ABS (bolum));
% участък (fftvolthis)
за I = 1:4096
bolum1 (I) = fftvolthis1 (I) / fftvolthis (I);
приключвампарцел (ABS (bolumson)); легенда ( '| S11 | ^ 2 | S12 | ^ 2');
xlabel ( 'честота (GHz)');
сет (GCA, "XTick", [0 52 103])
сет (GCA, "XTickLabel", ('0 '; '10 ";" 20'))
сет (GCA, "YLim", [0,3 1,1])
ос стегнат
и параметри са намерени чрез преобразувание на Фурие
------------------------
Г-н = 1, е = 1
М = г-н * 4 * пи * 10 ^ -7
E = е 8,854 * * 10 ^ -12
H = 0,075 * 10 ^ -2
DX = ч; DZ = ч;
C = 1 / (M * E) ^ 0,5
DT = 0,5 * H / C
Т = 20 * DT
а = 1,5 * 10 ^ -2
Ex = нули (370);
Ez = нули (370);
Хъ = нули (370);
C1 = 350; C2 = 350;
AB1 = (C * DT-DZ) / (C * DT DZ);
AB2 = (C * C * E * DZ * DT) / (2 * DX * (C * DT DZ));
volt_history1 = нули (1,3000);
volt_history2 = нули (1,3000);
Hy_eskilk1 = нули (20,3000);
Hy_eskilk2 = нули (20,3000);
Hy_eskison1 = нули (20,3000);
Hy_eskison2 = нули (20,3000);
син = нули (1511);
son1 = нули (1511);
fftvoltson = нули (1,4096);
fftvolthis = нули (1,4096);
bolum = нули (1,4096);
bolum1 = нули (1,4096);
bolumson = нули (1103);
за к = 1:3000
ако J <= 6 * Т / DT;
Ex (C2: C2 19,1) = Годен до (- (к * DT-3 * T) ^ 2 / T ^ 2);
% elseif J> 6 * Т / DT & & J <= 6 * Т / DT 10
% Ex (C2: C2 19,1) = 0;
% elseif J> 6 * Т / DT 10
% за B = C2: C2 19
% Хъ (B, 1) = Hy_eskilk2 (B-C2 1)-AB1 * (Hy_eskilk1 (B-C2 1)-Хъ (B, 2)) AB2 * (Ez ( 1,2 б) Ez (б 1,1)-Ez (B, 2)-Ez (B, 1));
% в края
приключвам
% ABC buraya
за х = 1:20
Hy_eskilk1 (X, J) = Хъ (C2 X-1, 1);
Hy_eskilk2 (X, J) = Хъ (C2 X-1, 2);
Hy_eskison1 (X, J) = Хъ (1, C2 X-1);
Hy_eskison2 (X, J) = Хъ (2, C2 X-1);
приключвам
за к = 1:370
за I = 1:369
ако к> 1
ако к == 350 & & I <350
Ex (I, K) = 0;
elseif к == 370
Ex (I, K) = 0;
в противен случай
Ex (I, K) = Ex (I, K) - (DT / (DZ * E)) * (Хъ (I, K)-Хъ (I, K-1));
приключвам
в противен случай
Ex (I, K) = Ex (I, K) - (DT / (DZ * E)) * Хъ (I, K);
приключвам
Ако K <370
ако к == 1 & & I> 349 & & J> 6 * Т / DT 10
Хъ (I, 1) = Hy_eskilk2 (I-C2 1, J-1)-AB1 * (Хъ (I, 1)-Хъ (I, 2)) AB2 * (Ez (I 1,2) Ez (I 1,1)-Ez (I, 2)-Ez (I, 1));
% elseif аз == 1 & & K> 349
% Хъ (1, I) = Hy_eskison2 (к 1- C1)-AB1 * (Hy_eskison1 (к 1- C1)-Хъ (2, I))-AB2 * (Ex (2, I 1) Ex (1, I 1)-Ex (2, I)-Ex (1, I));
elseif аз == 1 & & K> 349 & & J> 3
Хъ (1, к) = Hy_eskison2 (к 1- C1, J-1)-AB1 * (Хъ (1, к)-Хъ (2, к))-AB2 * (Ex (2, K 1) Ex (1, к 1)-Ex (2, K)-Ex (1, K));
elseif I ~ = 1
Хъ (I, K) = Хъ (I, K) - (DT / (DZ * m)) * (Ex (I, K 1)-Ex (I, K)) (DT / (M * DX) ) * (Ez (I 1, K)-Ez (I, K));
приключвам
ако аз == 349 & & K <350
Ez (I 1, к) = 0;
elseif аз == 369
Ez (I 1, к) = 0;
elseif I> 1
Ez (I, K) = Ez (I, K) (DT / (DX * E)) * (Хъ (I, K)-Хъ (I-1, K));
в противен случай
Ez (I, K) = Ez (I, K) (DT / (DX * E)) * Хъ (I, K);
приключвам
приключвам
приключвам
приключвам
voltage1 = 0;
за х = 1:20
voltage1 = voltage1 Ex (C2 X-1, C1-15);
приключвам
volt_history1 (к) = voltage1 * DX;
voltage2 = 0;
за х = 1:20
voltage2 = voltage2 Ez (C1-15, C2 X-1);
приключвам
volt_history2 (к) = voltage2 * DX;
% при = C1: C1 19
% Хъ (1, а) = Hy_eskison2 (а 1- C1)-AB1 * (Hy_eskison1 (а 1- C1)-Хъ (2, а))-AB2 * (Ex (2, 1) Ex (1, 1)-Ex (2, а)-Ex (1, а));
% в края
приключвам
% комплект (GCF, "Renderer", "zbuffer ');
% на окото (Ех);
% за I = 1700:2100
Volt_history1% (I) = volt_history1 (I);
Volt_history2% (I) = volt_history2 (I);
% в края
% второстепенна сюжетна линия (2,1,1);
% участък (volt_history1)
% второстепенна сюжетна линия (2,1,2);
% участък (volt_history2)
volt1% = FFT (volt_history1, 4096);
парцел% (ABS (volt1))
%------------------------------------------------- ----------------------
EX1 = нули (20,2900);
Ez1 = нули (20,2900);
Hy1 = нули (20,2900);
C1 = 350; C2 = 350;
AB1 = (C * DT-DZ) / (C * DT DZ);
AB2 = (C * C * E * DZ * DT) / (2 * DX * (C * DT DZ));
volt_history11 = нули (1,3000);
за к = 1:3000
ако J <= 6 * Т / DT;
EX1 (1:19,1) = Годен до (- (к * DT-3 * T) ^ 2 / T ^ 2);
приключвам
за к = 1:2000
за I = 1:19
ако к> 1
EX1 (I, K) = EX1 (I, K) - (DT / (DZ * E)) * (Hy1 (I, K)-Hy1 (I, K-1));
в противен случай
EX1 (I, K) = EX1 (I, K) - (DT / (DZ * E)) * Hy1 (I, K);
приключвам
Hy1 (I, K) = Hy1 (I, K) - (DT / (DZ * m)) * (EX1 (I, K 1)-EX1 (I, K)) (DT / (M * DX) ) * (Ez1 (I 1, K)-Ez1 (I, K));
ако аз == 1 | | I == 20
Ez1 (I, K) = 0;
в противен случай
Ez1 (I, K) = Ez1 (I, K) (DT / (DX * E)) * (Hy1 (I, K)-Hy1 (I-1, K));
приключвам
приключвам
приключвам
voltage11 = 0;
за х = 1:19
voltage11 = voltage11 EX1 (X, C1-15);
приключвам
volt_history11 (к) = voltage11 * DX;
приключвам
% участък (volt_history11)
% комплект (GCF, "Renderer", "zbuffer ');
% на окото (Ez);
за I = 1:3000
voltson (I) = volt_history1 (I)-volt_history11 (I);
приключвам
% участък (volt_history1);
fftvoltson = FFT (voltson, 4096);
fftvolthis = FFT (volt_history11, 4096);
fftvolthis1 = FFT (volt_history2, 4096);
за I = 1:4096
bolum (I) = fftvoltson (I) / fftvolthis (I);
приключвам
парцел% (ABS (bolum));
% участък (fftvolthis)
за I = 1:4096
bolum1 (I) = fftvolthis1 (I) / fftvolthis (I);
приключвампарцел (ABS (bolumson)); легенда ( '| S11 | ^ 2 | S12 | ^ 2');
xlabel ( 'честота (GHz)');
сет (GCA, "XTick", [0 52 103])
сет (GCA, "XTickLabel", ('0 '; '10 ";" 20'))
сет (GCA, "YLim", [0,3 1,1])
ос стегнат
S
Supaswing
Guest
Доста лесно:
с oAB се на FDTD отговор на пристанището А към extication, идващи от Б.
S11 = FFT (O11, N). / FFT (i1, N);
S21 = FFT (o21, N). / FFT (i1, N);
..
S34 = FFT (O34, N). / FFT (i4, N);
и така нататък.
Изберете N по отношение на честотата способност, която искате и времето за вземане на проби са, p.ex N = T / timestep
С уважение,
Дейвид
с oAB се на FDTD отговор на пристанището А към extication, идващи от Б.
S11 = FFT (O11, N). / FFT (i1, N);
S21 = FFT (o21, N). / FFT (i1, N);
..
S34 = FFT (O34, N). / FFT (i4, N);
и така нататък.
Изберете N по отношение на честотата способност, която искате и времето за вземане на проби са, p.ex N = T / timestep
С уважение,
Дейвид
R
ryfnjust
Guest
Това е лесно, Supaswing е прав!
S
smruti
Guest
ей supaswing,
може ли отново обясни как това твърдение се прилагат в кода публикувано от cakalhunter, за да се определи на S-параметри.
може ли отново обясни как това твърдение се прилагат в кода публикувано от cakalhunter, за да се определи на S-параметри.
S
SAJ25
Guest
благодарение friendds тези кодове са много полезни