|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
程序为:
syms w1 m
xj=9; q=1.6*10^(-19); delta=10^(-4);
d1=1.5; d2=0.19; d3=0.4; d4=0.4;
N1=2.5*10^15; N2=2*10^16; N3=6*10^16; N4=2*10^15;
v=0.12; x=9:0.01:11.45;
e01=(q*N4*((xj+w1)^2-xj^2)/(2*delta))/xj
g1=int(e01*xj/m-q*N4*(m^2-xj^2)/(m*2*delta),m,xj,xj+w1)+v
w1=solve(g1,'w1');
e01=(q*N4*((xj+w1)^2-xj^2)/(2*delta))/xj;
y=(e01*xj./x-q*N4.*(x.^2-xj^2)./(x.*2*delta)).*(xj<=x&x<=xj+w1);
plot(x,y)
错误提示:
Warning: Explicit solution could not be found.
> In solve at 81
In f1 at 17
??? Error using ==> mupadmex
Error in MuPAD command: Either base or exponent must be a scalar.
Error in ==> sym.sym>sym.mpower at 198
B = mupadmex('mllib::mpower',A.s,p.s);
Error in ==> f1 at 18
e01=(q*N4*((xj+w1)^2-xj^2)/(2*delta))/xj;
[ 本帖最后由 ChaChing 于 2010-7-13 23:51 编辑 ] |
|