Attempt to execute SCRIPT pdepe as a function:
3 次查看(过去 30 天)
显示 更早的评论
Hi dears
I have problem with executing this program. Can anyone explain me what is going on here? /
Attempt to execute SCRIPT pdepe as a function: /Applications/MATLAB_R2019b.app/toolbox/matlab/funfun/pdepe.m Error in kpz (line 7) sol = pdepe(m,@heat1D,@heatIC,@heatBC,x,t);
opts = odeset('RelTol',1e-8,'AbsTol',1e-8);
N=200; xi=0; xf=5; ti=0; tf=5; TIME=tf-ti;
x = linspace(xi,xf,N);
t = linspace(ti,tf,200);
m = 0;
tic
sol = pdepe(m,@heat1D,@heatIC,@heatBC,x,t);
toc
w = sol(:,:,1);
surf(x,t,w)
xlabel('x')
ylabel('t')
zlabel('w(x,t)')
view([40 55])
%plot(w(200,:))
function [c,f,s] = heat1D(x,t,w,dwdx)
D=4; %0.2;
c = 1;
f = D*dwdx;
s = 1.5*dwdx.^2; %% NO NOISE
%s = 1*dwdx.^2 - 0.5/t; %%White
%s = 0.5*dwdx.^2 - 1/(x*sqrt(t)); %%Pink
%s = -0.5*dwdx.^2 - 0.5/(x*x); %%Brown
%s = 10.5*dwdx.^2 - 0.1*exp(-x*x/t)/t; %%Gauss
end
function w0 = heatIC(x)
w0=2;
end
function [pl,ql,pr,qr] = heatBC(xl,wl,xr,wr,t)
pl = wl;
ql = 0;
pr = 0;
qr = 1;
end
2 个评论
madhan ravi
2020-3-21
which pdepe -all % what do you get ?
You have a script named "pdepe.m" probably rename or remove it.
采纳的回答
Cris LaPierre
2020-3-21
Your code runs as is for me. Perhaps you have overwritten the pdepe function? Try clearing it by typing
clear pdepe
Run the code again. If you still get the error, try typing
which pdepe
The result should be something like this
C:\Program Files\MATLAB\R2020a\toolbox\matlab\funfun\pdepe.m
14 个评论
Cris LaPierre
2020-3-21
编辑:Cris LaPierre
2020-3-21
This error suggests that MATLAB is finding multiple files called pdepe.m. Let's see what files MATLAB is seeing. What is the result of
which pdepe -all
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Eigenvalue Problems 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!