资源描述：
《MATLAB编程Zero-poleplane》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、Projects1:PartAProblem1forthefollowingsystem:y(n)+ay(n-1)=x(n);where(a)a=0.9;(b)a=-0.9(a)Zero-poleplane,impulseresponse,
2、H(w)
3、andangleH(w):infigure1.Figure1.DFToftheimpulseresponse(magnitudeandangularspectrum):infigure2.Figure2.ComparisonbetweenDFTandDTFT:Whe
4、nN=10:infigure3.Figure3.WhenN=30:infigure4Figure4.WhenN=50:infigure5.Figure5.(b)Zero-poleplane,impulseresponse,
5、H(w)
6、andangleH(w):infigure6.Figure6.DFToftheimpulseresponse(magnitudeandangularspectrum):infigure7.Figure7.ComparisonbetweenDFTandDTFT:WhenN=10:inf
7、igure8.Figure8.WhenN=30:infigure9.Figure9.WhenN=50:infigure10.Figure10.Problem2:repeatstepsaboveforthefollowingsystem:H(Z)=Where(a)ω0=π/8andr=0.95;(b)ω0=π/4andr=0.95;(c)ω0=π/8andr=0.99;(a)Zero-poleplane,impulseresponse,
8、H(w)
9、andangleH(w):infigure11.Figure11.D
10、FToftheimpulseresponse(magnitudeandangularspectrum):infigure12.Figure12.ComparisonbetweenDFTandDTFT:WhenN=10:infigure13.Figure13.WhenN=30:infigure14.Figure14.WhenN=50:infigure15.Figure15.(b)Zero-poleplane,impulseresponse,
11、H(w)
12、andangleH(w):infigure16.Figure16
13、.DFToftheimpulseresponse(magnitudeandangularspectrum):infigure17.Figure17.ComparisonbetweenDFTandDTFT:WhenN=10:infigure18.Figure18.WhenN=30:infigure19.Figure19.WhenN=50:infigure20.Figure20.(c)Zero-poleplane,impulseresponse,
14、H(w)
15、andangleH(w):infigure21.Figure
16、21.DFToftheimpulseresponse(magnitudeandangularspectrum):infigure22.Figure22.ComparisonbetweenDFTandDTFT:WhenN=10:infigure23.Figure23WhenN=30:infigure24.Figure24.WhenN=50:infigure25.Figure25.Problem1:(a)a=0.91（Z-Pplane&impulseresponse&DTFT）>>b=[1];>>a=[10.9];>
17、>pul=[1zeros(1,59)];>>subplot(221)>>zplane(b,a)>>title('zero-poleplane')>>subplot(222)>>h1=impz(b,a,60);>>stem(h1)>>title('impulseresponse')>>subplot(223)>>[H,w]=freqz(b,a,'whole');>>plot(w,abs(H))>>title('
18、H(w)
19、')>>xlabel('w')>>subplot(224)>>plot(w,angle(H))
20、>>title('angleH(w)')>>xlabel('w')2(DFT)>>b=[1];>>a=[10.9];>>L=100;>>NFFT=2^nextpow2(L);>>N=10;>>h1=impz(b,a,N);>>[H,w]=freqz(b,a,NFFT);>>y=fft(h1,NFFT);>>magY=abs(y);>>phaY=angle(y);>>sub
