非线性规划示例(fmincon)
clc,clear;
C=0.85;%流量系数
delta_p=60;%压强变化量
E=2786.4;%弹性模量
A=0.7^2*pi;%小孔的面积
rho1=0.85;%变化前的燃油密度
rho=E*rho1/(E-delta_p);%变化后的燃油密度
Q=C*A*sqrt(2*delta_p/rho);%进出高压流管的流量
V2=440;%出油量
x0=0;
[t,fval,exitflag,output]=fmincon(@(t)(1000/(t+10)*Q*t-V2)^2,x0,[],[],[],[],[],[],@mycon)
function[c,ceq]=mycon(t)
C=0.85;%流量系数
delta_p=60;%压强变化量
E=2786.4;%弹性模量
A=0.7^2*pi;%小孔的面积
rho1=100;%变化前的燃油密度
rho=E*rho1/(E-delta_p);%变化后的燃油密度
Q=C*A*sqrt(2*delta_p/rho);%进出高压流管的流量
V2=440;%出油量
c=1000/(t+10)*Q*t-V2-1;
c=-1000/(t+10)*Q*t-V2-1;
ceq=[];
end