版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (4157)
>虫友互识 (752)
>导师招生 (551)
>论文投稿 (181)
>招聘信息布告栏 (158)
>文献求助 (138)
>文学芳草园 (103)
>博后之家 (92)
>硕博家园 (88)
>学术会议 (75)
>考博 (59)
>休闲灌水 (44)
>考研 (41)
>公派出国 (35)
>找工作 (34)
>论文道贺祈福 (33)

返回列表

zju-hjm

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1014.3
帖子: 33
在线: 22.8小时
虫号: 758742
注册: 2009-04-27
专业: 高分子合成化学

[交流] 【求助】用matlab解含第二类边界条件的偏微分方程已有1人参与

一个简单的抛物线型偏微分方程，求u（x，t）
u_t 表示对t的一阶导数，u_x为对x的一阶导数，u_xx表示对x的二阶导数

方程 u_t =u_xx      （1）

边界条件
x=0,t>0，u_x=0；  （2）
x=1,t>0,  u(1,t)=0;
  （3）
初始条件
t=0, u(x,0)=1;          （4）
我用有限差分法，主要是第二类边界条件（2）式怎么在matlab里面利用
我大致解了一下，但觉得结果有问题，而且我认为问题是出在（2）式上。真心求教这个问题怎么解，谢谢虫友！

附：我写的程序--------------------------------------------------------------------------------------
function [u,x,t]=model2(xf,T)
% solve  u_xx=u_t for 0<=x<=xf,0<=t<=T
% Initial condition:u(x,0)=1
% boundary conditon:u(1,t)=0;u_x(0,t)=0

dx=xf/50;x=[0:50]'*dx;
dt=T/10000;t=[0:10000]*dt;
for i=1:51,u(i,1)=1;end % Initial condition:u(x,0)=1
for n=1:10001,u(51,n)=0;end % boundary conditon:u(1,t)=0
r=dt/dx/dx;r1=1-2*r;
for k=1:10000
u(1,k+1)=2*r*u(2,k)+r1*u(1,k);end % boundary conditon:u_x(0,t)=0 ?
for k=1:10000
for i=2:50
      u(i,k+1)=r1*u(i,k)+r*(u(i-1,k)+u(i+1,k));% solve  u_xx=u_t for 0<=x<=xf,0<=t<=T
end
end

--------------------------------------------------------------------------------------
命令窗口输入：

xf=1;
T=1;
[u,x,t]=model2(xf,T);
figure(1),clf,mesh(t,x,u);
==============================================

[ Last edited by zju-hjm on 2009-12-17 at 16:28 ]

回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2009-12-17 16:24:48

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

keventu