matlab 振幅谱函数,好用的画包络谱和幅值谱matlab函数

画包络函数：

function hua_baol (y,fs,style,varargin)

%画包络函数是hua_baoluo的代码简化版

%输入参数要么是3个，要么是5个

%当三个的时候

y_hht=hilbert(y);%希尔伯特变换

y_an=abs(y_hht);%包络信号

%y_an=y_an-mean(y_an);%去除直流分量

if nargin==3

hua_fft(y_an,fs,style);

elseif nargin==5

f1=varargin{1};

f2=varargin{2};

hua_fft(y_an,fs,style,f1,f2);

else

error('调用函数的输入参数数目不正确，输入参数只能是三个或者五个');

end

end

测试函数

x=@(t) (1+0.5*cos(9*pi*t)).*cos(200*pi*t+2*cos(10*pi*t))+sin(pi*t).*sin(30*pi*t);

t=0:0.01:9.99;

fs=1000;

y=x(t);

subplot(311);plot(t,y);

subplot(312);hua_fft(y,fs,1);%直接画出幅值谱

subplot(313);hua_fft(y,fs,1,0,250);%直接画出幅值谱，并且只画出0-250Hz范围的幅值谱

figure;

subplot(311);plot(t,y);

subplot(312);hua_baol(y,fs,1);%直接画出包络谱

subplot(313);hua_baol(y,fs,1,0,250);%直接画出包络谱，并且只画出0-250Hz范围的包络