三角级数法进行人工地震动合成matlab程序

请问一下，这是一个用三角级数法进行人工地震动合成matlab程序，关于这个w的取值是怎么计算的呀，她这个直接从a中提取，应该是直接先计算在读取数据，但不知道他是如何计算的，希望帮忙得到解答，感谢感谢

img

a = xlsread( 'C:\Users\13737\Desktop\地震动合成\M.xlsx',1,'C3: C1002') ;
w= a; % 频谱范围
WN= length( w) ;
TT= 40; % 持时
dltw= 2* pi /TT; % 求△w
ag = zeros( 40 /0.02+1,1) ;
sw= 0;
kist = 0;
for n = 0: 40 /0.02
    ckn = 0;
for i = 1: WN
    ck = sqrt(4* swl(w(i)) * dltw) ;
    ckn = ck* cos(w(i) * n* 0.02+rand( 1) * 2* pi) +ckn;
end
ag1(n+1) = ckn;
ag2(n+1) = ft( n* 0.02) ;
ag(n+1) = ft( n* 0.02) * ckn;
end