Question

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

 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

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

0 个答案: