MATLAB: How to draw evolution of concentration in 2D [solved]

concentrationplot2dshade

Hi all,

I have coded a small program in order to compute the evolution of concentration of a plume of smoke and would like to plot/shade it as follow (I did it with Ferret but would like to know how to do it with Matlab) :

My matrix is defined along x & z axis with concentration(x,z) defined as a sinusoidal function

Thanks for all,

Florian

Best Answer

Why not make save be a 3D array instead of a cell array? Why not just use imshow() or image() to display each 2D image of "save" right after you created it.

Related Solutions

MATLAB: Solving linear convection equation (wave equation) by Lax Wendroff Scheme

N=100; %No. of grid points
Tmax=1; % time period=1
alpha=1;  %given
h=0.01; %given value of delX
delt=0.005;
maxt=Tmax/delt; %Number of time steps
c=alpha*delt/h; %c=0.5
u=zeros(N+1,maxt+1);
x=zeros(N+1,1);
% Initial condition
for i=1:N+1
  x(i)=(i-1)*h;
  u(i,1)=sin(2*pi*x(i));
end
for k=1:maxt
  u0 = u(N,k);
  u(1,k+1)=u(1,k)-c/2*(u(2,k)-u0)+c^2/2*(u(2,k)-2*u(1,k)+u0);
  for i=2:N
    u(i,k+1) =u(i,k)-c/2*(u(i+1,k)-u(i-1,k))+(c^2)/2*(u(i+1,k)-2*u(i,k)+u(i-1,k));
  end
  uNp2 = u(2,k);
  u(N+1,k+1) =u(N+1,k)-c/2*(uNp2-u(N,k))+c^2/2*(uNp2-2*u(N+1,k)+u(N,k));  
end
plot(x,u(:,10),x,u(:,20),x,u(:,30),x,u(:,40),x,u(:,50),x,u(:,60),x,u(:,70))

MATLAB: The usage of [~,I]=min(abs(A+B)) on a specific example

I don't really see how min can help you. Basic logical indexing would be simpler, e.g.:

>> nnz(u<umin | u>umax)
ans =  4

Note you can get both values using one logical index:

>> idx = u<umin | u>umax;
>> nnz(idx)
ans =  4
>> nnz(~idx)
ans =  6

Best Answer

Related Solutions

MATLAB: Solving linear convection equation (wave equation) by Lax Wendroff Scheme

MATLAB: The usage of [~,I]=min(abs(A+B)) on a specific example

Related Question