版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

mathda

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 514.6
散金: 5
帖子: 55
在線: 52.6小時
蟲號: 1043825
注冊: 2010-06-19
性別: GG
專業(yè): 常微分方程與動力系統(tǒng)

[求助] matlab 求脈沖微分方程

程序代碼運行出錯，求教
function predprey
a=0.25;b=-0.01;c=-1.0;d=0.01; impcount=0;
Y=-1.2*c/d; p=0.2; tstar=0; tfinal=100;
while true
if (impcount ==0 )
options = ddeset('Events',@events,'AbsTol',1e-9,...
'RelTol',1e-9);
[t,y] = ode15s(@ddes,[tstar,tfinal],[80,30], options);
plot(y(:,1),y(:,2));
else
% Specify the new solution at impulse times...
options = ddeset(options, 'InitialY',...
[(1-p).*y(:,1); y(:,2)]);
[t,y] = ode15s(@ddes, [tstar,tfinal],...
[(1-p).*y(:,1);y(:,2)],options);
plot(y(:,1),y(:,2));
end
tstar = max(t);
impcount = impcount + 1;
if (tstar >= tfinal)
break;
end
end
%===Nested functions==============================
% Evaluate the DDES...
function dydt = ddes(t,y)
dydt = [a*y(1)+b*y(1)*y(2); ...
c*y(2)+d*y(1)*y(2)];
end
% Evaluate the event function residuals...
function [value,isterminal,direction] = events(t,y)
value = y(1)-Y;
isterminal = 1;
direction = 1;
end
end

回復(fù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有209人回復(fù)
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

求助：如何用MATLAB的偏微分方程工具箱畫出不同區(qū)域內(nèi)的溫度分布圖？已經(jīng)有5人回復(fù)
matlab求解偏微分方程組已經(jīng)有3人回復(fù)
求助matlab求解復(fù)雜微分方程已經(jīng)有7人回復(fù)
高手進(jìn)來用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較已經(jīng)有4人回復(fù)
matlab 線性常微分方程求解已經(jīng)有5人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復(fù)
高金求助matlab解微分方程組已經(jīng)有12人回復(fù)
matlab解微分方程組已經(jīng)有15人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)

1樓 2012-11-01 16:10:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強(qiáng)帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

else里面的初始y設(shè)置錯誤，第一次計算的結(jié)果y應(yīng)該是多行2列，在else里再整列調(diào)用賦給方程的初值就錯了。應(yīng)該選取一行吧？

贊一下(1人)

回復(fù)此樓

showmethemoney

2樓2012-11-01 16:30:44