исходный набор сэмплов децимицимируется до 44.1 KHz ( из 3.5MHz ) и записывается просто в wav ( важно что 16бит на канал )
однако появляется некий шум - вот и хочется разобраться в причине этого шума - исходные данные ( врядли ) или что то еще ?
пример my.m и test.m вместе с исходными данными ( осторожно - большой файл ) и получаемым .wav приатачен...
Полная версия этой страницы: fdesign.decimator и шум
Форум разработчиков электроники ELECTRONIX.ru > Cистемный уровень проектирования > Вопросы системного уровня проектирования
Цитата(DevL @ Dec 11 2011, 23:13) 
исходный набор сэмплов децимицимируется до 44.1 KHz ( из 3.5MHz ) и записывается просто в wav ( важно что 16бит на канал )
однако появляется некий шум - вот и хочется разобраться в причине этого шума - исходные данные ( врядли ) или что то еще ?
пример my.m и test.m вместе с исходными данными ( осторожно - большой файл ) и получаемым .wav приатачен...
однако появляется некий шум - вот и хочется разобраться в причине этого шума - исходные данные ( врядли ) или что то еще ?
пример my.m и test.m вместе с исходными данными ( осторожно - большой файл ) и получаемым .wav приатачен...
Ну а где собственно сам CIC фильтр? Я вижу только компенсирующий фильтр для CIC-а.
Или Вы думали, что h = fdesign.decimator(79, 'CIC Compensator',...) сделает еще и сам CIC? Увы, нет
filterbuilder генерирует очень схожее и для "чистого" CIC, приатаченно ниже
а так же вот кусок, который "с шумом тоже"
а так же вот кусок, который "с шумом тоже"
Код
function Hd = test_cic
%TEST_CIC Returns a discrete-time filter object.
%
% MATLAB Code
% Generated by MATLAB(R) 7.13 and the DSP System Toolbox 8.1.
%
% Generated on: 12-Dec-2011 16:01:07
%
D = 4; % Differential delay
Fpass = 1; % Passband Frequency
Astop = 1; % Stopband Attenuation (dB)
h = fdesign.decimator(79, 'CIC', D, 'Fp,Ast', Fpass, Astop);
Hd = design(h, 'multisection');
set(Hd, 'Arithmetic', 'fixed', ...
'InputWordLength', 16, ...
'InputFracLength', 15, ...
'FilterInternals', 'Specifyprecision', ...
'SectionWordLengths', 18, ...
'SectionFracLengths', 15, ...
'OutputWordLength', 16, ...
'OutputFracLength', 15);
% [EOF]
%TEST_CIC Returns a discrete-time filter object.
%
% MATLAB Code
% Generated by MATLAB(R) 7.13 and the DSP System Toolbox 8.1.
%
% Generated on: 12-Dec-2011 16:01:07
%
D = 4; % Differential delay
Fpass = 1; % Passband Frequency
Astop = 1; % Stopband Attenuation (dB)
h = fdesign.decimator(79, 'CIC', D, 'Fp,Ast', Fpass, Astop);
Hd = design(h, 'multisection');
set(Hd, 'Arithmetic', 'fixed', ...
'InputWordLength', 16, ...
'InputFracLength', 15, ...
'FilterInternals', 'Specifyprecision', ...
'SectionWordLengths', 18, ...
'SectionFracLengths', 15, ...
'OutputWordLength', 16, ...
'OutputFracLength', 15);
% [EOF]
Цитата(DevL @ Dec 11 2011, 23:13)
однако появляется некий шум - вот и хочется разобраться в причине этого шума - исходные данные ( врядли ) или что то еще ?
Мельком глянул входные данные, это что, ΔΣ АЦП? Странный он у Вас какой-то. И на кой ляд делать 3.48МГц, если данные с АЦП все равно меняются с частотой (примерно) 23кГц???? И для начала неплохо бы тестировать на банальной синусоиде, постепенно усложняя сигнал.
А шум...ну тут все просто, это шумит все то, что не задавилось фильтром и завернулось в полосу 22кГц, и завернуться там есть чему при таком исходном сигнале
.
все это - только промежуточная обработка,
к сожалению нельзя поменять входной сигнал и его характеристики, даже синусоида будет "непросто"
вот в этом и вся сложность
какие рекомендации насчет " задавилось фильтром и завернулось в полосу 22кГц" ?
к сожалению нельзя поменять входной сигнал и его характеристики, даже синусоида будет "непросто"
вот в этом и вся сложность
какие рекомендации насчет " задавилось фильтром и завернулось в полосу 22кГц" ?
Для просмотра полной версии этой страницы, пожалуйста, пройдите по ссылке.