请给出matlab完整程序，机载雷达辐射模型仿真，无源侦收机位置固定，机载雷达在侦收机斜上方匀速直线运动，需要考虑雷达天线方向图及波束指向的变化，波束指向为扇形扫描，无源侦收机的侦收结果作为输出，尽量贴近真实场景

在MATLAB中建立一个简单的机载雷达辐射模型仿真，我们需要考虑以下几个关键部分：

% 定义常量和变量
RadarPos = [0; 0; 150]; % 雷达的位置，假设在地面以上150米处
ReceiverPos = [0; 0; 0]; % 无源接收机位置，假设在地面上
Speed = 500; % 雷达的速度，单位是米每秒
ScanAngle = pi / 6; % 扇形扫描的角度，这里设置为30度
AntennaGain = 10; % 雷达天线增益

% 时间间隔和时间向量
dt = 0.1; % 每个时间步长，单位秒
timeSteps = 100; % 总时间步数
time = linspace(0, Speed * timeSteps, numel(timeSteps) + 1); % 时间范围

% 初始化接收信号矩阵
ReceivedSignal = zeros(size(time));

for t = 1:length(time)
    % 计算雷达与接收机之间的距离
    Range = norm(RadarPos - ReceiverPos);
    
    % 根据雷达速度计算雷达的位置
    RadarNewPos = RadarPos + Speed * [cosd(RadarPos(3) / Range) ; sind(RadarPos(3) / Range)] * dt;
    
    % 更新雷达位置
    RadarPos = RadarNewPos;
    
    % 计算天线的方向角
    Azimuth = atan2(RadarNewPos(1) - ReceiverPos(1), RadarNewPos(2) - ReceiverPos(2));
    
    % 考虑波束指向变化，仅在雷达上方区域接受信号
    if abs(Azimuth) <= ScanAngle
        ReceivedSignal(t) = AntennaGain / (Range^4); % 简单的瑞利衰落模型
    end
end

% 输出接收结果
figure;
plot(time, ReceivedSignal, 'LineWidth', 2);
xlabel('Time (s)');
ylabel('Received Signal Strength');
title('Machine Radar Scanning Simulation');

% 提供一些额外信息