楼主: 主啊上帝
2206 1

[问答] 求教用fmincon解方程组~~ [推广有奖]

  • 3关注
  • 0粉丝

大专生

53%

还不是VIP/贵宾

-

威望
0
论坛币
0 个
通用积分
0.0774
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
251 点
帖子
28
精华
0
在线时间
86 小时
注册时间
2014-12-3
最后登录
2017-4-2

楼主
主啊上帝 发表于 2015-7-22 19:49:06 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
最近在解这样一个方程组,如图片所示:
GD_$K04RJ@[Y2I`0H8H9YK5.png
限制条件是J个肉的和等于一个定值(就是数据中的一个变量sh(1)),每个肉都大于零,图片中橘黄色的是未知数,因为有很多方程组,每个的方程个数都不一样,我暂时只用了J等于2的一个方程组(数据就是n7.mat,只有两条观测值)编程,因为其他的方程组的J是不一样的,所以为了以后方便做,我没有直接写2,而用了s来代替,看能不能用fmincon做出来,我做的程序如下:
目标函数:

function  f=myfun(x,phi_fht,a1,b1)
    load n7.mat
    s=size(n7,1)
    phi_fht=n7(:,106)
    a=n7(:,108)
    a1=a(1,1)
    b=n7(:,107)
    b1=b(1,1)
    phi_ft=x(1,1)*ones(s,1)
     rou=x(2:s+1,1)
     D=rou*rou'
     d=diag(D)
     y=phi_fht-phi_ft-a1*rou-b1*d
    f=y'*y

其他的:
sh=n7(:,111)
  Aeq=[0,ones(1,s)]
  Beq=sh(1,1)
  lb=zeros(s+1,1)
  x0=0.5*ones(s+1,1)
    x= fmincon(@myfun,x0,[],[],Aeq,Beq,lb)

结果matlab报错:
WF`UUKCN%[{V}~%CFQRJ]MP.png
有没有大神能帮我看看这是怎么回事??是向量之间不能减吗,还是这个fmincon不能这么做?
万分感谢!!不好意思,我是菜鸟,并没有积蓄可以给予奖励,但求好心人帮帮忙!
二维码

扫码加我 拉你入群

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

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

关键词:fmincon 解方程组 Inc 解方程 FMI fmincon 方程组 未知数 橘黄色 图片 程序

沙发
主啊上帝 发表于 2015-7-23 08:47:03
有没有大神能帮帮我,到底是哪里错了~~~~自己顶下!

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2026-1-2 12:08