fseries.m

MATLAB 基础编程 包含一些很基础的程序

function A=fseries(S0,T,k)
%周期矩形脉冲信号的频谱。
%S0=zeros(1,100);
%S0(1:10)=1;S0(91:100)=1;
plot(S0);
%T=100;
omg1=2*pi/T;t=linspace(0,T,length(S0));
y=S0.*exp(-j.*(k-1).*omg1.*t);
A=trapz(t,y)./T;
figure;
stem(k-1,abs(A));