fseries3.m

来自「MATLAB 基础编程 包含一些很基础的程序」· M 代码 · 共 17 行

function A=fseries3(S0,T)
%周期矩形脉冲信号的频谱。
%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));
for k=1:30
    y=S0.*exp(-j.*(k-1).*omg1.*t);
    A=trapz(t,y)./T;
end
figure;
stem(k-1,abs(A));