syms Ct
syms G
function f=myfun(Ct,G);
sum=0;
for sr=1:14;e=1:14;p=1:14;gr=1:14;d=1:14;
GR=Grb(gr,1);
SR=Srb(sr,1);
E=EV(e,1);
P=PV(p,1);
D=DAM(d,1);
if P>E
sum1=(P-E)SR;
G=CtE+(sum1-DE);
s=diff(f); %% 一阶导数
s2=diff(f,2); %% 二阶导数
h=double(solve(s); %% 一阶导数为零点
for i=1:length(h)
if subs(s2,x,h(i))<0
disp(num2str(h(i)
disp(num2str(subs(G,x,h(i)
end
end
end
if E>P
sum2=(E-P)GR;
G=CtE-(sum2+DE);
s=diff(f); %% 一阶导数
s2=diff(f,2); %% 二阶导数
h=double(solve(s); %% 一阶导数为零点
for i=1:length(h)
if subs(s2,x,h(i))<0
disp(num2str(h(i)
disp(num2str(subs(G,x,h(i)
end
end
end
end
end
刚学 不要嫌弃
输出一下
disp(num2str(h(i)
disp(num2str(subs(G,x,h(i)