Matlab中用ode函数求二元微分方程，方程中有一个范围的常数怎么chu'li？

为什么Matlab中用ode函数求二元微分方程，方程中有一个范围的常数b 0，1会出现错误？由于b 在微分方程中要求是一个常数，那如何来处理b是一个向量数组，此时我们可以采用循环的方法，求出当b为某一个常数时的t，y,然后绘出其图形。

实现部分代码：

b0 = 0:05:1;

y0=[20,15];

for i=1:length(b0)

b=b0(i);

[t,y]=ode45(@(t,y)hopf(t,y,b),[0,20],y0); %ode45

plot(y(:,1),y(:,2));hold on

xlabel('x'),ylabel('y')

end

function dydt = hopf(t,y,b)

dydt(1) = y(1) - 01y(1)^2 - (by(1)y(2))/(1 + b04y(1));

dydt(2) = -y(2) - 01y(2)^2 + (by(1)y(2))/(1 + b04y(1));

dydt = dydt(:);

end

运行结果

“矗立”指高耸地立着，如:大街两旁矗立着高楼大厦;“屹立”指像山峰一样高耸而稳固地立着，常用来比喻坚定不可动摇，如:人民英雄纪念碑屹立在天安门广场上;“伫立”指长时间地站着，如:伫立窗前;“耸立”指高高地直立，如:群山耸立等。

屹立，汉语词汇。

拼音：yì lì

释义：高耸挺立。像山峰一样高耸而稳固地挺立着，

通常用来比喻人的意志十分坚定无法动摇。

中文名

屹立

外文名

stand erect

拼音

yì lì

近义词

矗立，耸立，伫立，挺立

典源

宋·张世南 《游宦纪闻》卷六

出力[chū lì] 拿出力量；尽力：～不讨好。他为人耿直，干工作又肯～。

造句：大家都在为集体、为社会出力，你总打个人小算盘，不惭愧吗？

出力不讨好：成语，指白费功夫，没有取得成效，相反带来坏的影响。

解释：指白费功夫，没有取得成效，相反带来坏的影响。示例：我们没必要去干这～的事情。用法：作宾语、定语；指白费功夫

出力近义词:效能、效用、效劳、效力、功效、着力、效率

形容出力的成语：

通力合作：通力：一起出力。不分彼此，一齐出力。出自《论语·颜渊》：“盍彻乎”朱熹注：“一夫受田百亩，而与同沟共井之人通力合作，计亩均收。”

志士仁人：原指仁爱而有节操，能为正义牺牲生命的人。现在泛指爱国而为革命事业出力的人。出自《论语·卫灵公》：“志士仁人，无求生以害仁，有杀身以成仁。”

葆力之士：葆力：勤劳任力。指愿意出力而勤于办事的人。