clc
clear?all
clf

Fs=5000;
n=4;%螺旋槳葉片數(shù)
s=4;%螺旋槳的轉(zhuǎn)速
%%%%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）%%%%隨機(jī)相位
%%%K為線譜數(shù)，AkfIk分別為第k條線譜對(duì)應(yīng)的幅度、頻率、相位
%%求取功率譜%%%
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（‘時(shí)間t/s‘）
ylabel（‘幅值/V‘）
title（‘線

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????5263??2015-06-17?10:03??jianchuanfushezaoshengmoxing.m