Matlab寻峰程序.doc
需积分: 38 125 浏览量
2020-04-20
16:59:14
上传
评论 2
收藏 37KB DOC 举报
function [output] = autopeak(x,y)
% INITIALIZATION 初始化
% set detection sensitivity parameters 设置灵敏度参数为黄金分割1.618
% the golden ratio seems to work exceptionally fine :)
nfwhms = 1.6180;%1.6180; 半高全宽
nstds = 1.6180; %1.6180; 标准偏差
% put signal to columns if not so yet 检查输入信号
shape = size(x);
if shape(2) > 1
x = x';
end
shape = size(y);
if shape(2) > 1
y = y';
end
% smooth signal, determine its baseline and detection threshold 平滑信号
y = smooth(y);
threshold = nstds * std(y); % 设定阈值
% plot the signal % 画图
plot(x,y);
hold on;
% PEAK DETECTION RUNDOWN 开始寻峰
% hunt the peaks down one by one and measure them 依次寻峰并测量，之后从信号
中删除，知道小于阈值
% then eliminate them from signal and repeat until threshold
for i = 1 : length(x)
% determine peak by global maximum 找出最大值，并计算基准1（平均值）。
[peaki, ipeak] = max(y);
base(1) = mean(y);
% check if big enough, else stop
if peaki < threshold
break
end
% save the peak and plot it 保存峰值并画出
peak(i) = peaki;
loc(i) = x(ipeak);
plot(loc(i),peak(i),'ro');
% determine fwhm
% by shifting signal twice and locating minimas on both sides
% do it in two steps: using global baseline (j=1) and fwhm-based
local baseline (j=2)
资源评论
qq_42721399
- 粉丝: 0
- 资源: 1
最新资源
- 计算机工程与应用期刊写作模块（最新）
- jsp 查看页显示图片(使用a标签的形式)
- n_objectracker.m
- gin-demo我自己写的gindemo
- 2023-加速度-DevOps-状态报告-中国DevOps社区版
- Framework-CoreKit-2023.12.07.unitypackage
- Rsync+Sersync
- rustdesk-1.2.3-aarch64-signed.apk.1
- ee240课程，基于斯坦福A Basic Introduction to the gm ID-Based Design
- R214-0762（III）SG-2-D-20231130连廊补三四层变更_t3.dwg
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
