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!