数学建模社区-数学中国

标题: 求解两道题，fmincon函数与lsqcurvefit函数 [打印本页]
作者: 子钩日月    时间: 2013-9-4 01:04
标题: 求解两道题，fmincon函数与lsqcurvefit函数
编程求解下列最小值问题, v- W* U) X# W8 V8 p! B6 Y
min 4*x1^3-a*x1-2*x2,
- d1 O5 M. Q; ]* B( g( t8 m+ m4 e6 r4 c: x3 S( U
                    s.t.     x1+x2<=4,* C6 _2 f! L: @: O& y5 [
                              2*x1+x2<=5,
# z3 U2 `# v- G3 u: s; J                              -x1+b*x2>=2,
" r% Z7 L( {9 A' n1 A8 ~) f                               x1,x2>=0,
4 Y( i  n. ^7 d% `. s; D                              a = 0,1,2,3,4;b=2,4,6,7; V9 s5 K9 v% W! P+ I, Y4 U$ L
用Matlab命令randint(5,2,[0,10])生成5X2的随机矩阵，其中矩阵第1列的数据作为x的观测值，矩阵第2列的数据作为y对应的观测值，来拟合二次曲线方程
  @8 U6 B  a; i
               a*x^2+b*x*y+c*y^2=3
; Y* \2 w* r9 }* j" C# q8 n; D6 w2 A3 K9 U
此外我想问的是fmincon函数与lsqcurvefit函数的初值一般如何估计。。。
( y* g( k" k5 C/ c9 M$ ?8 T
作者: cuiyi    时间: 2015-7-10 16:01
没有人回复。。。。。。。。- _: \# a3 ?+ Y6 V$ V7 R) B+ j
作者: 森之张卫东    时间: 2015-7-11 12:34
第一问，求解结果。
3 L9 t; N8 H% B, J' R8 [8 Y& q

部分结果.PNG (18.89 KB, 下载次数: 437)

部分结果.PNG

bang_ren_15_7_11.m

256 Bytes, 下载次数: 1, 下载积分: 体力 -2 点

程序

fun1.m

259 Bytes, 下载次数: 1, 下载积分: 体力 -2 点

程序
作者: 森之张卫东    时间: 2015-7-11 12:41
fminimax初值一般选取zeros（n，1）；
# o* M0 K0 |4 T2 P) b' [4 f$ sc = lsqcurvefit ('fun', c0, xdata, ydata) ;
0 Q$ \- ^/ }# B' M2 U/ \' dc0 = c; %一般以计算出的 c为初值进行迭代;5 j+ e* P7 i$ X+ I$ T+ m
作者: 森之张卫东    时间: 2015-7-11 17:55
再回一贴，二题。1 R2 R1 }& d) x8 Y7 d# J! N

结果.PNG (26.9 KB, 下载次数: 428)

结果.PNG

bang_ren_15_7_11_2.m

150 Bytes, 下载次数: 1, 下载积分: 体力 -2 点

fun2.m

144 Bytes, 下载次数: 1, 下载积分: 体力 -2 点

ti_qu_han_shu.m

96 Bytes, 下载次数: 1, 下载积分: 体力 -2 点

先将隐函数提取出来，即是y=f（x）的形式



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5