Matlab - Pluma Gaussiana

%Pluma Gaussiana Básica

clc

clear all

Dx=0.000625; %Coeficiente de difusion

V=0.1;       %Velocidad de Viento

M=1;         %Masa del Contaminante

xmin=-0.05;

xmax=0.2;

ymin=0; ymax=100;

t=0.1:0.005:1;  %tiempo

x=linspace(xmin,xmax,size(t,2));

x1=size(x);

tt=ones(1,x1);

y=t;   %para que y considere el contorno de t

C=[];  %Matriz

for i=1:size(t,2)

    xx=x-V*t(i);

    C=[C:[M/sqrt(4*pi*Dx*t(i))*tt.*exp(-(xx.*xx)./(4*Dx*t(i)))]];

end

figure

plot(C');

hold on

title('Pluma Gaussiana Puff 1D')

xlabel('Espacio')

ylabel('Concentracion')

hold on

[A B]=meshgrid(x,y);

figure

contourf(A,B,C,10);

colorbar;

xlabel('Eje X')

ylabel('Tiempo')

%Pluma Gaussiana Básica

clc

clear all

Dx=0.000625; %Coeficiente de difusion

V=0.1; %Velocidad de Viento

M=1; %Masa del Contaminante

xmin=-0.05;

xmax=0.2;

ymin=0; ymax=100;

t=0.1:0.005:1; %tiempo

x=linspace(xmin,xmax,size(t,2));

x1=size(x,2);

tt=ones(1,x1);

y=t; %para que y considere el contorno de t

C=[]; %Matriz

for i=1:size(t,2)

xx=x-V*t(i);

C=[C;[M/sqrt(4*pi*Dx*t(i))*tt.*exp(-(xx.*xx)./(4*Dx*t(i)))]];

end

figure

plot(C');

hold on

title('Pluma Gaussiana Puff 1D')

xlabel('Espacio')

ylabel('Concentracion')

hold on

[A B]=meshgrid(x,y);

figure

contourf(A,B,C,10);

colorbar;

xlabel('Eje X')

ylabel('Tiempo')