高手请帮忙修改

xujiameng

我这个命令流是用到谐波叠加的，但循环语言太多，哪位高手帮忙改一下，让运行时间缩短，小弟感激涕零！
Pi=3.1415926;
CX=16;
CY=10;
CZ=8;
V10=32.6;
K=0.22;
N=2048;
t=120;
t0=0.2;
t1=0.1;
wu=1;
wk=0.001;
w0= (wu-wk)/ N ;
m=30;
X(1)= -25;
X(2)= -21.131;
X(3)= -17.101;
X(4)=-12.941;
X(5)=-8.6824;
X(6)=-4.3578;
X(7)=0;
X(8)=2;
X(9)=5;
X(10)=8;
X(11)=11;
X(12)=14;
X(13)=17;
X(14)=20;
X(15)=23;
X(16)=26;
X(17)=29;
X(18)=32;
X(19)=35;
X(20)=38;
X(21)=41;
X(22)=44;
X(23)=47;
X(24)=50;
X(25)=54.358;
X(26)=58.682;
X(27)=62.941;
X(28)=67.101;
X(29)=71.131;
X(30)=75;

Z(1)=43.301;
Z(2)=45.315;
Z(3)=46.985;
Z(4)=48.296;
Z(5)=49.24;
Z(6)=49.81;
Z(7)=50;
Z(8)=50;
Z(9)=50;
Z(10)=50;
Z(11)=50;
Z(12)=50;
Z(13)=50;
Z(14)=50;
Z(15)=50;
Z(16)=50;
Z(17)=50;
Z(18)=50;
Z(19)=50;
Z(20)=50;
Z(21)=50;
Z(22)=50;
Z(23)=50;
Z(24)=50;
Z(25)=49.81;
Z(26)=49.24;
Z(27)=48.296;
Z(28)=46.985;
Z(29)=45.315;
Z(30)=43.301;


Y(1)=0;
Y(2)=0;
Y(3)=0;
Y(4)=0;
Y(5)=0;
Y(6)=0;
Y(7)=0;
Y(8)=0;
Y(9)=0;
Y(10)=0;
Y(11)=0;
Y(12)=0;
Y(13)=0;
Y(14)=0;
Y(15)=0;
Y(16)=0;
Y(17)=0;
Y(18)=0;
Y(19)=0;
Y(20)=0;
Y(21)=0;
Y(22)=0;
Y(23)=0;
Y(24)=0;
Y(25)=0;
Y(26)=0;
Y(27)=0;
Y(28)=0;
Y(29)=0;
Y(30)=0;



for a=1:m
  V(a)=V10*(Z(a)/10)^K
end


for l=1:N
    wl=wk+(l-1/2)* w0
for b=1:m
   for c=1:m
       L(b,c,l)=[sqrt((16^2)*(X(b)-X(c))^2+(10^2)*(Y(b)-Y(c))^2+(8^2)*(Z(b)-Z(c))^2)]
   
      Ch(b,c,l)=exp((-2)*wl*L(b,c,l)/(V(b)+V(c)))

   end
end
end

for l=1:N
wl=wk+(l-1/2)* w0
for b=1:m
for c=1:m
Z0=Z(b)-Z(c);
   W=(wl*Z0)/(2*Pi*V(c));
if 0.1>=W
Q(b,c,l)=Z0*wl/8*V(c)
  else if     0.125>=W>0.1
Q(b,c,l)=-5*Z0*wl/V(c)+1.25
else
Q(b,c,l) = -Pi+(Pi+Pi)*rand(1)
      end
end
end
end
end
for l=1:N
wl=0.001+(l-1/2)*((1- 0.001)/2048)
X=1200*wl/ V10

for d=1:m
  
   S(d,d)=4*K*(V10^2)*X^2/wl* (1+ (X^2))^(4/3)
  end
end



for l=1:N

for b=1:m
   for c=1:m
      
S(b,c,l)=Ch(b,c,l)*sqrt(S(b,b)* S(c,c))*exp(i*Q(b,c,l))

end
end
H(:,:,l)=chol(S(:,:,l))

   end