SURF SALE NEGRO - GRADIENTE COLORMAP ARGUMENTOS-
Publicado por Frikilab (1 intervención) el 25/04/2013 11:08:59
Estoy intentando hacer Surf con el siguiente código y la Figura sale negra. Se supone que mido concentraciones y difusividad por lo tanto debería aparecer un gradiente de colores.
Tengo 4 scripts:
pdex1pde.m
%dudx=du/dx
% en nuestro caso dudx= DcDr
function [c,f,s] = pdefun(x,t,u,dudx)
Dij=1e-4; % este valor lo vamos cambiando
c=1/Dij;
f=dudx;
s=0;
pdex1bc.m
function [pl,ql,pr,qr] = pdex1bc(rl,Cl,rr,Cr,t)
Co= 10; %constante que podemos ir cambiando
%que valor le ponemos a Cl y Cr?????
%%Condición C(0)= Co y C(L)=Co/2
pl = Cl - Co;
ql = 0.0;
pr = Cr - Co/2;
qr = 0.0;
% pl = 0;
% ql = 1.0;
% pr = Cr - Co;
% qr = 0.0;
end
pdex1ic.m
function u0 = pdex1ic(x)
Co = 10;
u0 = Co;
end
main.m
%% MAIN
%c=u(concentración) c(r,t)
%t=t time, variable temporal
%x=r radio, variable espacial
clear all
clc
m=2;
r=linspace(0.5,1,100);
t=linspace(0,100,1000);
%linspace(0,1,20) te crea un vector con 20 valores entre 0 y 1 (random)
sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,r,t);
C = sol(:,:,1);
surf(r,t,C)
title('Numerical solution ')
xlabel('Radius')
ylabel('Time')
zlabel('Concentration')
¡MUCHAS GRACIAS!
Tengo 4 scripts:
pdex1pde.m
%dudx=du/dx
% en nuestro caso dudx= DcDr
function [c,f,s] = pdefun(x,t,u,dudx)
Dij=1e-4; % este valor lo vamos cambiando
c=1/Dij;
f=dudx;
s=0;
pdex1bc.m
function [pl,ql,pr,qr] = pdex1bc(rl,Cl,rr,Cr,t)
Co= 10; %constante que podemos ir cambiando
%que valor le ponemos a Cl y Cr?????
%%Condición C(0)= Co y C(L)=Co/2
pl = Cl - Co;
ql = 0.0;
pr = Cr - Co/2;
qr = 0.0;
% pl = 0;
% ql = 1.0;
% pr = Cr - Co;
% qr = 0.0;
end
pdex1ic.m
function u0 = pdex1ic(x)
Co = 10;
u0 = Co;
end
main.m
%% MAIN
%c=u(concentración) c(r,t)
%t=t time, variable temporal
%x=r radio, variable espacial
clear all
clc
m=2;
r=linspace(0.5,1,100);
t=linspace(0,100,1000);
%linspace(0,1,20) te crea un vector con 20 valores entre 0 y 1 (random)
sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,r,t);
C = sol(:,:,1);
surf(r,t,C)
title('Numerical solution ')
xlabel('Radius')
ylabel('Time')
zlabel('Concentration')
¡MUCHAS GRACIAS!
Valora esta pregunta
0