Replace ex3_noisy.m

ex3_noisy.m

@@ -32,7 +32,7 @@ axis_length = 16;  % size of region of interest
 b = 3*epsilon;  % band-width of the ramp filter in filtered backprojection USE!!!
 
 RadonDeltaS = pi/b; 
-RadonDeltaOmega = sqrt(2)/axis_length * pi/b;
+RadonDeltaOmega = 360 / (ceil(360 / (sqrt(2)/axis_length*pi/b)));
 
 plot_flag = 0;
 
@@ -105,7 +105,7 @@ set(gca,'LooseInset',get(gca,'TightInset'));
 RadonRotAngles = (0:RadonDeltaOmega:360);
 RadonRotAngles = RadonRotAngles(1:end-1);
 
-SensorPos_max = 1.1*sqrt(2)*axis_length;
+SensorPos_max = ceil(1.1*sqrt(2)*axis_length/RadonDeltaS)*RadonDeltaS;
 RadonSensorPos = (-SensorPos_max:RadonDeltaS:SensorPos_max);
 
 RadonData = zeros(length(RadonSensorPos),length(RadonRotAngles));