全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
1476 2
2010-09-18
这是一个求极小值的问题
fun.m文件程序如下
    function f=fun(x)
    P1=102.6402;
    t1=3.6;
    c1=0.0352866;c2=0.377616;c3=-0.008852;c4=82.768353;c5=-3.594521;
   c6=0.0392739;c7=1;c8=1;
  
f=-(-844/2*log(2*Pi)-1/2*c7-1/2*c8-(844/2)*log(x(1)^2)-1/2*(x(1)^(-2))*(c1+(c2)*...
       x(2)*x(3)+(c3)*x(2)+c4*x(2)^2*x(3)^2+(c5)*x(2)^2*x(3)+c6*x(2)^2)-M/2*log(2*Pi)-M/2*log(x(4)^2)...
       -1/2*(x(4)^(-2))*((logP1-log(2*((x(2)+x(5))^2+2*x(1)^2)^(1/2)*exp((x(2)+x(5)+((x(2)+x(5))^2+...
       2*x(1)^2)^(1/2))*t1/2)/((x(2)+x(5)+((x(2)+x(5))^2+2*x(1)^2)^(1/2))*(exp(((x(2)+x(5))^2+2*x(1)^2)^(1/2)*t1)-1)...
       +2*((x(2)+x(5))^2+2*x(1)^2)^(1/2)))*2*x(2)*x(3)/(x(1)^2)-2*(exp(((x(2)+x(5))^2+2*x(1)^2)^(1/2)*t1)-1)/...
       ((x(2)+x(5)+((x(2)+x(5))^2+2*x(1)^2)^(1/2))*(exp(((x(2)+x(5))^2+2*x(1)^2)^(1/2)*t1)-1)+2*((x(2)+x(5))^2+2*x(1)^2)...
       ^(1/2))*t1)^2));
>> x=[1,1,1,1,1];
>> x=fminsearch('fun',x)
??? Capitalized internal function Pi; Caps Lock may be on.

Error in ==> E:\Program Files\Matlab\work\fun.m
On line 14   ==> f=-(-844/2*log(2*Pi)-1/2*c7-1/2*c8-(844/2)*log(x(1)^2)-1/2*(x(1)^(-2))*(c1+(c2)*...

Error in ==> E:\Program Files\Matlab\toolbox\matlab\funfun\fminsearch.m
On line 125   ==> fv(:,1) = feval(funfcn,x,varargin{:});
错在哪里哪位高手帮忙改一下
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2010-9-18 10:55:55
Pi-->pi
logP1-->log(P1)
M=?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-9-18 19:36:24
能拜托你把Pi换成pi吗?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群