版块导航
正在加载中...
客户端APP下载
论文辅导
调剂小程序
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
论坛更新日志
(1660)
>
考研
(276)
>
导师招生
(155)
>
虫友互识
(68)
>
招聘信息布告栏
(19)
>
休闲灌水
(15)
>
海外博后
(11)
>
硕博家园
(10)
>
考博
(10)
>
博后之家
(9)
>
基金申请
(9)
>
有机交流
(9)
>
教师之家
(6)
>
找工作
(6)
>
论文投稿
(6)
>
公派出国
(5)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
程序语言
»
MATLAB/Mathematica
»
谁能帮我编写个Matlab程序啊,急求。。。
北京石油化工学院2026年研究生招生接收调剂公告
18
2/2
返回列表
上一页
1
2
查看: 1952 | 回复: 17
只看楼主
@他人
存档
新回复提醒
(忽略)
收藏
在APP中查看
ybkooo
至尊木虫
(著名写手)
^^
应助: 32
(小学生)
金币: 11360
散金: 178
红花: 15
帖子: 2310
在线: 543.5小时
虫号: 2710519
注册: 2013-10-09
性别: GG
专业: 控制理论与方法
不共享一下吗?
回复此楼
abcd
11楼
2015-01-04 21:48:52
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
吃木头在虫子
新虫
(小有名气)
应助: 0
(幼儿园)
金币: 3444.9
帖子: 209
在线: 14小时
虫号: 3437721
注册: 2014-09-24
性别: GG
专业: 机械测试理论与技术
引用回帖:
11楼
:
Originally posted by
ybkooo
at 2015-01-04 21:48:52
不共享一下吗?
可以啊,怎么共享?据大婶说要matlab2012b以上版本才能运行,matlab7.0不行的
赞
一下
回复此楼
12楼
2015-01-05 10:30:06
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
NicoleLas
木虫
(正式写手)
应助: 66
(初中生)
金币: 4288.7
红花: 26
帖子: 514
在线: 177.2小时
虫号: 2872321
注册: 2013-12-15
性别: GG
专业: 计算机软件
吃木头在虫子: 回帖置顶
2015-01-06 14:04:57
我个人觉得先前的应助的结果有点问题,我没有MATLAB,使用Octave算了下,代码与MATLAB差不多。
代码:
CODE:
## Copyright (C) 2015 Leiz
##
## This program is free software; you can redistribute it and/or modify
## it under the terms of the GNU General Public License as published by
## the Free Software Foundation; either version 3 of the License, or
## (at your option) any later version.
##
## This program is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
## GNU General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with Octave; see the file COPYING. If not, see
## <http://www.gnu.org/licenses/>.
## demoSurf
## Author: Leiz <Leiz@LEIZ-PC>
## Created: 2015-01-05
clear all
clc
x1 = linspace( 0, 6, 50 );
x2 = linspace( 0, 6, 50 );
[ xx1, xx2 ] = meshgrid( x1, x2 );
tempx1 = xx1;
tempx2 = xx2;
tempx1( xx1 + xx2 > 6 ) = NaN;
tempx1( -xx1 + 2*xx2 > 2 ) = NaN;
tempx2( xx1 + xx2 > 6 ) = NaN;
tempx2( -xx1 + 2*xx2 > 2 ) = NaN;
fx = -2*tempx1 - 6*tempx2 + tempx1.^2 - 2*tempx1.*tempx2 + 2*tempx2.^2;
figure(1)
surf( tempx1, tempx2, fx );
h1 = get( gcf, "currentaxes" );
set( h1, "linewidth", 2 );
xlabel( 'x_{1}' );
ylabel( 'x_{2}' );
zlabel( 'f(x)' );
title( 'Surface of f(x) = -2x_{1} - 6x_{2} + x_{1}^{2} - 2x_{1}x_{2} + 2x_{2}^{2}' );
view( -60, 30 );
print( "demoSurf.png", "-dpng" );
figure(2)
contour( tempx1, tempx2, fx, 20 );
h2 = get( gcf, "currentaxes" );
set( h2, "linewidth", 2 );
axis( [ 0, 6, 0, 3 ] );
xlabel( 'x_{1}' );
ylabel( 'x_{2}' );
title( 'Contour of f(x) = -2x_{1} - 6x_{2} + x_{1}^{2} - 2x_{1}x_{2} + 2x_{2}^{2}' );
print( "demoContour.png", "-dpng" );
结果:
验证:
因为这是一个QP问题,Octave这方面求解的函数有一些不过我为了简便使用了LINGO来求解。
code:
CODE:
MIN = -2*X1 - 6*X2 + X1^2 - 2*X1*X2 + 2*X2^2;
X1 + X2 <= 6;
-X1 + 2*X2 <= 2;
solution:
CODE:
Global optimal solution found.
Objective value: -15.11111111082
Infeasibilities: 0.000000000000
Total solver iterations: 6
Elapsed runtime seconds: 0.02
Model is convex quadratic
Model Class: QP
Total variables: 2
Nonlinear variables: 2
Integer variables: 0
Total constraints: 3
Nonlinear constraints: 1
Total nonzeros: 6
Nonlinear nonzeros: 3
Variable Value Reduced Cost
X1 3.333333333380 0.3832140051194E-10
X2 2.666666666506 0.4779796153731E-10
Row Slack or Surplus Dual Price
1 -15.11111111082 -1.000000000000
2 0.1140998406868E-09 1.111111111122
3 0.3667064429891E-09 0.4444444448297
赞
一下
(2人)
回复此楼
» 本帖已获得的红花(最新10朵)
吃木头在虫子
Be water, my friend.
13楼
2015-01-05 14:31:41
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
ybkooo
至尊木虫
(著名写手)
^^
应助: 32
(小学生)
金币: 11360
散金: 178
红花: 15
帖子: 2310
在线: 543.5小时
虫号: 2710519
注册: 2013-10-09
性别: GG
专业: 控制理论与方法
引用回帖:
12楼
:
Originally posted by
吃木头在虫子
at 2015-01-05 10:30:06
可以啊,怎么共享?据大婶说要matlab2012b以上版本才能运行,matlab7.0不行的...
楼下的才是大神。我是菜鸟……
赞
一下
回复此楼
abcd
14楼
2015-01-05 15:55:28
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
吃木头在虫子
新虫
(小有名气)
应助: 0
(幼儿园)
金币: 3444.9
帖子: 209
在线: 14小时
虫号: 3437721
注册: 2014-09-24
性别: GG
专业: 机械测试理论与技术
送红花一朵
引用回帖:
13楼
:
Originally posted by
NicoleLas
at 2015-01-05 14:31:41
我个人觉得先前的应助的结果有点问题,我没有MATLAB,使用Octave算了下,代码与MATLAB差不多。
代码:
## Copyright (C) 2015 Leiz
##
## This program is free software; you can redistribute it and/o ...
我想知道这到底是什么玩意,太牛了,我都听不懂你在说什么,o(︶︿︶)o 唉,想哭啊
赞
一下
回复此楼
15楼
2015-01-06 14:04:01
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
NicoleLas
木虫
(正式写手)
应助: 66
(初中生)
金币: 4288.7
红花: 26
帖子: 514
在线: 177.2小时
虫号: 2872321
注册: 2013-12-15
性别: GG
专业: 计算机软件
内容已删除
回复此楼
Be water, my friend.
16楼
2015-01-06 14:38:50
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
feichedt
银虫
(小有名气)
应助: 3
(幼儿园)
金币: 709.1
散金: 10
红花: 3
帖子: 169
在线: 241.7小时
虫号: 1672940
注册: 2012-03-07
性别: GG
专业: 内流流体力学
楼主请问下能共享下代码么。想学习下,谢谢!
赞
一下
回复此楼
17楼
2015-01-06 17:58:17
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
--八云--
专家顾问
(知名作家)
博览众家之长,终集百家之短
专家经验: +143
应助: 216
(大学生)
金币: 4097.5
散金: 1408
红花: 70
帖子: 9216
在线: 1033.8小时
虫号: 1368331
注册: 2011-08-14
性别: GG
专业: 凝聚态物性 II :电子结构
管辖:
物理
看来又是一个被数学建模毒害的孩纸啊。。。。。。。
赞
一下
(1人)
回复此楼
人原本是不死的,后来他创造了死亡,故而人们歌颂他。
18楼
2015-01-06 18:16:16
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
相关版块跳转
第一性原理
量子化学
计算模拟
分子模拟
仿真模拟
程序语言
我要订阅楼主
吃木头在虫子
的主题更新
18
2/2
返回列表
上一页
1
2
如果回帖内容含有宣传信息,请如实选中。否则帐号将被全论坛禁言
普通表情
龙
兔
虎
猫
百度网盘
|
360云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站,将文件上传后,然后将下载链接复制到帖子内容中就可以了。
最具人气热帖推荐
[查看全部]
作者
回/看
最后发表
[
考研
]
081700,311,求调剂
+14
冬十三
2026-04-04
15/750
2026-04-06 01:57
by
BruceLiu320
[
考研
]
285求调剂
+7
AZMK
2026-04-04
9/450
2026-04-06 00:06
by
永字号
[
考研
]
085600,320分求调剂
+7
大馋小子
2026-04-01
8/400
2026-04-05 21:19
by
学员8dgXkO
[
考研
]
284求调剂
+7
徐同学_001
2026-04-04
13/650
2026-04-05 17:19
by
yulian1987
[
考研
]
312求调剂
+3
Say Never
2026-04-04
3/150
2026-04-05 11:19
by
guoweigw
[
考研
]
083200 333求调剂
+3
十二!!
2026-04-04
3/150
2026-04-05 08:28
by
barlinike
[
考研
]
302求调剂一志愿华中师范大学
+8
小江小江江江
2026-04-02
8/400
2026-04-04 19:50
by
蓝云思雨
[
考研
]
081200-11408-276学硕求调剂
+5
崔wj
2026-03-31
5/250
2026-04-04 19:45
by
1753564080
[
考研
]
363求调剂
+7
zh096
2026-04-04
7/350
2026-04-04 17:13
by
dongzh2009
[
考研
]
一志愿华中农业071010,总分320求调剂
+7
困困困困坤坤
2026-04-02
7/350
2026-04-03 17:26
by
Yuena_Wang
[
基金申请
]
请问共同通讯和共同一作的认可度问题
10
+4
psa1234
2026-04-01
10/500
2026-04-03 11:08
by
Kittylucky
[
考研
]
330求调剂
+3
白神呜呼呼
2026-04-02
3/150
2026-04-03 10:15
by
蓝云思雨
[
考研
]
279求调剂
+6
qazplm0852
2026-04-02
6/300
2026-04-03 10:03
by
蓝云思雨
[
考研
]
326求调剂
+3
9ahye
2026-04-02
4/200
2026-04-03 08:43
by
Jaylen.
[
考研
]
环境工程297分求调剂一志愿杭高院
+15
GENJIOW
2026-03-31
16/800
2026-04-02 17:56
by
cyh—315
[
考研
]
362求调剂
+14
西南交材料专硕3
2026-03-31
14/700
2026-04-02 17:50
by
yunlongyang
[
考研
]
286分调剂
+20
Faune
2026-03-30
22/1100
2026-04-02 13:24
by
clyblh
[
考研
]
复试调剂
+7
双马尾痞老板2
2026-03-31
7/350
2026-03-31 19:49
by
Dyhoer
[
考研
]
考研调剂求助
+7
13287130938
2026-03-31
7/350
2026-03-31 16:39
by
690616278
[
考研
]
262求调剂
+7
ZZ..000
2026-03-30
8/400
2026-03-31 10:05
by
cal0306
信息提示
关闭
请填处理意见
关闭
确定