- 1、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
南京工程学院通信工程学院数字信号处理第九章
1.(1)n1=0;n2=10;a1=0.5;
na1=n1:n2; x1=a1.^na1;
stem(na1,x1,k);
(2)n1=0;n2=10;a1=2;
na1=n1:n2; x1=a1.^na1;
stem(na1,x1,k);
(3)n=0:20;
x=sin(n*pi/5);
stem(n,x);
(4)n=0:20;
x=cos(n*pi/10-pi/5);
stem(n,x);
2.(1)a=[3,4,1];b=[1,1,0];n=32;
hn=impz(b,a,n);
stem(hn,k);
(2)a=[5/2,6,10];b=[1,0,0];n=32;
hn=impz(b,a,n);
stem(hn,k);
3. n1=0:15;n01=0; h1=[(n1-n01)=0];
n2=n1;n02=10; h2=[(n2-n02)=0];
h=(7/8).^n1.*(h1-h2);
x1=[(n1-n01)=0];
n2=n1;n03=5;x2=[(n2-n03)=0];
x=(x1-x2);
y=conv(x,h);
stem(y);
4. b=[2,16,44,56,32];a=[3,3,-15,18,-12];
[r p c]=residuez(b,a)
N=20;n=0:N-1;
x=r(1)*p(1).^n +r(2)*p(2).^n
stem(n,x,k);
title(用部分分式法求反变换x(n));
axis([0,N-1,1.1*min(x),1.1*max(x)]);
5.(1)b=[0,2,-1.6,-0.9];a=[1,-2.5,1.96,-0.48];
rz=roots(b);rp=roots(a);
subplot(1,2,1),zplane(b,a);title(系统的零极点分布图);
subplot(1,2,2),impz(b,a,20);title(系统的冲激响应);xlabel(n);ylabel(h(n));
不是因果稳定系统
(2)b=[0,0,1,-1];a=[1,-0.9,-0.65,0.873];
rz=roots(b);rp=roots(a);
subplot(1,2,1),zplane(b,a);title(系统的零极点分布图);
subplot(1,2,2),impz(b,a,20);title(系统的冲激响应);xlabel(n);ylabel(h(n));
因果稳定系统
6. b=[1,0,0];a=[1,-3/4,1/8];
n=(0:1000)*pi/1000;
[h,w]=freqz(b,a,n);
subplot(2,1,1),plot(n/pi,abs(h));grid
axis([0,1,1.1*min(abs(h)),1.1*max(abs(h))]);ylabel(幅度);
subplot(2,1,2),plot(n/pi,angle(h));grid
axis([0,1,1.1*min(angle(h)),1.1*max(angle(h))]);ylabel(相位);xlabel(以pi为单位的频率);
7. n1=0:10;n2=0:18;N=20;
x1=0.8.^n1;x2=0.6.^n2;
y=cconv(x1,x2,N)
8. n=0:20;N=20;
x=(0.9*exp(j*pi/3)).^n;
y=fft(x,20)
9. fp=3400;Omgp=2*pi*fp; fs=4000;Omgs=2*pi*fs; Rp=3;As=40;
[n,Omgc]=buttord(Omgp,Omgs,Rp,As,s);
[z0,p0,k0]=buttap(n);
b0=k0*real(poly(z0)); a0=real(poly(p0)); freqs(b0,a0);
10. T=1;wp=0.2*pi/T;ws=0.3*pi/T;rp=1;rs=15;
[N,wc]=cheb1ord(wp,ws,rp,rs,s);
[B,A]=cheby1(N,rp,wc,s);[Bz,Az]=impinvar(B,A);
[H,w]=freqz(Bz,Az);
plot(w/pi,abs(H));xlabel(频率/{\pi});ylabel(幅度);
11. T=1;Fs=1/T;wp=0.25*pi/T;ws=0.4*pi/T;rp=2;rs=20;
Omgp=(2/T)*tan(wp/2);Omgs=(2/T)*tan(ws/2);
[N,Omgc]=ellipord(Omgp,
文档评论(0)