sth galt i min Matlab kode ..

W

was29e

Guest
hej ..
im skrive en kode til at designe en flerlaget optisk filter, koden går sådan her:
Jeg har N lag, hvor hvert af dem er beskrevet ved hjælp af en karakteristisk matrix m, og den resulterende matrix for hele strukturen er M = Πm fra det første lag til den sidste ..

Jeg skrev koden, men theres sth.forkert det holder fortæller mig forkert brug vertcat!

kan man tage et kig på det og fortælle mig, hvad er det prob.'coz jeg kan ikke finde den ..og en ting mere: Hvordan kan jeg ændre udvidelse af en Matlab fil. m??*********************************************

for w = w% w = 1000:2000 bølgelængde
M = [1 0; 01];
d = dem (1, Nmax) * WC / 4% d: tykkelsen af lag
% n lag hver defineret ved m
for i = 1: Nmax
hvis det mod (i, 2) == 0
ni = NH;
anden
ni = nL;
ende

g = (2 * pi * ni * d1i) / w;% M: charactarestic matrix for alle
% g = fase lag = multiplikation af
m for hvert lag
M11 = cos (g);
M12 = (sin (g)) / ni;
M21 = ni * sin (g);
M22 = cos (g);
M = M * [M11 i * M12, M21 M22];
ende
t = 2 * N0 / (((M (1) i * M (3) * ns) * N0) (i * M (2) M (4) * ns));
% t = fractional transmission

tstar = conj (t);
T = t * tstar .* (ns/n0);
R = 1-T;
ende
plot (R)
axis ([Wmin wmax 0 1])

 
du har skrive problemet på den anden linje:

M = [1 0; 01];

bør

M = [1 0; 0 1];

du glemte rum mellem '0 'og '1'.så matlab Thoth var det kun ét nummer.

 
hej .. thanx, i didnt notcie at ..

en ting mere, kan u fortælle mig, hvordan jeg kan få R som en matrix i den store for-løkken, jeg kan ikke få hele resultaterne af løkken, det giver mig den sidste kun

<img src="http://www.edaboard.com/images/smiles/icon_sad.gif" alt="Ked af det" border="0" />

Lagt efter 4 timer 10 minutter:Så jeg gjorde noget i min kode, og det køre perfekt, da jeg forsøgte at køre efter 3 timer ikke sine arbejdsmetoder og igen jeg kan ikke finde, hvad der er galt ....N0 = 1
ns = 3
w = 1000:2000
wc = 1500
NL = 2
nH = 2,5
Nmax = 3

for w = w
R (W) = R (W);
M = [1 0; 0 1];
d = dem (1, Nmax) * WC / 4;% d tykkelse af lag

for i = 1: Nmax
hvis det mod (i, 2) == 0
ni = NH;
anden
ni = nL;
ende

g = (2 * pi * ni * d (1, i)) / w
% g = fase
M11 = cos (g);
M12 = (sin (g)) / ni;
M21 = ni * sin (g);
M22 = cos (g);
M = M * [M11 i * M12; i * M21 M22];
ende
t = 2 * N0 / (((M (1) i * M (3) * ns) * N0) (i * M (2) M (4) * ns));
% t = fractional transmission
tstar = conj (t);
T = t * tstar .* (ns/n0);
R (W) = 1-T;

ende

w = [Wmin: wmax];
plot (w, R (w))

 
hvad skal R butik?

hvad følgende linie betyder: R (W) = R (W);

 
R bør gemme resultatet af T-1
der er R = R (w)

Om R (W) = R (w) jeg skrev dette, så det holder give mig nye vales af R (W) i slutningen af løkken

 
Måske du mener dette:

N0 = 1
ns = 3
w = 1000:2000
wc = 1500
NL = 2
nH = 2,5
Nmax = 3
R = [];

for w = w
% R (W) = R (W);
M = [1 0; 0 1];
d = dem (1, Nmax) * WC / 4;% d tykkelse af lag

for i = 1: Nmax
hvis det mod (i, 2) == 0
ni = NH;
anden
ni = nL;
ende

g = (2 * pi * ni * d (1, i)) / w
% g = fase
M11 = cos (g);
M12 = (sin (g)) / ni;
M21 = ni * sin (g);
M22 = cos (g);
M = M * [M11 i * M12; i * M21 M22];
ende
t = 2 * N0 / (((M (1) i * M (3) * ns) * N0) (i * M (2) M (4) * ns));
% t = fractional transmission
tstar = conj (t);
T = t * tstar .* (ns/n0);
R = [R 1-T];

ende

% w = [Wmin: wmax];
% plot (w, R)
plot (R)

 

Welcome to EDABoard.com

Sponsor

Back
Top