程序41
最小二乘法消除多项式趋势项

clear
清除内存中所有变量和函数
clc
清除工作窗口中所显示的内容
closeallhidde

关闭所有隐藏的窗口

提示用键盘输入输入数据文件名
f
ii
put消除多项式趋势项输入数据文件名s
以只读方式打开数据文件
fidfope
f
ir
sffsca
ffidf1读入采样频率值
mfsca
ffidd1读入拟合多项式阶数
f
ofsca
ffids1读入输出数据文件名
xfsca
ffidfi
f读入时程数据存成列向量
关闭数据文件
statusfclosefid
取信号数据长度

le
gthx
建立离散时间列向量
t01sf
1sf
计算趋势项的多项式待定系数向量a
apolyfittxm
用x减去多项式系数a生成的趋势项
yxpolyvalat
将分成2行1列的图形窗口的第1列设为当前绘图区域
subplot211
绘制x对于t的时程曲线图形
plottx
在图幅上添加坐标网格
grido

将分成2行1列的图形窗口的第2列设为当前绘图区域
subplot212
绘制y对于t的时程曲线图形
plotty
在图幅上添加坐标网格
grido

以写的方式打开文件或建立一个新文件
fidfope
f
ow
进行
次循环将计算结果写到输出数据文件中
fork1

每行输出两个实型数据，t为时间，y为消除趋势项后的结果
fpri
tffidff
tkyk
循环体结束语句
fe
d关闭数据文件statusfclosefid
程序42五点滑动平均法平滑处理clearclccloseallhidde
f
ii
put五点滑动平均法平滑处理输入数据文件名sfidfope
f
irsffsca
ffidf1采样频率mfsca
ffidd1平滑次数f
ofsca
ffids1输出数据文件名xfsca
ffidfi
f输入数据存成列向量statusfclosefid取信号数据长度
le
gthx建立离散时间列向量t01sf
1sf将x赋值给aax循环m次进行平滑处理计算fork1m
b13a12a2a3a45b24a13a22a3a410
forj3
2bjaj2aj1ajaj1aj25
e
db
1a
32a
23a
14a
10b
a
3a
22a
13a
5abe
d将a赋值给yya将分成2行1列的图形窗口的第1列设为当前绘图区域subplot211绘制平滑前的时程曲线图形plottx
f添加横向坐标轴的标注xlabel时间s添加纵向坐标轴的标注ylabel加速度g在图幅上添加坐标网格grido
将分成2行1列的图形窗口的第2列设为当前绘图区域subplot212绘制平滑后的时程曲线图形plotty添加横向坐标轴的标注xlabel时间s添加纵向坐标轴的标注ylabel加速度g在图幅上r