信号的时域特征分析matlab程序

%信号的时域特征分析
N=length(x);
t=linspace(0,1,8192);
%时域图
figure;plot(t,x);xlabel(‘Time’);ylabel(‘Amplitude’);title(‘时域图’)
%信号的平均值
Xmean=sum(x)/N
%信号的有效值
Xrms=sqrt(sum(x.^2)/N)
%信号的方根幅值
Xr=(sum(abs(x))/N)^2
%信号的绝对平均幅值
abX=sum(abs(x))/N
%信号的歪度
alpha=sum(x.^3)/N
%信号的峭度
beta=sum(x.^4)/N
%无量纲的统计特征值 波形指标
Sf=Xrms/abX
%峰值指标
Xmax=max(x);
Cf=Xmax/Xrms
%脉冲指标
If=Xmax/abX
%裕度指标
CLf=Xmax/Xr
%峭度指标
Kv=beta/(Xrms^4)
%概率密度函数
[f1,xi]=ksdensity(x);figure;plot(xi,f1);title(‘概率密度函数’);
%信号的频域分析
fs=8192;
f=linspace(0,fs,N);
%功率谱
figure;y=abs(fft(x));plot(f(1:N/2),y(1:N/2));xlabel(‘Frequency’);ylabel(‘Amplitude’);title(‘频域图’);
Sx=(y(1:N/2)).^2;figure;plot(f(1:N/2),Sx);xlabel(‘Frequency’);ylabel(‘Amplitude’);title(‘功率谱’);