создание движущихся снежинок)

if (~isempty(daqfind)) 
    stop(daqfind) % якщо запущено програми збору данних, припиняємо їхню роботу 
end
AO = analogoutput('winsound', 0); % присвоюємо змінній аналоговий вихід
addchannel(AO, [1 2]); % додаемо функціонування в режимі стерео
load handel  %завантажуємо приклад аудіофайлу
% y- дані які будуть виводитися через звукову карту
% Fs - задавана частота дискретизації данних для оригінального аудіофайлу
set(AO, 'StandardSampleRates','Off') %виключаємо оригінальну частоту дискретизації
set(AO, 'SampleRate', Fs); % підставляємо значення задаваної Fs частоти дискретизіції
Fs
data = [y y];% завантажені на вихід дані повинні завантажуватись послідовно командою збору данних PUTDATA
putdata(AO, data); % дані які йдуть на вихід повинні містити стовпчик відповідних значень для кожного каналу аналогового виходу
startindex = 1; % стартове значення
increment = 500; % максимальне значення( значення збільшення)
start(AO); % запуск аудіо файлу
 
while isrunning(AO)  % цикл для графічного відображення зміни... в реальному часі
    while (AO.SamplesOutput < startindex + increment -1), end
    try
        x = AO.SamplesOutput;
        plot(y(x:x+increment-1));
        set(gca, 'YLim', [-0.8 0.8], 'XLim',[1 increment])
        drawnow;
        startindex =  startindex+increment;
    end
end
% delete(ao); % видалення данних з пам'яті

clear, clc
fps = 20; % кадров в сек.
dt = 1/fps; % шаг времени
T = 10; % длительность
t = 0:dt:T; % массив времени
n = 6; % кол-во снежинок
% начальные точки
X0 = 3*randn(1,n);
Y0 = 0*ones(1,n);
% радиусы траекторий
R = randint(1,n,[1 5]);
a = 3*rand(1,n); % скорость опускания
for i = 1:length(t)
   X = cos( t(1:i) )'*R + ones(i,1)*X0;
   Y = abs(sin( t(1:i)' ))*R + ones(i,1)*Y0 - ones(i,1)*a*t(i);
   plot(X,Y,'*')
   F(i) = getframe;
end
%movie(F,fps)

clear, clc
fps = 20; % кадров в сек.
dt = 1/fps; % шаг времени
T = 10; % длительность
t = 0:dt:T; % массив времени
n = 6; % кол-во снежинок
% начальные точки
X0 = 3*randn(1,n);
Y0 = 0*ones(1,n);
% радиусы траекторий
R = randint(1,n,[1 5]);
a = 3*rand(1,n)+0.5; % скорость опускания
w = 2*pi*(rand(1,n)+1);
for i = 1:length(t)
   X = X0 + cos( w*t(i) ).*R;
   Y = Y0+R.*abs(sin(w*t(i))) - a*t(i);
   plot(X,Y,'*')
   F(i) = getframe;
end
% movie(F,fps)

clear, clc
fps = 20; % кадров в сек.
dt = 1/fps; % шаг времени
T = 10; % длительность
t = 0:dt:T; % массив времени
n = 10; % кол-во снежинок
% начальные точки
X0 = 3*randn(1,n);
Y0 = 0*ones(1,n);
% радиусы траекторий
R = randint(1,n,[1 5]);
a = 3*rand(1,n)+0.5; % скорость опускания
w = 2*pi*rand(1,n);
for i = 1:length(t)
   X = X0 + cos( w*t(i) ).*R;
   Y = Y0-R.*abs(sin(w*t(i))) - a*t(i);
   plot(X,Y,'*')
   axis([-10 10 -35 0])
   pause(dt)
   %F(i) = getframe;
end