数学建模社区-数学中国

标题: 求解两道题,fmincon函数与lsqcurvefit函数 [打印本页]

作者: 子钩日月    时间: 2013-9-4 01:04
标题: 求解两道题,fmincon函数与lsqcurvefit函数
编程求解下列最小值问题
% e' e" n* k0 S+ H0 pmin 4*x1^3-a*x1-2*x2,3 M% c' e! ]9 N
" T/ ]- _9 \4 k+ N" ~; `$ A
                    s.t.     x1+x2<=4,
; R0 ^9 {5 l( q) A1 I* k/ w                              2*x1+x2<=5,9 {6 n: _. S2 `4 t2 n, D$ y
                              -x1+b*x2>=2,
" a9 n1 g2 b; E+ O9 b                               x1,x2>=0,% T- n$ d2 o% G( t( K9 n" B
                              a = 0,1,2,3,4;b=2,4,6,7
3 O, u% _; Z, G' s0 L, O* Z' U  c6 v
Matlab命令randint(5,2,[0,10])生成5X2的随机矩阵,其中矩阵第1列的数据作为x的观测值,矩阵第2列的数据作为y对应的观测值,来拟合二次曲线方程

: s# t9 U1 h! k               a*x^2+b*x*y+c*y^2=3& v$ E9 n6 L; C9 c: Z; q
/ I$ I+ j3 ^, H; F) K$ g/ ]
此外我想问的是fmincon函数与lsqcurvefit函数的初值一般如何估计。。。4 x5 B) w0 I, P# m

作者: cuiyi    时间: 2015-7-10 16:01
没有人回复。。。。。。。。5 m! ]$ W+ x5 U9 o) ?8 ~/ }& i/ ?7 y

作者: 森之张卫东    时间: 2015-7-11 12:34
第一问,求解结果。
9 B" H" p0 O2 b( i! j6 ]+ k

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

部分结果.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);
: N( g4 [  S) H$ l( Tc = lsqcurvefit ('fun', c0, xdata, ydata) ;) ^$ G" e8 `8 ], s/ w! u
c0 = c; %一般以计算出的 c为初值进行迭代;) }, m* S0 S7 T) h% t1 q* c

作者: 森之张卫东    时间: 2015-7-11 17:55
再回一贴,二题。. P: p- E. _! c$ y& K8 E

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

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