![]() = Īn easy way to get the same answer in scilab is to use their iir filter design function, which does both the analog design and bilinear transformation in one step, like so: -> hz=iir(N,'lp','butt',5./15.,)Ġ.3318051 +0.9954154z +0.9954154z² +0.3318051z³ It looks like you've missed a sign inversion when converting the output from that website to canonical form. ![]() Xtitle('The ''filter'' filtered gaussian noise','t','y') Xtitle('The ''csim'' filtered gaussian noise','t','y') Y_res = filter(Cnum, Cden, Input) // Filter the signal with filter Y_csim = csim(Input,t,elatf) // Filter the signal with csim T = t*0.01 // Convert sample index into time steps Input = rand(1,1000) // Produce a random gaussian noise ///////////////// plot an exemple to compare csim and filter Generate the equivalent linear system of the filterĭisp('coefficients : Num / Den : ',Cnum,Cden) conclusion : zpbutt et analpf donnent la même sortie disp("Gain : Zpbutt ",gainZP, "Analpf ",gainAna) disp("Pole : Zpbutt ",poleZP, "Analpf ",poleAna) ![]() ///////////////// compute different functions to compaire Butterworth //////////////////// variable declarationįcut = 5 //cut off frequency hz (delta 1)įsampl = 15 //sampling frequency hz (delta 2)ĭelta1_in_dB = -3 // attenuation value at fcut
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |