[求助]二维离散随机介质matlab如何实现

基质上随机分布球形粒子，粒子折射率为大于1，均相同，基质折射率为1，粒子的半径也为随机分布。高手请指点一二，有兴趣的大家一起讨论 .. 2K/t[.8  
yR{rje*  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

最近研究这个程序，就是基质上随机分散球形粒子，粒子的折射率都一样，基质的折射率为1，不知道该如何表达？希望高手可以指点一二，大家也可以与我讨论。

这是连续的二维随机介质的程序，离散的该如何下手呢？ wO!>kc<  
function b=randmedia(Z,X,A,B,mean,variance) v0&DD&mp  
mu=4000;dalta=0.05; EGv]K|  
A=20;%%%%%%%%%%%%%%%%%%%%%%%参数A Y cL((6A  
B=10;%%%%%%%%%%%%%%%%%%%%%%%参数B .i )K#82  
Z=2^9;%%%%%%%%%%%%%%%%%%%%%%%模型大小Z方向 :rufnmsP<U  
X=2^9;%%%%%%%%%%%%%%%%%%%%%%%模型大小Y方向 @JdeOL;  
dz=1;dx=1;%%%%%%%%%%%%%%%%%%%%%%%网格间距 l_04b];  
,'YKL",  
dalta=dalta*mu; [-W~o.`  
AUTO=zeros(Z,X); J<"Z6 '0v  
for i = 1:Z u#u/uS"  
    for j = 1:X OUN~7]OD%  
        AUTO(i,j)=rand(1,2); 46l*ui_  
    end $u,A/7\s  
end qD"~5vtLqQ  
AUTO = fft2(AUTO); @$p6w  
AUTO = real(AUTO)^2+imag(AUTO)^2; beYGP  
% AUTO = abs(AUTO); D=D.s)ns*  
N1y,~Z  
rand('twister', sum(100*clock)); nip*Y@-F  
RAND=rand(Z,X).*(2*pi); ]{|l4e4P  
% RAND=variance+2*variance.*rand(Z,X); aRn""3[  
% P9`CW  
`T%nGVl>\  
REAL=AUTO.*cos(RAND); 7?_gm>]a  
IMAGE=AUTO.*sin(RAND); 9"Dt3>Z  
Rand = complex(REAL,IMAGE); 0#8lg@e8  
Rand = fft2(Rand); %,S:^Rvv  
Rand = real(Rand)^2+imag(Rand)^2; ~(c<M>Q8  
% Rand = abs(Rand); -zt*C&)b  
tmoclK-  
% figure(1); .VmRk9Z  
% imagesc(Rand); &jnBDr  
    a=double(Rand);  % a表示原始图像灰度矩阵 MdkL_YP}.  
(Cti,g~  
    hw=0;   y^X]q[-?  
D&N3LH  
    [m n]=size(a); 2=7[r-*E  
]$UTMuOQl  
    b=zeros(Z-2*hw,X-2*hw); PK{FQ3b2{  
"K|':3n|  
    for i=1+hw:Z-hw  %求平滑矩阵 4/SltWU  
N/2WUp  
        for j=1+hw:X-hw NDe FY  
Qe9}%k6@E  
            wi=a(i-hw:i+hw,j-hw:j+hw); %6V=G5+W  
a9 S&n5  
            b(i-hw,j-hw)=mean2(wi); .",BLuce  
>*l2]3'`  
        end fgo3Gy*#  
cI7aTLC"s  
    end 9CIQRc  
   =zp{ ^mC  
b=zscore(b); aiw4J  
b=mu+dalta*b; n34d"l3  
figure(4); -=u9>S)!c  
mesh(b) 
0+_;6  
[Nz,Nx]=size(b);%randmedia ay{]Vqi9  
MIN=min(b);%randmedia /Vv)00  
MIN=min(MIN); MpJ3*$Dr  
MAX=max(b);%randmedia PUd/|Rc/}  
MAX=max(MAX); O15~\8#'  
N=50; *li5/=UC5*  
step=(MAX-MIN)/N; *TxR2pC}  
nV=zeros(1,N); S->S
p  
% for i = 1:Nz %~Ymb&ugg  
%     for j = 1:Nx O]25{L  
%         for n=1:N ${t$:0R,h  
%           if((MIN+step*(n-1))<b(i,j) & (MIN+step*(n))>=b(i,j))%randmedia UG
]5Dxk  
%             nV(1,n)=nV(1,n)+1; z`dnS]q9  
%             break; uE#i3( J  
%           end z3fU|*_c  
%         end JJ_KfnH  
%     end ; ?,'jI*1  
% end OtT*)8*c  
% figure(5); s"wz !{G4  
% plot([MIN:step:MIN+(N-1)*step],nV); b%j4W)Z  

你好，我想请教一下：加平滑矩阵的目的是什么啊？还有你是想画出随机介质的示意图还是其它什么啊？

我也挺感兴趣的，希望我们可以讨论一下哈。