matlab - 如何将不规则波变为傅里叶级数

Question

我制作了由 16 个具有随机相位差的规则波组成的不规则波。此外，我试图将不规则波分成几个波，它们是傅立叶级数的形式。所以我在第 11 行下面编写了代码。但是当我运行这段代码时，即使 20 分钟后我也无法得到结果。需要很长时间才能得到结果也让我担心，但更大的问题是我不确定结果是否正确。如何将“不规则波”转换为“傅立叶级数”？

 function regular
 t = (0 : 0.01 : 300); %% time
 w = (0.25:0.05:1);   %% wave frequency
 A = [0 0.14 0.45 0.8 0.95 0.95 0.88 0.77 0.67 0.6 0.53 0.45 0.4 0.36 0.32 0.28]; %% wave amplitude

  phase1 = 2*pi*rand(1,length(w));   %% phase difference
   for ii = 1:length(w)
     for jj = 1:length(t)       
     elev(ii,jj) = (A(ii)*cos(-w(ii)*t(jj)+phase1(ii))); %% wave elevation 
   end
 end
syms t
 T0 = 125;
 w0 = 2*pi/T0;
 a0 = (2/T0)*int(sum(elev),t,-T0/2,T0/2);


 for k = 1:5
     a(k)=(2/T0)*int(sum(elev)*cos(k*w0*t),t,-T0/2,T0/2);
     b(k)=(2/T0)*int(sum(elev)*sin(k*w0*t),t,-T0/2,T0/2);
     F1(k) = a(k)*cos(k*w0*t)+b(k)*sin(k*w0*t);
 end
figure;
plot(t,a0/2+sum(F1));
end