matlab pdefun,PDEFUN 的输出异常。，PDEFUN 必须返回长度为 3 的3哥lie'xiang

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

请问下各位大佬我这个的错误是怎么修改的，卡在这里了

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]