24小时热门版块排行榜

返回列表

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

chugezhuang

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 46
帖子: 14
在线: 9小时
虫号: 1407185
注册: 2011-09-19
性别: GG
专业: 岩土与基础工程

[求助] 哪位大仙帮我看下程序，我用的matlab7.1，运行提示du未定义

function [c,f,s]=pdefun(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]=pdebc(xa,ua,xb,ub,t)
pa=[0;ua(2)];
qa=[1;0];
pb=[ub(1)-1;0];
qb=[0;1];
function u0=pdeic(x)
u0=[1;0];
clc
x=0:0.05:1;
t=0:0.05:2;
m=0;
sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);
figure('numbertitle','off','name','PDE Demo——by Matlabsky')
subplot(211)
surf(x,t,sol(:,:,1))
title('The Solution of u_1')
xlabel('X')
ylabel('T')
zlabel('U')
subplot(212)
surf(x,t,sol(:,:,2))
title('The Solution of u_2')
xlabel('X')
ylabel('T')
zlabel('U')

回复此楼

» 猜你喜欢

售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有5人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有4人回复
上海工程技术大学激光智能制造课题组 2027级博士研究生招生已经有6人回复
课题组招2027级博士上海工程技术大学激光智能制造方向已经有6人回复
现代”学阀”该如何界定已经有14人回复
求合成方法已经有8人回复
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？已经有4人回复
有机合成以后会不会被AI改变？做科研的虫友怎么看已经有8人回复
我的奶奶已经有3人回复

Life isn't about waiting for the storm to pass, it's about learning to dance in the rain.

1楼 2012-12-03 10:11:23

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

Kuiyanlong

铁虫 (初入文坛)

应助: 4 (幼儿园)
金币: 1135.5
帖子: 25
在线: 27.1小时
虫号: 2166508
注册: 2012-12-04
性别: GG
专业: 机械动力学

【答案】应助回帖

★
xiegangmai: 金币+1, 谢谢参与 2012-12-06 19:45:17

我试了一下，while循环结果不是很理想，用for循环验算，发现u的结果太小
ep=1e-6;
p0=0.25;
t=0.02;
Cv=0.04;
Tv=Cv*t;
L=0.5;
z=3;
n=0;
c0=2.*L.*p0./pi;
u=0;
for n=0:100
m=2.*n+1
u=u+1/m*sin(m/(2*L))*pi*z*eps(-Tv*m^2*pi*pi/4)
end

赞一下

回复此楼

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

chugezhuang

Bebornthisway

8楼2012-12-06 18:50:32

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

查看全部 12 个回答

Kuiyanlong

铁虫 (初入文坛)

应助: 4 (幼儿园)
金币: 1135.5
帖子: 25
在线: 27.1小时
虫号: 2166508
注册: 2012-12-04
性别: GG
专业: 机械动力学

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与 2012-12-05 19:55:09

换一下表达就可以了，将 sol函数提前
function ABC
x=0:0.05:1;
t=0:0.05:2;
m=0;
sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);
figure('numbertitle','off','name','PDE Demo——by Matlabsky')
subplot(211)
surf(x,t,sol(:,:,1))
title('The Solution of u_1')
xlabel('X')
ylabel('T')
zlabel('U')
subplot(212)
surf(x,t,sol(:,:,2))
title('The Solution of u_2')
xlabel('X')
ylabel('T')
zlabel('U')

function [c,f,s]=pdefun(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]=pdebc(xa,ua,xb,ub,t)
pa=[0;ua(2)];
qa=[1;0];
pb=[ub(1)-1;0];
qb=[0;1];
function u0=pdeic(x)
u0=[1;0];
clc

赞一下

回复此楼

Bebornthisway

2楼2012-12-05 00:47:40

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