fseries2.m

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

function A=fseries2
%周期矩形脉冲信号的频谱。
S0=zeros(1,100);
S0(1:10)=1;S0(91:100)=1;
plot(S0);
T=100;
N=30;A=zeros(1,N);
[m,n]=size(S0);
for k=1:N
    omg1=2*pi/T;t=linspace(0,T,n);
    y=S0.*exp(-j.*(k-1).*omg1.*t);
    A(k)=trapz(t,y)./T;
end
figure;
tp=1:N;
stem(tp-1,abs(A));