MATLAB ADC 频谱分析问题

needforspeed

一段分析ADC频谱的MATLAB程序，fft输出很奇怪，出现了-inf，请大家看看有什么问题。多谢！
clear;

afs = 1;   % 1V full scale
fs = 1e6;    %  1MHz sample rate
N = 2^10;   % number of samples
cyc = 67;
fx = fs*cyc/N;

t = linspace(0,(N-1)/fs,N);

y = afs*sin(2*pi*fx*t);

B = 10;                        % internal ADC resolution
delta = afs/2^B;
y = cos(2*pi*fx/fs*[0:N-1]);
%quantize samples to delta=1LSB
y=round(y/delta)*delta;
s = 20*log10(abs(fft(y)/N/afs*2));
s = s(1:N/2);  % drop redundant half

f = (0:length(s)-1) / N;  % frequency vector (normalized to fs)

figure(1);
plot(f,s);
grid on;

victor0o0

其实就是计算fft的时候，matlab对数据存储的问题，你可以比较一下下面s0(理想正弦),s1(ADC 输出),sd(量化误差)。理论上s1=s0+sd，实际计算中s1有些数据丢失了。
clear;
afs = 1;   % 1V full scale
fs = 1e6;    %  1MHz sample rate
N = 2^10;   % number of samples
cyc = 67;
fx = fs*cyc/N;
t = linspace(0,(N-1)/fs,N);
B = 10;                        % internal ADC resolution
delta = afs/2^B;
y0 = cos(2*pi*fx/fs*[0:N-1]);
y1 =round(y0/delta)*delta;

s0 = 20*log10(abs(fft(y0)/N/afs*2));
s1 = 20*log10(abs(fft(y1)/N/afs*2));
sd = 20*log10(abs(fft(y1-y0)/N/afs*2));

plot([s, s1, sd])

a2744200000

回复 1# needforspeed


   顶了！下来看看。

needforspeed

回复 2# victor0o0


    我理解你说的，量化之后有一定的误差是正常的。但这是ADC设计和分析的基本问题吧，如何对量化之后的信号做比较准确的频谱分析？

victor0o0

回复  victor0o0


    我理解你说的。但如何对量化之后的信号做准确的频谱分析？
needforspeed 发表于 2015-1-16 09:43

                               
登录/注册后可看大图

   1. 做出的fft结果只是小量被忽略了，不会影响你计算SNR，SFDR之类。   2. 这种结果实际不会出现，随便加点非理想因素，都不会出现fft里有0,比如在输入信号里叠加一点随机噪声y=y+randn(size(y))*2^-16
   3. 如果你非得看那些小量，先把信号成分去了，再做fft

auch0311

谢谢