<?xml version="1.0" encoding="UTF-8"?>
<rss version="2.0">
  <channel>
    <link>http://www.mathworks.com/matlabcentral/newsreader/view_thread/261029</link>
    <title>MATLAB Central Newsreader - help for constrained optimization  -using fmincon</title>
    <description>Feed for thread: help for constrained optimization  -using fmincon</description>
    <language>en-us</language>
    <copyright>&amp;copy;1994-2012 by MathWorks, Inc.</copyright>
    <webmaster>webmaster@mathworks.com</webmaster>
    <generator>MATLAB Central Newsreader</generator>
    <docs>http://blogs.law.harvard.edu/tech/rss</docs>
    <ttl>60</ttl>
    <image>
      <title>MathWorks</title>
      <url>http://www.mathworks.com/images/membrane_icon.gif</url>
    </image>
    <item>
      <pubDate>Thu, 17 Sep 2009 05:28:02 -0400</pubDate>
      <title>help for constrained optimization  -using fmincon</title>
      <link>http://www.mathworks.com/matlabcentral/newsreader/view_thread/261029#680734</link>
      <author>Ramana Murthy</author>
      <description>Hi &lt;br&gt;
I am trying to minimize a function (y) with linearity constraint. &lt;br&gt;
I was using the fmincon and I'm not getting desired results. Terminating message is coming.&lt;br&gt;
So I also passed on gradient also. But now program is showing error. Can someone tell me where is the problem.&lt;br&gt;
&lt;br&gt;
with regards,&lt;br&gt;
murthy&lt;br&gt;
&lt;br&gt;
******************&lt;br&gt;
clear all&lt;br&gt;
clc&lt;br&gt;
global lk sk;&lt;br&gt;
lk = [52,81,77,76,16,106,48,99,66,79];&lt;br&gt;
sk = [0.39,0.56,1.03,0.47,0.93,1.20,1.74,1.86,1.23,0.54];&lt;br&gt;
K = 10;&lt;br&gt;
R = 350;&lt;br&gt;
x0 = lk/4;&lt;br&gt;
% f = @(x)parameterfun(x,sk,lk);&lt;br&gt;
A = [];&lt;br&gt;
b = [];&lt;br&gt;
Aeq = ones(1,K);&lt;br&gt;
beq = R;&lt;br&gt;
lb = zeros(1,K);&lt;br&gt;
ub = lk;&lt;br&gt;
options = optimset('Algorithm','active-set');&lt;br&gt;
options = optimset(options,'GradObj','on');&lt;br&gt;
tic&lt;br&gt;
[x,fval]=fmincon(@parameterfun,x0,A,b,Aeq,beq,lb,ub,[],options);&lt;br&gt;
[x' rk']&lt;br&gt;
toc&lt;br&gt;
**********&lt;br&gt;
function [y,g] = parameterfun(x)&lt;br&gt;
global lk sk;&lt;br&gt;
y = 0;&lt;br&gt;
for i = 1:numel(x)&lt;br&gt;
&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;y = y + sk(i)*exp(-2*pi/sk(i))*(exp(2*pi*x(i)/(lk(i)*sk(i)))-1);&lt;br&gt;
end&lt;br&gt;
if nargout  &amp;gt; 1&lt;br&gt;
&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;g = 0;&lt;br&gt;
&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;for i = 1:numel(x)&lt;br&gt;
&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;g = g + 2*pi/lk(i)*exp(-2*pi/sk(i))*exp(2*pi*x(i)/(lk(i)*sk(i)));&lt;br&gt;
&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;end&lt;br&gt;
end&lt;br&gt;
**************</description>
    </item>
    <item>
      <pubDate>Thu, 17 Sep 2009 06:29:04 -0400</pubDate>
      <title>Re: help for constrained optimization  -using fmincon</title>
      <link>http://www.mathworks.com/matlabcentral/newsreader/view_thread/261029#680738</link>
      <author>Bruno Luong</author>
      <description>&quot;Ramana murthy&quot; &amp;lt;omurthy@yahoo.com&amp;gt; wrote in message &amp;lt;h8shd2$sb6$1@fred.mathworks.com&amp;gt;...&lt;br&gt;
&amp;gt; Hi &lt;br&gt;
&amp;gt; I am trying to minimize a function (y) with linearity constraint. &lt;br&gt;
&amp;gt; I was using the fmincon and I'm not getting desired results. Terminating message is coming.&lt;br&gt;
&amp;gt; So I also passed on gradient also. But now program is showing error. Can someone tell me where is the problem.&lt;br&gt;
&amp;gt; &lt;br&gt;
&lt;br&gt;
The gradient is a VECTOR. I see you compute some sort of scalar, so it must NOT be a gradient.&lt;br&gt;
&lt;br&gt;
Bruno</description>
    </item>
  </channel>
</rss>

