matlab基础——在工作区内显示函数内的变量
注意：在通常情况下，matlab的工作区中不会显示函数内的变量。
以如下的SEIR模型为例：
function SEIR(xN,xE,xI,xR,xr,xr2,xB,xB2,xa,xy)
%第一阶段fun1(837300,0,1,0,10,10,0.045,0.045,0.1,0.9)
N = xN; %总人数
E = xE; %潜伏人数
I = xI; %感染人数
R = xR; %康复人数
S = N - I-R; %易感人数
r = xr; %一个感染者的接触人数
r2 = xr2 ; %一个潜伏者的接触人数
B = xB;%x1;%0.045; %接触感染者，易感到潜伏的传染率
B2 = xB2; %接触潜伏者，易感到潜伏的传染率
a = xa; %潜伏到感染的传染率
y = xy; %感染者的治愈率
T = 1:150;%天数
for i = 1:length(T)-1 %也就是1：149
S(i+1) = S(i) - r*B*S(i)*I(i)/N-r2*B2*S(i)*E(i)/N ;
E(i+1) = E(i) + r*B*S(i)*I(i)/N(1)+r2*B2*S(i)*E(i)/N(1)-a*E(i);
I(i+1) = I(i) + a*E(i) - y*I(i);
R(i+1) = R(i) + y*I(i);
Y(i+1)=I(i+1)+I(i)+ R(i+1);
end
subplot(1,2,1);
plot(T,Y,'r');grid on;%绘制多条曲线并显示网格
hold on;
subplot(1,2,2);
plot(T,I,'r');grid on;%绘制多条曲线并显示网格
hold on;
end
执行代码发现工作区不显示数组Y,I,T的值：
解决方法：
以函数输出的形式定义希望输出的变量：
function [A1,B1,C1]=SEIR(xN,xE,xI,xR,xr,xr2,xB,xB2,xa,xy)
%
% xr,xr2大于0 xB,xB2,xa,xy都属于[0，1]
%第一阶段fun1(837300,0,1,0,10,10,0.045,0.045,0.1,0.9)
N = xN; %总人数
E = xE; %潜伏人数
I = xI; %感染人数
R = xR; %康复人数
S = N - I-R; %易感人数
r = xr; %一个感染者的接触人数
r2 = xr2 ; %一个潜伏者的接触人数
B = xB;%x1;%0.045; %接触感染者，易感到潜伏的传染率
B2 = xB2; %接触潜伏者，易感到潜伏的传染率
a = xa; %潜伏到感染的传染率
y = xy; %感染者的治愈率
T = 1:150;%天数
for i = 1:length(T)-1 %也就是1：149
S(i+1) = S(i) - r*B*S(i)*I(i)/N-r2*B2*S(i)*E(i)/N ;
E(i+1) = E(i) + r*B*S(i)*I(i)/N(1)+r2*B2*S(i)*E(i)/N(1)-a*E(i);
I(i+1) = I(i) + a*E(i) - y*I(i);
R(i+1) = R(i) + y*I(i);
Y(i+1)=I(i+1)+I(i)+ R(i+1);
end
%
%经典SEIR模型的差分方程（间隔为1天）
A1 = Y;
B1 = I;
C1 = T;
subplot(1,2,1);
plot(T,Y,'r');grid on;%绘制多条曲线并显示网格
hold on;
subplot(1,2,2);
plot(T,I,'r');grid on;%绘制多条曲线并显示网格
hold on;
end
改变的地方在这：
参考文章