| 5 | 1/1 | 返回列表 |
| 查看: 362 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
yangnan26金蟲 (小有名氣)
|
[求助]
請(qǐng)幫忙看看下面的程序有沒有問題! 已有2人參與
|
|
|
編了一段matlab程序但運(yùn)行出不來想要的效果,不知道問題在哪?請(qǐng)各位蟲友幫忙看看,謝謝! close all; clear all; delta_t=0.05;%離散時(shí)間間隔 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兩個(gè)方程 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)') |
至尊木蟲 (著名寫手)
驃騎將軍

新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9 | Ymlll 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研]
|
.6lL 2026-03-18 | 6/300 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|