| 查看: 461 | 回复: 1 | ||
zpzgraduate铁杆木虫 (著名写手)
|
[求助]
matlab求解偏微分方程组遇到的问题求助
|
|
求解偏微分方程组时遇到的，请高手指点。 下面是文件和错误提示 function pdex4 m = 0; x = [0 0.005 0.01 0.05 0.1 0.2 0.5 0.7 0.9 0.95 0.99 0.995 1]; t = [0 0.005 0.01 0.05 0.1 0.5 1 1.5 2]; sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t); u1 = sol(:,:,1); u2 = sol(:,:,2); u3=sol(:,:,2); figure surf(x,t,u1) title('u1(x,t)') xlabel('Distance x') ylabel('Time t') figure surf(x,t,u2) title('u2(x,t)') xlabel('Distance x') ylabel('Time t') figure surf(x,t,u3) title('u3(x,t)') xlabel('Distance x') ylabel('Time t') % -------------------------------------------------------------- function [c,f,s] = pdex4pde(x,t,u,DuDx) c = [1; 1;1]; f = [0; 0;0] .* DuDx; F = 6*u(3); s = [-F*u(1); u(2)*F;-20000*u(1)-2-4*u(2)]; % -------------------------------------------------------------- function u0 = pdex4ic(x); u0 = [0.03;0;0]; % -------------------------------------------------------------- function [pl,ql,pr,qr] = pdex4bc(xl,ul,xr,ur,t) pl = [0; ul(2)-1;0]; ql = [1; 1;1]; pr = [ur(1)-1; 0;0]; qr = [1;1;1]; 错误信息 Warning: Failure at t=1.497828e-002. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (2.775558e-017) at time t. > In ode15s at 809 In pdepe at 317 In pdefun001 at 6 Warning: Time integration has failed. Solution is available at requested time points up to t=1.000000e-002. > In pdepe at 323 In pdefun001 at 6 ??? Error using ==> surf at 78 Data dimensions must agree. Error in ==> pdefun001 at 12 surf(x,t,u1) |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有48人回复
-
FileCode能看出啥？
已经有17人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有6人回复
-
filecode
已经有9人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有23人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有6人回复
-
为什么网上很多人说本周 12号出结果
已经有5人回复
-
帮忙看看fileCode
已经有11人回复
-
基金中了
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- matlab解偏微分方程求助
已经有12人回复
- 求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？
已经有5人回复
- matlab求解偏微分方程组
已经有3人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab解微分方程组
已经有15人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】如何求解下列偏微分方程
已经有13人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
zpzgraduate
铁杆木虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 14109.8
- 散金: 40
- 帖子: 1152
- 在线: 202.1小时
- 虫号: 116996
- 注册: 2005-11-25
- 专业: 聚合物共混与复合材料
2楼2013-07-01 05:49:22
100