%五点滑动平均法平滑处理.DOCVIP

%五点滑动平均法平滑处理 %%%%%%%%%%%%%%%%%%%%%% clear clc close all hidden %%%%%%%%%%%%%%%%%%%%%% fni=input(五点滑动平均法平滑处理-输入数据文件名:,s); fid=fopen(fni,r); sf = fscanf(fid,%f,1); %采样频率 m = fscanf(fid,%d,1); %平滑次数 fno = fscanf(fid,%s,1);%输出数据文件名 x = fscanf(fid,%f,inf);%输入数据存成列向量 status=fclose(fid); %取信号数据长度 n=length(x); %建立离散时间列向量 t=(0:1/sf:(n-1)/sf); %将x赋值给a a=x; %循环m次进行平滑处理计算 for k=1:m b(1)=(3*a(1)+2*a(2)+a(3)-a(4))/5; b(2)=(4*a(1)+3*a(2)+2*a(3)+a(4))/10; for j=3:n-2 b(j)=(a(j-2)+a(j-1)+a(j)+a(j+1)+a(j+2))/5; end b(n-1)=(a(n-3)+2*a(n-2)+3*a(n-1)+4*a(n))/10; b(n)=(-a(n-3)+a(n-2)+2*a(n-1)+3*a(n))/5; a=b; end %将a赋值给y y=a; %将分成2行1列的图形窗口的第1列设为当前绘图区域 subplot(2,1,1); %绘制平滑前的时程曲线图形 plot(t,x); %添加横向坐标轴的标注 xlabel(时间 (s)); %添加纵向坐标轴的标注 ylabel(加速度(g)); %在图幅上添加坐标网格 grid on; %将分成2行1列的图形窗口的第2列设为当前绘图区域 subplot(2,1,2); %绘制平滑后的时程曲线图形 plot(t,y); %添加横向坐标轴的标注 xlabel(时间 (s)); %添加纵向坐标轴的标注 ylabel(加速度(g)); %在图幅上添加坐标网格 grid on; %打开文件输出平滑后的数据 fid=fopen(fno,w); for k=1:n %每行写两个实型数据，t为时间，y为平滑后的数据 fprintf(fid,%f %f\n,t(k),y(k)); end status=fclose(fid);