数学建模
CAMlive
炼丹术士
展开
-
生产任务分配问题 matlab+lingo
A题 生产任务分配问题有一个生产管理人员决定下个时期生产10批产品,他现在需要决定怎样将这些生产任务分配到其管理的车间中的不同机器上。此车间内的5台可用机器都可以加工这些产品,但由于它们是在不同时期生产的,因此其加工速度各不相同。此外,由于存在维护和调整期,在规划的时期内每台机器只能工作有限时间(小时数)。表1列出了所有机器加工各批产品所需的时间以及此机器可用的工时数。各批次产品的生产成本取决于处理此批次的机器。每台机器每小时运行成本取决于其采用的技术,此机器的制造年代以及其他耗费(如电能,机油等),原创 2021-07-03 22:59:48 · 4947 阅读 · 4 评论 -
管住嘴迈开腿
1.问题重述将模型假设修改为“每周体重减少不要超过 1kg”,检查所制定的减肥计划是否满足,如不满足,重新制定计划。2.模型假设体重增加正比于吸收的热量,平均每8000 ????????????????(????????????????为非国际单位制单位,1 ???????????????? = 4.2????????)增加体重1 ????????。身体正常代谢引起的体重减少正比于体重,每周每千克体重消耗热量一般在200 ????????????????至320 ????????????????之原创 2021-07-03 22:46:06 · 7520 阅读 · 2 评论 -
药物中毒急救 matlab
药物中毒急救1. 问题重述2. 模型假设1.胃肠道中药物向血液系统的转移率与药量成正比,且(???? > 0),总剂量????的药物在???? = 0瞬间进入胃肠道。2.血液系统中药物排除率与药量????(????)成正比,比例系数????(???? > 0),假设???? = 0时血液中无药物。3.氨茶碱被吸收的半衰期为5ℎ,排除的半衰期为6ℎ。4.孩子的血液总量为2000????????。3.符号说明4.模型建立与解决根据假设对胃肠道中药量????(????)和血液系统原创 2021-07-03 22:34:55 · 9706 阅读 · 13 评论 -
运输问题 Lingo
1. 问题重述某公司有 6 个建筑工地要开工,每个工地的位置(用平面坐标系 a,b 表示,距离单位:千米 )及水泥日用量 d(吨)由下表给出.目前有两个临时料场位于A(5,1),B(2,7),日储量各有 20 吨.假设从料场到工地之间均有直线道路相连.(1)试制定每天的供应计划,即从 A,B 两料场分别向各工地运送多少吨水泥,使总的吨千米数最小.2.模型假设1.不考虑实际运输成本问题,只考虑吨千米最低。2.为了便于计算,不考虑水泥的黏附损耗问题。3.不考虑运输过程其他因素对水泥质量的改变,仅考原创 2021-07-03 22:24:16 · 23673 阅读 · 19 评论 -
倾倒的啤酒杯
倾倒的啤酒杯1. 问题重述本文就如何调整啤酒杯中液面的高度,来使得酒与酒杯整体的质心达到最低点,讨论同时杯子的几何形状等其他因素,是否对该方面有一定影响。2.模型假设1.忽略酒杯底部厚度,不考虑其在重心计算内,且就要求而言对结果无影响。2.为了便于计算,不考虑杯底与杯沿接触时的形变。3.不考虑其他因素对溶液密度的影响,仅考虑溶液均匀分布的情况。3.符号说明4.问题分析为了使得杯与其内部的酒最不容易倾倒,首先需要考虑杯子保持静态平衡的条件,即质心的竖直投影应在物体与其支撑面的接触范围内。原创 2021-07-03 22:03:48 · 2905 阅读 · 2 评论 -
拯救小白鼠
拯救小白鼠1.问题重述假设有999个桶,有且仅有一个桶里面装了毒药,其他的都装了水。这些桶从外面看上去完全相同,如果一只小白鼠喝了毒药,它将在15分钟内死去,在一个小时内至少需要多少只小白鼠才能判断出哪一个桶里装的是毒药?2.模型假设1.假设所有的小白鼠都是相同的,且只考虑极限死亡时间15分钟。(待定)2.为了便于计算,不考虑老鼠喝水的时间。3.不考虑其他因素对装了毒药的桶的位置的影响,仅考虑随机作用。3.符号说明4.问题分析极端条件下,如果有足够精确的设备与条件(15min),一只小白原创 2021-07-03 21:50:21 · 254 阅读 · 0 评论 -
数学建模 员工指派问题 Lingo求解优化问题
第七次作业——员工指派问题1.问题重述2.模型假设1.不考虑实际通讯成本问题,只考虑给定的计算表达式最低。2.为了便于计算,不考虑通讯时间的变化问题。3.不考虑通讯单价随通话时间的增加而增加,仅考虑给定价格不变。3.符号说明4.问题分析本题可以处理为一个01规划问题,这个优化问题的目标是通话费用最小,即最小化通讯成本,同时又满足各地区人员的通话需求,要做的决策是工作人员安排方案,即十个不同的工作人员在十个工作地区的安排情况。决策受到三个条件的约束:通话时间,通话单价,同时决策变量默认为非负原创 2021-04-22 10:11:32 · 12198 阅读 · 4 评论 -
数学建模之倾倒的啤酒杯
数学建模之倾倒的啤酒杯代码clear;clca=10;k=0.3w1=200;w2=230;w3=50;r=15;s3=1;e=0.3;Pi=3.1415926;for i=4:7k=i/10;f=@(x)(w1*3*(x^4-k^4)*x^3/(4*(1-k^3))+2*r*e*Pi*a*(1-k^3)/3+w3*s3)/(w1*(x^3-k^3)/(1-k^3)+w2+w3)fplot(f,[0,1]);hold on;endylabel("质心高度");原创 2021-04-11 21:53:23 · 3632 阅读 · 0 评论 -
Leetcode 458. 可怜的小猪or拯救小白鼠
Leetcode 458. 可怜的小猪or拯救小白鼠代码clear;clcnum=999;mouse_num=144;T=5000;Time=zeros(1,T);i=1;while i<T box=zeros(1,num); loc=randi(num); box(loc)=-1; flag=0; sum=0; while flag==0 mouse=randi(num,1,mouse_num); for j=1:mouse_num原创 2021-04-11 21:47:53 · 121 阅读 · 0 评论 -
数学建模之出租车数目估计
数学建模之出租车数目估计代码m=100;%样本数目n=20;%样本容量x=zeros(5,m);y=zeros(5,20);for w=1:20 n=w*20;for k=1:m pe=randperm(1000); temp=pe(1:n); x(1,k)=2*mean(temp)-1; x(2,k)=2*median(temp)-1; x(3,k)=max(temp)+min(temp)原创 2021-04-11 21:41:17 · 10183 阅读 · 9 评论 -
雇员与雇主的关系之满意度曲线
1.问题重述2.模型假设假设雇员所得的总工资可以近似于其满意度。假设雇员对相同单位的工时与工资变化的喜恶程度相同。3.符号说明4.问题分析在实际情况中,雇主总是优先考虑自己的利益,即使利益最大化。此种情况下,雇主会尽量以最低的工资获得最大的收益,然而在市场的调节作用和人的主观判断下,不会妥协于同样的工资的更长时间的工作而无任何的加班费用,劳动人民会以各种形式联合起来打倒资本家。因此在雇员与雇主的相互交涉中,工资与工时的关系总会趋于稳定。对雇主来说,他只会着眼于总收益,在成熟的市场下,考原创 2021-04-11 21:36:18 · 1738 阅读 · 0 评论 -
matlab BP神经网络 用matlab库实现
BP神经网络 用matlab库实现先把代码存在这里,以后用了方便原理BP神经网络是一种多层的前馈神经网络,其主要的特点是:信号是前向传播的,而误差是反向传播的。具体来说,对于如下的只含一个隐层的神经网络模型:BP神经网络的过程主要分为两个阶段,第一阶段是信号的前向传播,从输入层经过隐含层,最后到达输出层;第二阶段是误差的反向传播,从输出层到隐含层,最后到输入层,依次调节隐含层到输出层的权重和偏置,输入层到隐含层的权重和偏置。计算流程神经网络的基本组成单元是神经元。神经元的通用模型如图所示,其原创 2021-04-11 20:57:31 · 1310 阅读 · 1 评论 -
模拟退火 matlab
模拟退火,求解函数的最小值。函数function result=func(x)result=3*cos(x(1)*x(2))+x(1)+x(2)^2;end算法部分clear; clc; close all;Xmax=5; Xmin=-5;Ymax=5; Ymin=-5;L=200; %马尔科夫链长度K=0.999; %衰减参数S=0.02; %步长因子T=100; %初始温度YZ=1e-8; %容差P=0;Prex=rand*(Xmax-Xmin)+Xmin;原创 2021-04-01 15:06:29 · 218 阅读 · 0 评论 -
蚁群算法求解TSP matlab
蚁群算法求解TSPclear; clc; close all;Ant_Num=50;Alpha=1; %信息素权重Beta=5; %启发式因子权重Rho=0.1; %信息素消失系数G=200; %迭代次数Q=100;%信息素增加强度系数City_Num=20; %城市数量border=100; %城市边界City=[rand(City_Num,1)*border,rand(City_Num,1)*border]; %随机城市坐标% City=[9.442061706065957e+01原创 2021-04-01 15:03:20 · 258 阅读 · 0 评论 -
粒子群算法matlab以求解函数最优解为例
粒子群算法matlab以求解函数最优解为例clear; clc; close all;N=100; %粒子个数D=2; %粒子维数MaxIter=500; %最大迭代次数C1max=1.8; %权重参数,自适应C2max=1.8; C1=1.2; C2=1.2;w=0.79;Wmax=0.8; %对自己速度的记忆Wmin=0.4;Xmax=4; Xmin=-4; %自变量的范围Vmax=1; Vmin=-1; %速度的范围x=rand(N,D)*(Xmax-Xmin原创 2021-04-01 15:01:43 · 1162 阅读 · 0 评论 -
遗传算法 matlab ---以tsp为例
美赛的时候准备的代码,没用上,不过现代优化算法挺有意思的,之后再发些其他的。close all clear allclc%读入初始数据N=20;border=100;%C=[rand(N,1)*border,rand(N,1)*border];C=[9.442061706065957e+01 5.554527461175734e+01 7.033585927645424e+01 8.210485000479132e+01 8.51388516080158原创 2021-04-01 14:59:14 · 225 阅读 · 0 评论 -
最短路之迪杰斯特拉算法
这个代码忘了是哪来的,最短路,迪杰斯特拉算法clc,clear alla=zeros(6); %输入最大的横坐标或纵坐标a(1,2)=50;a(1,4)=40;a(1,5)=25;a(1,6)=10;%把已知的权值输入进去 a(2,3)=15;a(2,4)=20;a(2,6)=25;a(3,4)=10;a(3,5)=20;a(4,5)=10;a(4,6)=25;a(5,6)=55;a=a+a'原创 2021-04-01 10:05:20 · 89 阅读 · 0 评论 -
PCA 主成分分析 Matlab
记得有些东西要改%% 数据导入处理clcclear allA = xlsread('数据绝对地址.xlsx','B2:I16');%% 数据标准化处理a = size(A,1);b = size(A,2);for i = 1:b SA(:,i) = (A(:,i) - mean(A(:,i)))/std(A(:,i));end%% 计算相关系数矩阵的特征值和特征向量CM = corrcoef(SA); %计算相关系数矩阵[V,原创 2021-03-31 09:53:15 · 285 阅读 · 0 评论 -
阻滞增长模型Logistic model Matlab
%Part1:Logistic Model % dP/dt=rP*dt,P(t0)=P0 % P:t时刻人口数 r:人口自然增长率%fplot(@(t)1000000*ones(1,2));K=1000;%设定人口上限r=0.18;syms P tP=dsolve('DP/((1-P/1000)*P)=0.18','P(0)=200',t);%K<P0P=matlabFunction(P);x=0.5:0.5:200;ft=[];for i=0.5:0.5:200.原创 2021-03-31 09:51:01 · 4134 阅读 · 1 评论 -
人工神经网络demo matlab
clear;clc;close all;x=[1,2,3,4,5,6,7];y=[1,4,9,16,25,36,49];net=fitnet(5);%创建网络net=train(net,x,y);%根据数据训练y=net(x);plot(x,y,'bx-');原创 2021-03-31 09:48:30 · 178 阅读 · 0 评论 -
聚类分析 matlab
%聚类分析代码:注意A为聚类矩阵,聚类时是以列为聚类单位。A=[1,2,3;2,3,4;3,4,5];BXA=zscore(A); % 标准化数据矩阵YA=pdist(A) % 用欧氏距离计算两两之间的距离DA=squareform(YA) % 欧氏距离矩阵ZA = linkage(YA) % 最短距离法TA = cluster(ZA,3) %等价于 { TA=clusterdata(A,3) }find(TA==3) % 第3类集合中的元素[HA,TA]=dend原创 2021-03-31 09:46:16 · 1780 阅读 · 0 评论 -
灰色预测法 matlab
最近会把自己总结的数学建模代码一并上传,欢迎交流clc;clear;%建立符号变量a(发展系数)和b(灰作用量)syms a b;c = [a b]';%原始数列 AA = [174, 179, 183, 189, 207, 234, 220.5, 256, 270, 285]; %输入待预测的数据n = length(A);%对原始数列 A 做累加得到数列 BB = cumsum(A);%对数列 B 做紧邻均值生成for i = 2:n C(i) = (B(i) + B(i原创 2021-03-31 09:27:55 · 345 阅读 · 3 评论 -
最短路之 Floyd弗洛伊德算法
三重循环很简单,不多说。欢迎交流% floyd.m% 采用floyd算法计算图a中每对顶点最短路% d是矩离矩阵% r是路由矩阵function [d,r]=floyd(a)n=size(a,1);% 初始化距离矩阵d=a;% 初始化路由矩阵for i=1:n for j=1:n r(i,j)=j; end end % Floyd算法开始for k=1:n for i=1:n for j=1:n if原创 2021-03-31 09:15:10 · 236 阅读 · 0 评论