zheshiyige消除信号趋势项的程序

				clear;
				clc;
				close all hidden;
				fni=input('消除多项式趋势项-输入文件名：','s')
				fidin=fopen(fni,'r')
				fs=100
				%实际振动中常采用m=1~3来去趋势项
				m=2
				x=fscanf(fidin,'%f',inf)
				status=fclose(fidin)
				%取数据长度
				n=length(x)
				t=(0:1/fs:(n-1)/fs)'
				%求趋势项多项式待定系数向量
				a=polyfit(t,x,m)
				y=x-polyval(a,t)
				subplot(3,1,1)
				plot(t,x)
				xlabel('t/s')
				ylabel('转矩/N.m')
				title('消趋势项前')
				grid on
				subplot(3,1,2)
				plot(t,y)
				xlabel('t/s')
				ylabel('转矩/N.m')
				title('消趋势项后')
				grid on
				z=polyval(a,t)
				subplot(3,1,3)
				plot(t,z)
				title('趋势项')
				grid on