?? exa030201.m
字號:
%----------------------------------------------------------------------------
% exa030201, for example 3.2.1 and fig 3.2.3;
% to show Gibbs phenomena;
%----------------------------------------------------------------------------
wc=0.6*pi;a=[1];Nf=256;
Hd1=ones(1,0.6*Nf+1);
Hd2=zeros(1,0.4*Nf);
Hd=[Hd1,Hd2]';
% 改變窗函數的長度,觀察Gibbs現象
%
M=10; % for M=10
N=2*M+1;
for n=1:N
if (n-M)==0;
h(n)=wc/pi;
else
h(n)=sin(wc*(n-M))/(pi*(n-M));
end
end
[H,w]=freqz(h,a,Nf,1);
subplot(221)
plot(w,Hd,'-',w,abs(H));grid on;
%
M=30; % for M=30
N=2*M+1;
for n=1:N
if (n-M)==0;
h(n)=wc/pi;
else
h(n)=sin(wc*(n-M))/(pi*(n-M));
end
end
[H,w]=freqz(h,a,Nf,1);
subplot(222)
%
M=50; %for M=50
N=2*M+1;
for n=1:N
if (n-M)==0;
h(n)=wc/pi;
else
h(n)=sin(wc*(n-M))/(pi*(n-M));
end
end
[H,w]=freqz(h,a,Nf,1);
subplot(223)
plot(w,Hd,'-',w,abs(H));grid on;
%
M=100; % for M=100
N=2*M+1;
for n=1:N
if (n-M)==0;
h(n)=wc/pi;
else
h(n)=sin(wc*(n-M))/(pi*(n-M));
end
end
[H,w]=freqz(h,a,Nf,1);
subplot(224)
plot(w,Hd,'-',w,abs(H));grid on;
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -