| 查看: 1795 | 回复: 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];  |
回复此楼» 猜你喜欢
有没有大神帮我看看基金代码
已经有10人回复
-
今天务委会开完了，明天出结果吗
已经有13人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助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
