我利用MATLAB中的fmincon函数求非线性规划求电机转矩的最大值

Question

但无论我的转速变为何值，所得的最优解都不变，想了很久改变了初值等都好像没有用，自己感觉程序也是对的，大家看一下我的程序是不是哪里错了？
z=[1;2;3;110];
for i=1:1:4
we=z(i)

%%%%%%
x0=[-219,339];%在极限点时的dq电流值

[x,fev]=fmincon(@Tc,x0,[],[],[],[],[-285*sqrt(2),0],[0,285*sqrt(2)],@(x) yueshutiaojian(x,we))

end
%%%%%%%%
function [g,ceq]=yueshutiaojian(x,we)

g=[(0.33e-3*x(2))^2-(0.13e-3*x(1)+0.062)^2+(185/200)^2;x(1)^2+x(2)^2-2*285^2];
ceq=[];
end

%%%%%%%%%%%%%%
function f = Tc( x )
f=-6*(0.062+(0.13e-3-0.33e-3)*x(1))*x(2);

end
其中we就是转速，但不论转速给多少最优解算出来的转矩都不变，这是为什么？

Answer 1

使用匿名函数句柄就可以了
fmincon 需要的函数fun只有一个输入参数 fun(x)

而你的函数需要多一个能改变的变量 那么形式变成 myfun(p,x)
p就是你要传递的参数，x还是原来的输入

当然首先，你得写个m文件，假设为myfun.m，其内容大致如下
function [out]=myfun,(p,x)
...........
end
就是根据参数p和输入x，求的输出，当然名字可以自己取，然后就是主程序了

for p=p0:pn %参数p的取值
[x,fval]=fmincon(@(x) myfun(p,x),Bt0,[],[],[],[],lb,ub,@constr);

.......
end

其中关键就是@(x) myfun(p,x)
这会返回一个匿名的函数句柄，会使用当时的p值代入myfun
得到一个只有一个输入x的函数，这样就能把额外的参数传递给需要优化的函数
如果你的函数有等多的参数例如
myfun(a,b,c,x)
a=..
b=..
c=..
同样用@(x) myfun(a,b,c,x)把a，b，c参数的值传给要优化的函数

在以上的基础上，加几句，存储每次的参数p和对应的优化结果
然后画图就可以了追问

这个我知道呢，你看一下我的程序就是这样做的啊，可是如果求最小值就没有问题，一旦求转矩最大值，运行出来的结果不变