应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 请帮忙看看下面的程序有没有问题！
31/1返回列表
查看: 324  |  回复: 2
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yangnan26

金虫 (小有名气)

[求助] 请帮忙看看下面的程序有没有问题！已有2人参与

编了一段matlab程序但运行出不来想要的效果，不知道问题在哪？请各位虫友帮忙看看，谢谢！
close all;
clear all;
delta_t=0.05;%离散时间间隔
N=10000;M=4000;
S=zeros(N+1,1);
P=zeros(N+1,1);
I=zeros(N+1,M+1);
dI=zeros(M+1,1);
beta=zeros(M+1,1);

tao=0.15;
S(1)=9.2;P(1)=0.8;
a=0:M;I(1,：)=1.472*exp(-0.5*a*delta_t);%I(0,a)
dI(tao/delta_t:end)=0.5;
beta(tao/delta_t:end)=0.3;


% 求解1,3两个方程
for t=1:N
    S(t+1)=S(t)+(0.2*S(t)-0.2*P(t)*S(t))*delta_t;
    P(t+1)=P(t)+(0.5-0.7*P(t)-0.1*P(t)*S(t)+I(t,: )*beta*delta_t)*delta_t;
    I(t+1,1)=0.2*P(t)*S(t);
    for a=2:M+1
        I(t+1,a)=I(t,a-1)-dI(a)*I(t,a)*delta_t;
    end
end

X=1:M+1;
Y=1:N+1;
[X,Y]=meshgrid(X,Y);
surf(X,Y,I)
title('I(t,a)')


figure;
plot(S)
title('S(t)')
figure;
plot(P)
title('P(t)')
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-04-09 13:39:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qq2537742567

新虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
yangnan26: 金币+2 2014-04-09 21:48:59

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
yangnan26: 金币+2 2014-04-09 21:49:08
相关版块跳转 我要订阅楼主 yangnan26 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