舰船辐射噪声的生成和建模仿真

主要介绍舰船辐射噪声的生成，包括宽带连续谱和线谱的建模仿真，得到了不错的效果

clc
clear all
clf

Fs=5000;
n=4;%螺旋桨叶片数
s=4;%螺旋桨的转速
%%%%100Hz以下线谱频率可表示为如下%%%%
Ts=0:1/Fs:2.5;
ni=8;%八条线谱
f1=zeros（1ni）;
C=[];
for m1=1:ni
    f1=m1*n*s;
    C（:m1）=f1;
  %%%线谱频率
end
f=C
Ak=0.5*rand（110）%%线谱幅度
Ik=0.25*pi*rand（110）%%%%随机相位
%%%K为线谱数，AkfIk分别为第k条线谱对应的幅度、频率、相位
%%求取功率谱%%%
fi=80;
L=4;%线谱最高频率不能超过100Hz最高为ni*fi/L
fait=0.25*pi;
temp=1:1:ni;
G=0;
for i=1:1:ni              
   G1=Ak（i）*sin（temp（i）*（2*pi*fi/L*Ts+fait）+Ik（i））;
   G=G+G1;
end

G=G+7.5*sin（2*pi*350*Ts）+6.5*sin（2*pi*700*Ts）;
% G7=Ak（17）*sin（2*pi*350*n*Ts‘+Ik（17））;
% G8=Ak（18）*sin（2*pi*500*n*Ts‘+Ik（18））;
% G9=Ak（19）*sin（2*pi*700*n*Ts‘+Ik（19））;
% G10=Ak（110）*sin（2*pi*900*n*Ts‘+Ik（110））;
% G=G2‘+G7+G8+G9+G10;
figure（1）
plot（TsG‘k‘）
xlabel（‘时间t/s‘）
ylabel（‘幅值/V‘）
title（‘线

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        5263  2015-06-17 10:03  jianchuanfushezaoshengmoxing.m