24小时热门版块排行榜

>论坛更新日志 (4397)
>文献求助 (274)
>导师招生 (245)
>博后之家 (225)
>论文投稿 (215)
>虫友互识 (160)
>硕博家园 (84)
>有机交流 (80)
>考博 (73)
>招聘信息布告栏 (45)
>学术会议 (37)
>休闲灌水 (31)
>基金申请 (21)
>找工作 (18)
>考研 (16)
>教师之家 (15)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

土老九

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1357.3
散金: 236
红花: 1
帖子: 126
在线: 39小时
虫号: 841055
注册: 2009-09-05
专业: 作物栽培与耕作学

[求助] 求助：偏微分方程求解

期盼高人指点

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 求助：偏微分方程.doc

2013-01-23 11:43:17, 20 K

» 本帖已获得的红花（最新10朵）

土老九

» 猜你喜欢

散金币祈福已经有187人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有185人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有114人回复
深圳大学2027级光电信息工程/物理学研究生推免已经有70人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有4人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有43人回复

耕者有其田

1楼 2013-01-23 11:43:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

土老九

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1357.3
散金: 236
红花: 1
帖子: 126
在线: 39小时
虫号: 841055
注册: 2009-09-05
专业: 作物栽培与耕作学

★
fegg7502: 金币+1, 鼓励交流 2013-01-29 09:03:57

引用回帖:

14楼: Originally posted by southpark007 at 2013-01-27 21:32:44
用PDEPE来做就可以了，把t设定成无穷大，可以绕过t无穷大时候的h情况。

function pdex1
m=0;
x=linspace(0,0.14,14);
t=linspace(0,100000,100);
sol=pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t);
% Extract the first solution component as u.
u=sol(:,:,1);
% A surface plot is often a good way to study a solution.
surf(x,t,u)
title('Numerical solution computed with 20 mesh points.')
xlabel('Distance x')
ylabel('Time t')
% A solution profile can also be illuminating.
figure
plot(x,u(end,

)
title('Solution at t=2')
xlabel('Distance x')
ylabel('u(x,2)')
% --------------------------------------------------------------
function [c,f,s]=pdex1pde(x,t,u,DuDx)
c=u^5;
f=0.0000001*DuDx;
s=0;
% -------------为保证收敛，初始u设置成反正切函数（三角形）增加-------------------------
function u0=pdex1ic(x)
u0=10*atan(100000*x)/(pi/2);
% --------------------------------------------------------------
function [pl,ql,pr,qr]=pdex1bc(xl,ul,xr,ur,t)
pl=u;
q1=0;
pr=0;
qr=1;