代码搜索结果

/****************************************Copyright (c)************************************************** ** 广州周立功单片机发展有限公司 ** 研 究

function t=domen0(a,b) %[a,b]为初始区间 t:极小值点 e=0.001; a0=a; b0=b; F0=(b0-a0)/e; F(1)=1; F(2)=1; for i=3:inf F(i)=F(i-1)+F(i-2); if(F(i)>=F0) n=i;

function t=domen0(a,b) %[a,b]为初始区间 t:极小值点 clc; e=0.001; a0=a; b0=b; lamda=0.618033989; rafa1=a0+(1-lamda)*(b0-a0); rafa2=a0+lamda*(b0-a0); y1=fai(rafa1); y2=fai(rafa2)

function a = setx0(a) global Bus DAE PV Syn Line jay if ~a.n, return, end % reset transmission line reactance and admittance matrix Line = addxl(Line,a.line,a.u.*a.xcs); if sum(a.u), Line = build_y

function a = setx0(a) global DAE PQ if ~a.n, return, end V1 = DAE.y(a.vbus); for k = 1:a.n idx = findbus(PQ,a.bus(k)); if isempty(idx) warn(a,k,'No PQ load found.') a.con(k,6) = 0;

function a = setx0(a) global Bus DAE Settings Wind if ~a.n, return, end check = 1; Pc = Bus.Pg(a.bus); Qc = Bus.Qg(a.bus); Vc = DAE.y(a.vbus); ac = DAE.y(a.bus); vds = -Vc.*sin(ac); vqs = Vc.*co

function a = setx0(a) global DAE if ~a.n, return, end % variable initialization DAE.x(a.Vs) = 0; Vs = DAE.x(a.Vs); % check time constants idx = find(a.con(:,4) == 0); if idx warn(a,idx,'Time con