Forum studentów IV roku EiT Politechniki Śląskiej

michasia

Mógłby ktoś wrzucić to co mieliśmy zrobić na laborce przed świętami? chodzi o grupę E1

Stefan99

Lesiu chyba gdzies wrzucil ten programik.

jc

programik nadesłany przez kumpla

clear,clf,clc;
fs=10000;
fp1=2*1000/fs;
fz1=2*1500/fs;
fz2=2*3000/fs;
fp2=2*3500/fs;
ap=3;
az=50;
Rp=3;
Rs=50;
%[N,fn]=ellipord(fp1,fz1,ap,az);
%[B,A]=ellip(N,Rp,Rs,fn,'low');
%[N,fn]=ellipord([fp1 fp2],[fz1 fz2],ap,az);
%[B,A]=ellip(N,Rp,Rs,fn,'stop');
%[N,fn]=buttord([fp1 fp2],[fz1 fz2],ap,az);
%[B,A]=butter(N,fn,'stop');
%[B,A]=butter(N,fn);
%[N,fn]=cheb1ord([fp1 fp2],[fz1 fz2],ap,az);
%[B,A]=cheby1(N,Rp,fn,'bandpass');
[N,fn]=cheb2ord([fp1 fp2],[fz1 fz2],ap,az);
[B,A]=cheby2(N,Rs,fn,'stop');
dirac=[1 zeros(1,100000)];
h=filter(B,A,dirac);
test=[ones(1,20)-ones(1,20) ones(1,20) -ones(1,20) ones(1,20) -ones(1,20)];
out=filter(B,A,test);
widmo=abs(fft(h,2^20));
widmo=20*log10(widmo);
subplot(311),plot(test);
subplot(312),plot(out);
%axis([0 100 -1 1]);
[z s]=size(widmo);
f=(1:s);
f=fs*f./s;
subplot(313),plot(f,widmo);grid;
axis([0 4000 -60 10]);
hold on;
a=-3*ones(1,10000);
b=-50*ones(1,10000);
plot(a,'r');
hold on;
plot(b,'r');

GREY

Tutaj troche przejrzysciej rozpisane formulki dla Buttera:

rudolf

Co to jest az i ap, a co to jest Rz i Rp? nie moze byc po prostu jedno lub drugie?
i jeszcze dlaczego przy normalizowaniu czestotliwosci wzgledem fs mnozymy wszystko przez 2?

GREY

az i Rz to to samo, ap i Rp tak samo. Po prostu sa to oznaczenia zaczerpniete z opisow funkcji. A mnozymy razy 2 gdyz wszystkie funkcje przyjmuja wspolczynniki znormalizowane wzgledem polowy Fs.