| 查看: 1768 | 回复: 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];  |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有54人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有29人回复
-
关于Filecode分析方法
已经有12人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有11人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助matlab问题
已经有5人回复
- 请教关于matlab矩阵转换问题
已经有6人回复
- 关于Matlab作图的问题
已经有14人回复
- matlab调试问题
已经有3人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab非线性参数拟合问题
已经有7人回复
- 【求助】麻烦大家看看我的邀请信有没有问题
已经有9人回复
- 【求助】大家看看，matlab安装过程的截图，找出错
已经有14人回复
- 【求助】MATLAB多个mat中 数据画图问题
已经有7人回复
- 【求助】大家帮我看看这个MATLAB程序？
已经有6人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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]; |
亮亮1283
4楼2012-05-28 14:35:31
2楼2012-05-26 21:43:32
3楼2012-05-27 12:52:04
5楼2012-05-28 15:20:24
