суббота, 25 октября 2008 г.

пятница, 17 октября 2008 г.

суббота, 11 октября 2008 г.

пятница, 10 октября 2008 г.

Метод Быкова в Maxima (вроде работает).

После плясания с бубном итог:
http://stream.ifolder.ru/8504762
Качаете Maxima
Запускаете wxMaxima, и открываете byk.wxm
Забыл в файле написать комментарий, что после того как набрали перевычислить ввод, надо в пункте главного меню Maxima нажать Reevaluate all (или Shift+Ctr+R).

среда, 1 октября 2008 г.

Пример использования autocorr















function y = noise
Fs = 4000; % Sampling frequency ???????????? ??????????? ?? 2
T = 1/Fs; % Sample time
L = 1096; % Length of signal
t = (0:L-1)*T; % Time vector

% Все параметры прикидочные НЕ ПО ЗАДАНИЮ
sigma = 1
w = 1;

% Подбором задается в зависимости от ширины спектра шума
gamma = 0.5;

t = 0:1:50;

% Коррелируются 4 отсчета
k = 0:1:3;

% Весовые коэфициенты №8 в таблице Быкова стр 105.
c = (2 * sigma) .* sqrt(gamma ./ pi) .* (1./(1 + 4 .* (gamma.^2) .* k.^2));


%a = (2.*sigma.*sqrt(gamma^3));
%c = a .*((1 - ( - 1) .^ abs(k) .* exp(-gamma/2))./(gamma^2 + 4*pi^2.*k.^2));

%dt = t(4) - t(3)

% Аналитическая автокоррел. функция
R = sigma ./ (1 + w.^2 .* t.^2);
%R = R/50;

% generates white noise
x = wgn(L,1,0);

% на стр 56 Быкова формула очень похожа на ту, котороя в MATLAB Function Reference
y = conv(c, x);
%autocorr(y, [], 2);
[ACF, Lags, Bounds] = autocorr(y, [], 2);
[ACF, Lags]

%[Lags, R] = eqtflength(ACF,R);

subplot(2, 1, 1)
plot(y);

subplot(2, 1, 2);

plot(Lags, ACF, t, R)
%axis([0 1 0 1]);