matlab 0-1规划求解指派问题
指派问题可以用匈牙利算法，嫌代码太长，用0-1整数规划，matlab求解很方便，另外lingo也特别方便
参照问题如下：
用0-1整数规划可建立如下模型，截的图排版有些，，，，但不影响：
matlab代码如下，我码了这么多注释纯属心情太好，业界良心！！！（对于参加数模懂编程但又不太熟练matlab的同学非常好）：
function [y,fval]=zhipai(C)%C为指派n*n系数矩阵
[m,n]=size(C);
Aeq=zeros(2*n,n*n);
Aeq(1:n,1+(i-1)*n:i*n)=eye(n,n);
Aeq(i+n,1+(i-1)*n:i*n)=ones(1,n);
beq=ones(2*n,1);
lb=zeros(n*n,1);
ub=ones(n*n,1);
x=linprog(f',[],[],Aeq,beq,lb,ub);
y=reshape(x,n,n);
y=round(y);
sol=zeros(n,n);
if y(i,j)==1
sol(i,j)=C(j,i);
fval=sum(sol(:));
一会上传一个匈牙利算法的M文件，想要的可以下。妈的要睡觉了，刚歇了一天，明儿又要为期四天的数模训练了，，，，，，
没有更多推荐了，
加入CSDN，享受更精准的内容推荐，与500万程序员共同成长！ 上传我的文档
 下载
 收藏
粉丝量：19
该文档贡献者很忙，什么也没留下。
 下载此文档
MATLAB习题
下载积分：30
内容提示：MATLAB习题
文档格式：PDF|
浏览次数：119|
上传日期： 07:18:04|
文档星级：
全文阅读已结束，如果下载本文需要使用
 30 积分
下载此文档
该用户还上传了这些文档
MATLAB习题
关注微信公众号只需一步，快速开始
扫一扫，访问微社区
查看: 3987|回复: 16|关注: 0
含参数一元高次方程求解问题
<h1 style="color:# 麦片财富积分
新手, 积分 10, 距离下一级还需 40 积分
请各位高手赐教！
论坛优秀回答者
<h1 style="color:#72 麦片财富积分
关注者： 353
你试过solve了吗？
<h1 style="color:# 麦片财富积分
你试过solve了吗？
& &solve('r^3-18*a*r^2+54*a^2*r-36*a^3=0','r')
& && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && &6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
6*a - ((-648*a^6)^(1/2) + 72*a^3)^(1/3)/2 - (9*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (3^(1/2)*(((-648*a^6)^(1/2) + 72*a^3)^(1/3) - (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3))*i)/2
6*a - ((-648*a^6)^(1/2) + 72*a^3)^(1/3)/2 - (9*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3) - (3^(1/2)*(((-648*a^6)^(1/2) + 72*a^3)^(1/3) - (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3))*i)/2
&& simple(6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3))
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
simplify(Steps = 100):
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
combine(sincos):
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
combine(sinhcosh):
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
combine(ln):
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
(6*a*((-648*a^6)^(1/2) + 72*a^3)^(1/3) + ((-648*a^6)^(1/2) + 72*a^3)^(2/3) + 18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
rewrite(exp):
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
rewrite(sincos):
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
rewrite(sinhcosh):
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
rewrite(tan):
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
mwcos2sin:
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
collect(a):
(18/((-648*a^6)^(1/2) + 72*a^3)^(1/3))*a^2 + 6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3)
6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
<h1 style="color:# 麦片财富积分
& &solve('r^3-18*a*r^2+54*a^2*r-36*a^3=0','r')
solve&&表示没有做出上述结果啊
论坛优秀回答者
<h1 style="color:#72 麦片财富积分
关注者： 353
没看懂你贴的这些代码，这是我求的结果：
syms a r
solve('r^3-18*a*r^2+54*a^2*r-36*a^3=0','r')复制代码
RootOf(z^3 - 18*a*z^2 + 54*a^2*z - 36*a^3, z)[1]
RootOf(z^3 - 18*a*z^2 + 54*a^2*z - 36*a^3, z)[2]
RootOf(z^3 - 18*a*z^2 + 54*a^2*z - 36*a^3, z)[3]
论坛优秀回答者
<h1 style="color:#72 麦片财富积分
关注者： 353
把MaxDegree设置为3可以解出具体的结果：
syms a r
sol=solve(r^3-18*a*r^2+54*a^2*r-36*a^3==0,r,'MaxDegree',3)复制代码
& && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && && & 6*a + ((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3)
6*a - ((-648*a^6)^(1/2) + 72*a^3)^(1/3)/2 - (9*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3) - (3^(1/2)*(((-648*a^6)^(1/2) + 72*a^3)^(1/3) - (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3))*1i)/2
6*a - ((-648*a^6)^(1/2) + 72*a^3)^(1/3)/2 - (9*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3) + (3^(1/2)*(((-648*a^6)^(1/2) + 72*a^3)^(1/3) - (18*a^2)/((-648*a^6)^(1/2) + 72*a^3)^(1/3))*1i)/2
至于你后面贴的那一些东西应该还要用到其他的操作吧
<h1 style="color:# 麦片财富积分
没看懂你贴的这些代码，这是我求的结果：
能把您的代码附上吗？拜托了
论坛优秀回答者
<h1 style="color:#72 麦片财富积分
关注者： 353
能把您的代码附上吗？拜托了
我已经把代码贴上了啊
<h1 style="color:# 麦片财富积分
把MaxDegree设置为3可以解出具体的结果：
但是还是没有求出题目中的结果啊& &我算的结果也是这样子的
<h1 style="color:# 麦片财富积分
我已经把代码贴上了啊
不好意思刚看到啊& &&&我们仔细分析一下。确实可以解得题目中的结果& && &将结果带入后次数是一致的& &&&能够保持证方程的成立& &可是该如何操作呢？
站长推荐 /3
筑起功能安全的堡垒 - 基于模型设计的软件开发
MATLAB中文论坛是全球最大的 MATLAB & Simulink 中文社区。用户免费注册会员后，即可下载代码，讨论问题，请教资深用户及结识书籍作者。立即注册加入我们吧！
MATLAB官方社交平台
MATLAB中文论坛微社区与《》相关：
- Copyright & 2017 www.xue63.com All Rights Reserved24小时热门版块排行榜&&&&
(小有名气)
在线: 76.3小时
虫号: 2705277
注册: 性别: GG专业: 电化学
Matlab习题求解已有1人参与
小弟不才，Matlab的程序代码出现问题，求大神给出代码模板或者此题程序代码。期末试题啊，求给力！感激涕零！
已知水在不同温度下的密度如下表所示,请用一元三点拉格朗日插值法计算水在12,33和66℃时的密度。
T（℃）& & & && & 0& & & && &20& & & && &40& & & && &60& & & && &80& & & &&&100
ρ（kg/m3）999.9& & & & 998.2& & & & 992.2& & & & 983.2& & & & 971.8& & & & 958.4
& 猜你喜欢
已经有1人回复
已经有10人回复
已经有80人回复
已经有4人回复
已经有5人回复
已经有2人回复
已经有0人回复
已经有0人回复
已经有133人回复
已经有116人回复
已经有0人回复
& 本主题相关价值贴推荐，对您同样有帮助:
已经有5人回复
已经有19人回复
已经有6人回复
已经有8人回复
已经有12人回复
已经有7人回复
已经有10人回复
已经有11人回复
已经有7人回复
已经有6人回复
已经有7人回复
已经有7人回复
已经有17人回复
已经有8人回复
已经有9人回复
已经有14人回复
已经有15人回复
已经有5人回复
已经有23人回复
已经有115人回复
已经有6人回复
已经有6人回复
已经有3人回复
已经有5人回复
已经有7人回复
已经有20人回复
已经有18人回复
已经有4人回复
已经有11人回复
翅膀的痕迹
(职业作家)
散金: 1853沙发: 4
在线: 2787.4小时
虫号: 1122189
注册: 专业: 宇宙学管辖:
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ 感谢参与，应助指数 +1杨幸遇: 金币+15, ★★★★★最佳答案
CODE:function Lag_interp
Ti=[12 33 66];
T=[ 0& && && & 20& && && & 40& && && & 60& && && & 80& && && &100];
p=[999.9& && &&&998.2& && &&&992.2& && &&&983.2& && &&&971.8& && &&&958.4];
for i=1:length(Ti);
& & p_answer(i)=polyinterp(T,p,Ti(i));
Result=[Ti' p_answer']
function v=polyinterp(x,y,u)
n=length(x);
v=zeros(size(u));
for k=1:n;
& & w=ones(size(u));
& & for j=[1:k-1 k+1:n];
& && &&&w=(u-x(j))./(x(k)-x(j)).*w;
v=v+w*y(k);
end计算结果：
(小有名气)
在线: 76.3小时
虫号: 2705277
注册: 性别: GG专业: 电化学
引用回帖:: Originally posted by 月只蓝 at
function Lag_interp
for i=1:length(Ti);
& & p_answer(i)=polyinterp(T,p,Ti(i));
function v=polyinterp(x,y,u)
n=length(x);
v=zeros(size(u)) ... 万分感谢！金币送上！
& 本帖已获得的红花（最新10朵）
翅膀的痕迹
(小有名气)
在线: 76.3小时
虫号: 2705277
注册: 性别: GG专业: 电化学
送红花一朵
引用回帖:: Originally posted by 杨幸遇 at
万分感谢！金币送上！... 大神 我可否再问你一道题啊，实在是对MATLAB编程无感，只为完成作业。
翅膀的痕迹
相关版块跳转
第一性原理
我要订阅楼主
的主题更新
小木虫,学术科研互动社区,为中国学术科研免费提供动力
违规贴举报删除请发送邮件至：
广告投放与宣传请联系 李想 QQ：
QQ：&&邮箱：
Copyright &
MuChong.com, All Rights Reserved. 小木虫 版权所有}