Changeset 256 in ETALON for simulations

Timestamp:

Aug 7, 2015, 1:29:55 PM (9 years ago)

Author:

malovyts

Message:

deleted colormap| scripts are too old

Location:

simulations

Files:

• FarZoneCriterium/pw_PlotCrit.m (modified) (1 diff)
• FarZoneCriterium/pw_farZoneCrit.m (modified) (1 diff)
• colormap (deleted)

simulations/FarZoneCriterium/pw_PlotCrit.m

@@ +1 @@
+%{
 pw_pitch = input('Please enter the pitch(period) of the grating (in mm):\n');
 pw_stripnum = input('Please enter the number of strips in grating:\n');
 pw_plotname = input('Please enter Name of the plot:\n', 's');
-pw_farZoneCrit(pw_pitch,pw_stripnum,pw_plotname);
+pw_r0 = input('Please enter distance from detector:\n');
+pw_th_min = input('Please enter minimum theta:\n');
+pw_th_max = input('Please enter maximum theta:\n');
+pw_farZoneCrit(pw_pitch,pw_stripnum,pw_plotname,pw_r0,pw_th_min,pw_th_max);
+%}
+pw_farZoneCrit(0.05,800,'Plot4',220,40,140);
+

simulations/FarZoneCriterium/pw_farZoneCrit.m

@@ -1 @@
-function result=pw_farZoneCrit(d,N,Name)
+function result=pw_farZoneCrit(d,N,Name,r0,theta_min,theta_max)
 y=@(x)N^2*d*(1+cos(x));
 theta=0:180;
 theta=theta*pi/180;
+grating_size=[-(d*N/2),(d*N/2)];
+zone_X=[0,(r0).*cos((theta_min:(theta_max-theta_min)/100:theta_max).*pi/180),0];
+zone_Y=[0,(r0).*sin((theta_min:(theta_max-theta_min)/100:theta_max).*pi/180),0];
 
 fig = figure('visible', 'off');
-plot(y(theta).*cos(theta),y(theta).*sin(theta), 'LineWidth', 2);
+plot(y(theta).*cos(theta),y(theta).*sin(theta),...
+y(theta).*cos(theta)./5,y(theta).*sin(theta)./5,...
+y(theta).*cos(theta)./20,y(theta).*sin(theta)./20,...
+'LineWidth', 2);
+    set(gca, 'FontSize', 18, 'LineWidth', 2);
 xlabel('Z , mm');
 ylabel('Y , mm');
 title({'Far zone criterium in ZY plane';['d =' num2str(d) 'mm N = ' num2str(N)]});
+legend('R far','R far/5','R Far/20');
 print( fig , '-dpng', [Name '_ZY.png']);
 close(fig);
 
 fig = figure('visible', 'off');
-plot(theta*180/pi,y(theta), 'LineWidth', 2);
+plot(y(theta).*cos(theta)/20,y(theta).*sin(theta)/20,zone_X,zone_Y,grating_size,[0 0] ,'LineWidth', 2);
+    set(gca, 'FontSize', 18, 'LineWidth', 2);
+xlabel('Z , mm');
+ylabel('Y , mm');
+title({'Far zone criterium in ZY plane, ';['d =' num2str(d) 'mm N = ' num2str(N)]});
+legend('Far zone criterium/20','Observation zone','Grating')
+print( fig , '-dpng', [Name 'log_ZY.png']);
+close(fig);
+
+
+fig = figure('visible', 'off');
+semilogy(theta*180/pi,y(theta),[theta(1) theta(end)]*180/pi,[r0 r0], 'LineWidth', 2);
+    set(gca, 'FontSize', 18, 'LineWidth', 2);
 xlabel('\theta , deg');
 ylabel('Criterium , mm');
 grid on;
+legend('Far zone criterium','Observation zone');
 %grid minor;
 title({'Far zone criterium related to theta';['d =' num2str(d) 'mm N = ' num2str(N)]});