写一段生成舰船辐射噪声的程序

生成舰船辐射噪声的程序需要考虑到舰船在不同工况下的噪声特征，并且需要使用到信号处理和数学模型等技术。下面是一个简单的程序示例，可以生成一个基本的舰船辐射噪声信号：

% 设置采样参数
fs = 44100; % 采样率
T = 1; % 采样时间

% 生成白噪声
noise = randn(1, fs*T);

% 设计滤波器
fc = 1000; % 截止频率
[b, a] = butter(4, fc/(fs/2));

% 滤波处理
filtered_noise = filter(b, a, noise);

% 加入舰船噪声模型
ship_noise = 0.5*sin(2*pi*1000*(0:1/fs:(T-1/fs))) + 0.3*sin(2*pi*2000*(0:1/fs:(T-1/fs))) + 0.2*sin(2*pi*3000*(0:1/fs:(T-1/fs)));

% 加权合成
weighted_noise = 0.7*filtered_noise + 0.3*ship_noise;

% 播放
soundsc(weighted_noise, fs);

这个程序首先生成了一个长度为`fs*T`的白噪声信号，然后使用Butterworth滤波器对其进行滤波处理，截止频率为1kHz。接着加入了一个简单的舰船噪声模型，包括3个正弦波，分别频率为1kHz、2kHz和3kHz，振幅权重分别为0.5、0.3和0.2。最后对滤波后的白噪声和舰船噪声进行加权合成，最终得到舰船辐射噪声信号并播放。

需要注意的是，这只是一个简单的示例程序，实际的舰船辐射噪声信号可能更加复杂，生成程序需要根据实际情况进行调整和优化。