• Author / Uploaded
• Rene Llipita Garcia

Citation preview

MODELO GAUSSIANO PARA UNA FUENTE PUNTUAL clear all clc %%% parametros Dx=0.000625; %difusividad v=0.1; %velocidad M=1; %masa %intervalos xmin=-0.05 ; xmax=2.015; %intervalo en eje x t=1:4:20; %...........proceso de ejecución……………….. x=linspace(xmin,xmax,10); c=[]; for i=1:size(t,2) xx=x-v*t(i); c=(M/sqrt(4*pi*Dx*t(i)))*ones(1,size(x)).*exp(-(xx.*xx)/4*Dx*t(i)); end %........................salida……………. plot(x,c') hold on; xlabel('space') ylabel('concentracion') title ('modelo gausiano en 1D') grid on

clear all clc %%% parametros Dx=0.000625; %difusividad v=0.1; %velocidad M=1; %masa %intervalos xmin=-0.05 ; xmax=2.015; %intervalo en eje x t=1:4:20; %...........proceso de ejecución……………….. x=linspace(xmin,xmax,100); c=[]; for i=1:size(t,2) xx=x-v*t(i); c=(M/sqrt(4*pi*Dx*t(i)))*ones(1,size(x)).*exp(-(xx.*xx)/4*Dx*t(i)); end %........................salida……………. plot(x,c','ro') hold on; xlabel('space') ylabel('concentracion') title ('modelo gausiano en 1D') grid on

clear all clc %%% parametros Dx=0.000625; %difusividad v=0.1; %velocidad M=1; %masa %intervalos xmin=-0.05 ; xmax=2.015; %intervalo en eje x t=1:4:20; %...........proceso de ejecución……………….. x=linspace(xmin,xmax,100); [x,t] = meshgrid(x,t); for i=1:size(t,2) xx=x-v*t(i); c=(M/sqrt(4*pi*Dx*t(i))).*exp(-(xx.*xx)/4*Dx*t(i)); end %........................salida……………. surf(x,t,c) shading interp xlabel('distancia') ylabel('tiempo') zlabel('concentración') title ('modelo gausiano en 3D') colorbar