передискретизация, как? Опции

[TigerSHARC]

Здравствуйте.
Есть сигнал:

Fs = 6400;
N = 128;
Ts = 1/Fs;
T = Ts*N;
f = 45;
t = 0:Ts:T-Ts;
A = 10;
Y = A*sin(2*pi*f*t)+ A*sin(2*pi*2*f*t)+ A*sin(2*pi*3*f*t);

как и какой функцией в MATLAB пересчитать частоту дискретизации с 6400 на 5760.
Причём с тем же числом отсчётов (128).

Это нужно для корректировки спектра сигнала.
Дело в том, что когда меняю частоту "руками":

Fs = 5760;
N = 128;
Ts = 1/Fs;
T = Ts*N;
f = 45;
t = 0:Ts:T-Ts;
A = 10;
df = Fs/N;
f1 = 0:df:Fs-(Fs/N);
Y = A*sin(2*pi*f*t)+ A*sin(2*pi*2*f*t)+ A*sin(2*pi*3*f*t);
X = fft(Y)/length(t);
>> stem(f1, abs(X)), grid

получается красивый спектр, отражающий верную информацию об амплитудах сигнала.
Пробовал для передискреизации функцию resample, но она выдаёт меньшее число отсчётов и спектр искажён...

да и вообще, не толкьо в MATLAB, как из 128 отсчётов получается 128 передискретизированных отсчётов... полюбому время выборки должно как-то варьироваться...

[TigerSHARC]

с точки зрения такой науки как численные методы возможно провести интерполяцию сигнала.

P.S. не использую оконное взвешиване здесь так как ГОСТ это запрещает (для моего случая) вот и приходится частоту менять программно.

Сообщение отредактировал TigerSHARC - Feb 22 2010, 17:15

[bahurin]

P.S. не использую оконное взвешиване здесь так как ГОСТ это запрещает (для моего случая)

В штате Индиана (США) в 1897 году был выпущен билль (см.: en:Indiana Pi Bill), законодательно устанавливающий значение числа Пи равным 3,2.

Вобще странная ситуация, когда гост запрещает чтото использовать. Почему бы не использовать окна, если они позволяют получить установленные гостом харакетристики? Чудно однако.