版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (2288)
>考研 (265)
>導師招生 (213)
>蟲友互識 (56)
>休閑灌水 (32)
>基金申請 (31)
>考博 (21)
>文獻求助 (18)
>博后之家 (15)
>論文道賀祈福 (15)
>教師之家 (15)
>論文投稿 (15)
>碩博家園 (9)
>找工作 (9)
>公派出國 (9)
>生物科學 (5)

北京石油化工學院2026年研究生招生接收調劑公告

返回列表

mathda

新蟲 (小有名氣)

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

[求助] 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

回復此樓

» 猜你喜歡

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

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

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

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

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

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

【答案】應助回帖

感謝參與，應助指數 +1

else里面的初始y設置錯誤，第一次計算的結果y應該是多行2列，在else里再整列調用賦給方程的初值就錯了。應該選取一行吧？

贊一下(1人)

回復此樓

showmethemoney

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

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

相關版塊跳轉我要訂閱樓主 mathda 的主題更新