Matlab - Grafico animado en matlab

D=0.05;

M=20;

N=20;

m=M+1;

n=N;

alf=5;

bet=4;

figure(1)

T=1;% aqui regula o tempo da figura animada

dtt=0.0125;

tt=0;

while tt<T

    tt=tt+dtt; % evolucion del tiempo

    for itr=1:m-1

    for i=1:n

    uu(itr,i)=(2*D*bet*pi*sin(pi*x(itr))*exp(-D*pi*pi*t(i)))./(alf+bet*cos(pi*x(itr))*exp(-D*pi*pi*t(i)));

    end

end

x=xl+(0:m-1)*h;t=tb+(0:n)*k;

    surfc(x, t, uu'), colorbar; shg

    caxis([-0.15 0.25])   % Pseudocolor axis scaling.

     title([' Simulation results after t = '  num2str(tt) '  [s]'])

     pause(dtt)

     axis tight

end

 xlabel('Width, [m] ')

 ylabel('Length, [m] ')

zlabel('Displacement, [m]')

% %% Contour plot

 figure(2)

 contourf(x, t, uu'), colorbar; shg

 title([' Simulation results after t = '  num2str(t) '  [s]'])

 axis tight

 xlabel('Width, [m] ')

 ylabel('Length, [m] ')

La variable  x donde se define.

Undefined function or variable 'x'

xl=0;

xr=1;

tb=0;

te=1;

M=20;

N=20;

m=M+1;

n=N;

h=(xr-xl)/M;

k=(te-tb)/N;

x=xl+(0:m-1)*h;t=tb+(0:n)*k;

>> animado1

  Name       Size            Bytes  Class     Attributes

  t          1x21              168  double

  uu        20x20             3200  double

  x          1x21              168  double

xl=0;

xr=1;

te=1;

tb=0;

M=20;

N=20;

alf=5;bet=4;D=0.05;

h=(xr-xl)/M; k=(te-tb)/N; m=M+1; n=N;

x=xl+(0:m-1)*h;t=tb+(0:n)*k;

T=1;% aqui regula o tempo da figura animada

dtt=0.0125;

tt=0;

uu=zeros(m,m);

while tt<T

    % condicion de frontera:

    uu(:, [1,end])  = 0;

    uu([1, end], :) = 0;

    %  evolucion del tiempo:

    tt=tt+dtt;

    for itr=2:m-1

    for i=2:n-1

    uu(itr,i)=(2*D*bet*pi*sin(pi*x(itr))*exp(-D*pi*pi*t(i)))./(alf+bet*cos(pi*x(itr))*exp(-D*pi*pi*t(i)));

    end

end

    surfc(x, t, uu'), colorbar; shg

    caxis([-0.15 0.25])   % Pseudocolor axis scaling.

     title([' Simulation results after t = '  num2str(tt) '  [s]'])

     pause(dtt)

     axis tight

end

for itr=2:m-1

for i=2:n-1

uu(itr,i)=(2*D*bet*pi*sin(pi*x(itr))*exp(-D*pi*pi*t(i)))./(alf+bet*cos(pi*x(itr))*exp(-D*pi*pi*t(i)));

end

surfc(x, t, uu'), colorbar; shg

end

clear

close all

xl=0;

xr=1;

te=1;

tb=0;

M=20;

N=20;

alf=5;bet=4;D=0.05;

h=(xr-xl)/M; k=(te-tb)/N; m=M+1; n=N;

x=xl+(0:m-1)*h;t=tb+(0:n)*k;

T=1;% aqui regula o tempo da figura animada

dtt=0.0125;

tt=0;

uu=zeros(m,m);

% condicion de frontera:

uu(:, [1,end]) = 0;

uu([1, end], :) = 0;

for itr=2:m-1

    for i=2:n-1

    uu(itr,i)=(2*D*bet*pi*sin(pi*x(itr))*exp(-D*pi*pi*t(i)))./(alf+bet*cos(pi*x(itr))*exp(-D*pi*pi*t(i)));

    end

    axis([0 1 0 1  0 0.5])

    tt=tt+dtt;

    surf(x, t, uu)

    view(86,14)

    axis([0 1 0 1  0 0.5])

    t1=sprintf('%8.6f',tt);

    title([' Simulation results after t = ', t1(1:6),'[s]'])

    pause(0.1)

end

>> animado1