250 likes | 462 Views
OT3OS1. 26.10.2010. Diskretizacija i rekonstrukcija signala - realni model. Signal, trajanje 5s, frekvencija se menja od 20Hz do 10kHz. Diskretizacija, fs=20kHz. Diskretizacija, fs=10kHz, bez filtra. Diskretizacija, fs=10kHz, sa filtrom, fg=5kHz.
E N D
OT3OS1 26.10.2010.
A/D konverzija fs=100000; dt=1/fs; t=(0:dt:0.01)'; x=-1+200*t; x_16=round(x*2^15)/2^15; x_8=round(x*2^7)/2^7; x_4=round(x*2^3)/2^3; plot(t,x,t,x_16,'o',t,x_8,'*',t,x_4,'d'),xlabel('t'),ylabel('x(n\DeltaT),x_q(n\DeltaT)'), legend('x(t)','x_{16b}(n\DeltaT)','x_{8b}(n\DeltaT)','x_{4b}(n\DeltaT)'); s_g_16=sum((x-x_16).^2)/length(x); disp([s_g_16 (1/2^15)^2/12]); s_g_8=sum((x-x_8).^2)/length(x); disp([s_g_8 (1/2^7)^2/12]); s_g_4=sum((x-x_4).^2)/length(x); disp([s_g_4 (1/2^3)^2/12]); figure,plot(t,x-x_4,t,x-x_8,t,x-x_16),xlabel('t'),ylabel('x(n\DeltaT)-x_q(n\DeltaT)'), legend('x_{4b}(n\DeltaT)','x_{8b}(n\DeltaT)','x_{16b}(n\DeltaT)'); figure,hist(x); figure,hist(x-x_16); figure,hist(x-x_8); figure,hist(x-x_4);
A/D konverzija greška
A/D konverzija Snaga šuma kvantizacije s_g_16=sum((x-x_16).^2)/length(x) s_g_8=sum((x-x_8).^2)/length(x) s_g_4=sum((x-x_4).^2)/length(x) s_g_16 = 7.7528e-011 s_g_8 = 5.0809e-006 s_g_4 = 0.0013
Komplement dvojke, x[-1, 1) 0.75+0.75+0.75 +(-0.5)+(-0.5)+(-0.5) 011+ 011= 110 110+ 011= 1001 001+ 110= 111 111+ 110= 1101 101+ 110= 1011
1=xms~xmq=1 Karakteristika kvantizacije, opisana s xmq, podrazumeva da je ulazni signal u opsegu [-1 1) Amplituda ulaznog signala xmsje 1
0.01=xms<<xmq=1 Karakteristika kvantizacije, opisana s xmq, podrazumeva da je ulazni signal u opsegu [-1 1) Amplituda ulaznog signala xmsje 0.01
Signal greške U oba slučaja A/D konvertor je isti pa je max greška približno ista
Rezultati 1=xms~xmq=1 0.01=xms~xmq=1
Signal govora – uniformna kvantizacija Uniformna kvantizacija nije dobro rešenje
Signal govora – neuniformna kvantizacija Kompresija menja raspodelu signala