clc;clear all;
Fp = 0.4;
Fst = 0.5;
Fs = 2;
R = (10^(0.02/20) - 1)/(10^(0.02/20) + 1); %Ripple
A = 10^(-93/20); %attenuation
Filt = firgr('minorder',[0 Fp/(Fs/2) Fst/(Fs/2) 1/(Fs/2)], [1 1 0 0],[R A]);
if (mod(length(Filt),2) == 0)
N = length(Filt);
else
N = length(Filt)+1;
end
Filt = firgr(N,[0 Fp/(Fs/2) Fst/(Fs/2) 1/(Fs/2)], [1 1 0 0],[R A]);
Filt = Filt.*blackman(length(Filt))';
figure;zplane(Filt);title('Pole-Zero Plot');
[h,w]=freqz(Filt);
H = abs(h);
H = 20*log10(H/max(H));
figure;
subplot(2,1,1);plot(w*(Fs/(2*pi)),H,'--k','Linewidth',2);
xlabel('Frequency(Hz)');ylabel('Magnitude');
title('1D Filter Magnitude Response');
subplot(2,1,2);plot(w*(Fs/(2*pi)),angle(h),'--k','Linewidth',2);
xlabel('Frequency(Hz)');ylabel('Angle (radians)');
title('1D Filter Phase Response');
H3 = ftrans2(Filt);
figure;freqz2(H3);
title('2D Filter Response');

clc;clear all;
Fs1 = 0.4;
Fp1 = 0.45;
Fp2 = 0.6;
Fs2 = 0.65;
Fs = 2;
R1 = (10^(0.04/20) - 1)/(10^(0.04/20) + 1); %Ripple1
R2 = (10^(0.04/20) - 1)/(10^(0.04/20) + 1); %Ripple2
Filt = firgr('minorder',[0,Fs1/(Fs/2),Fp1/(Fs/2),Fp2/(Fs/2),Fs2/(Fs/2),1/(Fs/2)],[0 0 1 1 0 0],[R2 R1 R2]);
if (mod(length(Filt),2) == 0)
N = length(Filt);
else
N = length(Filt)+1;
end
Filt = firgr(N,[0,Fs1/(Fs/2),Fp1/(Fs/2),Fp2/(Fs/2),Fs2/(Fs/2),1/(Fs/2)],[0 0 1 1 0 0],[R2 R1 R2]);
Filt = Filt.*hamming(length(Filt))';
figure;zplane(Filt);title('Pole-Zero Plot');
[h,w]=freqz(Filt);
H = abs(h);
H = 20*log10(H/max(H));
figure

1. Fields and Potentials in MATLAB Use the table 6.2, p.373, &quo...