matlab上机完美版全解.doc

  1. 1、本文档共19页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
matlab上机完美版全解

已知3阶椭圆IIR数字低通滤波器的性能指标为:通带截止频率0.4π,通带波纹为0.6dB,最小阻带衰减为32dB。设计一个6阶全通滤波器对其通带的群延时进行均衡。绘制低通滤波器和级联滤波器的群延时。 %Q1_solution %ellip(N,Ap,Ast,Wp)//双线性法设计低通滤波器 %N---The order of the filter %Ap--ripple in the passband %Ast-a stopband Rs dB down from the peak value in the passband %Wp--the passband width [be,ae]=ellip(3,0.6,32,0.4); hellip=dfilt.df2(be,ae); f=0:0.001:0.4; g=grpdelay(hellip,f,2); g1=max(g)-g; [b,a,tau]=iirgrpdelay(6,f,[0 0.4],g1); hallpass=dfilt.df2(b,a);//级联 hoverall=cascade(hallpass,hellip); hFVT=fvtool([hellip,hoverall]); set(hFVT,Filter,[hellip,hoverall]); legend(hFVT,Lowpass Elliptic filter,Compensated filter);//添加图例的标注 clear; [num1,den1]=ellip(3,0.6,32,0.4); [GdH,w]=grpdelay(num1,den1,512); plot(w/pi,GdH); grid xlabel(\omega/\pi); ylabel(Group delay, samples); F=0:0.001:0.4; g=grpdelay(num1,den1,F,2); % Equalize the passband Gd=max(g)-g; % Design the allpass delay equalizer [num2,den2]=iirgrpdelay(6,F,[0,0.4],Gd); [GdA,w] = grpdelay(num2,den2,512); hold on; plot(w/pi,GdH+GdA,r); legend(Original Filter,Compensated filter); 设计巴特沃兹模拟低通滤波器,其滤波器的阶数和3-dB截止频率由键盘输入,程序能根据输入的参数,绘制滤波器的增益响应。 clear; N=input(Type in the order N = ); Wn=input(Type in the 3-dB cutoff frequency Wn = ); %模拟频率 [num,den]=butter(N,Wn,s); w=0:2*Wn; h=freqs(num,den,w); plot(w,20*log(abs(h))),grid; 3.已知系统的系统函数为: 用MATLAB进行部分分式展开,并写出展开后的表达式。 % Partial-Fraction Expansion of Rational z-Transform num = [0 0 1 -0.2 0.5]; den = [1 3.2 1.5 -0.8 1.4]; [r,p,k] = residuez(num,den); disp(Residues);disp(r) disp(Poles);disp(p) disp(Constants);disp(k) 设计切比雪夫I型IIR数字高通滤波器,其性能指标为:通带波纹,最小阻带衰减,通带和阻带边缘频率和绘制所设计的滤波器增益响应。 %a4 disp(prewapping is done,and T=2); Wp = tan(0.75*pi/2); Ws = tan(0.5*pi/2); Rp = 0.5; Rs = 43; [N,Wn] = cheb1ord(Ws,Wp,Rp,Rs,s); [b,a] = cheby1(N,Rp,Wn,s); [bt,at]=lp2hp(b,a,Wp); [num,den]=bilinear(bt,at,0.5); [h,omega] = freqz(num,den); plot (omega/pi,20*log10(abs(h)));grid; xlabel(\omega/\pi); ylabel(Gain); title(Type I Chebyshev Highpass Filter); clear;%预畸变Rp=0.5; Rs=43; Wp=0.75

文档评论(0)

ee88870 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档