alpha.m

matlab中用于解决高等动力学hamilton系统的相关程序

clc
h=2;
r=1;
x=[-50:.01:50];
y = atan2((r-1/2*h./(x.^2+h^2).*(2*h*r - 2*sign(x).*(x.^4-x.^2*r^2+h^2*x.^2).^(1/2)))./x, 1/2./(x.^2+h^2).*(2*h*r - 2*sign(x).*(x.^4-x.^2*r^2+h^2*x.^2).^(1/2)))*180/pi

plot(x,y)
grid