4. 四方演化博弈代碼

論文標題：重大疫情期醫(yī)藥研究報道質(zhì)量監(jiān)管四方演化博弈分析

關(guān)注b站：譚小飛同學(xué)

%sifang.m

% 特殊的希臘字符換成其讀音字母表示挡闰，文中的r,m,g,p依次表示為y(1),y(2),y(3),y(4)
% 論文標題為《重大疫情期醫(yī)藥研究報道質(zhì)量監(jiān)管四方演化博弈分析》
% 由于公式較長，本人雖然幾番認真仔細核對，難免也可能出錯，建議讀者自行錄入公式烁涌，看看是否有出入
function dydt=sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp)
dydt=zeros(4,1);
dydt(1)=a*y(1)*(1-y(1))*((1-(1-y(2))*(1-y(3))*(1-y(4)))*phi+(1-y(2))*(1-y(3))*(y(4)+(1-y(4))*b)*Rr+(1-y(2))*y(3)*Fr-Crh+Crl);
dydt(2)=a*y(2)*(1-y(2))*((1-y(1))*(1-y(3))*(y(4)+(1-y(4))*b)*Rm+(1-y(1))*y(3)*Fm-(1-y(1))*(1-y(3))*(1-y(4))*psi+Cml-Cmh);
dydt(3)=a*y(3)*(1-y(3))*((y(1)+(1-y(1))*(1-y(2)))*(Cgl-Cgh)+(1-y(1))*(1-y(2))*(Fr+Fm+(1-y(4))*(Ng-b*Tg)));
dydt(4)=a*y(4)*(1-y(4))*((1-y(1))*(1-y(2))*(1-y(3))*(Np-b*Tp-Cp)-y(1)*Cp);
end

