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?

3 次查看（过去 30 天）
MOHD
MOHD 2013-4-3
  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,

请先登录，再回答此问题。

回答（1 个）

Image Analyst
Image Analyst 2013-4-3
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.

类别

在 Help Center 和 File Exchange 中查找有关 Startup and Shutdown 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by