对一个圆盘进行rcs计算，给出代码，需要的参数设置合理即可

对一个圆盘进行雷达截面积（Radar Cross Section, RCS）的计算，通常采用瑞利散射理论。在这个简单场景中，我们可以假设圆盘是完美的电导体，并且忽略边缘效应。以下是一个基础的MATLAB代码示例，用于计算一个圆形目标的RCS：

function rcs = disk_rcs(freq, wavelength, radius, incidence_angle)
% 定义常量
pi = pi;
sigma_0 = (4 * pi) / (wavelength^2); % 平方厘米每赫兹
cos_theta = cos(incidence_angle); % 入射角余弦值

% 圆盘半径到直径转换
diameter = 2 * radius;

% 计算阴影区域（减去圆盘部分）
shadow_area = (pi * diameter^2) / 4; % 圆盘面积

% 根据瑞利散射公式计算RCS
if cos_theta == 0 % 如果入射垂直于圆盘面
    rcs = sigma_0;
else
    rcs = sigma_0 * (4 - 2*cos_theta^2) / (1 + cos_theta^2); % 主瓣最大值，瑞利散射公式
    rcs = max(rcs, 0); % 因为主瓣范围内RCS为正值
    rcs = rcs - shadow_area; % 减去圆盘自身阴影部分
end
end

% 示例使用
freq_hz = 1e9; % 频率（赫兹）
lambda_m = freq_hz / 3e8; % 波长（米）
radius_cm = 5; % 直径半径（厘米）
incidence_degrees = 30; % 入射角（度）

wavelength_rad = lambda_m * pi / 180; % 转换为弧度
incidence_rad = incidence_degrees * pi / 180;

[rcs] = disk_rcs(freq_hz, wavelength_rad, radius_cm, incidence_rad);
disp(['RCS at ' num2str(freq_hz) ' Hz and ' num2str(incidence_degrees) ' degrees is: ' num2str(rcs) ' cm^2']);

%