λ 的值没给,假设为1,记=1.
首先编辑如下的函数文件并保存:
function y=hanshu(t,x)
u=0.25;w0=0.4;w=0.5;la=1;
y=zeros(2,1);
y(1)=x(2);
y(2)=la*sin(w*t)-2*u*x(2)-w0^2*x(1);
之后在命令窗口输入:
[t,y]=ode45(@hanshu,[0 300],[1 -1]);
plot(t,y(:,1));
xlabel('t')
ylabel('xd')
即可得到结果,这里只计算了图t=0到300的值了.