代码搜索结果

% make the contour plot for wftest % Copyright 1999 by Todd K. Moon x1 = [0:.05:1]; x2 = [0:.05:1]; [X1,X2] = meshgrid(x1,x2); J = 10*(cos(2*pi*(X1))+1).^(1) .* (2*X2-.7).^(2)+10; clf subplot

% Test the eigenfilter stuff % Copyright 1999 by Todd K. Moon Nplot = 256 ; wp = 0.2; ws = 0.25; N = 45; alpha = 0.2; h1 = eigfil(wp*pi,ws*pi,N,alpha); [mh1,w] = freqz(h1,1,Nplot); mh1 = mh

% Test the eigenfilter stuff % Copyright 1999 by Todd K. Moon Nplot = 256 ; wp = 0.2*pi; ws = 0.25*pi; N = 45; alpha = 0.2; W0 = [0 0.4*pi .5*pi .8*pi]'; d = [1 1 0 0]'; %$$$ W0 = [0 0.

function [P,Q] = eigmakePQ(wp,ws,m) % % Make the P and Q matrices for eigenfiltering % % function [P,Q] = eigmakePQ(wp,ws,N) % % wp = passband frequency % ws = stopband frequency % N = number

function [db,mag,pha,grd,w] = freqz_m(b,a); % Modified version of freqz subroutine % ------------------------------------ % [db,mag,pha,grd,w] = freqz_m(b,a); % db = Relative magnitude in dB comp

% Chapter 7 Example 7.25 % Highpass filter design using PM algorithm % ws = 0.6*pi; wp = 0.75*pi; Rp = 0.5; As = 50; delta1 = (10^(Rp/20)-1)/(10^(Rp/20)+1); delta2 = (1+delta1)*(10^(-As

% Chapter 7: Example 7.14 % Freq. Samp. Tech.: Lowpass,Naive method % M = 20; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l; Hrs = [1,1,1,zeros(1,15),1,1]; Hdr = [1,1,0,0]; wdl = [0,0.25

% Chapter 7: Example 7.18 % Freq. Samp. Tech.: Highpass, Optimum method T1 % ws=0.6pi, wp=0.8pi, Rp=1dB, As=50dB % M=33, T1 = 0.1095; T2 = 0.598; M = 33; alpha = (

% Chapter 7: Example 7.19 % Freq. Samp. Tech.: Differentiator % M = 33; alpha = (M-1)/2; Dw = 2*pi/M; l = 0:M-1; wl = Dw*l; k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1; Hrs = [j*

% Chapter 7: Example 7.15 % Freq. Samp. Tech.: Lowpass,Optimum method % % (a) T1 = 0.39 M = 40; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l; Hrs = [ones(1,5),0.3904,zeros(1,29),0.3904,