24小时热门版块排行榜    

查看: 357  |  回复: 2

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)')
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qq2537742567

新虫 (小有名气)


【答案】应助回帖

★ ★
感谢参与,应助指数 +1
yangnan26: 金币+2 2014-04-09 21:48:59
建议你最好加上注释和程序背景。不然很不容易看懂,也不理解你要达到的效果。
2楼2014-04-09 15:17:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
yangnan26: 金币+2 2014-04-09 21:49:08
把报错信息贴出来看看才知道问题在哪里啊
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
3楼2014-04-09 20:27:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangnan26 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿武汉理工材料工程专硕调剂 +5 Doleres 2026-03-19 5/250 2026-03-19 20:14 by 制度的
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 12/600 2026-03-19 16:17 by 脱颖而出
[考研] 304求调剂 +3 曼殊2266 2026-03-18 3/150 2026-03-19 14:42 by peike
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 346求调剂[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[考研] 一志愿天大材料与化工(085600)总分338 +5 蔡大美女 2026-03-13 5/250 2026-03-19 10:44 by 是小刘呀~
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 265求调剂 +3 梁梁校校 2026-03-17 3/150 2026-03-18 09:12 by zhukairuo
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
信息提示
请填处理意见