movievikings.m

% 
% VIKINGASKEPP -- THE _ANIMATED_ REVENGE
%
%    Huvudprogram, nu med animerade ror, vimpel och vatten.
%
%    verarbetning av laboration 3:10, Numeriska Metoder gkII.
%    Pontus Axelsson (d94-pax) & Fredrik Sandberg (d94-fsa).

clear,clf,format compact

global rpos botten nx ny wpos wneg % What can I say... they're all needed.

global oars         % En vektor med OBJECT HANDLES till de sex rorna.
global spunkter     % Vattnets skrningspunkter med skrovet.
global vimpel       % Ett OBJECT HANDLE till vimpeln.
noofframes=40;      % Totalt antal frames i hela animationen.

vikingplot;              % Plotta vikingaskepp utan ror.
water=watersurface;      % Berkna vattenlinjen.
waterplot(0,water,'b');  % Plotta vattnet i ursprungslge.
watervelvec=[];          % Initiera vattnets positionsvektor.
oarplot;                 % Plotta ror med konstanta x-vrden. Stt 'oars'.
newamplitude=0.5;        % Utgngsamplitud fr vimpeln.
vimpelplot(newamplitude);% Plotta vimpeln med amplitud 0.5.

% Utgngspositioner fr rorna.
for i=1.5:0.5:2.5
  rotate(oars(i/0.5-2),[0 0 1],-2*1/6*noofframes,[i 0.8 0.6-40*0.02])
  rotate(oars(i/0.5+1),[0 0 1],2*1/6*noofframes,[i -0.8 0.6-40*0.02])
end

disp('Frame no.:')

% Uppbyggnad av movie matrix M.
M=moviein(noofframes);
for i=1:noofframes
  title('Generating frames...')
  disp(i)

  delete(vimpel)
  newamplitude=newamplitude-sign(20-i-0.1)*1/20;
  vimpelplot(newamplitude);
  
  delete(wpos), delete(wneg)
%  watervelocity=(25-i)/24;
%  watervelvec=[watervelvec watervelocity*i/noofframes];

  if i<=5
    for j=1.5:0.5:2.5
      % Rotation kring x-axeln.
      rotate(oars(j/0.5-2),[1 0 0],-1,[j 0.8 0.6-40*0.02])
      rotate(oars(j/0.5+1),[1 0 0],1,[j -0.8 0.6-40*0.02])
      % Rotation kring y-axeln.
%      rotate(oars(j/0.5-2),[0 1 0],90/12.5,[j 0.8 0.6-40*0.02])
%      rotate(oars(j/0.5+1),[0 1 0],90/12.5,[j -0.8 0.6-40*0.02])
    end
    waterplot(i*2*pi/noofframes,water,'b');

  elseif i<=20
    for j=1.5:0.5:2.5
      % Rotation kring z-axeln.
      rotate(oars(j/0.5-2),[0 0 1],2,[j 0.8 0.6-40*0.02])
      rotate(oars(j/0.5+1),[0 0 1],-2,[j -0.8 0.6-40*0.02])
      % Rotation kring y-axeln.
%      rotate(oars(j/0.5-2),[0 1 0],sign(12.5-i)*90/12.5,[j 0.8 0.6-40*0.02])
%      rotate(oars(j/0.5+1),[0 1 0],sign(12.5-i)*90/12.5,[j -0.8 0.6-40*0.02])
    end
    waterplot(i*2*pi/noofframes,water,'b');

  elseif i<=25
    for j=1.5:0.5:2.5
      % Rotation kring x-axeln.
      rotate(oars(j/0.5-2),[1 0 0],1,[j 0.8 0.6-40*0.02])
      rotate(oars(j/0.5+1),[1 0 0],-1,[j -0.8 0.6-40*0.02])
      % Rotation kring y-axeln.
%      rotate(oars(j/0.5-2),[0 1 0],-90/12.5,[j 0.8 0.6-40*0.02])
%      rotate(oars(j/0.5+1),[0 1 0],-90/12.5,[j -0.8 0.6-40*0.02])
    end
    waterplot(i*2*pi/noofframes,water,'b');

  else
    for j=1.5:0.5:2.5
      % Rotation kring z-axeln.
      rotate(oars(j/0.5-2),[0 0 1],-2,[j 0.8 0.6-40*0.02])
      rotate(oars(j/0.5+1),[0 0 1],2,[j -0.8 0.6-40*0.02])
    end
    % waterplot((1-sum(watervelvec))/15*(i-25)*2*pi,water,'b');
    waterplot(i*2*pi/noofframes,water,'b');

  end

  M(:,i)=getframe;
end

% Uppspelning.
title('Roddbt DeLuxe av d94-fsa och d94-pax')

timestoplay=20;
movie(M,timestoplay);

% terinstllningar.
format
return