![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
matlab
CAMlive
炼丹术士
展开
-
生产任务分配问题 matlab+lingo
A题 生产任务分配问题有一个生产管理人员决定下个时期生产10批产品,他现在需要决定怎样将这些生产任务分配到其管理的车间中的不同机器上。此车间内的5台可用机器都可以加工这些产品,但由于它们是在不同时期生产的,因此其加工速度各不相同。此外,由于存在维护和调整期,在规划的时期内每台机器只能工作有限时间(小时数)。表1列出了所有机器加工各批产品所需的时间以及此机器可用的工时数。各批次产品的生产成本取决于处理此批次的机器。每台机器每小时运行成本取决于其采用的技术,此机器的制造年代以及其他耗费(如电能,机油等),原创 2021-07-03 22:59:48 · 5356 阅读 · 4 评论 -
管住嘴迈开腿
1.问题重述将模型假设修改为“每周体重减少不要超过 1kg”,检查所制定的减肥计划是否满足,如不满足,重新制定计划。2.模型假设体重增加正比于吸收的热量,平均每8000 ????????????????(????????????????为非国际单位制单位,1 ???????????????? = 4.2????????)增加体重1 ????????。身体正常代谢引起的体重减少正比于体重,每周每千克体重消耗热量一般在200 ????????????????至320 ????????????????之原创 2021-07-03 22:46:06 · 8866 阅读 · 3 评论 -
药物中毒急救 matlab
药物中毒急救1. 问题重述2. 模型假设1.胃肠道中药物向血液系统的转移率与药量成正比,且(???? > 0),总剂量????的药物在???? = 0瞬间进入胃肠道。2.血液系统中药物排除率与药量????(????)成正比,比例系数????(???? > 0),假设???? = 0时血液中无药物。3.氨茶碱被吸收的半衰期为5ℎ,排除的半衰期为6ℎ。4.孩子的血液总量为2000????????。3.符号说明4.模型建立与解决根据假设对胃肠道中药量????(????)和血液系统原创 2021-07-03 22:34:55 · 10633 阅读 · 13 评论 -
倾倒的啤酒杯
倾倒的啤酒杯1. 问题重述本文就如何调整啤酒杯中液面的高度,来使得酒与酒杯整体的质心达到最低点,讨论同时杯子的几何形状等其他因素,是否对该方面有一定影响。2.模型假设1.忽略酒杯底部厚度,不考虑其在重心计算内,且就要求而言对结果无影响。2.为了便于计算,不考虑杯底与杯沿接触时的形变。3.不考虑其他因素对溶液密度的影响,仅考虑溶液均匀分布的情况。3.符号说明4.问题分析为了使得杯与其内部的酒最不容易倾倒,首先需要考虑杯子保持静态平衡的条件,即质心的竖直投影应在物体与其支撑面的接触范围内。原创 2021-07-03 22:03:48 · 3492 阅读 · 2 评论 -
拯救小白鼠
拯救小白鼠1.问题重述假设有999个桶,有且仅有一个桶里面装了毒药,其他的都装了水。这些桶从外面看上去完全相同,如果一只小白鼠喝了毒药,它将在15分钟内死去,在一个小时内至少需要多少只小白鼠才能判断出哪一个桶里装的是毒药?2.模型假设1.假设所有的小白鼠都是相同的,且只考虑极限死亡时间15分钟。(待定)2.为了便于计算,不考虑老鼠喝水的时间。3.不考虑其他因素对装了毒药的桶的位置的影响,仅考虑随机作用。3.符号说明4.问题分析极端条件下,如果有足够精确的设备与条件(15min),一只小白原创 2021-07-03 21:50:21 · 265 阅读 · 0 评论 -
数学建模之倾倒的啤酒杯
数学建模之倾倒的啤酒杯代码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 · 3854 阅读 · 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 · 123 阅读 · 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 · 10350 阅读 · 9 评论 -
matlab BP神经网络 用matlab库实现
BP神经网络 用matlab库实现先把代码存在这里,以后用了方便原理BP神经网络是一种多层的前馈神经网络,其主要的特点是:信号是前向传播的,而误差是反向传播的。具体来说,对于如下的只含一个隐层的神经网络模型:BP神经网络的过程主要分为两个阶段,第一阶段是信号的前向传播,从输入层经过隐含层,最后到达输出层;第二阶段是误差的反向传播,从输出层到隐含层,最后到输入层,依次调节隐含层到输出层的权重和偏置,输入层到隐含层的权重和偏置。计算流程神经网络的基本组成单元是神经元。神经元的通用模型如图所示,其原创 2021-04-11 20:57:31 · 1331 阅读 · 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 · 230 阅读 · 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 · 261 阅读 · 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 · 1219 阅读 · 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 · 230 阅读 · 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 · 92 阅读 · 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 · 292 阅读 · 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 · 4241 阅读 · 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 · 182 阅读 · 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 · 1838 阅读 · 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 · 367 阅读 · 3 评论 -
最小二乘法拟合 matlab
一般先验证得到最优次项,再来拟合。%最小二乘拟合法%已知m个点的值,求表达式%容易存在过拟合的情况,此方法可以改进%对于1~m-1每个次方程,提出训练集与验证集%做出其验证集损失曲线,%其曲线第一个导数为零的点对应的训练集大小即为标准解clc;clear;x=[1,2,3];y=[2,8,18];%y=2*x^3m=size(x);m=m(1,2);num=2;%拟合多少次方的多项式<=m-1A=[];for i=0:num A=[A,(x').^i];endB=原创 2021-03-31 09:25:23 · 301 阅读 · 0 评论 -
复合辛普森公式求解定积分 matlab
复合辛普森积分已知函数表达式与积分区间精度esp正相关与1/num%复合辛普森积分%已知函数表达式与积分区间clc;clear;a=0;b=1;%积分范围num=1000;%积分准确度h=(b-a)/(2*num);f=@(x)exp(-x);%积分表达式I=0;%积分结果I=f(a)-f(b);for i=1:num I=I+(2*f(a+2*i*h)+4*f(a+(2*i-1)*h));endI=I*h/3...原创 2021-03-31 09:22:05 · 6332 阅读 · 0 评论 -
改进欧拉法求解常微分方程 matlab
最近开始将平时的作业整理上传,欢迎交流。%改进欧拉法求解常微分方程数值解%初始条件%dy/dx=f(x)%y(x0)=y0clc;clear;f=@(x,y)1/x+1/y;%方程x=[];y=[];%初值x(1)=1;y(1)=1;num=100;%迭代次数a=1;b=50;%自变量区间h=(b-a)/(num-1);for i=2:num x(i)=x(i-1)+h;endfor i=2:num y(i)=y(i-1)+h*f(x(i-1),y(i-1));原创 2021-03-31 09:18:55 · 3018 阅读 · 0 评论 -
最短路之 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 · 244 阅读 · 0 评论