| 查看: 2575 | 回复: 4 | ||||
Math露珠木虫 (知名作家)
|
[交流]
pde求解的MATLAB语言求助 已有3人参与
|
|
如下网址上的MATLAB语言， http://www.matlabsky.com/thread-981-1-1.html 我想按照他的代码写在M文件中，运行一下，但是一直运行不出来。 第四行，他就说我没有定义 ??? Input argument "ua" is undefined. Error in ==> pde at 4 pa=[0;ua(2)]; 以下是他的代码。求助大神帮我改一下。谢谢 %% 目标PDE函数 %%syms x t u function [c,f,s]=pdefun (x,t,u,du); c=[1;1]; %%dudx=diff(u,'x'); f=[0.024*du(1);0.17*du(2)]; temp=u(1)-u(2); s=[-1;1].*(exp(5.73*temp)-exp(-11.46*temp)); %% 边界条件函数 function [pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t) %a表示下边界，b表示上边界 pa=[0;ua(2)]; qa=[1;0]; pb=[ub(1)-1;0]; qb=[0;1]; %% 初值条件函数 function u0=pdeic(x) u0=[1;0]; %% 主调函数 clc x=0:0.05:1; t=0:0.05:2; m=0; sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t); figure('numbertitle','off','name','PDE Demo——by Matlabsky') subplot(211) surf(x,t,sol(:,:,1)) title('The Solution of u_1') xlabel('X') ylabel('T') zlabel('U') subplot(212) surf(x,t,sol(:,:,2)) title('The Solution of u_2') xlabel('X') ylabel('T') zlabel('U') |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有3人回复
-
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.2小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
程序没有问题! 可能你对Matlab的调用不清吧!注意函数顺序,下面是M文件及运行结果! function pdefcqj %% 主调函数 clc x=0:0.05:1; t=0:0.05:2; m=0; sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t); figure('numbertitle','off','name','PDE Demo——by Matlabsky') subplot(211) surf(x,t,sol(:,:,1)) title('The Solution of u_1') xlabel('X') ylabel('T') zlabel('U') subplot(212) surf(x,t,sol(:,:,2)) title('The Solution of u_2') xlabel('X') ylabel('T') zlabel('U') function [c,f,s]=pdefun (x,t,u,du); c=[1;1]; %%dudx=diff(u,'x'); f=[0.024*du(1);0.17*du(2)]; temp=u(1)-u(2); s=[-1;1].*(exp(5.73*temp)-exp(-11.46*temp)); %% 边界条件函数 function [pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t) %a表示下边界，b表示上边界 pa=[0;ua(2)]; qa=[1;0]; pb=[ub(1)-1;0]; qb=[0;1]; %% 初值条件函数 function u0=pdeic(x) u0=[1;0];  pdetx.jpg |
Math露珠
2楼2017-04-22 08:57:24
Math露珠
木虫 (知名作家)
- 应助: 10 (幼儿园)
- 金币: 7841.6
- 散金: 2548
- 红花: 123
- 沙发: 33
- 帖子: 8870
- 在线: 885.7小时
- 虫号: 4010226
- 注册: 2015-08-07
- 专业: 常微分方程与动力系统
3楼2017-04-23 16:00:51
4楼2017-04-24 11:26:09
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
刚接触matlab求解偏微分方程，按照pdepe帮助里的例子写的程序，一运行就崩溃，请求帮忙 function pdex2 clc clear m = 0; x = linspace(0,1,20); t = linspace(0,2,5); sol = pdepe(m,@pdex2pde,@pdex2ic,@pdex2bc,x,t); u = sol(:,:,1); figure surf(x,t,u) title('Numerical solution computed with 20 mesh points.') xlabel('Distance x') ylabel('Time t') figure plot(x,u(end,  ) title('Solution at t = 2') xlabel('Distance x') ylabel('u(x,2)') function [c,f,s] = pdex2pde(x,t,u,DuDx) c = pi^2; f = DuDx; s = 0; function u0 = pdex2ic(x) u0 = sin(pi*x); function [pl,ql,pr,qr] = pdex2bc(xl,ul,xr,ur,t) pl = ul; ql = 0; pr = pi * exp(-t); qr = 1; |
5楼2017-10-23 09:31:37