数学建模社区-数学中国

标题: 求解两道题，fmincon函数与lsqcurvefit函数 [打印本页]
作者: 子钩日月    时间: 2013-9-4 01:04
标题: 求解两道题，fmincon函数与lsqcurvefit函数
编程求解下列最小值问题1 M4 n( ?7 x7 F( {/ M0 r: F7 Q5 h+ \
min 4*x1^3-a*x1-2*x2,7 K6 e" }# y. ^" d0 P+ B
% z  Z. m1 ?# \/ u1 @5 f) i
                    s.t.     x1+x2<=4,
/ H) V9 D! P# G8 y3 ~0 b# E6 G; W                              2*x1+x2<=5,: U/ y4 a0 _/ f5 A: m8 A* p+ \
                              -x1+b*x2>=2,* }+ d, B& T6 Q( x+ E
                               x1,x2>=0,
' V/ l6 ?, g8 k9 C                              a = 0,1,2,3,4;b=2,4,6,7
$ o, i3 g  `* W' s2 K0 o/ d1 Y
用Matlab命令randint(5,2,[0,10])生成5X2的随机矩阵，其中矩阵第1列的数据作为x的观测值，矩阵第2列的数据作为y对应的观测值，来拟合二次曲线方程

8 {6 P) U' B4 J               a*x^2+b*x*y+c*y^2=3
, m6 F* c' M' A  T) F# j) _0 ~3 T9 z1 L* _
此外我想问的是fmincon函数与lsqcurvefit函数的初值一般如何估计。。。
5 X+ B. E8 T: I
作者: cuiyi    时间: 2015-7-10 16:01
没有人回复。。。。。。。。, B* [9 z8 q5 O" g
作者: 森之张卫东    时间: 2015-7-11 12:34
第一问，求解结果。& O5 v3 g' |# [9 F9 j0 B' y9 k

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

部分结果.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）；
' s& [" W( y# n9 _c = lsqcurvefit ('fun', c0, xdata, ydata) ;
" d1 t. Z" j- f% R  _( j3 Vc0 = c; %一般以计算出的 c为初值进行迭代;  L0 T* F/ j7 w) v/ V: E
作者: 森之张卫东    时间: 2015-7-11 17:55
再回一贴，二题。. y+ h+ S8 C. E/ u) u& q

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

结果.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