%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')