clc;
clear;
x=1:10:2560;
z=1:10:2560;
%a=zeros(256,256);
b=x'*x/200/200+z'*z/100/100;
c=exp(-sqrt(b));
d=fft2(c,256,256);
e=unifrnd(0,2*pi,256,256);
%normplot(e);
f=sqrt(d).*exp(-i*e);
g=ifft2(f,256,256);
h=real(g);
k=mean2(h(:));
l=var(h(:));
m=0.05/sqrt(l)*(h-k);
v=2000*(1+m);
imagesc(v);
- 1
- 2
前往页