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

loveeatmore

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

clc,clear,close all
warning off
1 ^( N5 |. w7 h" y) p+ m$ efeature jit off
" R, y: m; m; N. \" DD0 = 20; % 阻止的频率点与频域中心的距离
$ w& R1 ?+ j+ A  s4 P& dW = 20;  % 带宽
x = 0:.5:80;
# S6 H' T% h+ }+ Z/ A9 Fy = 0:.5:80;
" f2 t2 ^/ C# h( z3 I% T) i8 b) t[X,Y] = meshgrid(x,y);
for i=1:size(X,1)
    for j=1:size(X,2)
+ t9 V5 ~$ D/ _        if sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2) <D0-W/2
( A0 N( E! C+ Q$ W$ s6 ?& L            Z(i,j)=1;
        elseif sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2)<D0+W/2 ...
                && sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2)>=D0-W/2
- l- v5 a; W/ i2 ]1 u            Z(i,j)=0;
        elseif sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2)>D0+W/2
            Z(i,j)=1;
        end
# ~- |4 K+ x) ~- V2 U  \: a    end
end
figure('color',[1,1,1])
mesh(X,Y,Z)
% q' D3 ^3 A2 S3 p0 I
8 k+ }" v  ^7 n! s1 Q2 _请问为什么采用X(i,j)-40,Y(i,j)-40?
$ H0 o0 [0 ~. D$ k$ Q3 Y' L) u


$ t& D- v9 i0 x9 f. X# r( E

regngfpcb

顶一下