LC-метр из звуковой карты

DVK: чем больше гармоник, тем короче импульс и, соответственно, меньше энергия сигнала. Амплитуда, подразумевается, фиксирована и ограничена возможностями звуковой карты. Нельзя ли как-нибудь размазать эту энергию по времени, чтобы постоянно какая-нибудь колебалка шла?

rial_man: надо калибровку проводить чаще, желательно даже наверное раз в секунду.

rial_man: Беру случайный сигнал (белый шум)....,

rial_man: надо добавить в коробочку буферные усилители

>все довольно стабильно
1.ну и? насколько стабильно? за минуту на сколько уплывает? (мой расчет для моей звуковухи после 2 часов прогрева)
2. дифусилитель мне облом, надеюсь, что БПФ отфильтровывает 50 гц
3. емкости опорные тоже возможны, не только резисторы.

Первые 15мин КЗ убежало на ~ 30мОм, потом скачет +/-2.5мОм при усреднении 0.3сек, опорник 81Ом т.е.+/-30ппм от него, основная нестабильность вызвана нч шумом, если полосу снизу поставить 1кГц то болтанка получается меньше +/-1мОм (это на бортовом Реалтеке)

rial_man: . емкости опорные тоже возможны, не только резисторы.

Сейчас можно скомпенсировать Свх, Сл/р(эквивалентна прониканию между каналами)R и Lпроводов они задаются врукопашную и не полностью отражают частотные зависимости, но точность повышают прилично

Хочу попробовать в матлабе посмотреть на спектр и напряжения гармоник качающейся синусоиды, но не знаю как математикой сформировать такую синусоиду.
sin(2.0*pi*F*i) - это некачающаяся. За каждый отсчет дискредитации высчитывать прирост частоты и сразу F увеличивать? Правильно?

t=0:0.001:1;
sig=sin(2*pi*(5*t+20*t.*t));
plot(sig);

>точность повышают прилично
все ж надо б испытания хоть какие провести. на сколько врет, на каком диапазоне.

rial_man
Что вы там пишете? Откуда взялись 5 и 20 в ваших формулах? Как они связаны с минимальной и максимальной частотами качания?
Вот вырезка кода из моей программы, где я качаю от aMin до aMax гармоники (целые числа).
aA, aB - константы, зависящие от минимальной и максимальной гармоники качания,
aCount - количество точек выборки (размер буфера кратный 2 в целой степени, чтобы можно было БПФ использовать)
Pi - число Пи,
aWave - результирующий сигнал.
aA:= Pi/(aCount);
aB:= aA*(aMax-aMin)/aCount;
aA:= 2*aA*aMin;
for i:= 0 to aCount - 1 do begin
aWave[i]:= Sin(i*(aA + aB*i));
end;
Сгенерите такую синусоиду и возьмите от неё Фурье - посмотрите что будет.

Получается такое, точность измерения (вычисления) нечем проверить.
Качание за одну секунду.
-----------------------------------------------------------------------------------------------------
Качание частоты 20Гц - 20000Гц - V гармоник приблизительно в 136 раз меньше 1.
Качание частоты 200Гц - 20000Гц - V гармоник приблизительно в 135 раз меньше 1.
Качание частоты 1000Гц - 10000Гц - V гармоник приблизительно в 91 раз меньше 1.
Качание частоты 1000Гц - 2000Гц - V гармоник приблизительно в 30.5 раза меньше 1.
Качание частоты 1100Гц - 1600Гц - V гармоник приблизительно в 21.5 раза меньше 1.
------------------------------------------------------------------------------------------------------
В цикле матлабу нельзя синус задавать вычислять. Минутами результата приходится ждать. Нужно предварительно вычислить F и t, а потом
y =sin(2*pi*F.*t);
Позже надо у дельта импульса гармоники измерить.

Михалыч А
Ну вот так и должно быть - амплитуда падает как корень числа гармоник.
У вас ведь шаг по частоте где-то 1,08 Гц получается?

To DVK Я не могу ответить на счет 1.08Гц.
Вот мой матлабовский "вольметр" -

Fd = 192000; % частота дискредитации
Nbpf = 65536*2; % степень двойки
F= zeros (1,Nbpf) ; % качающаяся частота
t = 01/Fd):3;
t = t(1:Nbpf)';
F1 = 1100; % первая меньшая частота
F(1) = F1;
F2 = 1600; % вторая большая частота
dF = (F2 - F1)/Fd;
for i = 2:1:Nbpf
F(i) = F(i-1)+dF;
end
F = F';
y =sin(2*pi*F.*t);
Fft_Y=fft(y,Nbpf );
Pyy = 2*sqrt(Fft_Y.*conj(Fft_Y))/Nbpf ;
n=(1:1:Nbpf );
f = Fd *(n-1)/Nbpf ;
plot (f(1:20000),Pyy(1:20000))