Pluma Gaussiana
Publicado por Beto (4 intervenciones) el 02/02/2018 22:33:19
saludos! alguien podría ayudarme AL CORRER ESTE PROGRAMA EN LA MATRIZ TT DE ONES, ME GENERA UN ERROR, lo he cambiado pero no me da el resultado deseado, en el matlab 2014 en adelante
que se puede hacer?
gracias!
Error using ones
Size inputs must be scalar.
Error in PLUMA_GAUSSIANA_2018 (line 17)
tt=ones(1,x1);
que se puede hacer?
gracias!
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
%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')
Error using ones
Size inputs must be scalar.
Error in PLUMA_GAUSSIANA_2018 (line 17)
tt=ones(1,x1);
Valora esta pregunta
![Me gusta: Está pregunta es útil y esta clara Me gusta: Está pregunta es útil y esta clara](/img/img.png?11.51)
![NO me gusta: Está pregunta no esta clara o no es útil No me gusta: Está pregunta no esta clara o no es útil](/img/img.png?11.51)
0