EDA365电子论坛网

标题: MATLAB源程序代码分享：MATLAB实现三次样条插值法求信号的包络线 [打印本页]

---

作者: House    时间: 2020-3-12 10:43
标题: MATLAB源程序代码分享：MATLAB实现三次样条插值法求信号的包络线
MATLAB源程序代码分享：MATLAB实现三次样条插值法求信号的包络线
* y1 t) }  [- x; z2 J
4 m$ \0 f  Z# \+ J% Q
! n5 \3 f" @, }+ }  x' [/ j) J  z

%% 定义一个信号

clear;clc;close all

t=0:0.001:1.5;

y=sin(11*pi*2*t)+sin(7*pi*4*t);

%% 对信号求导, 查找信号里的极小值点和极大值点

d=diff(y);       % 对信号求导

n=length(d);

d1=d(1:n-1);

d2=d(2:n);

indmin = find(d1.*d2<0 & d1<0)+1;     % 找到极小值点的索引位置

indmax = find(d1.*d2<0 & d1>0)+1;     % 找到极大值点的索引位置

; l2 E& o0 @4 f

%% 使用三次样条插值的方法, 求信号的包络线

envmin=spline([t(1) t(indmin) t(end)],[y(1) y(indmin) y(end)],t);  % 所有极小值的点 + 信号的两个端点, 进行三次样条插值

envmax=spline([t(1) t(indmax) t(end)],[y(1) y(indmax) y(end)],t);  % 所有极大值的点 + 信号的两个端点, 进行三次样条插值

%% 绘制信号的包络线

figure

set(gcf,'units','normalized','position',[0.2 0.2 0.6 0.6]);  % 设置 figure 窗口的位置和尺寸

hold on

plot(t,y)               % 绘制原信号

plot(t,envmin,'r--');   % 绘制下包络线

plot(t,envmax,'m--');   % 绘制上包络线

xlabel('t')

ylabel('y')

legend({'原始信号','下包络线','上包络线'})

0 r- R/ \9 _' G) I( {! `  N1 H

---

作者: gaoxings    时间: 2020-3-12 18:16
MATLAB实现三次样条插值法求信号的包络线。

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2