matlab 三维 物体 运动 仿真

n=1; % 阻尼系数
w=6; % 圆频率
ww=sqrt（w^2-n^2）; % 等效的圆频率
B1=80; % 系数
B2=-20; % 系数
axis（[-105105-105105]）;hold on; % 设置坐标轴范围
set（gcf‘doublebuffer‘‘on‘）; % 设置图形窗口的渲染效果
G=plot（1000‘r‘）; % 轨迹曲线
x=100; % 初始位置的横坐标
y=0; % 初始位置的纵坐标
r=2; % 小球半径
t1=r*sin（0:.1:2*pi）; % 计算小球的边界数据
t2=r*cos（0:.1:2*pi）; % 计算小球的边界数据
X=x+r*t1; % 得到小球的横坐标数据
Y=y+r*t2; % 得到小球的纵坐标数据
Q=fill（XYrand（13））; % 画出小球
t=0; % 时间的初始数值
plot（[-100100][00]）; % 画出格线
plot（[00][-100100]）; % 画出格线
k=1; % 控制循环的参数
xlabel（‘Please press “space“ key and stop this program!‘...
    ‘fontsize‘14‘color‘‘r‘）; % 添加提示字符
title（[‘《damp turbination》 ‘‘B1=