matlab求解二次方程组clearformatsyms Ls Lw Ts Tb Te Tc Tns Tn Tp Q P

matlab求解二次方程组

clear

format

syms Ls Lw Ts Tb Te Tc Tns Tn Tp Q Ps1 P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 Ps2 S

k=1.38e-23;%波尔兹曼常数

T=300;%温度

e=1.6e-19;

m=0.067*9.108e-31;%电子有效质量

%m=0.5*9.108e-31;%空穴有效质量

E=0.1*1.6e-19;%势垒高度

Lb=10e-9;%垒层高度

Ls=100e-9;%SCH层高度

Lw=5e-9;%阱层厚度

Ts=45e-12;%载流子穿越SCH层时间 10ps 10ps

Tb=0.1e-12;%载流子扩散势垒层时间

Te=(2*pi*m*Lw*Lw/(k*T))^0.5*exp(E/k/T);

Tc=20e-12;%载流子阱间输运时间 1e-10

Tns=10e-9;%载流子SCH层复合寿命 10ns

Tn=10e-9;%载流子量子阱层复合寿命 10ns

M=10;

Tw=0.03;

vg=3e18;

eslo=1e-17;

T=M*Tw;

G0=150000;A=1e8;B=1e-16;C=1e-41;N=1e19;N0=1e18;

g=G0*log((A*N+B*N^2+C*N^3)/(A*N0+B*N0^2+C*N0^3));

w=1;%元频率=1

%S=1e22;%e22量级

%G=M*Tw*vg*g*(1-eslo*T*S)*S

%3e23*(1-3e-18*S)*S

f1=-i*w*Ps1+Q-(1/Ts+1/Tns)*Ps1+(P1/Te)*(Lw/Ls)-3e23*(1-3e-18*S)*S;

f2=-i*w*P1+(Ps1/Ts)*(Ls/Lw)-(P1-P2)/Tc-(1/Tn+1/Te)*P1-3e23*(1-3e-18*S)*S;

f3=-i*w*P2+(P1-P2)/Tc-(P2-P3)/Tc-P2/Tn-3e23*(1-3e-18*S)*S;

f4=-i*w*P3+(P2-P3)/Tc-(P3-P4)/Tc-P3/Tn-3e23*(1-3e-18*S)*S;

f5=-i*w*P4+(P3-P4)/Tc-(P4-P5)/Tc-P4/Tn-3e23*(1-3e-18*S)*S;

f6=-i*w*P5+(P4-P5)/Tc-(P5-P6)/Tc-P5/Tn-3e23*(1-3e-18*S)*S;

f7=-i*w*P6+(P5-P6)/Tc-(P6-P7)/Tc-P6/Tn-3e23*(1-3e-18*S)*S;

f8=-i*w*P7+(P6-P7)/Tc-(P7-P8)/Tc-P7/Tn-3e23*(1-3e-18*S)*S;

f9=-i*w*P8+(P7-P8)/Tc-(P8-P9)/Tc-P8/Tn-3e23*(1-3e-18*S)*S;

f10=-i*w*P9+(P8-P9)/Tc-(P9-P10)/Tc-P9/Tn-3e23*(1-3e-18*S)*S;

f11=-i*w*P10+(Ps2/Ts)*(Ls/Lw)+(P9-P10)/Tc-(1/Tn+1/Te)*P10-3e23*(1-3e-18*S)*S;

f12=-i*w*Ps2+-Ps2/Ts-Ps2/Tns+(P10/Te)*(Lw/Ls)-3e23*(1-3e-18*S)*S;

f13=-i*w*S+3e23*(1-3e-18*S)*S-S/Tp;

[Ps1,P1,P2,P3,P4,P5,P6,P7,P8,P9,P10,Ps2,S]=solve(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,f13,'Ps1','P1','P2','P3','P4','P5','P6','P7','P8','P9','P10','Ps2','S')

ps1=eval(Ps1);%

p1=eval(P1);%

p2=eval(P2);%

p3=eval(P3);%

p4=eval(P4);%

p5=eval(P5);%

p6=eval(P6);%

p7=eval(P7);%

p8=eval(P8);%

p9=eval(P9);%

p10=eval(P10);%

ps2=eval(Ps2);%

P=[ps1,p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,ps2]/1e6;

%P=[Ps1/1e6,P1/1e6,P2/1e6,P3/1e6,P4/1e6,P5/1e6,P6/1e6,P7/1e6,P8/1e6,P9/1e6,P10/1e6,P11/1e6,P12/1e6,Ps2/1e6]*1e6;

N=[1,2,3,4,5,6,7,8,9,10,11,12];

bar(N,P)

P

提示错误：

Error using ==> strcmp

Function 'strcmp' is not defined for values of class 'sym'.

Error in ==> D:MATLAB6p5toolboxmatlabspecgraphmakebars.m

On line 49 ==> if strcmp(varargin{nin},'3'),

Error in ==> D:MATLAB6p5toolboxmatlabspecgraphbar.m

On line 36 ==> [msg,x,y,xx,yy,linetype,plottype,barwidth,equal] = makebars(varargin{:});

Error in ==> D:MATLAB6p5work物构所量子阱tenqwdongtaiequ.m

On line 82 ==> bar(N,P)

该方程组是二次方程组,貌似超出matlab计算长度,谁能帮忙修改下,

fsolve也试过了,还是不行,谁能给个可运行程序啊

clear

P0 = [1e26;1e26;1e26;1e26;1e26;1e26;1e26;1e26;1e26;1e26;1e25;1e25;1e22]; % 初始点

options=optimset('Display','iter'); % 显示输出信息

[P,fval] = fsolve(@myfun2,P0,options)

P=[P(11);P(1);P(2);P(3);P(4);P(5);P(6);P(7);P(8);P(9);P(10);P(12)]

N=[1,2,3,4,5,6,7,8,9,10,11,12];

p=P/1e6;

bar(N,p);

1年前

悬赏10滴雨露 已收到1个回答