|
|
大龙展翅(金币+5): 2011-01-01 08:14:59
楼主试试下面这个函数行不行,不知道你要哪种随机平面
function phi=randomphase(Rin,N,wx,wy,a)
%Rin=165;N=1024;wx=20;wy=20;a=0.05;
x_in=linspace(-Rin,Rin,N);
[xx_in,yy_in]=meshgrid(x_in);
ran=unifrnd(-1,1,N);
gauss=exp(-(xx_in/wx).^2-(yy_in/wy).^2);
Ran=fft2(ran);
Gauss=fft2(gauss);
Fi=Ran.*Gauss;
phi=a*ifft2(Fi); |
|