| 1 | 1/1 | 返回列表 |
| 查看: 378 | 回復(fù): 0 | ||
mayaya26木蟲 (小有名氣)
|
[求助]
混沌仿真
|
|
子程序 function dy=Bao2(t,y) dy=zeros(3,1); global a b c a=1.812;c=0.5; dy(1)=(b/(a*(1-c)))*(-exp(-y(2))+1-y(3)); dy(2)=b/(a*c)*y(3); dy(3)=-(a*c*(1-c)/b)*(y(1)+y(2))-y(3)/a; 主程序 clear format long global a b c b1 dt=0.001; tspan=(0:dt:50); Y0=[5 5 5]; for i=1:1:1001 b=0.015*(i-1); sol=ode45(@Bao2,[tspan(1) tspan(end)],Y0); ss=deval(sol,tspan)'; s=ss(20000:end, ;aa=1; bb=0; cc=0; p=-1; func=((aa*s(:,1))+(bb*s(:,2))+(cc*s(:,3))+p)>0; func_n=(func(1:end-1)>0)&(xor(func(1:end-1),func(2:end))); m=find(func_n==1); t_intsec = -(aa*s(m,1)+bb*s(m,2)+cc*s(m,3)+p)./(aa*(s(m+1,1)-s(m,1))+bb*(s(m+1,2)-s(m,2))+cc*(s(m+1,3)-s(m,3))); z_map=s(m,3)+(s(m+1,3)-s(m,3)).*t_intsec; z1_map=log10(z_map); b1=log10(b*ones(size(s(m,3)))) H=plot(b1,z1_map); set(H,'linestyle','none','marker','.','markersize',1) hold on drawnow end 但是出現(xiàn)錯(cuò)誤 b1 = Empty matrix: 0-by-1 為什么?是初始值設(shè)錯(cuò)了?還是...... 有誰幫忙,急急急,謝謝 |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601 329分調(diào)劑 +3 | yzsa12 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +7 | JKSOIID 2026-03-26 | 7/350 |
|
|
[考研]
|
鐘llll 2026-03-26 | 5/250 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
[考研]
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
[考研] 292求調(diào)劑 +13 | 是妍子也是研子 2026-03-30 | 13/650 |
|
|
[考研] 310求調(diào)劑 +10 | 爭取九點(diǎn)睡 2026-03-30 | 10/500 |
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 化學(xué)308分調(diào)劑 +7 | 你好明天你好 2026-03-30 | 8/400 |
|
|
[考研] 375求調(diào)劑 +6 | 雨夏整夜 2026-03-29 | 6/300 |
|
|
[考研] 317求調(diào)劑 +10 | 蛋黃咸肉粽 2026-03-26 | 10/500 |
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研] 296求調(diào)劑 +4 | 汪。! 2026-03-25 | 7/350 |
|