Типы данных в Matlab

@Антов Феоктистов · Регистрация: 14.07.2009

Author24 — интернет-сервис помощи студентам

Здравствуйте,

пишу программу в Matlab, которая активно использует свертки методом overlap and add. Проблема вся в том, что нужно свернуть последовательность из 2*10^8 с последовательностью 5*10^6. За основу был взят алгоритм fftfilt и слегка подрихтован под свои задачи. Но все равно остается одна большая проблема - слишком долго выполняется и слишком много памяти требуется. Я пробовал перевести данные в формат int8, но светрка, написанная через преобразование Фурье, напрочь отказывается считать это преобразование в int8.

Тогда я попытался хотя бы понять как правильно работают типы данных в Матлабе и будет ли быстрее работать с интом, нежели с родным даблом. Запускал вот этот скрипт:

Matlab M

y = 0;
x = 0;
 
for h=1:20
 
    N = 1000*h;
 
    a = randint(1,N,[0, 10]);
    c = randint(1,N,[0, 10]);
    aa = int8(a);
    cc = int8(c);
 
    b = rand(1,N);
    bb = double(b);
    e = rand(1,N).*10;
    ee = double(e);
 
    sumint = 0;
    sumint = sumint
    sumfloat = 0;
    timeintsum = 0;
    timefloatsum = 0;
 
    for i=1:N
        tic
            sumint = sumint + aa(i)*cc(i);
        timeint = toc;
        timeintsum = timeintsum + timeint;
    sumint = sumint
       
        tic
            sumfloat = sumfloat + bb(i)*ee(i);
        timefloat = toc;
        timefloatsum = timefloatsum + timefloat;
    end;
 
    y = [y timeintsum];
    x = [x timefloatsum];
   
end;
 
figure
hold on;
plot(1:length(y),y,'r');
plot(1:length(x),x,'b');
grid;

Увы, результаты, мягко сказать шокировали - во-первых, так просто оно не работает, и sumint = sumint + aa(i)*cc(i); не преобразует автоматически инт8 в инт16\32. Поэтому, меняем ее на sumint = int32(sumint) + int32(aa(i))*int32(cc(i));

Но и это дает проблемы - в десять раз длиннее вычисления :\ в чем подвох? В переводе типов? возможно ли уменьшить количество памяти и сделать программу быстрее, если использовать целочисленные типы данных вместо плавающей точки?

p.s. почитав литературу по Си было сказано - испольуйте те типы, которые приведут вас к меньшей погрешности. Затрата времени на преобразование типов чаще больше, нежели на перевод их из одного типа в другой. так ли это?

Заранее благодарен

	16.09.2010, 12:36