应用matlab求下式积分：9、一物体做平抛运动从高度为200m位置水平抛出抛出点的位置坐标为0200m在位置90m90m画圆圈将圆圈设置为合适的大小手动调节抛出时的水平速度参数使得抛出的物体的轨迹穿过圆圈只要轨迹穿过圆圈即可要求图中添加横坐标注释、纵坐标注释然后保存程序和运行得到的图。重力加速度取10 ms2不考虑空气的阻力作用用matlab编程来实现。请给出完整的MATLAB具体程序

下面是一个MATLAB程序来实现上述要求：

% 物体平抛运动轨迹
function projectile_motion()

    % 初始化参数
    g = 10; % 重力加速度
    h0 = 200; % 初始高度
    x0 = 0; % 初始水平位置
    y0 = h0; % 初始垂直位置
    x1 = 90; % 圆圈的水平位置
    y1 = 90; % 圆圈的垂直位置

    % 手动调节水平速度参数
    vx = 10; % 初始水平速度

    % 设置时间步长
    dt = 0.01;

    % 计算物体的运动轨迹
    t = 0;
    x = x0;
    y = y0;
    while y >= 0
        % 计算加速度
        ax = 0;
        ay = -g;

        % 更新速度
        vx = vx + ax * dt;
        vy = 0;

        % 更新位置
        x = x + vx * dt;
        y = y + vy * dt;

        % 更新时间
        t = t + dt;

        % 判断物体是否穿过圆圈
        if (x - x1)^2 + (y - y1)^2 <= 10^2
            break;
        end
    end

    % 绘制轨迹
    figure;
    plot([x0, x], [y0, y]);
    hold on;
    rectangle('Position',[x1-10, y1-10, 20, 20],'Curvature',[1,1],'EdgeColor','r');
    xlabel('横坐标');
    ylabel('纵坐标');
    title('物体平抛运动轨迹');
    legend('轨迹','圆圈');
    grid on;

    % 保存程序和运行得到的图
    saveas(gcf, 'projectile_motion.png');
end

将以上代码保存为一个.m文件，例如projectile_motion.m，然后在MATLAB命令窗口中运行projectile_motion函数即可得到运行结果，并保存为projectile_motion.png图像文件