像上面我举的例子
常微分方程为5x”+x’+5x=0 初速度v0=0,初位置x0=10
先将其化为一阶方程组
x'=v x0=10
v'=0.2v-x v0=0
在程序中x(1)相当于x,x(2)相当于v
%jswlx_8_4_3f.m 阻尼振动函数
function y=jswlx_8_4_3f(t,x)
y=[x(2);-0.2*x(2)-x(1)] %一阶方程组,列向量
%jswlx_8_4_3.m
Hf_1=figure;
set(Hf_1,'color','w')
tDelt=0.05;T=70.0;%变量t的间隔和最终T
x10=10;x20=0.0;x0=[x10,x20];%初始条件
ts=0:tDelt:T;
[t,x]=ode45('jswlx_8_4_3f',ts,x0)
plot(t,x(:,1)),%grid,
xlabel('T'),ylabel('X')
像上面我举的例子
常微分方程为5x”+x’+5x=0 初速度v0=0,初位置x0=10
先将其化为一阶方程组
x'=v x0=10
v'=0.2v-x v0=0
在程序中x(1)相当于x,x(2)相当于v
%jswlx_8_4_3f.m 阻尼振动函数
function y=jswlx_8_4_3f(t,x)
y=[x(2);-0.2*x(2)-x(1)] %一阶方程组,列向量
%jswlx_8_4_3.m
Hf_1=figure;
set(Hf_1,'color','w')
tDelt=0.05;T=70.0;%变量t的间隔和最终T
x10=10;x20=0.0;x0=[x10,x20];%初始条件
ts=0:tDelt:T;
[t,x]=ode45('jswlx_8_4_3f',ts,x0)
plot(t,x(:,1)),%grid,
xlabel('T'),ylabel('X')