使用Matlab求解一列离散数据的功率谱密度,通过实验采集到了随时间变化的加速度数据,其中时间存在txt文档的第一列,加速度数据存在第三列,采样率为5000Hz,采样时间为100秒,下面是我的一个Matl
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/06 00:25:46
![使用Matlab求解一列离散数据的功率谱密度,通过实验采集到了随时间变化的加速度数据,其中时间存在txt文档的第一列,加速度数据存在第三列,采样率为5000Hz,采样时间为100秒,下面是我的一个Matl](/uploads/image/z/1488479-23-9.jpg?t=%E4%BD%BF%E7%94%A8Matlab%E6%B1%82%E8%A7%A3%E4%B8%80%E5%88%97%E7%A6%BB%E6%95%A3%E6%95%B0%E6%8D%AE%E7%9A%84%E5%8A%9F%E7%8E%87%E8%B0%B1%E5%AF%86%E5%BA%A6%2C%E9%80%9A%E8%BF%87%E5%AE%9E%E9%AA%8C%E9%87%87%E9%9B%86%E5%88%B0%E4%BA%86%E9%9A%8F%E6%97%B6%E9%97%B4%E5%8F%98%E5%8C%96%E7%9A%84%E5%8A%A0%E9%80%9F%E5%BA%A6%E6%95%B0%E6%8D%AE%2C%E5%85%B6%E4%B8%AD%E6%97%B6%E9%97%B4%E5%AD%98%E5%9C%A8txt%E6%96%87%E6%A1%A3%E7%9A%84%E7%AC%AC%E4%B8%80%E5%88%97%2C%E5%8A%A0%E9%80%9F%E5%BA%A6%E6%95%B0%E6%8D%AE%E5%AD%98%E5%9C%A8%E7%AC%AC%E4%B8%89%E5%88%97%2C%E9%87%87%E6%A0%B7%E7%8E%87%E4%B8%BA5000Hz%2C%E9%87%87%E6%A0%B7%E6%97%B6%E9%97%B4%E4%B8%BA100%E7%A7%92%2C%E4%B8%8B%E9%9D%A2%E6%98%AF%E6%88%91%E7%9A%84%E4%B8%80%E4%B8%AAMatl)
使用Matlab求解一列离散数据的功率谱密度,通过实验采集到了随时间变化的加速度数据,其中时间存在txt文档的第一列,加速度数据存在第三列,采样率为5000Hz,采样时间为100秒,下面是我的一个Matl
使用Matlab求解一列离散数据的功率谱密度,
通过实验采集到了随时间变化的加速度数据,其中时间存在txt文档的第一列,加速度数据存在第三列,采样率为5000Hz,采样时间为100秒,下面是我的一个Matlab程序,首先画出了原始时域信号,然后画出傅里叶变化后幅值的模,然后画出功率谱密度随频率的变化,其应该是横轴为频率,纵轴为功率谱密度值.可是我通过程序得出的图和自带软件画出的图差别挺大,
M =load('data.txt');
T=M(:,1);
A=M(:,3);
figure(1);
plot(T,A);
y=fft(A,50000);
m=abs(y);
f1=(0:length(y)/2-1)'*5000/length(y);%计算变换后不同点对应的幅值
figure(2);
plot(f1,m(1:length(y)/2));
p=y.*conj(y)/50000; %计算功率谱密度
ff=5000*(0:length(y)/2-1)/50000; %计算变换后不同点对应的频率值
figure(3);
plot(ff,p(1:length(y)/2));
使用Matlab求解一列离散数据的功率谱密度,通过实验采集到了随时间变化的加速度数据,其中时间存在txt文档的第一列,加速度数据存在第三列,采样率为5000Hz,采样时间为100秒,下面是我的一个Matl
完全没有问题
y是一个复数,
m=abs(y);
p=y.*conj(y)/50000;
这两个都能把y转化成实数,但p的每一个值都是m对应值的平方.一般p用的多点