| 查看: 346 | 回复: 0 | ||
[求助]
求解二元微分方程组
|
|
请看我写的主函数： clc; clear; tp=0.1; n=5; ttp=tp/n; tspan=[0:ttp:tp]; y0=[0,4.5,0,4.5]'; [t y]=ode45(@myfunn,tspan,y0); plot(4.5*t/0.075,y(:,4)/4.5); 下面是myfunn子函数的代码： function dy=myfunn(t,y) dy=zeros(4,1); Hf=0.075; M2=80; k=353.33*353.33*M2; c=0.2*2*(k*M2)^0.5; Mf=M2/50; A=Mf/(2700*Hf); M1=250*Mf; ee=y(1)/Hf; Y=0； if y(1)>=Y p=3.5e6+0.6e6*ee/(1-ee)^2; else p=0 ; end dy(1)=y(2); dy(2)=c/M1*(y(4)-y(2))+k/M1*(y(3)-y(1))-A*p/M1+9.8; dy(3)=y(4); dy(4)=-c/M2*(y(4)-y(2))-k/M2*(y(3)-y(1))+9.8; end if y(1)>=Y Y=y(1); else Y=Y; end 我的问题来了，如下： 在这个微分方程中，p值是随y(1)的变化而变化的。假如在第t步计算中，y(1)>Y则p=3.5e6+0.6e6*ee/(1-ee)^2;当y(1)<Y时，p=0。然后计算(t+1)步的y(1)值。把(t+1)步的y(1)和Y对比。当y(1)>Y,则用y(1)替换Y。我写的代码中每步的Y初始值都是0开始，但是每一次y(1)与Y对比之后，Y的值有可能发生变化。而且我的代码中y(1)在每一次调用子函数的过程中都是不变的。 请问怎么解决？ |
回复此楼» 猜你喜欢
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有18人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有119人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有0人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复