| 查看: 959 | 回复: 1 | ||
[求助]
帮忙看一下程序解方程部分老说越界，不能正常运行 已有1人参与
|
|
这是主程序，想把sctrt子程序的结果注入sctrr中求出信号功率等画图，但是出现问题，??? Attempted to access tspan(1); index out of bounds because numel(tspan)=0. Error in ==> dde23 at 100 t0 = tspan(1); Error in ==> ctr at 82 solt=dde23(@sctrt,xt0,[ts,tfs],[],p); 下面附上子程序：function xdot=sctrt(x,p) i=p.it;%DFB的偏置电流； %********transmitter laser with all-optical feedback********* dnt=i/(p.V*p.q)-1/p.tn*x(1)-x(2)*p.Gn*(x(1)-p.n0)*(1+p.eps*x(2))^-1; dst=1/p.tn*p.gm*p.yt*x(1)-1/p.tpt*x(2)+x(2)*p.Gn*(x(1)-p.n0)*((1+p.eps*x(2))^-1)*p.yt; dft=1/2*p.af*(p.Gn*(x(1)-p.n0)*((1+p.eps*x(2))^-1)*p.yt-1/p.tpt); xdot=[dnt; dst; dft]; 另一个子程序：function xdot=sctrr(x,p,st) kcj=p.kinj/p.tin; if t-p.tj>=0 [xt,xpt]=deval(st,t-p.tj); else xt=[0,0,0]; end %*receiver laser with all-optical feedback and transmitter light injection* g2=(1-p.epss*x(2)^2-p.epc*x(4)^2); g4=(1-p.epss*x(4)^2-p.epc*x(2)^2); dnr=I/(p.q*p.Vs)-x(1)/p.tns-p.vg*p.Gs*(x(1)-p.n0s)*(g2*x(2)^2+g4*x(4)^2); dsr=1/2*(p.yts*p.vg*p.Gs*(x(1)-p.n0s)*g2-1/p.tps)*x(2)+kcj*xt(2)*cos((x(3)-xt(3))-p.shixie*t); dfr=1/2*p.as*(p.yts*p.vg*p.Gs*(x(1)-p.n0s)*g2-1/p.tps)-2*pi*(p.miucw-p.vs)-kcj*xt(2)/x(2)*sin((x(3)-xt(3))-p.shixie*t); dsr2=1/2*(p.yts*p.vg*p.Gs*(x(1)-p.n0s)*g4-1/p.tps)*x(4); dfr2=1/2*p.as*(p.yts*p.vg*p.Gs*(x(1)-p.n0s)*g4-1/p.tps)-2*pi*(p.miuccw-p.vs); xdot=[dnr; dsr; dfr; dsr2; dfr2]; 请大家帮我看看如何正确得出结果。如果能出来正确的结果愿意全部的金币奉上 [ Last edited by xiegangmai on 2013-5-9 at 15:25 ] |
回复此楼» 猜你喜欢
289求调剂
已经有4人回复
-
292求调剂
已经有4人回复
-
080500求调剂
已经有3人回复
-
291求调剂
已经有3人回复
-
材料292调剂
已经有8人回复
-
一志愿211 初试270分 求调剂
已经有6人回复
-
求调剂
已经有3人回复
-
274求调剂
已经有4人回复
-
300求调剂，材料科学英一数二
已经有5人回复
-
招08考数学
已经有15人回复
2楼2017-03-23 14:32:14
