随笔分类 - 数模
摘要:单目标优化 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(
阅读全文
摘要:代码: function GaTSPChen CityNum = 30; % 城市数目,可以选 10, 30, 50, 75 [dislist, Clist] = tsp(CityNum); % dislist 为城市之间相互的距离,Clist 为各城市的坐标 inn = 30; % 初始种群大小
阅读全文
摘要:非刚性 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
阅读全文
摘要:文件资源: https://github.com/Academic-monster/map-of-China 对于这种图可以找地图模板,然后在上面画,或者在excel三维地图上画
阅读全文
摘要:    ![
阅读全文
摘要:#模型 #代码1 fun1.m function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 beta = 0.1; % 易感染者与已感染者接触且被传染的强度 sigma = 0.2; % 潜伏者转换为感染者的速率 gamma = 0.02; % 康复率 dx = zero
阅读全文
摘要:#SIRS模型 #代码1 fun1.m function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 global TOTAL_N % 定义总人数为全局变量 beta = 0.1; % 易感染者与已感染者接触且被传染的强度 gamma = 0.03; % 康复率 alpha
阅读全文
摘要:#SIR模型 #代码 fun1.m function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 beta = 0.1; % 易感染者与已感染者接触且被传染的强度 gamma = 0.02; % 康复率 dx = zeros(3,1); % x(1)表示S x(2)表示I
阅读全文
摘要:#模型 #代码 fun1.m % 情况一:甲可以独自生存,乙不能独自生存 function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 r1=0.5; r2=0.5; % 甲的增长率和乙的死亡率 N1=300; N2=500; % 甲乙的最大数量 % sigma1: 单位数
阅读全文
摘要:#SI模型 #代码 fun1.m function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 global TOTAL_N % 定义总人数为全局变量 beta = 0.1; % 易感染者与已感染者接触且被传染的强度 dx = zeros(2,1); % x(1)表示S x
阅读全文
摘要:#SIS模型 #代码 fun1.m function dx=fun1(t,x) % 大家可以修改里面的参数,来看结果的变化 global TOTAL_N % 定义总人数为全局变量 beta = 0.1; % 易感染者与已感染者接触且被传染的强度 alpha = 0.06; % 由感染状态I恢复为易感
阅读全文
摘要:#模型 #代码 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
阅读全文
摘要:#模型 #代码 fun.m function dx=fun(t,x) % 大家可以修改里面的参数,来看结果的变化 r1=0.5; r2=0.5; % 甲乙的增长率 % r1=0.8; r2=1; % 甲乙的增长率 N1=300; N2=500; % 甲乙的最大数量 % sigma1: 单位数量的乙种
阅读全文
摘要:#模型: #代码: %% 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替
阅读全文
摘要:#微分方程的框架 #导弹引例问题 #基本概念和怎样建立微分方程 #MATLAB求微分方程解析解 %% 例1 clear;clc dsolve('y-Dy=2*x','x') % 这里要指定自变量为x % 2*x + C1*exp(x) + 2 (这里的C1表示任意常数,有时候也会出现C2 C3等)
阅读全文
摘要:#小技巧:MATLAB鼠标放在变量上面能够显示数值 这样就能显示啦 #Matlab求一阶微分方程数值解 数值解就是求一个一个点的值,然后连接在一起 #例一 df1.m function dy = df1(x,y) % 微分方程:y-y'=2x(函数名称可以任意取) dy = y - 2*x; % 写
阅读全文
摘要:更多关于Matlab求方程的介绍可看这个博客:https://blog.csdn.net/weixin_30724853/article/details/99004382 code.m %% matlab求解方程和方程组 % 不同的MATLAB版本之间的语法存在不兼容的情况:https://www.
阅读全文