咨询一下关于理想带阻滤波器参数的问题，请大神帮忙啊

loveeatmore

咨询一下关于理想带阻滤波器参数的问题，请大神帮忙啊

clc,clear,close all
warning off
# k0 n: a/ g/ k1 z$ Kfeature jit off
D0 = 20; % 阻止的频率点与频域中心的距离
8 W; e* p: o9 ^0 eW = 20;  % 带宽
  |  m4 r! l9 M3 w9 ~) |' [0 ]* rx = 0:.5:80;
: T. n, r4 x2 V- Uy = 0:.5:80;
9 q' Y) ?. i/ q3 z[X,Y] = meshgrid(x,y);
for i=1:size(X,1)
    for j=1:size(X,2)
! C$ P) E$ P) A* d4 ?- u        if sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2) <D0-W/2
* I: _6 @& a/ M+ y) `            Z(i,j)=1;
        elseif sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2)<D0+W/2 ...
1 p1 g  x; H, |; O4 r                && sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2)>=D0-W/2
  P0 l; w+ m) E* F            Z(i,j)=0;
        elseif sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2)>D0+W/2
            Z(i,j)=1;
* t, v1 q" o1 W9 ~' }3 b" t8 H1 l; x        end
: m! O2 o: O% R/ n4 ~    end
/ j( Q) N! e. I; G6 a! fend
figure('color',[1,1,1])
mesh(X,Y,Z)
7 _* k* p; R9 q( o7 H0 `
请问为什么采用X(i,j)-40,Y(i,j)-40?
  V$ u. J$ h2 Y: v

regngfpcb

顶一下