sns_fur_ser.m

来自「radar equation for target estimation」· M 代码 · 共 38 行

%periodic signal

t=[-8 -8 -6 -6 0 0 2 2 8 8 10 10];
x=[0 1 1 0 0 1 1 0 0 1 1 0];
% figure(1);
subplot(2,1,1),plot (t,x);
axis ([-10 12 -1.5 1.5]);
xlabel('timt (t) -t');
ylabel('x(t)');
title('rectangular pulse train');

T=8;
W=2*pi*1/8;
% t=linspace(-8,12,1000);
t=-8:.00001:16;
a0=.25;
n=1:50;
an=(1./(pi*n)).*sin(n*pi/2);
bn=(1./(pi*n)).*(1-cos(n*pi/2));

%loop
X1=a0;
for i=1:50
    X1=X1+an(i)*cos(i*W*t)+bn(i)*sin(i*W*t);
end
% 
% figure(2);
subplot(2,1,2),plot(t,X1);
xlabel('time -> t');
ylabel('amp of X(t)');
title('rectangular pulse using fourier series');