jiaoantu3_5.m

经典《信号与系统》教程的matlab例程,对深入理解信号与系统相关概念有很大帮助

%========================================================================
%               Line spectrum of periodic signal
%    Program name:jiaoantu3_5
%========================================================================
clear
syms t  n   k  x 
T=input('input the period T=');
Nf=5*T;
x=sym('Heaviside(t+.5)-Heaviside(t-.5)');
%==========================================================================
%                下面的程序段计算傅里叶级数的系数
%==========================================================================
A0=int(x,t,-T/2,T/2);                     %求出三角函数展开系数A0
Fn=int(x*exp(-j*2*pi*n*t/T),t,-T/2,T/2);
for k=1:Nf;
    F_sym(k+1)=double(vpa(subs(Fn,n,k)));
end
Fn=F_sym;
Fn(1)=A0(1);
%====================================
%      下段程序绘制信号的频谱图
%====================================
fn1=fliplr(Fn);
fn=ones(1,length(fn1)-1);
for i=1:Nf;
    fn(i)=fn1(i);
end
Fn=([fn,Fn]);
n=-Nf:Nf;
subplot(2,1,1);
plot(n,Fn,'r')
%stem(n,(Fn),'.r');
ylabel('T*Fn');
xlabel('nW1');