xcqsx.m

来自「zheshiyige消除信号趋势项的程序」· M 代码 · 共 34 行

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