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

登录注册

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

24小时热门版块排行榜

返回列表

zhangzhiguang

至尊木虫 (文坛精英)

应助: 121 (高中生)
金币: 53088.8
红花: 24
沙发: 222
帖子: 26519
在线: 881.4小时
虫号: 492270
注册: 2008-01-08
性别: GG
专业: 基础物理学

别太着急,就这几天,突然有点事儿要办

赞一下

回复此楼

好好学习,天天向上!

11楼2013-01-26 11:46:51

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

zhangzhiguang

至尊木虫 (文坛精英)

应助: 121 (高中生)
金币: 53088.8
红花: 24
沙发: 222
帖子: 26519
在线: 881.4小时
虫号: 492270
注册: 2008-01-08
性别: GG
专业: 基础物理学

【答案】应助回帖

★
csgt0: 金币+1, 谢谢 2013-01-28 09:35:01

clear;clc;
h=zeros(100,100);
h(:,1)=10;
h(1,

=0.1;
h(:,100)=0.1;
omega=1.1;

for d=1:100
for n=1:98
      for i=2:99

         h(100,n+1)=1e-5*(h(100,n))^-5*(h(99,n)-2*h(100,n)+h(99,n))+h(100,n);
         h(i,n+1)=(1-omega)* h(i,n+1)+omega*(1e-5*(h(i,n))^-5*(h(i+1,n)-2*h(i,n)+h(i-1,n))+h(i,n));
      end
end
end

mesh(h)

赞一下

回复此楼

好好学习,天天向上!

12楼2013-01-27 13:22:00

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

zhangzhiguang

至尊木虫 (文坛精英)

应助: 121 (高中生)
金币: 53088.8
红花: 24
沙发: 222
帖子: 26519
在线: 881.4小时
虫号: 492270
注册: 2008-01-08
性别: GG
专业: 基础物理学

【答案】应助回帖

★
csgt0: 金币+1, 谢谢 2013-01-28 09:34:43

那个笑脸是:加),说明,空间步长是0.01H,时间是0.01T,取得时间长度为100个步长,这里你回去自己修改吧,如果取得较大把mesh(h)去掉，不知你计算的是什么，感觉算出来不太对劲，你那公式里h的5次方h是系数吗？我是按函数算得

赞一下(1人)

回复此楼

好好学习,天天向上!

13楼2013-01-27 13:26:08

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

southpark007

禁虫 (著名写手)

★
csgt0: 金币+1, 应助指数+1, 谢谢 2013-01-28 09:35:13

本帖内容被屏蔽

14楼2013-01-27 21:32:44

已阅回复此楼关注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;

赞一下

回复此楼

耕者有其田

15楼2013-01-28 08:42:44

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