matlaB函数
连续周期三角波信号频谱图
N=10;
n1=-N:-1;
C1=-4*j*sin(n1*pi/2)/pi^2./n1.^2;
C0=0;
n2=1:N;
C2=-4*j*sin(n2*pi/2)/pi^2./n2.^2;
Cn=[C1 C0 C2];
n=-N:N;
subplot(2,1,1);stem(n,abs(Cn));
subplot(2,1,2);stem(n,angle(Cn));
周期性方波频谱
a=1;tao=0.1;t=0.5;
n0=t/tao;
n=0:2*n0;
fn_p=a*tao/t*(sin(n*pi*tao/t+eps*(n==0)))./(n*pi*tao/t+eps*(n==0)); fn_pabs=abs(fn_p);
fn_pang=angle(fn_p);
fn_mabs=fliplr(fn_pabs(2:(2*n0)+1));
fn_mang=-fliplr(fn_pang(2:(2*n0)+1));
fnabs=[fn_mabs fn_pabs];
fnang=[fn_mang fn_pang];
subplot(2,1,1);stem((-2*n0:2*n0),fnabs);
text(4,0.11,'amplitude spectrum');
subplot(2,1,2);stem((-2*n0:2*n0),fnang);
text(-2,2,'phase spectrum');