Eg:画出x(t)=cos(πt/8),|t|<=4 信号的x(t)的离散幅度和相位谱
main.m
echo on
fnct='Fcos';
a=-4;
b=4;
n=24;
tol=0.1;
xx=Fseries(fnct,a,b,n,tol);
xx1=xx(n+1:-1:2);
xx1=[conj(xx1),xx];
absxx1=abs(xx1);
pause % Press any key to see a plot of the magnitude spectrum.
n1=[-n:n];
stem(n1,absxx1)
title('The Discrete Magnitude Spectrum')
phasexx1=angle(xx1);
pause % Press any key to see a plot of the phase.
stem(n1,phasexx1)
title('The Discrete Phase Spectrum')
Fcos.m
function y=Fcos(x)
% y=Fcos(x) generates a cos signal
y=cos(pi*x/8);
end
Fserious.m
function xx=Fseries(funfcn,a,b,n,tol,p1,p2,p3)
j=sqrt(-1);
args0=[];
for nn=1:nargin-5
args0=[args0,',p',int2str(nn)];
end
args=[args0,')'];