代码搜索结果

% 第七章: 例 7.16 % 频率采样技术: 低通, 最优法T1 & T2 % wp = 0.2pi, ws=0.3pi, Rp=0.25dB, As=50dB % T1 = 0.5925, T2=0.1099 M = 60; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l; Hrs = [ones(1

% 第七章 例 7.25 % 用PM算法进行的高通滤波器设计 % 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/20)); deltaH = max(delta1,delta2); deltaL = m

% 第七章: 例 7.19 % 频率采样技术: 差分器 % 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*Dw*k1,-j*Dw*(M-k2)]; angH = [-alpha*Dw*k1, a

%第三章: 例3.4: % 用矩阵-向量乘法求有限长序列的DTFT % n = -1:3; x = 1:5; % x(n)序列 k = 0:500; w = (pi/500)*k; % [0, pi] 轴分为501点. X = x * (exp(-j*pi/500)) .^ (n'*k); % 用矩阵-向量乘法求DTFT magX

% 第三章: 例 3.6: % subplot(1,1,1) n = -5:5; x = (-0.9).^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 ax

% 第三章: 例 3.12: % n = -5:10; x = sin(pi*n/2); k = -100:100; w = (pi/100)*k; % -pi和+pi之间的频率 X = x * (exp(-j*pi/100)).^(n'*k); % x的DTFT % 信号分解 [xe,xo,m] = evenodd(x,n); % 偶部和奇部

% 第三章: 例 3.9: % n = 0:100; x = cos(pi*n/2); k = -100:100; w = (pi/100)*k; % 把 -pi 到 +pi 分为201点 X = x * (exp(-j*pi/100)).^(n'*k); % x 的 DTFT % y = exp(j*pi*n/4).*x; % 把信

% 第三章: 例 3.3: % w = [0:1:500]*pi/500; % 把[0, pi] 轴分为501点. 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(w/

% 第三章: 例 3.99: % w = [-100:100]*pi/100; % -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+ep

% 第三章: 例 3.5: % 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,mag