| 查看: 1232 | 回复: 5 | ||
独孤神宇版主 (知名作家)
|
[求助]
一阶 偏微分 方程怎么求解 ？？？？ 已有1人参与
|
|
请问个路大神，如图所示的 一阶 偏微分 方程怎么求解？ 除了 特征曲线法外，能不能直接用命令求解？？？ 比如 MATLAB中的 pdepe 函数， 或者 maple 里面的 pdsolve ???  IMG_1043.JPG |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab龙格库塔法求解初值为零高阶微分方程编程数学问题请教
已经有3人回复
- 求解四阶偏微分方程的解析解
已经有1人回复
- 求大神帮忙求解偏微分方程
已经有0人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18157.7
- 红花: 97
- 帖子: 8891
- 在线: 1486小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
【答案】应助回帖
|
function pde %主函数 x=0:0.05:1; t=0:0.05:2; m=0; sol=pdepe( m,@mpde,@mpic,@mpbc,x,t ); surf( x,t,sol( :,:,1 ) ) figure; surf( x,t,sol( :,:,2 ) ) function [c,f,s]=mpde( x,t,u,du ) %给出偏微分方程的函数描述 c=[1;1]; 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]=mpbc( xa,ua,xb,ub,t ) %边界条件描述 pa=[0;ua( 2 )]; qa=[1;0]; pb=[ub( 1 )-1;0]; qb=[0;1]; function u0=mpic( x ) %初值描述 u0=[1;0]; Matlab求解微分方程（2）——偏微分方程的求解 https://blog.sina.com.cn/s/blog_b0ecb51d0101bogk.html |
2楼2016-08-07 22:32:32
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
3楼2016-08-08 09:28:31
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18157.7
- 红花: 97
- 帖子: 8891
- 在线: 1486小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
|
PDEPE Automatic Solution of 1D Initial Boundary Value Problems (IBVP's) https://people.sc.fsu.edu/~jburkardt/m_src/pdepe/pdepe.html |
4楼2016-08-08 09:33:50
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
5楼2016-08-08 15:23:11
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18157.7
- 红花: 97
- 帖子: 8891
- 在线: 1486小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
★ ★ ★ ★ ★
独孤神宇(月只蓝代发): 金币+5, 感谢指导 2018-05-19 10:55:33
独孤神宇(月只蓝代发): 金币+5, 感谢指导 2018-05-19 10:55:33
|
參考 https://scicomp.stackexchange.co ... e-peclet-number-for https://scicomp.stackexchange.co ... boundary-conditions |
6楼2016-08-08 15:51:58