sifang2.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%政府大圖
figure(1)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%子圖1
clc;clear;
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=12,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
%subplot(3,1,1)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,200],[0.4,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
% hold on放在函數(shù)ode45前面的話，生成的圖像不是封閉的方框酒觅，而是坐標系的第一象限
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:50:200],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 200 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$C_{gh}=12$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖2
clc;clear;
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
figure(2)
%subplot(3,1,2)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,200],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:50:200],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 200 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$C_{gh}=6$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖3
clc;clear;
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=2,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
figure(3)
%subplot(3,1,3)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,200],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:50:200],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 200 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$C_{gh}=2$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%群眾大圖
figure(4)
%%%%%%%%%%%%%%%子圖1
clc;clear;
a=0.7,phi=22,b=0.3,Rr=12,Fr=10,Crh=10,Crl=2,Rm=6,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
%subplot(3,1,1)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.4,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$R_{r}=12,R_{m}=6$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖2
clc;clear;
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
figure(5)
%subplot(3,1,2)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$R_{r}=25,R_{m}=15$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖3
clc;clear;
a=0.7,phi=22,b=0.3,Rr=50,Fr=10,Crh=10,Crl=2,Rm=35,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
figure(6)
%subplot(3,1,3)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
%下面一行是加方框的撮执，不過也可以在成圖界面去手動制圖
%rectangle('position',[0 0 2 1],'edgecolor','g','linewidth',1)
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$R_{r}=50,R_{m}=35$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%%%% 群眾第三張圖的小圖
%下面一行決定小圖的位置及其大小
axes('position',[0.35 0.2 0.2 0.2]);
clc;clear;
a=0.7,phi=22,b=0.3,Rr=50,Fr=10,Crh=10,Crl=2,Rm=35,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:1:3],'YTick',[0:0.2:1])
hold on
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 3 -0.05 1.05])
hold on 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%民間組織大圖
figure(7)
%%%%%%%%%%%%%%%子圖1
clc;clear;
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
%subplot(3,1,1)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$\beta=0.3$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖2
clc;clear;
a=0.7,phi=22,b=0.31,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
figure(8)
%subplot(3,1,2)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$\beta=0.31$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖3
clc;clear;
a=0.7,phi=22,b=0.5,Rr=50,Fr=10,Crh=10,Crl=2,Rm=35,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
figure(9)
%subplot(3,1,3)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
%又到了加方框的時候了
%rectangle('position',[0 0 2 1],'edgecolor','g','linewidth',1)
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$\beta=0.5$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%%%% 民間組織第三張圖的小圖
axes('position',[0.35 0.2 0.2 0.2]);
clc;clear;
a=0.7,phi=22,b=0.5,Rr=50,Fr=10,Crh=10,Crl=2,Rm=35,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:1:3],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 3 -0.05 1.05])
hold on 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3d figure
%%%%%%%%%%%%%%%%%%%%%%%%%g=0的子圖
clc;clear; 
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;  
figure(10)
%subplot(2,1,1)
for i=0.1:0.2:1
    for j=0.1:0.2:1
        for k=0
            for l=0.1:0.2:1
            [t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0 50],[i j k l]);
            grid on
            %搞清楚你畫的是哪三維的關(guān)系
            plot3(y(:,1),y(:,2),y(:,4),'linewidth',1);
            set(gca,'XTick',[0:0.2:1],'YTick',[0:0.2:1],'ZTick',[0:0.2:1])
            set(gca,'XTickLabel',num2str(get(gca,'XTick')','%.1f'));
            set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
            set(gca,'ZTickLabel',num2str(get(gca,'ZTick')','%.1f'));
            hold on
            axis([0 1 0 1 0 1])
            end
        end
    end
end
xlabel('$r$','interpreter','latex');
ylabel('$m$','interpreter','latex');
zlabel('$p$','interpreter','latex','Rotation',360);
title('g=0','interpreter','latex');
%%%%%%%%%%%%%%%%%%%%%%%%%g=0.8的子圖
clc;clear; 
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;  
figure(11)
%subplot(2,1,2)
for i=0.1:0.2:1
    for j=0.1:0.2:1
        for k=0.8
            for l=0.1:0.2:1
            [t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0 50],[i j k l]);
            grid on
            %搞清楚你畫的是哪三維的關(guān)系
            plot3(y(:,1),y(:,2),y(:,4),'linewidth',1);
            set(gca,'XTick',[0:0.2:1],'YTick',[0:0.2:1],'ZTick',[0:0.2:1])
            set(gca,'XTickLabel',num2str(get(gca,'XTick')','%.1f'));
            set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
            set(gca,'ZTickLabel',num2str(get(gca,'ZTick')','%.1f'));
            hold on
            axis([0 1 0 1 0 1])
            end
        end
    end
end
xlabel('$r$','interpreter','latex');
ylabel('$m$','interpreter','latex');
zlabel('$p$','interpreter','latex','Rotation',360);
title('{\it\fontname{Bodoni MT}g}=0.8');

image.png

請注意，上面圖片與原論文中的樣圖還是有很大的差別舷丹，下面我們使用subplot()函數(shù)及MATLAB的圖像編輯器來逼近原樣圖（具體操作看上方鏈接視頻）抒钱。

稍微改變下代碼如下：

% 特殊的希臘字符換成其讀音字母表示，文中的r,m,g,p依次表示為y(1),y(2),y(3),y(4)
% 論文標題為《重大疫情期醫(yī)藥研究報道質(zhì)量監(jiān)管四方演化博弈分析》
% 由于公式較長颜凯，本人雖然幾番認真仔細核對谋币，難免也可能出錯，建議讀者自行錄入公式症概，看看是否有出入
function dydt=sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp)
dydt=zeros(4,1);
dydt(1)=a*y(1)*(1-y(1))*((1-(1-y(2))*(1-y(3))*(1-y(4)))*phi+(1-y(2))*(1-y(3))*(y(4)+(1-y(4))*b)*Rr+(1-y(2))*y(3)*Fr-Crh+Crl);
dydt(2)=a*y(2)*(1-y(2))*((1-y(1))*(1-y(3))*(y(4)+(1-y(4))*b)*Rm+(1-y(1))*y(3)*Fm-(1-y(1))*(1-y(3))*(1-y(4))*psi+Cml-Cmh);
dydt(3)=a*y(3)*(1-y(3))*((y(1)+(1-y(1))*(1-y(2)))*(Cgl-Cgh)+(1-y(1))*(1-y(2))*(Fr+Fm+(1-y(4))*(Ng-b*Tg)));
dydt(4)=a*y(4)*(1-y(4))*((1-y(1))*(1-y(2))*(1-y(3))*(Np-b*Tp-Cp)-y(1)*Cp);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%政府大圖
figure(1)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%子圖1
clc;clear;
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=12,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
subplot(3,1,1)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,200],[0.4,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
% hold on放在函數(shù)ode45前面的話蕾额，生成的圖像不是封閉的方框，而是坐標系的第一象限
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:50:200],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 200 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$C_{gh}=12$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖2
clc;clear;
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
%figure(2)
subplot(3,1,2)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,200],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:50:200],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 200 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$C_{gh}=6$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖3
clc;clear;
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=2,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
%figure(3)
subplot(3,1,3)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,200],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:50:200],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 200 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$C_{gh}=2$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%群眾大圖
figure(4)
%%%%%%%%%%%%%%%子圖1
clc;clear;
a=0.7,phi=22,b=0.3,Rr=12,Fr=10,Crh=10,Crl=2,Rm=6,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
subplot(3,1,1)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.4,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$R_{r}=12,R_{m}=6$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖2
clc;clear;
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
%figure(5)
subplot(3,1,2)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$R_{r}=25,R_{m}=15$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖3
clc;clear;
a=0.7,phi=22,b=0.3,Rr=50,Fr=10,Crh=10,Crl=2,Rm=35,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
%figure(6)
subplot(3,1,3)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
%下面一行是加方框的彼城，不過也可以在成圖界面去手動制圖
%rectangle('position',[0 0 2 1],'edgecolor','g','linewidth',1)
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$R_{r}=50,R_{m}=35$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%%%% 群眾第三張圖的小圖
%下面一行決定小圖的位置及其大小
axes('position',[0.35 0.2 0.1 0.1]);
clc;clear;
a=0.7,phi=22,b=0.3,Rr=50,Fr=10,Crh=10,Crl=2,Rm=35,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:1:3],'YTick',[0:0.2:1])
hold on
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 3 -0.05 1.05])
hold on 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%民間組織大圖
figure(7)
%%%%%%%%%%%%%%%子圖1
clc;clear;
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
subplot(3,1,1)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$\beta=0.3$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖2
clc;clear;
a=0.7,phi=22,b=0.31,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
%figure(8)
subplot(3,1,2)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$\beta=0.31$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%子圖3
clc;clear;
a=0.7,phi=22,b=0.5,Rr=50,Fr=10,Crh=10,Crl=2,Rm=35,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
%figure(9)
subplot(3,1,3)
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
%又到了加方框的時候了
%rectangle('position',[0 0 2 1],'edgecolor','g','linewidth',1)
hold on
set(gca,'XTick',[0:20:100],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 100 -0.05 1.05])
xlabel('$t$','interpreter','latex');
ylabel('概率');
zhuti=title('$\beta=0.5$');
set(zhuti,'interpreter','latex')
legend('醫(yī)藥研究機構(gòu)({\it\fontname{Bodoni MT}r})','媒體平臺({\it\fontname{Bodoni MT}m})','政府監(jiān)管部門({\it\fontname{Bodoni MT}g})','群眾({\it\fontname{Bodoni MT}p})');
%%%%%%%%%%%%%%%%%% 民間組織第三張圖的小圖
axes('position',[0.35 0.2 0.1 0.1]);
clc;clear;
a=0.7,phi=22,b=0.5,Rr=50,Fr=10,Crh=10,Crl=2,Rm=35,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;
set(0,'defaultfigurecolor','w')
[t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0,100],[0.5,0.3,0.2,0.3]);
points=1:1:length(t);
plot(t,y(:,1),'r^-','linewidth',1,'markersize',3,'markerfacecolor','r','markerindices',points);
hold on
plot(t,y(:,2),'b-','linewidth',1);
hold on
plot(t,y(:,3),'y-.','linewidth',1);
hold on
plot(t,y(:,4),'g--','linewidth',1);
hold on
set(gca,'XTick',[0:1:3],'YTick',[0.0:0.2:1.0])
set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
axis([0 3 -0.05 1.05])
hold on 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3d figure
%%%%%%%%%%%%%%%%%%%%%%%%%g=0的子圖
clc;clear; 
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;  
figure(10)
subplot(2,1,1)
for i=0.1:0.2:1
    for j=0.1:0.2:1
        for k=0
            for l=0.1:0.2:1
            [t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0 50],[i j k l]);
            grid on
            %搞清楚你畫的是哪三維的關(guān)系
            plot3(y(:,1),y(:,2),y(:,4),'linewidth',1);
            set(gca,'XTick',[0:0.2:1],'YTick',[0:0.2:1],'ZTick',[0:0.2:1])
            set(gca,'XTickLabel',num2str(get(gca,'XTick')','%.1f'));
            set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
            set(gca,'ZTickLabel',num2str(get(gca,'ZTick')','%.1f'));
            hold on
            axis([0 1 0 1 0 1])
            end
        end
    end
end
xlabel('$r$','interpreter','latex');
ylabel('$m$','interpreter','latex');
zlabel('$p$','interpreter','latex','Rotation',360);
title('g=0','interpreter','latex');
%%%%%%%%%%%%%%%%%%%%%%%%%g=0.8的子圖
clc;clear; 
a=0.7,phi=22,b=0.3,Rr=25,Fr=10,Crh=10,Crl=2,Rm=15,Fm=6,psi=12,Cml=1,Cmh=5,Cgl=1,Cgh=6,Ng=20,Tg=10,Np=12,Tp=5,Cp=10;  
%figure(11)
subplot(2,1,2)
for i=0.1:0.2:1
    for j=0.1:0.2:1
        for k=0.8
            for l=0.1:0.2:1
            [t,y]=ode45(@(t,y) sifang(t,y,a,phi,b,Rr,Fr,Crh,Crl,Rm,Fm,psi,Cml,Cmh,Cgl,Cgh,Ng,Tg,Np,Tp,Cp),[0 50],[i j k l]);
            grid on
            %搞清楚你畫的是哪三維的關(guān)系
            plot3(y(:,1),y(:,2),y(:,4),'linewidth',1);
            set(gca,'XTick',[0:0.2:1],'YTick',[0:0.2:1],'ZTick',[0:0.2:1])
            set(gca,'XTickLabel',num2str(get(gca,'XTick')','%.1f'));
            set(gca,'YTickLabel',num2str(get(gca,'YTick')','%.1f'));
            set(gca,'ZTickLabel',num2str(get(gca,'ZTick')','%.1f'));
            hold on
            axis([0 1 0 1 0 1])
            end
        end
    end
end
xlabel('$r$','interpreter','latex');
ylabel('$m$','interpreter','latex');
zlabel('$p$','interpreter','latex','Rotation',360);
title('{\it\fontname{Bodoni MT}g}=0.8');

image.png

出圖效果還是有差距诅蝶，在編輯器中修改。

image.png

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者

人面猴
序言：七十年代末募壕，一起剝皮案震驚了整個濱河市调炬，隨后出現(xiàn)的幾起案子，更是在濱河造成了極大的恐慌舱馅，老刑警劉巖缰泡，帶你破解...
沈念sama閱讀 218,546評論 6贊 507
死咒
序言：濱河連續(xù)發(fā)生了三起死亡事件，死亡現(xiàn)場離奇詭異代嗤，居然都是意外死亡棘钞，警方通過查閱死者的電腦和手機缠借，發(fā)現(xiàn)死者居然都...
沈念sama閱讀 93,224評論 3贊 395
救了他兩次的神仙讓他今天三更去死
文/潘曉璐我一進店門，熙熙樓的掌柜王于貴愁眉苦臉地迎上來宜猜，“玉大人泼返，你說我怎么就攤上這事”Χ瘢” “怎么了符隙？”我有些...
開封第一講書人閱讀 164,911評論 0贊 354
道士緝兇錄：失蹤的賣姜人
文/不壞的土叔我叫張陵趴捅，是天一觀的道長垫毙。經(jīng)常有香客問我，道長拱绑，這世上最難降的妖魔是什么综芥？我笑而不...
開封第一講書人閱讀 58,737評論 1贊 294
?港島之戀（遺憾婚禮）
正文為了忘掉前任，我火速辦了婚禮猎拨，結(jié)果婚禮上膀藐，老公的妹妹穿的比我還像新娘。我一直安慰自己红省，他們只是感情好额各，可當我...
茶點故事閱讀 67,753評論 6贊 392
惡毒庶女頂嫁案：這布局不是一般人想出來的
文/花漫我一把揭開白布。她就那樣靜靜地躺著吧恃，像睡著了一般虾啦。火紅的嫁衣襯著肌膚如雪。梳的紋絲不亂的頭發(fā)上痕寓，一...
開封第一講書人閱讀 51,598評論 1贊 305
城市分裂傳說
那天傲醉，我揣著相機與錄音，去河邊找鬼呻率。笑死硬毕，一個胖子當著我的面吹牛，可吹牛的內(nèi)容都是我干的礼仗。我是一名探鬼主播吐咳，決...
沈念sama閱讀 40,338評論 3贊 418
雙鴛鴦連環(huán)套：你想象不到人心有多黑
文/蒼蘭香墨我猛地睜開眼，長吁一口氣：“原來是場噩夢啊……” “哼元践！你這毒婦竟也來了挪丢？” 一聲冷哼從身側(cè)響起，我...
開封第一講書人閱讀 39,249評論 0贊 276
萬榮殺人案實錄
序言：老撾萬榮一對情侶失蹤卢厂，失蹤者是張志新（化名）和其女友劉穎乾蓬，沒想到半個月后，有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體慎恒，經(jīng)...
沈念sama閱讀 45,696評論 1贊 314
?護林員之死
正文獨居荒郊野嶺守林人離奇死亡任内，尸身上長有42處帶血的膿包…… 初始之章·張勛以下內(nèi)容為張勛視角年9月15日...
茶點故事閱讀 37,888評論 3贊 336
?白月光啟示錄
正文我和宋清朗相戀三年撵渡，在試婚紗的時候發(fā)現(xiàn)自己被綠了。大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片死嗦。...
茶點故事閱讀 40,013評論 1贊 348
活死人
序言：一個原本活蹦亂跳的男人離奇死亡趋距，死狀恐怖，靈堂內(nèi)的尸體忽然破棺而出越除，到底是詐尸還是另有隱情节腐，我是刑警寧澤，帶...
沈念sama閱讀 35,731評論 5贊 346
?日本核電站爆炸內(nèi)幕
正文年R本政府宣布摘盆，位于F島的核電站翼雀，受9級特大地震影響，放射性物質(zhì)發(fā)生泄漏孩擂。R本人自食惡果不足惜狼渊，卻給世界環(huán)境...
茶點故事閱讀 41,348評論 3贊 330
男人毒藥：我在死后第九天來索命
文/蒙蒙一、第九天我趴在偏房一處隱蔽的房頂上張望类垦。院中可真熱鬧狈邑，春花似錦、人聲如沸蚤认。這莊子的主人今日做“春日...
開封第一講書人閱讀 31,929評論 0贊 22
一樁弒父案，背后竟有這般陰謀
文/蒼蘭香墨我抬頭看了看天上的太陽砰琢。三九已至蘸嘶，卻和暖如春，著一層夾襖步出監(jiān)牢的瞬間氯析，已是汗流浹背亏较。一陣腳步聲響...
開封第一講書人閱讀 33,048評論 1贊 270
情欲美人皮
我被黑心中介騙來泰國打工，沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留掩缓，地道東北人雪情。一個月前我還...
沈念sama閱讀 48,203評論 3贊 370
代替公主和親
正文我出身青樓，卻偏偏與公主長得像你辣，于是被迫代替她去往敵國和親巡通。傳聞我的和親對象是個殘疾皇子，可洞房花燭夜當晚...
茶點故事閱讀 44,960評論 2贊 355

4. 四方演化博弈代碼

推薦閱讀更多精彩內(nèi)容