该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
请问下各位大佬我这个的错误是怎么修改的,卡在这里了
function newbuild_jiaru_N
clear
format long
m =0;
x=0:10e-10:300e-9;
t=0:10e-15:10e-12;
sol = pdepe(m,@pdex5pde,@pdex5ic,@pdex5bc,x,t);
n = sol(:,:,1);
c = sol(:,:,2);
%v = sol(:,:,3);
subplot(2,2,1)
mesh(x,t,n);
title('Distribution of Electron Temperature--Te(x,t)');
xlabel('Distance x/m');
ylabel('Time t/s');
subplot(2,2,3)
mesh(x,t,c);
title('Distribution of Lattice Temperature--Ti(x,t)');
xlabel('Distance x/m');
ylabel('Time t/s');
subplot(2,2,2)
plot(t,n(:,1),'r',t,c(:,1),'b');
text(0.35e-12,2500,' \leftarrow Te','FontSize',12)
text(1e-12,380,' \uparrow Ti','FontSize',12)
title('Distribution of Electron and Lattice Temperature--Te(0,t),Ti(0,t).');
xlabel('Time t/s');
ylabel('Temperature T/K');
subplot(2,2,4)
plot(x,n(50,:),'r',x,c(50,:),'b');
title('Distribution of Electron and Lattice Temperature--Te(x,1ps),Ti(x,1ps).');
xlabel('Distance x/m');
ylabel('Temperature T/K');
function [c,f,s] = pdex5pde(x,t,u,DuDx)
Ke=0.08*u(1);
kb=1.38*10^-23;
g=1.5*kb*u(1)/0.5*10^-12;
Eg=1.12;
D0=1.8*10^-3*300/u(1);
h=6.63*10^-34;
omg=2.4*10^15; %rad/s
delta=3.6*10^10*exp(-1.5*Eg/(kb*u(1)));
alpha=502000*exp(300/430);
beta=2e-11;
xita=5*10^-22*u(1)/300;
Ce=1.5*kb*u(1);
Ci=1978000+354*u(2)-3680000*u(2)^-2;
R=0.6713;
tao=120e-15;
J=23.89e4; %%%%%有问题
miu=3.8*10^-43;
I=(1-R)/tao*J*(exp(-x/delta-4*log(2)*(t/tao).^2)); %%%%%这个delta没确定
c = [Ce;Ci;1];
f = [ Ke*DuDx(1);0;D0*DuDx]; %这里的舍弃还得斟酌
s1 = -g*(u(1)-u(2))+(alpha+u(3)*xita)*I+beta*I^2;
s2 =g*(u(1)-u(2));
s = [s1; s2;(delta*u(3)-miu*u(3)^3+alpha*I/(h*omg)+beta*I^2/(2*h*omg))];
function u0 = pdex5ic(x)
u0 = [300; 300;10^16];
function [pl,ql,pr,qr] = pdex5bc(xl,ul,xr,ur,t)
pl = [0;0;0];
ql = [1;1;1];
pr = [0;0;0];
qr = [1;1;1]