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

loveeatmore

咨询一下关于理想带阻滤波器参数的问题，请大神帮忙啊
* Y0 `* |0 n2 x6 f& I
3 B" `9 R# ]4 Z( Wclc,clear,close all
warning off
feature jit off
9 C2 p; ]) F$ r/ pD0 = 20; % 阻止的频率点与频域中心的距离
W = 20;  % 带宽
' @9 y! q$ n* p/ |x = 0:.5:80;
y = 0:.5:80;
! X3 |, V- v4 M: u[X,Y] = meshgrid(x,y);
7 m4 `! v4 ]. n+ U9 a  S# U( Afor i=1:size(X,1)
    for j=1:size(X,2)
) c8 Z5 `5 d: h3 K" d        if sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2) <D0-W/2
            Z(i,j)=1;
, P( C/ Z  m1 G9 K, P/ n        elseif sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2)<D0+W/2 ...
. h3 l" \3 S: q8 n                && sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2)>=D0-W/2
; }" C* G$ `! S: z9 L% L            Z(i,j)=0;
& _3 V3 I. L+ K" v1 H9 I' A, Z# ?        elseif sqrt( (X(i,j)-40).^2 + (Y(i,j) -40).^2)>D0+W/2
            Z(i,j)=1;
        end
    end
& E0 ?, |8 z6 mend
, G- ]$ y5 l# x2 {figure('color',[1,1,1])
2 e' j; T/ b! X( emesh(X,Y,Z)

+ w2 a8 a3 t" V7 P/ g* {) \) S请问为什么采用X(i,j)-40,Y(i,j)-40?

regngfpcb

顶一下