【转帖】振动信号趋势项消除的Matlab程序实现

振动信号消除趋势项matlab程序%用最小二乘法消除振动信号的趋势项
clear;
clc;
fni=input('数据文件名：','s');
fid=fopen(fni,'r');
fs=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/fs:(n-1)/fs);
a=polyfit(t,x,m); %计算多项式待定系数向量a
y=x-polyval(a,t); %用x减去多项式系数a生成的趋势项
subplot(2,1,1);
plot(t,x);
title('original signal');
grid on;
subplot(2,1,2);
plot(t,y);
title('redifined signal');
grid on;
fno=input('输出数据的文件名：','s');
save(fno,'t','y');
本贴来自： http://www.chinavib.com/forum/thread-47811-1-1.html