matlab解微分方程组（自己编4阶龙格库塔法），步长真么选取？
其中：1.r=0.001,a=2*104,s=0.015,b=10-4;x(0)=1200,y(0)=6002.r=0.001,a=2*10-6,s=0.01,b=10-...
其中：
1. r=0.001,a=2*104, s=0.015, b=10-4; x(0)=1200, y(0)=600
2. r=0.001,a=2*10-6, s=0.01, b=10-6; x(0)=12000, y(0)=6000 展开
1. r=0.001,a=2*104, s=0.015, b=10-4; x(0)=1200, y(0)=600
2. r=0.001,a=2*10-6, s=0.01, b=10-6; x(0)=12000, y(0)=6000 展开
2个回答
展开全部
步长是选时间t的间隔，看你要计算多长的时间段t：[t0，tn]，一般t0=0，tn由题意或自己选取最多20s就可以了，步长选取h=0.01。
本回答被提问者采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
富港检测技术（东莞）有限公司_
2024-04-02 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大（共振点）；正弦振动在任一瞬间只包含一种频率的振动，而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率，所以样品上的共振点会同时激发...
点击进入详情页
本回答由富港检测技术（东莞）有限公司_提供
展开全部
function R=rk4(f,a,b,ya,N)
%y'=f(x,y)
%a,b左右端点
%N为迭代步长
%h为步长
%ya为初值
h=(b-a)/N;
T=zeros(1,N+1);
Y=zeros(1,N+1);
T=a:h:b;
Y(1)=ya;
for j=1:N
k1=h*feval(f,T(j),Y(j));
k2=h*feval(f,T(j)+h/2,Y(j)+k1/2);
k3=h*feval(f,T(j)+h/2,Y(j)+k2/2);
k4=h*feval(f,T(j)+h/2,Y(j)+k3);
Y(j+1)=Y(j)+(k1+2*k2+2*k3+k4)/6;
end
R=[T' Y']
%y'=f(x,y)
%a,b左右端点
%N为迭代步长
%h为步长
%ya为初值
h=(b-a)/N;
T=zeros(1,N+1);
Y=zeros(1,N+1);
T=a:h:b;
Y(1)=ya;
for j=1:N
k1=h*feval(f,T(j),Y(j));
k2=h*feval(f,T(j)+h/2,Y(j)+k1/2);
k3=h*feval(f,T(j)+h/2,Y(j)+k2/2);
k4=h*feval(f,T(j)+h/2,Y(j)+k3);
Y(j+1)=Y(j)+(k1+2*k2+2*k3+k4)/6;
end
R=[T' Y']
本回答被网友采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
