代码搜索结果

%例程3-12 求解单位冲击响应的幅频 N=20; a=[0.4 0.3 -1]; b=[1 -0.2 0.6 0.3]; x=[1 zeros(1,N-1)]; k=0:N-1; y=filter(a,b,x); Fs=1000; w=0:pi/Fs:pi; h=freqz(a,b,w); subplot(2,1,1); stem(k,y);

%例程3-1 计算序列抽样和与序列能量 t=(0:0.1:10)'; x=sin(2*pi*t)+2*cos(pi*t); y=x+randn(size(t)); s=sum(x(1:100)); E=sum(abs(y).^2); subplot(2,2,1); plot(t,x); title('X'); subplot(2,2,2); plot(t,y); ti

%例程3-3 产生方波波形 t=0:0.1*pi:4*pi; x=square(t); y=square(t,30); subplot(2,1,1); plot(t,x); axis([0 5*pi -1.5 1.5]); subplot(2,1,2); plot(t,y); axis([0 5*pi -1.5 1.5]);

%例程3-6 产生Dirichlet函数波形 t=[-4*pi:0.1:4*pi]; x=diric(t,7); y=diric(t,6); subplot(2,1,1); plot(t,x); subplot(2,1,2); plot(t,y);

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

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

%例程4－25 基于频率抽样方法2实现Ⅲ型FIR数字带通滤波器 clear all; N=41; T1=0.1095; alpha= (N-1)/2; l=0:N-1; wl= (2*pi/N)*l; Hrs=[zeros(1,6),T1,ones(1,7),T1,zeros(1,11),T1,ones(1,7),T1,zeros(1,6)];

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

%例程4－23 利用汉宁窗设计Ⅰ型数字带阻滤波器 clear all; Wpl=0.2*pi; Wph=0.8*pi; Wsl=0.4*pi; Wsh=0.6*pi; tr_width=min((Wsl-Wpl),(Wph-Wsh)); %过渡带宽度 N=ceil(6.2*pi/tr_width)

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