clear; close all; clc; N=2^5; w_kaiser=kaiser(N,1); amp=2.0*fft(w_kaiser,2048)/N; freq = linspace(-0.5,0.5,length(amp)); mag = 20*log10(abs(fftshift(amp/max(abs(amp))))); plot(freq,mag); title('Kaiser window: beta=1.0'); figure(); w_kaiser…

clear; close all; clc; N=2^5; w_hamming = hamming(N); amp = 2.0*fft(w_hamming,2048)/N; freq = linspace(-0.5,0.5,length(amp)); mag = 20*log10(abs(fftshift(amp/max(abs(amp))))); plot(freq,mag);

clear; close all; clc; N=2^5; wave = boxcar(N); amp =2.0*fft(wave, 2048)/N; freq = linspace(-0.5,0.5,length(amp)); mag = 20*log10(abs(fftshift(amp/max(abs(amp))))); plot(freq,mag);

信号 FFT clear; close all; clc; N=2^8; fs = 100; T=1/fs; f1=10;f2=20;f3=40; A1=1;A2=0.5;A3=0.8; t=0:T:(N-1)*T; scale = 0.3; noise = rand(1,length(t))*scale; x=A1*sin(2*pi*f1*t)+A2*sin(2*pi*f2*t)+A3*sin(2*pi*f3*t)+noise; plot(t,x); xlim([0,…

信号 フーリエ変換 逆フーリエ変換 clear; close all; clc; N=2^8; fs = 100; T=1/fs; f1=10;f2=20;f3=40; A1=1;A2=0.5;A3=0.8; t=0:T:(N-1)*T; scale = 0.3; noise = rand(1,length(t))*scale; x=A1*sin(2*pi*f1*t)+A2*sin(2*pi*f2*t)+A3*sin(2*pi*f3*t)+no…

clear; close all; clc; % fs = 3; time = linspace(0,1,100); plot(time,SF(time,3),'DisplayName','Partial Sum (fs=3)'); hold on ; plot(time,SF(time,4),'DisplayName','Partial Sum (fs=4)'); plot(time,SF(time,5),'DisplayName','Partial Sum (fs=5)…