alpha.m

来自「matlab中用于解决高等动力学hamilton系统的相关程序」· M 代码 · 共 8 行

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