sg滤波 matlab,Matlab，SG滤波器

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

function Y=fpeaks                                                          %  fpeaks函数，生成高斯峰信号，返回一个总的模拟信号向量

Peaks=load('peaks.txt');                                                   %  载入高斯峰文件，每个高斯峰占一行三列，按照峰高、峰中心位置、半峰宽写入

[row,col]=size(Peaks);                                                     %  读取矩阵大小参数

X=0:1:255;                                                                 %  生成位置自变量x

YM=zeros(row,256);                                                         %  将在YM中形成各峰模拟信号，YM初始化为0

for t=1:1:row                                                              %  每行载入一个对应于peaks.txt行的峰信号，t是循环变量

YM(t,:)=Peaks(t,1)*exp((-(X-Peaks(t,2)).^2)/((Peaks(t,3)^2)/(2*log(2))));

end;

Y=zeros(size(X));                                                          %  初始化返回向量Y为0

for t=1:1:row;                                                             %  各行相加合成总信号函数，t为循环变量

Y=Y+YM(t,:);

end;

YN=rand(size(X))-0.5;                                                      %  用伪随机数生成噪声信号

Y=Y+YN;

plot(Y)