MATLAB复现——基于演化博弈的南方林地森林装备网络平台模式研究

今天主要是介绍对一个中文期刊“农业与技术”的一篇演化博弈的论文《基于演化博弈的南方林地森林装备网络平台模式研究》的数值仿真部分进行MATLAB编程复现。

摘要:本文基于农户与经营平台双方关键主体，运用演化博弈理论的方法，构建了理性经济人假设下的农户与经营平台的合作演化博弈模型，分析了博弈双方合作以及不合作条件下采取惩罚机制的动态博弈关系，并利用MATLAB进行模拟仿真，对比了在不同影响因素下的博弈策略选择。研究结果表明:农户与经营平台选择合作策略与所获得的超额收益、经济损失密切相关，单方不合作带来的处罚罚金也会影响到博弈结果。因此需要将农户与经营平台的经济利益紧密结合,建立良好的处罚机制对双方合作进行监督。

1.内容

2.数值仿真复现结果

3.程序

需要完整程序，可以进行打赏后截图（50元及以上），进行联系，或者公众号“联系掌门”进行联系，或者在公众号内回复截图，几小时内会回复。论文复现编程不易，还请见谅！

主程序

clc;
clear all;
close all;
%参数取值
alpha = 7/12;
V = 60;
P3 = (1-alpha)*V;
P4 = alpha*V ;
C1 = 8;
C2 = 12;
L2 = 7;
L1 = 10;
F1 = 20;
F2 = 30;
C3 = 15;
m = 2;
h = 0.2;
tspan = [0:h*0.1:1];
x0 = [0.6;0.6]; %初值
P3_I = [15 35 45 55];
result = [];
for i = 1:length(P3_I)
    P3 = P3_I(i);
    [t,x] = ode45(@(t,Y) myfunY1(t,Y,P3,C1,L2,F1,C3,m,P4,C2,L1,F2),tspan,x0);
    result = [result;x'];
end
figure(1);
plot(t',result(1,:),'b-o');
hold on;
plot(t',result(3,:),'b-p');
plot(t',result(5,:),'b-x');
plot(t',result(7,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('X:P3 = 15','X:P3 = 35','X:P3 = 45','X:P3 = 55');
grid on;
ylim([ 0 1]);
P4_I = [15 25 45 55];
P3 = (1-alpha)*V;
result1 = [];
for i = 1:length(P4_I)
    P4 = P4_I(i);
    [t1,x1] = ode45(@(t,Y) myfunY1(t,Y,P3,C1,L2,F1,C3,m,P4,C2,L1,F2),tspan,x0);
    result1 = [result1;x1'];
end
figure(2);
plot(t1',result1(2,:),'b-o');
hold on;
plot(t1',result1(4,:),'b-p');
plot(t1',result1(6,:),'b-x');
plot(t1',result1(8,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('Y:P4 = 15','Y:P4 = 25','Y:P4 = 45','Y:P4 = 55');
grid on;
ylim([ 0 1]);


P3 = (1-alpha)*V;
P4 = alpha*V ;
F1_I = [10 15 25 30];
result2 = [];
for i = 1:length(F1_I)
   F1 = F1_I(i);
    [t2,x2] = ode45(@(t,Y) myfunY1(t,Y,P3,C1,L2,F1,C3,m,P4,C2,L1,F2),tspan,x0);
    result2 = [result2;x2'];
end
figure(3);
plot(t2',result2(1,:),'b-o');
hold on;
plot(t2',result2(3,:),'b-p');
plot(t2',result2(5,:),'b-x');
plot(t2',result2(7,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('X:F1 = 10','X:F1 = 15','X:F1 = 25','X:F1 = 30');
grid on;
ylim([ 0 1]);
P3 = (1-alpha)*V;
P4 = alpha*V ;
F1 = 20;
F2_I = [20 25 35 40];
result3 = [];
for i = 1:length(F2_I)
   F2 = F2_I(i);
    [t3,x3] = ode45(@(t,Y) myfunY1(t,Y,P3,C1,L2,F1,C3,m,P4,C2,L1,F2),tspan,x0);
    result3 = [result3;x3'];
end
figure(4);
plot(t3',result3(2,:),'b-o');
hold on;
plot(t3',result3(4,:),'b-p');
plot(t3',result3(6,:),'b-x');
plot(t3',result3(8,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('Y:F2 = 20','Y:F2 = 25','Y:F2 = 35','Y:F2 = 40');
grid on;
ylim([ 0 1]);
P3 = (1-alpha)*V;
P4 = alpha*V ;
F1 = 20;
F2 = 30;
L1_I = [5 8 15 20];
result4 = [];
for i = 1:length(L1_I)
   L1 = L1_I(i);
    [t4,x4] = ode45(@(t,Y) myfunY1(t,Y,P3,C1,L2,F1,C3,m,P4,C2,L1,F2),tspan,x0);
    result4 = [result4;x4'];
end
figure(5);
plot(t4',result4(2,:),'b-o');
hold on;
plot(t4',result4(4,:),'b-p');
plot(t4',result4(6,:),'b-x');
plot(t4',result4(8,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('Y:L1 = 5','Y:L1 = 8','Y:L1 = 15','Y:L1 = 20');
grid on;
ylim([ 0 1]);
P3 = (1-alpha)*V;
P4 = alpha*V ;
F1 = 20;
F2 = 30;
L1 = 10;
L2_I = [1 3 9 12];
result5 = [];
for i = 1:length(L2_I )
   L2 = L1_I(i);
    [t5,x5] = ode45(@(t,Y) myfunY1(t,Y,P3,C1,L2,F1,C3,m,P4,C2,L1,F2),tspan,x0);
    result5 = [result5;x5'];
end
figure(6);
plot(t5',result5(1,:),'b-o');
hold on;
plot(t5',result5(3,:),'b-p');
plot(t5',result5(5,:),'b-x');
plot(t5',result5(7,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('X:L2 = 1','X:L2 = 3','X:L2 = 9','X:L2 = 12');
grid on;
ylim([ 0 1]);
P3 = (1-alpha)*V;
P4 = alpha*V ;
F1 = 20;
F2 = 30;
L1 = 10;
L2 = 7;
m_I= [0 6 8 10];
result6 = [];
for i = 1:length(m_I)
   m = m_I(i);
    [t6,x6] = ode45(@(t,Y) myfunY1(t,Y,P3,C1,L2,F1,C3,m,P4,C2,L1,F2),tspan,x0);
    result6 = [result6;x6'];
end
figure(7);
plot(t6',result6(1,:),'b-o');
hold on;
plot(t6',result6(3,:),'b-p');
plot(t6',result6(5,:),'b-x');
plot(t6',result6(7,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('X:m = 0','X:m = 6','X:m = 8','X:m = 10');
grid on;
ylim([ 0 1]);
figure(8);
plot(t6',result6(2,:),'b-o');
hold on;
plot(t6',result6(4,:),'b-p');
plot(t6',result6(6,:),'b-x');
plot(t6',result6(8,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('Y:m = 0','Y:m = 6','Y:m = 8','Y:m = 10');
grid on;
ylim([ 0 1]);
F1 = 20;
F2 = 30;
L1 = 10;
L2 = 7;
alpha_I = [0.1 0.3 0.5 0.7];
V = 60;
m = 2;


result7 = [];
for i = 1:length(alpha_I)
    alpha = alpha_I(i);
    P3 = (1-alpha)*V;
    P4 = alpha*V ;
    [t7,x7] = ode45(@(t,Y) myfunY1(t,Y,P3,C1,L2,F1,C3,m,P4,C2,L1,F2),tspan,x0);
    result7 = [result7;x7'];
end
figure(9);
plot(t7',result7(1,:),'b-o');
hold on;
plot(t7',result7(3,:),'b-p');
plot(t7',result7(5,:),'b-x');
plot(t7',result7(7,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('X:\alpha = 0.1','X:\alpha = 0.3','X:\alpha = 0.5','X:\alpha = 0.7');
grid on;
ylim([ 0 1]);
figure(10);
plot(t7',result7(2,:),'b-o');
hold on;
plot(t7',result7(4,:),'b-p');
plot(t7',result7(6,:),'b-x');
plot(t7',result7(8,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('Y:\alpha = 0.1','Y:\alpha = 0.3','Y:\alpha = 0.5','Y:\alpha = 0.7');
grid on;
ylim([ 0 1]);


F1 = 20;
F2 = 30;
L1 = 10;
L2 = 7;
alpha = 7/12;
V_I = [40 50 70 80];
m = 2;


result8 = [];
for i = 1:length(V_I )
    V = V_I (i);
    P3 = (1-alpha)*V;
    P4 = alpha*V ;
    [t8,x8] = ode45(@(t,Y) myfunY1(t,Y,P3,C1,L2,F1,C3,m,P4,C2,L1,F2),tspan,x0);
    result8 = [result8;x8'];
end
figure(11);
plot(t8',result8(1,:),'b-o');
hold on;
plot(t8',result8(3,:),'b-p');
plot(t8',result8(5,:),'b-x');
plot(t8',result8(7,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('X:V = 40','X:V = 50','X:V = 70','X:V = 80');
grid on;
ylim([ 0 1]);


figure(12);
plot(t8',result8(2,:),'b-o');
hold on;
plot(t8',result8(4,:),'b-p');
plot(t8',result8(6,:),'b-x');
plot(t8',result8(8,:),'b-s');
xlabel('Time');
ylabel('Solution');
legend('Y:V = 40','Y:V = 50','Y:V = 70','Y:V = 80');
grid on;
ylim([ 0 1]);

本文内容来源于网络，仅供参考学习，如内容、图片有任何版权问题，请联系处理，24小时内删除。

作 者 | 郭志龙
编 辑 | 郭志龙
校 对 | 郭志龙