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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 哪位大仙帮我看下程序，我用的matlab7.1，运行提示du未定义
121/2返回列表12下一页
查看: 1186  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

chugezhuang

铁虫 (初入文坛)

[求助] 哪位大仙帮我看下程序，我用的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')
回复此楼

» 猜你喜欢
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

铁虫 (初入文坛)

【答案】应助回帖

★ ★
感谢参与，应助指数 +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的回帖

Kuiyanlong

铁虫 (初入文坛)

【答案】应助回帖

太可惜了运行出来的图片放不进去
赞一下 回复此楼
Bebornthisway
3楼2012-12-05 00:52:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chugezhuang

铁虫 (初入文坛)
引用回帖:
3楼: Originally posted by Kuiyanlong at 2012-12-05 00:52:47
太可惜了运行出来的图片放不进去

已经解决了，不过还是非常感谢
赞一下 回复此楼
Life isn't about waiting for the storm to pass, it's about learning to dance in the rain.
4楼2012-12-05 16:47:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chugezhuang

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by Kuiyanlong at 2012-12-05 00:47:40
换一下表达就可以了，将 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 ...

你好，再麻烦你下，能不能帮我把这个解的程序编出来，我想要解的图，可我不太会弄。L,Tv,p0都是常数

TC[OK9ZXP3BTFH{S825$[49.jpg
赞一下 回复此楼
Life isn't about waiting for the storm to pass, it's about learning to dance in the rain.
5楼2012-12-05 16:58:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Kuiyanlong

铁虫 (初入文坛)

【答案】应助回帖

t在哪里？？？？
赞一下 回复此楼

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

chugezhuang
Bebornthisway
6楼2012-12-05 21:06:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chugezhuang

铁虫 (初入文坛)
送鲜花一朵
引用回帖:
6楼: Originally posted by Kuiyanlong at 2012-12-05 21:06:53
t在哪里？？？？

Tv=Cv*t，Cv是常数。非常感谢。我看的头大了这几天，呵呵
赞一下 回复此楼
Life isn't about waiting for the storm to pass, it's about learning to dance in the rain.
7楼2012-12-06 10:38:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Kuiyanlong

铁虫 (初入文坛)

【答案】应助回帖

★
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的回帖

Kuiyanlong

铁虫 (初入文坛)

【答案】应助回帖

★
xiegangmai: 金币+1, 谢谢参与 2012-12-06 19:45:07
eps=1e-6;
p0=0.25;
t=0.02;
Cv=0.04;
Tv=Cv*t;
L=0.5;
z=3;
n=1;
c0=2.*L.*p0./pi;
u=0;
u=inline('1/m*sin(m/(2*L))*pi*z*eps(-Tv*m^2*pi*pi/4)','m','L','Tv','z');
u=u(n,L,Tv,z);
eu=0;
while eu>=eps
    m=2*n+1
    eu=u(m,L,Tv,z)
u=u+eu;
eu=co*eu;
end
u=c0*u-p0

还是上面说的问题，那个求和的部分太小
赞一下 回复此楼
Bebornthisway
9楼2012-12-06 19:12:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Kuiyanlong

铁虫 (初入文坛)

【答案】应助回帖

★
xiegangmai: 金币+1, 谢谢参与 2012-12-06 19:45:00
eps=1e-6;
p0=0.25;
t=0.02;
Cv=0.04;
Tv=Cv*t;
L=0.5;
z=3;
n=1;
c0=2.*L.*p0./pi;
c=inline('1/m*sin(m/(2*L))*pi*z*eps(-Tv*m^2*pi*pi/4)','m','L','Tv','z');
u=c(n,L,Tv,z);
eu=1;
while abs(eu)>=eps
    m=2*n+1;
    eu=c(m,L,Tv,z);
u=u+eu;
n=n+1;
eu=c0*eu;
end
u=c0*u-p0;
赞一下 回复此楼
Bebornthisway
10楼2012-12-06 19:27:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chugezhuang 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +18 大米饭！ 2026-03-15 22/1100 2026-03-21 20:20 by HH领袖
[考研] 考研调剂 +3 呼呼？~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 311求调剂 +12 冬十三 2026-03-15 13/650 2026-03-21 19:07 by ColorlessPI
[考研] 一志愿南大，0703化学，分数336，求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 15/750 2026-03-21 10:16 by 脱颖而出
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 一志愿南昌大学，327分，材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
信息提示
关闭
请填处理意见
关闭