Hello everyone, I am doing graduation design need to use the matlab simulation, because never used before, there is no code, so self-study stage now and ready to start from the simplest scene simulation, I try to describe the problem clearly, so that you browse,
I is the ultimate goal of this topic simulation through iterative optimization of unmanned aerial vehicle (uav) weight distribution of horizontal position with users and make a uav communication model in the network to achieve maximum utility, the first thing I simulate 1 unmanned aerial vehicle (uav) 2 user scenario, except for unmanned aerial vehicle (uav) horizontal position (the independent variable x) and users (random) unknown outside the horizontal position, the rest of the variables (user weight, bandwidth, unmanned aerial vehicle (uav) height, etc.) all fixed, the following functions are the scene communication model:
F=0.5 * log10 (20 * 0.5 * 10 ^ 6 * log2 (1 + (10 ^ (4) * 1.42 * 10 ^ (4) * 2.9088)/(10 ^ (9.6) * (dist (rand (1, 2), x ') + 300 ^ 2)))) - 0.5 * log10 (20 * 0.5 * 10 ^ 6 * log2 (1 + (10 ^ (4) * 1.42 * 10 ^ (4) * 2.9088)/(10 ^ (9.6) * (dist (rand (1, 2), x ') + 300 ^ 2))));
Thinking is the minimum worth by fmincon toolbox to calculate the negative to a maximum of network utility, but finally ran out of the result is: the initial point is a local minimum that satisfies the constraints, no matter how my initial point set (0; 0 to 10 w; 10 w), the result is the initial point is the local optimal solution,
My own analysis possibility is that maybe too little number of users or choose network utility function is not suitable, so I add users to 10, and modify the network utility function (from the proportional fairness to the minimum potential delay fairness, avoid log10 logs base),
The following functions are an unmanned aerial vehicle (uav) with arbitrary channel model between 10 users:
The function f=fun (x)
F=symsum (0.5 * (- (1/(1/(n) * 20 * 10 ^ 6 * log2 (1 + (10 ^ (4) * 1.42 * 10 ^ (4) * 2.9088)/(10 ^ (9.6) * (dist (rand (1, 2), x) + 300 ^ 2)))))), n, 0, 10);
End

F network utility, n number of users, symsum this command I am query help manual, want to implement a network utility and calculate the; X is a horizontal coordinate vector (2 x1 matrix), fmincon each parameter Settings in the two problems as follows:
X0=[0, 0];
A=(1, 1]; B=[1000].
Aeq=[]; beq=[];
VLB=(1000-1000); VUB=(1000-1000);
[x, fval]=fmincon (@ fun, x0, A, b, Aeq, beq, VLB, VUB);

At the same time, I'm not sure whether it will set the qualification can be said has been one of the horizontal ordinate range within 1000,
After modifying function, the system told me undefined function or variable 'n', do it here I don't know the way of thinking is correct, and I don't know should how to turn this function with n and x as fmincon order to solve the optimization of parameters, so come first on the BBS for help once,

Thank you very much!