Low SOI
Low SOI
wav');
N1=length(y1);
widmo1=fft(y1);
widmo1=abs(fft(y1));
widmo1(1)=widmo1(1)/N1;
widmo1(2:end)=widmo1(2:end)./(N1/2);
df=Fs/N1;
f=0:df:(Fs-df);
figure;
subplot(211);
stem(f,widmo1,'-*r')
N=200;
wn=9500/(Fs/2);
fo=[0 9500 10000 (Fs/2)]/(Fs/2);
Ho=[1 1 0 0];
b=firpm(N,fo,Ho);
fy=filter(b,1,y1)
a=1;
N2=length(fy);
widmo2=fft(fy);
widmo2=abs(fft(fy));
widmo2(1)=widmo2(1)/N2;
widmo2(2:end)=widmo2(2:end)./(N2/2);
df=Fs/N2;
f=0:df:(Fs-df);
subplot(212)
stem(f,widmo2,'-*r')
figure
freqz(b,a,512,Fs)
figure
subplot(211)
plot(y1)
subplot(212)
plot(fy)
sound(fy,Fs)