如何使用matlab对一段离散数据进行频谱分析.docx
《如何使用matlab对一段离散数据进行频谱分析.docx》由会员分享,可在线阅读,更多相关《如何使用matlab对一段离散数据进行频谱分析.docx(14页珍藏版)》请在冰豆网上搜索。
如何使用matlab对一段离散数据进行频谱分析
相信大家做实验或者做工程的时候,都会遇到这样的情况,导出的数据是一组离散数据,那如何对这组离散数据进行频谱分析呢,下面我用MATLAB对一组离散数据做一个频谱分析案例。
0.560
0.5604-0.04
0.56080.02
0.56120.02
0.5616-0.02
0.5620.02
0.5624-0.04
0.56280.06
0.5632-0.02
0.56360.08
0.5640
0.5644-0.02
0.5648-0.06
0.56520.04
0.5656-0.02
0.566-0.02
0.5664-0.02
0.5668-0.04
0.5672-0.04
0.56760
0.568-0.02
0.56840
0.5688-0.02
0.56920
0.5696-0.04
0.57-0.06
0.57040.06
0.5708-0.06
0.5712-0.02
0.5716-0.04
0.572-0.08
0.5724-0.02
0.5728-0.1
0.5732-0.06
0.5736-0.08
0.574-0.02
0.5744-0.08
0.5748-0.02
0.5752-0.1
0.5756-0.08
0.576-0.08
0.5764-0.06
0.5768-0.08
0.5772-0.08
0.5776-0.06
0.578-0.08
0.5784-0.1
0.5788-0.08
0.5792-0.14
0.5796-0.1
0.58-0.12
0.5804-0.16
0.5808-0.12
0.5812-0.14
0.5816-0.16
0.582-0.1
0.5824-0.18
0.5828-0.12
0.5832-0.18
0.5836-0.12
0.584-0.18
0.5844-0.18
0.5848-0.28
0.5852-0.16
0.5856-0.22
0.586-0.2
0.5864-0.22
0.5868-0.28
0.5872-0.28
0.5876-0.28
0.588-0.3
0.5884-0.3
0.5888-0.3
0.5892-0.3
0.5896-0.34
0.59-0.38
0.5904-0.28
0.5908-0.42
0.5912-0.38
0.5916-0.42
0.592-0.46
0.5924-0.42
0.5928-0.56
0.5932-0.46
0.5936-0.56
0.594-0.52
0.5944-0.56
0.5948-0.5
0.5952-0.62
0.5956-0.58
0.596-0.64
0.5964-0.62
0.5968-0.64
0.5972-0.66
0.5976-0.7
0.598-0.72
0.5984-0.74
0.5988-0.72
0.5992-0.84
0.5996-0.8
0.6-0.82
0.6004-0.88
0.6008-0.86
0.6012-0.9
0.6016-0.94
0.602-0.9
0.6024-0.98
0.6028-0.94
0.6032-1
0.6036-0.94
0.604-1.06
0.6044-1.08
0.6048-1.14
0.6052-1.02
0.6056-1.1
0.606-1.1
0.6064-1.12
0.6068-1.14
0.6072-1.16
0.6076-1.12
0.608-1.16
0.6084-1.16
0.6088-1.18
0.6092-1.14
0.6096-1.14
0.61-1.16
0.6104-1.06
0.6108-1.16
0.6112-1.1
0.6116-1.06
0.612-1.1
0.6124-1
0.6128-1.08
0.6132-0.96
0.6136-1
0.614-0.9
0.6144-0.9
0.6148-0.8
0.6152-0.86
0.6156-0.74
0.616-0.7
0.6164-0.66
0.6168-0.62
0.6172-0.54
0.6176-0.52
0.618-0.42
0.6184-0.36
0.6188-0.3
0.6192-0.34
0.6196-0.16
0.62-0.16
0.6204-0.12
0.62080
0.62120.06
0.62160.08
0.6220.2
0.62240.2
0.62280.34
0.62320.34
0.62360.44
0.6240.42
0.62440.5
0.62480.5
0.62520.64
0.62560.64
0.6260.7
0.62640.7
0.62680.76
0.62720.78
0.62760.84
0.6280.86
0.62840.92
0.62880.92
0.62920.96
0.62960.98
0.630.96
0.63041.08
0.63080.98
0.63121.06
0.63161.06
0.6321.02
0.63241.04
0.63281
0.63321.08
0.63361.02
0.6341.1
0.63441.02
0.63481.1
0.63520.98
0.63561.04
0.6361.02
0.63641.02
0.63681
0.63721.02
0.63760.98
0.6380.96
0.63840.92
0.63880.96
0.63920.86
0.63960.94
0.640.86
0.64040.82
0.64080.88
0.64120.84
0.64160.78
0.6420.84
0.64240.72
0.64280.8
0.64320.72
0.64360.76
0.6440.66
0.64440.68
0.64480.6
0.64520.7
0.64560.62
0.6460.64
0.64640.54
0.64680.58
0.64720.54
0.64760.52
0.6480.54
0.64840.5
0.64880.5
0.64920.5
0.64960.48
0.650.44
0.65040.46
0.65080.4
0.65120.42
0.65160.42
0.6520.34
0.65240.4
0.65280.3
0.65320.36
0.65360.28
0.6540.36
0.65440.32
0.65480.36
0.65520.24
0.65560.28
0.6560.26
0.65640.26
0.65680.26
0.65720.26
0.65760.26
0.6580.24
0.65840.24
0.65880.22
0.65920.18
0.65960.26
0.660.2
0.66040.2
0.66080.22
0.66120.18
0.66160.18
0.6620.2
0.66240.16
0.66280.22
0.66320.16
0.66360.22
0.6640.1
0.66440.18
0.66480.1
0.66520.22
0.66560.12
0.6660.14
0.66640.1
0.66680.14
0.66720.1
0.66760.12
0.6680.12
0.66840.1
0.66880.08
0.66920.12
0.66960.1
0.670.08
0.67040.14
0.67080.06
0.67120.08
0.67160.08
0.6720.02
0.67240.1
0.67280
0.67320.08
0.6736-0.02
0.6740.1
0.67440.04
0.67480.1
0.67520
0.67560.08
0.6760.02
0.67640.06
0.67680.06
0.67720.08
0.67760.02
0.6780.06
0.67840.04
0.67880.04
0.67920.02
0.67960.08
0.680.04
0.68040
0.68080.06
0.68120.02
0.68160.02
0.6820.06
0.68240
0.68280.1
0.68320.04
0.68360.06
0.6840.02
0.68440.06
0.6848-0.02
0.68520.08
0.68560
0.6860.02
0.68640
0.68680.06
0.6872-0.04
0.6876-0.02
0.6880.04
0.68840
0.68880
0.68920.02
0.68960.02
0.690
0.69040.04
0.6908-0.04
0.6912-0.02
0.69160
0.692-0.02
0.69240.06
0.6928-0.04
0.69320
0.6936-0.04
0.6940.06
0.6944-0.04
0.69480.06
0.6952-0.08
0.6956-0.02
0.696-0.02
0.69640
0.69680
0.69720
0.6976-0.02
0.6980.02
0.69840
0.6988-0.02
0.6992-0.02
0.69960.02
0.70.02
以上一组离散数据第一列是时间,第二列对应的是该时间点线圈传感器检测磁场变化所输出的电压值。
现在我写一段matlab程序对该数据进行分析。
fs=2500;//因为以上数据是每隔0.0004秒采集一次数据,所以采样率应该是2500
n=351;//因为以上离散数据一共有351组数据,所以n取值351
t=0.56:
0.0004:
0.7;//以上离散数据第一列时间点数据是从0.56秒开始,每隔0.0004秒采集一次,0.7秒结束
baxxxxse=importdata('C:
\Users\Administrator\Desktop\AGM\111.txt');//这句作用是导入存储该离散数据的TXT文件
y1=fft(baxxxxse(:
2));
y2=fftshift(y1);
f=(0:
n-1)*fs/n-fs/2;
figure
(1)
plot(tbaxxxxse(:
2)'r');//打印出离散数据原始图
figure
(2)
plot(fabs(y2)'b');//打印出离散数据的频谱分析线
axis([015002000]);//横坐标、纵坐标的显示范围,可根据实际情况而修改。
以下两个图形是对该组离散数据的频谱分析结果。