%------------------------------------------------------------------------
% for example 13.3.2,for fig13.3.2 (c),(d)
%------------------------------------------------------------------------
clear;
%设计IIR滤波器;
N=100;
wp=.4*pi;ws=.6*pi;Fs=1;rp=1.;rs=40;
wap=tan(wp/2);was=tan(ws/2);
[n,wn]=cheb1ord(wap,was,rp,rs,'s');
[z,p,k]=cheb1ap(n,rp);
[bp,ap]=zp2tf(z,p,k);
[bs,as]=lp2lp(bp,ap,wap);
[b,a]=bilinear(bs,as, Fs/2);
[h,w] = freqz(b,a,512,1); g = 20*log10(abs(h));
%对应并联实现;
[r,p,k] = residuez(b,a);
[b1,a1]=residuez(r(1:2),p(1:2),0);
[b2,a2]=residuez(r(3:4),p(3:4),0);
[b3,a3]=zp2tf(0,p(5),r(5));
% 对滤波器系数按4bit量化;
bq1 = truncation(b1,4);
aq1 = truncation(a1,4);
bq2 = truncation(b2,4);
aq2 = truncation(a2,4);
bq3 = truncation(b3,4);
aq3 = truncation(a3,4);
kq = truncation(k,4);
[h1,w] = freqz(bq1,aq1,512,1);
[h2,w] = freqz(bq2,aq2,512,1);
[h3,w] = freqz(bq3,aq3,512,1);
[h4,w] = freqz(kq,1,512,1);
hq=h1+h2+h3+h4;
gq = 20*log10(abs(hq));
subplot(221)
plot(w,g,'b',w,gq,'r:');grid
Ylabel('Gain, dB');
% 对滤波器系数按5bit量化;
bq1 = truncation(b1,5);
aq1 = truncation(a1,5);
bq2 = truncation(b2,5);
aq2 = truncation(a2,5);
bq3 = truncation(b3,5);
aq3 = truncation(a3,5);
kq = truncation(k,5);
[h1,w] = freqz(bq1,aq1,512,1);
[h2,w] = freqz(bq2,aq2,512,1);
[h3,w] = freqz(bq3,aq3,512,1);
[h4,w] = freqz(kq,1,512,1);
hq=h1+h2+h3+h4;
gq = 20*log10(abs(hq));
subplot(222)
plot(w,g,'b',w,gq,'r:');grid
Ylabel('Gain, dB');
没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
收起资源包目录
DSP_MATLAB.rar (99个子文件)
DSP_MATLAB
exa060702_1.m 988B
exa100800_xcov.m 567B
exa100800_corrcoef.m 390B
exa020807.m 593B
exa130302_ab.m 1KB
exa020802_impz.m 389B
exa011003_randn.m 589B
exa030102.m 567B
exa070902_fir2.m 647B
fast_LOT.m 1KB
exa130302_cd.m 2KB
rounding.m 568B
exa090803_deconv.m 524B
exa060701_1.m 747B
exa020804_zplane.m 511B
exa060701_2.m 729B
exa050403.m 1KB
exa030101.m 558B
exa130201.m 696B
exa090300_2.m 586B
exa050902_grpdelay.m 718B
exa070903_remez_2.m 600B
exa070903_remez_1.m 594B
exa011001_rand.m 605B
exa121201_pmcov.m 449B
exa020801_filter.m 500B
exa010101.m 634B
exa020805_residuez.m 635B
exa030203.m 496B
exa070101.m 819B
exa121201_pburg.m 444B
exa121201_pmem.m 445B
exa050903_tf2latc.m 643B
exa020806.m 626B
exa031101_fftfilt.m 814B
exa060603.m 941B
exa080901_dct1.m 642B
exa020803_freqz.m 560B
exa070500_1.m 425B
exa060702_2.m 593B
exa031101_conv.m 718B
exa011005_chirp.m 534B
exa030701.m 694B
exa130100.m 486B
exa100801_cohere.m 1KB
exa070500_3.m 546B
exa110702_specgram.m 665B
exa050501.m 1KB
exa100800_cov.m 389B
exa050401_2.m 604B
exa121201_pcov.m 441B
exa020802.m 974B
exa030202.m 620B
exa130301.m 1KB
exa050402.m 855B
exa070904_remezord_2.m 687B
exa060701_3.m 524B
exa070904_remezord_1.m 602B
truncation.m 575B
exa090802_modulate.m 823B
exa110701_pwelch.m 605B
exa090801_in_de_re.m 844B
exa050905_tf2ss.m 527B
exa050401_1.m 609B
exa121201_peig.m 564B
exa011002_rand.m 540B
exa050101.m 627B
exa080901_dct1_test.m 352B
exa070900.m 925B
exa080902_dct2.m 1KB
exa100802_tfe.m 608B
exa030903.m 957B
exa040801_fft.m 557B
exa030201.m 1KB
test.mat 2KB
exa090501.m 1009B
exa070901_fir1.m 721B
exa060601.m 947B
exa060703_2.m 611B
exa020503.m 992B
exa121202.m 865B
exa050901_filtfilt.m 565B
exa011006_conv.m 561B
exa030702.m 973B
exa040802_czt.m 765B
exa011004_sinc.m 405B
exa090300_1.m 733B
exa050904_latcfilt.m 745B
exa031102_hilbert.m 550B
exa020502.m 494B
exa130401.m 880B
exa121201_pmusic.m 508B
GIRL.BMP 17KB
exa011007_xcorr.m 895B
exa070500_2.m 449B
exa080903_LOT_1D.m 1KB
exa090804_svd.m 478B
exa060703_1.m 835B
exa121201_pyulear.m 451B
共 99 条
- 1
资源评论
shujianyaoxin
- 粉丝: 0
- 资源: 4
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功