| 查看: 436 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
请帮忙看看下面的程序有没有问题！ 已有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)') |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 符号方程组 4个方程 3个参数 求哥哥姐姐帮忙看一下哈 急~
已经有7人回复
- 请大家帮帮忙，关于CSC及申请流程的事，谢谢大家
已经有26人回复
- 求助大神帮我看看这段程序是干什么的？（就是定义这个函数算什么）
已经有6人回复
- 求助帮帮我编写一个MATLAB程序求解下面一个方程组！！多谢了
已经有5人回复
- 计算出错求帮忙
已经有6人回复
- 求帮忙计算 本人不会
已经有12人回复
- 帮忙看一下这条语句是什么意思
已经有7人回复
- linux终端运行程序出错，求大神帮忙
已经有6人回复
- 请帮忙编写一个perl程序处理文本
已经有26人回复
- 帮忙看看这个平面图用matlab怎么画
已经有19人回复
- 大家帮我看一个fortran的程序，我总是计算不出正确的结果
已经有20人回复
- 有会看fortran程序的，能帮忙看看这个程序吗？
已经有10人回复
- 求帮忙，关于含有bessel函数的方程的编程求解问题
已经有11人回复
- 【求助】50金币求大家帮忙改一段程序
已经有13人回复
- 【求助】关于向量投影的问题
已经有3人回复
- 【求助】帮忙看一个程序
已经有16人回复
- 【求助】帮忙改写展宽程序
已经有11人回复
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
3楼2014-04-09 20:27:54
2楼2014-04-09 15:17:24
