clear all;
close all;
set(0,'defaultAxesFontSize',20)
set(0,'defaultTextFontSize',20)
set(0,'defaultAxesFontName','Helvetica')
set(0,'defaultTextFontName','Helvetica')


# 問3 振幅1V, 時間幅 1/4 sec, 周期 1 secの方形パルス列の、複素型フーリエ級数の両側スペクトル成分を図示せよ。
#
A=1 #volt
tau = 1/4 #sec
T = 1 #sec
N=199 #N次高調波まで計算

t= -2*T: T/50: 2*T;

V=zeros(1,N+1);
f=zeros(1,N+1);
V0 = A*tau/T;
v=V0;#DC成分
for n = 1:N
  f(n)=n/T;#周波数
  x= n*pi*tau/T;
  V(n) = V0*sin(x)/x; #=Vn
  v += V(n)*(exp(j*2*pi*n*t/T)+exp(j*2*pi*-n*t/T));#振幅は偶対象、位相は奇対象
endfor

  subplot(2,1,1);
    plot (t, v,'linewidth', 2);
    xlabel('t [sec]')
    ylabel('Amplitude [V]')
  subplot(2,1,2);
    freq(N+1)=0;
    VV(N+1)=V0;
    for n=1:N
      freq(n)=-f(N-n+1);
      freq(N+1+n)=f(n);
      VV(n)=V(N-n+1);
      VV(N+1+n)=V(n);
    endfor  
    if N<11
      stem(freq(1:2*N+1),VV(1:2*N+1),'linewidth', 2);
    else  
      stem(freq(N-8:N+10),VV(N-8:N+10),'linewidth', 2);
    endif
    xlabel('f [Hz]')#,'fontsize',20)
    ylabel('Amplitude [V]')#,'fontsize',20)

#print -deps foo.eps
print eq_1_23.pdf