ip_01_02.m
上传用户:loeagle
上传日期:2013-03-02
资源大小:1236k
文件大小:0k
- % MATLAB script for Illustrative Problem 2, Chapter 1.
- echo on
- fnct='lambda';
- 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')