ode_weight.asv

matlab从入门到精通第6章书后源码对初学者很有帮助

% ode_weight.m
% 求解加权函数微分方程组
m1=0.1;m2=0.1;L=1;g=9.81;
% 微分方程的参数
tspan=linspace(0,4,25)l
y0=[0;4;2;20;-pi/2;2];
options=odeset('Mass',@mass);
% 求解加权函数微分方程组
theta=y(1,5);X=y(1,1);Y=y(1,3);
xvals=[X X+L*cos(theta)];
yvals=[Y Y+L*sin(theta)];
% 绘制变量的图形
figure;
plot(xvals,yvals,xvals(1),yvals(1),'ro',xvals(2),yvals(2),'bo'
ti)