随笔分类 -  数模

摘要:单目标优化 function y=fun1(x); %x 为行向量 y=x(1)*2+x(1).^2*3+x(2)*3+x(2).^2+x(3); y=-y; function [f,g]=fun2(x); f=[x(1)+2*x(1)^2+x(2)+2*x(2)^2+x(3)-10 x(1)+x( 阅读全文
posted @ 2022-05-09 10:50 司砚章 阅读(80) 评论(0) 推荐(0) 编辑
摘要:代码: function GaTSPChen CityNum = 30; % 城市数目,可以选 10, 30, 50, 75 [dislist, Clist] = tsp(CityNum); % dislist 为城市之间相互的距离,Clist 为各城市的坐标 inn = 30; % 初始种群大小 阅读全文
posted @ 2022-05-09 10:48 司砚章 阅读(180) 评论(0) 推荐(0) 编辑
摘要:非刚性 ode45 function dy=vdp1(t,y); dy=[y(2);(1-y(1)^2)*y(2)-y(1)]; [T,Y]=ode45('vdp1',[0 20],[2;0]); plot(T,Y(:,1),'-',T,Y(:,2),'--') title('Solution of 阅读全文
posted @ 2022-05-09 10:48 司砚章 阅读(103) 评论(0) 推荐(0) 编辑
摘要:文件资源: https://github.com/Academic-monster/map-of-China 对于这种图可以找地图模板,然后在上面画,或者在excel三维地图上画 阅读全文
posted @ 2022-05-09 10:48 司砚章 阅读(72) 评论(0) 推荐(0) 编辑
摘要:![image](https://img2020.cnblogs.com/blog/2140721/202109/2140721-20210909211504107-354693971.png) ![image](https://img2020.cnblogs.com/blog/2140721/202109/2140721-20210909213441231-1127858178.png) ![i 阅读全文
posted @ 2022-05-09 10:46 司砚章 阅读(216) 评论(0) 推荐(0) 编辑
摘要:采样时间 上面的采集总时长计算:拿文件一为例->因为数据总量为185725行 采样频率为1HZ 也就是1s采集一次 所以185725/2600=51.59个小时 另外用Python计算两个事件相差的小时: from dateutil.parser import parse a = parse('20 阅读全文
posted @ 2022-05-09 10:46 司砚章 阅读(239) 评论(0) 推荐(0) 编辑
摘要:#题目 问题1: 首先分析密度和压力传播的变化特点,在高压环境下,将压力的非均匀变化问题简化为均匀变化问题,然后进行模型的建立与求解:第一步是用最小二乘法对附件3中弹性模量与压力的对应关系进行拟合,得出二次多项式关系;第二步是根据注1中给出的密度变化量和压力变化量的正比例关系计算出附件2中压力值对应 阅读全文
posted @ 2022-05-09 10:46 司砚章 阅读(204) 评论(0) 推荐(0) 编辑
摘要:![image](https://img2020.cnblogs.com/blog/2140721/202108/2140721-20210825164138653-1195078028.png) ![image](https://img2020.cnblogs.com/blog/2140721/202108/2140721-20210825164145744-1769734731.png) ![ 阅读全文
posted @ 2022-05-09 10:40 司砚章 阅读(418) 评论(0) 推荐(0) 编辑
摘要:#模型 #代码1 fun1.m function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 beta = 0.1; % 易感染者与已感染者接触且被传染的强度 sigma = 0.2; % 潜伏者转换为感染者的速率 gamma = 0.02; % 康复率 dx = zero 阅读全文
posted @ 2022-05-09 10:40 司砚章 阅读(1805) 评论(0) 推荐(0) 编辑
摘要:#SIRS模型 #代码1 fun1.m function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 global TOTAL_N % 定义总人数为全局变量 beta = 0.1; % 易感染者与已感染者接触且被传染的强度 gamma = 0.03; % 康复率 alpha 阅读全文
posted @ 2022-05-09 10:39 司砚章 阅读(2092) 评论(0) 推荐(0) 编辑
摘要:#SIR模型 #代码 fun1.m function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 beta = 0.1; % 易感染者与已感染者接触且被传染的强度 gamma = 0.02; % 康复率 dx = zeros(3,1); % x(1)表示S x(2)表示I 阅读全文
posted @ 2022-05-08 10:21 司砚章 阅读(1072) 评论(0) 推荐(0) 编辑
摘要:#模型 #代码 fun1.m % 情况一:甲可以独自生存,乙不能独自生存 function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 r1=0.5; r2=0.5; % 甲的增长率和乙的死亡率 N1=300; N2=500; % 甲乙的最大数量 % sigma1: 单位数 阅读全文
posted @ 2022-05-08 10:19 司砚章 阅读(2170) 评论(0) 推荐(0) 编辑
摘要:#SI模型 #代码 fun1.m function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 global TOTAL_N % 定义总人数为全局变量 beta = 0.1; % 易感染者与已感染者接触且被传染的强度 dx = zeros(2,1); % x(1)表示S x 阅读全文
posted @ 2022-05-08 10:19 司砚章 阅读(867) 评论(0) 推荐(0) 编辑
摘要:#SIS模型 #代码 fun1.m function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 global TOTAL_N % 定义总人数为全局变量 beta = 0.1; % 易感染者与已感染者接触且被传染的强度 alpha = 0.06; % 由感染状态I恢复为易感 阅读全文
posted @ 2022-05-08 10:19 司砚章 阅读(924) 评论(0) 推荐(0) 编辑
摘要:#模型 #代码 post_war.m % 战后 function dx=post_war(t,x) dx=zeros(2,1); dx(1)=x(1)*(0.9-0.1*x(2)); dx(2)=x(2)*(-0.6+0.02*x(1)); end pre_war.m % սǰ function d 阅读全文
posted @ 2022-05-08 10:17 司砚章 阅读(993) 评论(0) 推荐(0) 编辑
摘要:#模型 #代码 fun.m function dx=fun(t,x) % 大家可以修改里面的参数,来看结果的变化 r1=0.5; r2=0.5; % 甲乙的增长率 % r1=0.8; r2=1; % 甲乙的增长率 N1=300; N2=500; % 甲乙的最大数量 % sigma1: 单位数量的乙种 阅读全文
posted @ 2022-05-08 10:17 司砚章 阅读(518) 评论(0) 推荐(0) 编辑
摘要:#模型: #代码: %% Malthus模型(马尔萨斯模型) clear;clc x = dsolve('Dx=r*x','x(0)=x0','t') % x = dsolve('Dx=r*x','x(t0)=x0','t') % x = x0*exp(r*t) % 怎么把上面这个式子中的x0和r替 阅读全文
posted @ 2022-05-08 10:16 司砚章 阅读(2455) 评论(0) 推荐(0) 编辑
摘要:#微分方程的框架 #导弹引例问题 #基本概念和怎样建立微分方程 #MATLAB求微分方程解析解 %% 例1 clear;clc dsolve('y-Dy=2*x','x') % 这里要指定自变量为x % 2*x + C1*exp(x) + 2 (这里的C1表示任意常数,有时候也会出现C2 C3等) 阅读全文
posted @ 2022-05-08 10:08 司砚章 阅读(228) 评论(0) 推荐(0) 编辑
摘要:#小技巧:MATLAB鼠标放在变量上面能够显示数值 这样就能显示啦 #Matlab求一阶微分方程数值解 数值解就是求一个一个点的值,然后连接在一起 #例一 df1.m function dy = df1(x,y) % 微分方程:y-y'=2x(函数名称可以任意取) dy = y - 2*x; % 写 阅读全文
posted @ 2022-05-08 10:08 司砚章 阅读(963) 评论(0) 推荐(0) 编辑
摘要:更多关于Matlab求方程的介绍可看这个博客:https://blog.csdn.net/weixin_30724853/article/details/99004382 code.m %% matlab求解方程和方程组 % 不同的MATLAB版本之间的语法存在不兼容的情况:https://www. 阅读全文
posted @ 2022-05-08 00:19 司砚章 阅读(222) 评论(0) 推荐(0) 编辑

点击右上角即可分享
微信分享提示