How we can display images at different times by using for loop instead of following program and how we can make movie from following program?

1. clear all;
2. clc;
3. close all;
4. dim=32;
5. [kx,ky]=meshgrid(-3:6/(dim-1):3);
6. circ=sqrt(kx.^2+ky.^2)<1;
7. r=sqrt(kx.^2+ky.^2);
8. alp=asin(0.95);
9. k0=1/sin(alp);
10. k=256;
11. kz=sqrt(k0^2-(kx.^2+ky.^2));
12. Gx=-1i*sqrt(k0./kz).*((k0*ky.^2+kz.*kx.^2)./(k0*(kx.^2+ky.^2)));
13. Gy=1i*sqrt(k0./kz).*((kz-k0).*kx.*ky)./(k0*(kx.^2+ky.^2));
14. Gz=1i*sqrt(k0./kz).*(kx./k0);
16. Fx=1i*sqrt(k0./kz).*((kz-k0).*kx.*ky)./(k0*(kx.^2+ky.^2));
17. Fy=-1i*sqrt(k0./kz).*((k0*kx.^2+kz.*ky.^2)./(k0*(kx.^2+ky.^2)));
18. Fz=1i*sqrt(k0./kz).*(ky./k0);
20. ph=(atan2(kx,ky)+pi);
21. ph=rot90(ph,3);
23. la=-sin(ph);
24. lb=rot90(cos(ph),3);
25. z=-128:1:127;
26. c = numel(z);
27. mz=(pi/(3*8*2*tan(alp/2)));
28. n = size(kx);
29. Ex = zeros([k k c]);
30. Ey = zeros([k k c]);
31. Ez = zeros([k k c]);
33. for jj=1:c
35. Ex(:,:,jj)=fftshift(fft2(exp(1i*kz*mz*z(jj)).*(Gx.*la.*circ)-rot90(exp(1i*kz*mz*z(jj)).*Gy.*lb.*circ,3),k,k));
36. Ey(:,:,jj)=fftshift(fft2(exp(1i*kz*mz*z(jj)).*(Gy.*la.*circ)+rot90(exp(1i*kz*mz*z(jj)).*Gx.*lb.*circ,3),k,k));
37. Ez(:,:,jj)=fftshift(fft2(exp(1i*kz*mz*z(jj)).*(Gz.*la.*circ)+rot90(exp(1i*kz*mz*z(jj)).*Gz.*lb.*circ,3),k,k));
38. end
41. Ex=Ex(k/2,:,:);
42. Ey=Ey(k/2,:,:);
43. Ez=Ez(k/2,:,:);
44. Ex=squeeze(Ex);
45. Ey=squeeze(Ey);
46. Ez=squeeze(Ez);
47. Ix=Ex.*conj(Ex);
48. Iy=Ey.*conj(Ey);
49. Iz=Ez.*conj(Ez);
51. I=Ix+Iy+Iz;
53. t=0;
54. I1=((I).*exp(1i*pi*t));
55. t=(2/7);
56. I2=((I).*exp(1i*pi*t));
58. t=(2*2/7);
59. I3=((I).*exp(1i*pi*t));
61. t=(3*2/7);
62. I4=((I).*exp(1i*pi*t));
64. t=(4*2/7);
65. I5=((I).*exp(1i*pi*t));
67. t=(5*2/7);
68. I6=((I).*exp(1i*pi*t));
70. t=(6*2/7);
71. I7=((I).*exp(1i*pi*t));
73. t=7*2/7;
74. I8=((I).*exp(1i*pi*t));
77. figure(7)
78. subplot(2,4,1);imagesc(I1.*conj(I1)),colormap gray;axis image;axis off
79. subplot(2,4,2);imagesc(I2.*conj(I2)),colormap gray;axis image;axis off
80. subplot(2,4,3);imagesc(I3.*conj(I3)),colormap gray;axis image;axis off
81. subplot(2,4,4);imagesc(I4.*conj(I4)),colormap gray;axis image;axis off
82. subplot(2,4,5);imagesc(I5.*conj(I5)),colormap gray;axis image;axis off
83. subplot(2,4,6);imagesc(I6.*conj(I6)),colormap gray;axis image;axis off
84. subplot(2,4,7);imagesc(I7.*conj(I7)),colormap gray;axis image;axis off
85. subplot(2,4,8);imagesc(I8.*conj(I8)),colormap gray;axis image;axis o
86. % % movie(F(t))

 # I what to two display time dependence image of I by using for loop and how we can make movie from this time dependence of I.
#  how we get all images at all times;
# by using for loop t=0:2/7:2;
#   if anyone help me in this regard it is great help for me
#  thanks  in advance,

[Walter Roberson]

Those decimal numbers at the beginning of every line are going to cause syntax errors.

[Image Analyst]

Yeah, makes it so much of a pain to fix that I'm not even going to attempt to copy and paste this into MATLAB. I suggest the author read this tutorial on how to put code in so that we can copy and paste it, as well as read it: http://www.mathworks.com/matlabcentral/answers/13205-tutorial-how-to-format-your-question-with-markup

Answer 1

[Image Analyst]

I can't run your code, so all I can suggest is the FAQ: http://matlab.wikia.com/wiki/FAQ#How_can_I_create_a_movie_from_my_MATLAB_figures.3F and to put in a "drawnow" in the loop to force it to update/refresh the display immediately instead of waiting until it gets around to it or finishes the loop.