1 % radiale Gauss-Filter unterschiedlicher Breite bauen
   2 
   3 % gewünschte Bildschirmauflösung
   4 X_Pixel = 1280;
   5 Y_Pixel = 1024;
   6 % X_Pixel = 1024;
   7 % Y_Pixel = 768;
   8 % X_Pixel = 800;
   9 % Y_Pixel = 600;
  10 
  11 for groesse=10:20:190
  12     maxx=X_Pixel;
  13     maxy=Y_Pixel;
  14     [xi,yi] = meshgrid(1:1:maxx,1:1:maxy);
  15     xc=xi-round(maxx/2);
  16     yc=yi-round(maxy/2);
  17     zi=exp(-(xc.*xc+yc.*yc)/(groesse*maxy));
  18 %     [xn,yn] = meshgrid(1:20:X_Pixel,1:20:Y_Pixel);
  19 %     zn=interp2(xi,yi,zi,xn,yn,'linear',0);
  20     fname=sprintf('GaussHole_%d_%dx%d.dat', (groesse/10+1)/2, X_Pixel, Y_Pixel);
  21     f=fopen(fname,'w','ieee-le');
  22     fwrite(f,zi','float32');
  23     fclose(f);
  24 end