代码搜索结果

% Chapter 8: Example 8.12 % Impulse Invariance Transformation % Chebyshev-1 Lowpass Filter Design % % Digital Filter Specifications: wp = 0.2*pi;

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

/* whattime.c For Free Chat By Bill Kendrick kendrick@zippy.sonoma.edu http://zippy.sonoma.edu/kendrick/ September 29, 1996 - June 10, 1997 */ #include #incl

function [B,A] = adsgn(Fs); % ADSGN Design of a A-weighting filter. % [B,A] = ADSGN(Fs) designs a digital A-weighting filter for % sampling frequency Fs. Usage: Y = FILTER(B,A,X). % Warni

subplot(1,1,1) n = 0:10; x = (0.9*exp(j*pi/3)).^n; k = -200:200; w = (pi/100)*k; X = x * (exp(-j*pi/100)) .^ (n'*k); magX = abs(X); angX =angle(X); subplot(2,1,1); plot(w/pi,magX);grid axis([-2,

w = [0:1:500]*pi/500; % [0, pi] axis divided into 501 points. X = exp(j*w) ./ (exp(j*w) - 0.9*ones(1,501)); magX = abs(X); angX = angle(X); subplot(2,1,1); plot(w/pi,magX); grid; axis([0,1,0,10])

n = 0:100; x = cos(pi*n/2); k = -100:100; w = (pi/100)*k; % frequency between -pi and +pi X = x * (exp(-j*pi/100)).^(n'*k); % DTFT of x % y = exp(j*pi*n/4).*x; % signal mu

w = [-100:100]*pi/100; % frequency between -pi and +pi N = 5; X1 = sin((w+eps)*N/2)./sin((w+eps)/2); N = 15; X2 = sin((w+eps)*N/2)./sin((w+eps)/2); N = 25; X3 = sin((w+eps)*N/2)./sin((w+eps)/2);

w = [0:1:500]*pi/500; % [0, pi] axis divided into 501 points. X = exp(j*w) ./ (exp(j*w) - 0.5*ones(1,501)); magX = abs(X); angX = angle(X); realX = real(X); imagX = imag(X); subplot(2,2,1); plot(

% Chapter 3: Example 3.4: % DTFT of finite duration seq using mat-vec operation % n = -1:3; x = 1:5; % sequence x(n) k = 0:500; w = (pi/500)*k; % [0, pi] axis divided in