代码搜索结果

%例程4－6 带通椭圆滤波器的幅频函数 wp1=10; wp2=20; Fs=100; rp=0.5; rs=50; wp1=2*pi*wp1; wp2=2*pi*wp2; Bw=wp2-wp1; Wo=sqrt(wp2*wp1); [z,p,k]=ellipap(7,rp,rs); [A,B,C,D]=zp2ss(z,p,k); [At,Bt,Ct,Dt]=lp2b

%例程4－21 利用三角窗设计Ⅲ型数字带通滤波器 clear all; Wpl=0.4*pi; Wph=0.6*pi; Wsl=0.2*pi; Wsh=0.8*pi; tr_width=min((Wpl-Wsl),(Wsh-Wph)); %过渡带宽度 N=ceil(6.1*pi/tr_width)

%例程4－4 基于巴特沃斯模拟滤波器设计数字带通滤波器 Wp=[0.3*pi,0.4*pi]; Ws=[0.2*pi,0.5*pi]; Ap=3; As=18; [N,wn]=buttord(Wp/pi,Ws/pi,Ap,As); %计算巴特沃斯滤波器阶次和截止频率 [b,a]=butter(N,wn,'bandpass'); %频率变换法设计巴特沃斯带通滤波器

%例程4－3 基于巴特沃斯模拟滤波器设计数字高通滤波器 Wp=0.6*pi; Ws=0.4*pi; Ap=1; As=15; [N,wn]=buttord(Wp/pi,Ws/pi,Ap,As); %计算巴特沃斯滤波器阶次和截止频率 [b,a]=butter(N,wn,'high'); %频率变换法设计巴特沃斯高通滤波器 [db,mag,pha,grd,w]

%例程4－24 基于频率抽样方法1实现Ⅰ型FIR数字高通滤波器 clear all; N=61; T1=0.1095; T2=0.598; alpha= (N-1)/2; l=0:N-1; wl= (2*pi/N)*l; Hrs=[zeros(1,22),T1,T2,ones(1,14),T2,T1,zeros(1,21)]; %理想振幅采样响应 Hdr

%例程4－17 利用汉宁窗设计Ⅳ型数字高通滤波器 clear all; Wp=0.6*pi; Ws=0.4*pi; tr_width=Wp-Ws; %过渡带宽度 N=ceil(6.2*pi/tr_width)+1 %滤波器长度 n=0:1:N-1; Wc=(Ws+Wp)/2;

%例程4－20 利用布拉克曼窗设计Ⅱ型数字带通滤波器 clear all; Wpl=0.4*pi; Wph=0.6*pi; Wsl=0.2*pi; Wsh=0.8*pi; tr_width=min((Wpl-Wsl),(Wsh-Wph)); %过渡带宽度 N=ceil(11*pi/tr_width)+1

%例程4－1 利用冲激响应不变法设计数字低通滤波器 %利用模拟巴特沃斯滤波器设计数字低通滤波器 %冲激响应不变法 wp=0.2*pi; ws=0.3*pi; Rp=1; As=15; T=1; %性能指标 Rip=10^(-Rp/20); Atn=10^(-As/20); OmgP=wp*T; OmgS=ws*T; [N,OmgC]=buttord(OmgP,Omg

%例程7－13自适应陷波滤波器 %自适应陷波器仿真 clear all clc N=400; %总采样长度 t=0:N-1; s=sin(2*pi*t/20); %原始正弦信号 A=0.5; %干扰信号的幅值 fai=pi/3; %干扰

clear all N=10; amf=[-(N-1)/2:(N-1)/2]'; b=[0:0.01:pi]; a=exp(j*pi*amf*cos(b)); wh=[1 1 1 1 1 1 1 1 1 1]; p0=abs(wh*a); p=10*log(p0/max(p0)); plot(b,p); axis([0 pi -100 0]);