求助高手帮忙改一下matlab程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

这是一个求极小值的问题
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{:});
错在哪里哪位高手帮忙改一下

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MATLAB程序 MATLAB matla atlab 求助高手 求助 MATLAB 程序 高手 帮忙

Pi-->pi
logP1-->log(P1)
M=?

能拜托你把Pi换成pi吗？