How to solve Index out of bounds because numel(A)=4

1 view (last 30 days)
My aim is to generate the intensity profile for the output light passing through a polariser and an analyser for different angles of analyser . I am cosidering a gaussian light beam.
gaussian function-
function H = my_gaussian(n, sigma)
two_sigma_sq = 2*sigma*sigma;
H = zeros(n);
for i = 1:n
for j = 1:n
x = i - n/2;
y = j - n/2;
H(i,j) = exp(-(x*x + y*y)/two_sigma_sq);
Intensity profile-
H = my_gaussian(256, 50) ;
Ph = 0; E1 = cosd(Ph); E2 = sind(Ph);
E = [E1 ; E2];
Th = [0 30 90];
for i = 1:length(Th)
A11 = cosd(Th(i))*cosd(Th(i));
A12 = cosd(Th(i))*sind(Th(i));
A21 = cosd(Th(i))*sind(Th(i));
A22 = sind(Th(i))*sind(Th(i));
A = [A11, A12; A21, A22]
O = A(i)*E
I = O'*O
IP = I*H;
figure, imshow(IP);axis image;colorbar;
I can put maximum 3 values of Th, it gives me correct answer but as soon as I put 4 values, the 4th value is Incorrect. If I try putting more than 4 values error comes up saying
" Attempted to access A(5); index out of bounds because numel(A)=4.
Error in Untitled2 (line 31)
O = A(i)*E "
What should I change in this program?
  1 Comment
KSSV on 12 Feb 2021
What inputs you have tried? The given code now workd fine.

Sign in to comment.

Accepted Answer

KSSV on 12 Feb 2021
This line:
O = A(i)*E ;
should be changed to:
O = A*E

More Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!