www.pudn.com > Eliminate-the-trend-method.zip > 正弦信号去趋势项.m, change:2012-05-09,size:658b


clear; 
 %fs为采样频率; 
 fs=1000; 
 %建立离散时间列向量 
 t=1/fs:1/fs:10; 
 %一个带有趋势项的正弦信号 
 x=2*sin(pi*t)-3*t; 
 %绘制x对于t的时程曲线图形 
 original=plot(t,x); 
  %横坐标的标签为时间t 
 xlabel('时间(t)') 
 %在图幅上添加坐标网格 
 grid on; 
 %当前轴及图形保持而不被刷新,准备接受此后将绘制 
 hold on; 
 %拟合多项式的阶数 
 m=2; 
 %计算多项式待定系数向量a 
 a=polyfit(t,x,m); 
 %用x减去多项式系数a生成的趋势项 
 y=x-polyval(a,t);  
 %在图幅上添加坐标网格 
 grid on; 
 %把去除趋势项后的时程曲线绘成红色 
 Leastsquaremethod=plot(t,y,':'); 
 %把纵坐标标签为位移 
 ylabel('位移(mm)') 
 %绘制曲线所用线型.色彩或数据点形图例 
 legend('输入','输出')