以下のフーリエ変換を考える。 フーリエ変換 フーリエ変換 振幅スペクトル clear; close all; clc; tau = 1; f0=5; dt=0.01; T=1/f0; t=-T/2:dt:T/2; y=f(tau,f0,t); plot(t,y,'DisplayName','original'); hold on; xf = @(f) tau*sinc(2*tau*(f-f0))+tau*si…

y(t)=, y(t+2π)=y(t) のフーリエ級数展開。（周期T=2π） 固有周波数f0=1/T=1/2π。 clear; close all; clc; T=2*pi; f0=1/T; dt=0.01; t=-T/2:dt:T/2; f=@(t) exp(cos(t)).*cos(sin(t)); y=f(t); plot(t,y,'DisplayName','original'); hold on; syms x; x2 = …

y(t)=, y(t+2π)=y(t) のフーリエ級数展開。（周期T=2π） 固有周波数f0=1/T=1/2π。以下のように式変形する。 積和公式 したがって N=m+1で完全に元の波を再現できる。 以下はm=20の例。N=19では再現できないが、N=21で完全に復元できる。 clear; close all; c…