| 查看: 1415 | 回复: 6 | |||
[求助]
matlab程序问题，总是出错，麻烦大家看看问题在那里？
|
|
function pdex1() clc clear Up =400; Un =8000; Kb =1.38*10^(-23); h =6.63*10^(-34); e =1.6*10^(-19); r =0.35; y =0.6; n0 =1*10^16; p0 =1*10^16; tn =1*10^(-9); tp =1*10^(-9); p1 =p0; n1 =n0; T =300; a =1.65*10^4; v =(3/532)*10^17; X =10; P=10; %主函数---------------------------------------------------------------- m=0; x=0:0.1:10; t=0:20:1000; sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t); figure('numbertitle','off','name','PDE Demo-by matlab') subplot(211) surf(x,t,sol(:,:,1)) title('The Solution of p') xlable('x') ylable('t') zlable('p') subplot(212) surf(x,t,sol(:,:,2)) title('The Solution of n') xlable('x') ylable('t') zlable('n') %目标函数---------------------------------------------------------------- function [c,f,s] = pdex1pde(x,t,u,DuDx) c= [1;1]; f= [(Up*Kb*T/e)*DuDx(1); (Un*Kb*T/e)*DuDx(2)]; G= a*y*(1-r)*P*exp(-a*x)/(h*v); R1= u(1)*u(2)-p0*n0; R2=tn*(u(1)+p1)+tp*(u(2)+n1); R=R1/R2; s= [1;1].*(G-R); %初值条件函数------------------------------------------------------------ function u0 = pdex1ic(x) u0= [p0;n0]; %边界条件函数------------------------------------------------------------ function [pa,qa,pb,qb] = pdex1bc(xa,ua,xb,ub,t) %a为下边界， b为上边界 pa =[0; 0]; qa =[1; 1]; pb =[0; 0]; qb =[1; 1];  |
回复此楼» 猜你喜欢
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助matlab问题
已经有5人回复
- 请教关于matlab矩阵转换问题
已经有6人回复
- 关于Matlab作图的问题
已经有14人回复
- matlab调试问题
已经有3人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab非线性参数拟合问题
已经有7人回复
- 【求助】麻烦大家看看我的邀请信有没有问题
已经有9人回复
- 【求助】大家看看，matlab安装过程的截图，找出错
已经有14人回复
- 【求助】MATLAB多个mat中 数据画图问题
已经有7人回复
- 【求助】大家帮我看看这个MATLAB程序？
已经有6人回复
2楼2012-05-26 21:43:32
3楼2012-05-27 12:52:04
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
亮亮1283: 金币+20, ★★★★★最佳答案, 哈哈，就是这个原因。太谢谢你了！ 2012-05-28 15:17:02
感谢参与，应助指数 +1
亮亮1283: 金币+20, ★★★★★最佳答案, 哈哈，就是这个原因。太谢谢你了！ 2012-05-28 15:17:02
|
不知道对不对，画出来的结果好像都是1*10^16 原程序主要有两个问题，你的那些参数需要用全局变量，子函数不能直接用主函数的值的 还有xlabel，ylabel，zlabel，不是xlable等 function pdex1() clc clear global Up Kb T P a y r h v Un e tp p0 n0 n1 p1 G tn Up =400; Un =8000; Kb =1.38*10^(-23); h =6.63*10^(-34); e =1.6*10^(-19); r =0.35; y =0.6; n0 =1*10^16; p0 =1*10^16; tn =1*10^(-9); tp =1*10^(-9); p1 =p0; n1 =n0; T =300; a =1.65*10^4; v =(3/532)*10^17; X =10; P=10; %主函数---------------------------------------------------------------- m=0; x=0:0.1:10; t=0:20:1000; sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t); % u =[ sol(:,:,1), sol(:,:,2)]; figure('numbertitle','off','name','PDE Demo-by matlab') subplot(211) surf(x,t,sol(:,:,1)) title('The Solution of p') xlabel('x') ylabel('t') zlabel('p') subplot(212) surf(x,t,sol(:,:,2)) title('The Solution of n') xlabel('x') ylabel('t') zlabel('n') %目标函数---------------------------------------------------------------- function [c,f,s] = pdex1pde(x,t,u,DuDx) global Up Kb T P a y r h v Un e tp p0 n0 n1 p1 G tn c= [1;1]; f= [(Up*Kb*T/e)*DuDx(1); (Un*Kb*T/e)*DuDx(2)]; G= a*y*(1-r)*P*exp(-a*x)/(h*v); R1= u(1)*u(2)-p0*n0; R2=tn*(u(1)+p1)+tp*(u(2)+n1); R=R1/R2; s= [1;1].*(G-R); %初值条件函数------------------------------------------------------------ function u0 = pdex1ic(x) global p0 n0 u0= [p0;n0]; %边界条件函数------------------------------------------------------------ function [pa,qa,pb,qb] = pdex1bc(xa,ua,xb,ub,t) %a为下边界， b为上边界 pa =[0; 0]; qa =[1; 1]; pb =[0; 0]; qb =[1; 1]; |
» 本帖已获得的红花（最新10朵）
亮亮1283
4楼2012-05-28 14:35:31
5楼2012-05-28 15:20:24
6楼2012-05-28 15:35:59
7楼2012-05-28 16:15:44