×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Are you an
Engineering professional?
Join Eng-Tips Forums!
• Talk With Other Members
• Be Notified Of Responses
• Keyword Search
Favorite Forums
• Automated Signatures
• Best Of All, It's Free!

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

#### Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

# Is it RIGHT?

## Is it RIGHT?

(OP)
n0 = 1;
n1 = 2.52;
n2 = 1.9;
n3 = 2.11;
n4 = 2.39;
n5 = 2.46;
lamda = 330*(10^-3);
pi = 3.14;
p = (n0)/(n2);
q = n0;
r = 1/(n2);
loop = 3;
tsap = 430;
taln = 900;
talgan1 = 107;
talgan2 = 128;
talgant = 583;
nrows = tsap + taln + (talgan1 + talgan2)*loop + talgant;

R = zeros(1,nrows);

for  t0 = 1:1:tsap;              % for saphhire
Gr0 = 0;
Ph0 = (2*(pi)*2*(n0)*(Gr0)*(t0))/lamda;
X = [cos((Ph0)/2) (j/n0)*sin((Ph0)/2); (j*n0)*sin((Ph0)/2) cos((Ph0)/2)];
Ro0 = ((p)*(X(1,1))+(q)*(X(1,2))-(r)*(X(2,1))-(X(2,2)));
Ao0 = ((p)*(X(1,1))+(q)*(X(1,2))+(r)*(X(2,1))+(X(2,2)));
Re0 = (Ro0)./(Ao0);
Rc0 = conj(Re0);
R(t0) = ((Re0).*(Rc0))*100;
tl = t0;
end

for t1 = 1:1:taln;        % for AlN

Gr1 = 0.0003;
t0 = tl + t1;
Ph1 = (2*(pi)*2*(n3)*(Gr1)*(t1))/lamda;
Y = [cos((Ph1)/2) (j/n3)*sin((Ph1)/2); (j*n3)*sin((Ph1)/2) cos((Ph1)/2)];
Ro1 = ((p)*(Y(1,1))+(q)*(Y(1,2))-(r)*(Y(2,1))-(Y(2,2)));
Ao1 = ((p)*(Y(1,1))+(q)*(Y(1,2))+(r)*(Y(2,1))+(Y(2,2)));
Re1 = (Ro1)./(Ao1);
Rc1 = conj(Re1);
R(t0) = ((Re1).*(Rc1))*100;
end
tl = tl + t1;

for k = 1:1:loop;
for t2 = 1:1:talgan1;     % for 20% AlGaN
t0 = tl + (k-1)*(talgan1 + talgan2)+t2;
Gr2 = 0.0003;
Ph2 = (2*(pi)*2*(n1)*(Gr2)*(t2))/lamda;
Z = [cos((Ph2)/2) (j/n1)*sin((Ph2)/2); (j*n1)*sin((Ph2)/2) cos((Ph2)/2)];

if k == 1
M = (Y)*(Z);
elseif  (k > 1 && k<= 20)
M = (M1)*(Z);
end

Ro2 = ((p)*(M(1,1))+(q)*(M(1,2))-(r)*(M(2,1))-(M(2,2)));
Ao2 = ((p)*(M(1,1))+(q)*(M(1,2))+(r)*(M(2,1))+(M(2,2)));
Re2 = (Ro2)./(Ao2);
Rc2 = conj(Re2);
R(t0) = ((Re2).*(Rc2))*100;
end

for t3 = 1:1:talgan2;     % for 70% AlGaN
t0 = tl + k * talgan1 + (k-1) * talgan2 + t3;
Gr3 = 0.0003;
Ph3 = (2*(pi)*2*(n4)*(Gr3)*(t3))/lamda;
W = [cos((Ph3)/2) (j/n4)*sin((Ph3)/2); (j*n4)*sin((Ph3)/2) cos((Ph3)/2)];

M1 = (W)*(M);

Ro3 = ((p)*(M1(1,1))+(q)*(M1(1,2))-(r)*(M1(2,1))-(M1(2,2)));
Ao3 = ((p)*(M1(1,1))+(q)*(M1(1,2))+(r)*(M1(2,1))+(M1(2,2)));
Re3 = (Ro3)./(Ao3);
Rc3 = conj(Re3);
R(t0) = ((Re3).*(Rc3))*100;
end
end
tl = tl + loop * (talgan1 + talgan2);

for t4 = 1:1:talgant;     % for 15% AlGaN
t0 = tl + t4;
Gr4 = 0.0003;
Ph4 = (2*(pi)*2*(n5)*(Gr4)*(t4))/lamda;
V = [cos((Ph4)/2) (j/n5)*sin((Ph4)/2); (j*n5)*sin((Ph4)/2) cos((Ph4)/2)];

M2 = (V)*(M1);

Ro4 = ((p)*(M2(1,1))+(q)*(M2(1,2))-(r)*(M2(2,1))-(M2(2,2)));
Ao4 = ((p)*(M2(1,1))+(q)*(M2(1,2))+(r)*(M2(2,1))+(M2(2,2)));
Re4 = (Ro4)./(Ao4);
Rc4 = conj(Re4);
R(t0) = ((Re4).*(Rc4))*100;
end

plot(R);
grid on; xlabel('Time (s)'); ylabel('Reflectivity (%)');

I am trying to model a Bragg Reflector with Sapphire Substrate and AlGaN with varying Al molar fractions at growth temp and 330 nm wavelength.

Is the graph characteristics for the loop that includes the 20 iterations right

#### Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

#### Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Close Box

# Join Eng-Tips® Today!

Join your peers on the Internet's largest technical engineering professional community.
It's easy to join and it's free.

Here's Why Members Love Eng-Tips Forums:

• Talk To Other Members
• Notification Of Responses To Questions
• Favorite Forums One Click Access
• Keyword Search Of All Posts, And More...

Register now while it's still free!