Witam,
Chciałem zrobić filtr K-krotnie interpolujący o długości N=2M+1=31 i zastosować
w nim okno prostokątne.
Zrobiłem coś takiego
clear all;
M = 15; %połowa długości filtra (N=2*M+1=31)
typ=1;
N=2*M+1; n= 1 : M; m = -M : 1 :M;
if (typ==1)
K=5;
fc=1/(2*K);
wc=2*pi*fc;
h = 2*fc*sin(wc*m)./(wc*m);
end
w = boxcar(N); %okno prostokątne
w=w';
hw = h .* w; %wymnożenie odpowiedzi impulsowej z oknem
NF=500; %obliczanie widma Fouriera
fn=0.5*(1:NF)/(NF+1);
for k=1:NF
H(k)=sum(h .* exp(-j*pi*fn(k)*m) );
HW(k)=sum(hw .* exp(-j*2*pi*fn(k)*m) );
end
chciałem wyświetlić do tego charakterystyki amplitudową i fazaową. Jednak nie wiem jak. Jak wykorzystuje funkcję fvtool to pojawia się komunikat że parametr musi być skończony. Nie wiem w jaki inny sposób mam je uzyskać. Proszę o pomoc,
Dzięki
Chciałem zrobić filtr K-krotnie interpolujący o długości N=2M+1=31 i zastosować
w nim okno prostokątne.
Zrobiłem coś takiego
clear all;
M = 15; %połowa długości filtra (N=2*M+1=31)
typ=1;
N=2*M+1; n= 1 : M; m = -M : 1 :M;
if (typ==1)
K=5;
fc=1/(2*K);
wc=2*pi*fc;
h = 2*fc*sin(wc*m)./(wc*m);
end
w = boxcar(N); %okno prostokątne
w=w';
hw = h .* w; %wymnożenie odpowiedzi impulsowej z oknem
NF=500; %obliczanie widma Fouriera
fn=0.5*(1:NF)/(NF+1);
for k=1:NF
H(k)=sum(h .* exp(-j*pi*fn(k)*m) );
HW(k)=sum(hw .* exp(-j*2*pi*fn(k)*m) );
end
chciałem wyświetlić do tego charakterystyki amplitudową i fazaową. Jednak nie wiem jak. Jak wykorzystuje funkcję fvtool to pojawia się komunikat że parametr musi być skończony. Nie wiem w jaki inny sposób mam je uzyskać. Proszę o pomoc,
Dzięki