Lx=0.1;
Ly=0.2;
Lz=1;
st=0;
wc=2*pi*600;%低頻范圍最高頻率為200Hz
c0=344;
k=wc/c0;
lmax=fix(k*Lx/pi);
mmax=fix(k*Ly/pi);
nmax=fix(k*Lz/pi);
vol=(lmax+1)*(mmax+1)*(nmax+1);
%ymax=max([lmax,mmax,nmax]);
%ymin=min([lmax,mmax,nmax]);
%ymid=lmax+mmax+nmax-ymax-ymin;
%t3,t2,t1,
y=rand(3,1);
y(1)=fix((vol-1)/((mmax+1)*(nmax+1)))
y(2)=fix((vol-1-y(1)*(mmax+1)*(nmax+1))/(nmax+1))
y(3)=vol-1-y(1)*(mmax+1)*(nmax+1)-y(2)*(nmax+1)
lmax
mmax
nmax