代码搜索结果

%例程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

t=0:pi/100:4*pi; s=sin(t+3*pi/4); plot(s); title('原始信号s'); %==================================== %对s进行小波分解：db1 3层 [c,l]=wavedec(s,3,'db1'); %==================================== %提取小波分解的低频系数a

# include # include # include void main() { //program d2r9 //driver for routine polin2 int n,i,j; double pi,x1,x2,f,dy,y,x1a[6],x2a[6],ya[6][6];

/*巴特沃斯滤波器； ln:为滤波器的阶次N；k: 第k个极点；n:分子分母多项式长度；d[]为分子系数；c[]为分母系数；*/ void main() { int ln,k,n,i; double d[],c[]; double pi,tmp; ln= pi=4.0*atan(1.0); /*初始化d[],c[]，分子分母常数项必为1。*/ d