王济matlab在振动信号处理中的应用代码.docVIP

程序4-1 %最小二乘法消除多项式趋势项 %%%%%%%%%%%%%%%%%%%%%%%% 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); %计算趋势项的多项式待定系数向量a a=polyfit(t,x,m); %用x减去多项式系