matlab对频谱傅里叶逆变换,基于功率谱的傅里叶逆变换问题
各位大神好,
我最近在做路面不平度。就是用傅里叶逆变换的方法利用该方法将公路路面功率谱密度进行合理的离散化得到功率谱密度的数据,然后对这些数据进行计算得到路面不平度的离散傅立叶变换,对离散傅立叶变换的数据按照一定规则补齐后再进行傅立叶逆变换,便得到路面不平度数据,其功率谱密度与所给定的应该是准确一致。但是我逆变换出来的功率谱与原来的功率谱之间总是差5倍,求解答??
程序如下:
clear
clc
v=10;
fmin=0.5;
fmax=30;
nmin=0.01; %fmin/v;
nmax=3; %fmax/v;
Gd0=16;
fs=10; %采样频率fs(frequency of sampling>2fmax
deltt=1/fs; %采样时间间隔
N=2^12; %采样点数
deltl=0.1; %v*deltt; %采样距离间隔
L=deltl*N; %采样距离
deltf=1/L; %采样空间频率间隔
%deltf=1/deltt/N; %采样时间频率间隔
Neff=(nmax-nmin)/deltf; %有效频段内的采样点数
N0=nmin/deltf; %需要满足条件(N0+Neff
k=0;
while(k*deltf
k=k+1; %判断下截止频率所在位置
end
l=0;
while(l*deltf
l=l+1; %判断上截止频率所在位置
end
Gd=[];
x=[];
for ii=1:k
Gd(ii)=0;
x(ii)=0;
end
for ii=k+1:1:l+1
Gd(ii)=Gd0*((ii-1)*deltf/0.1)^(-2);
x(ii)=sqrt(N*Gd(ii)/2/deltl);
%x(ii)=N*sqrt(Gd(ii)*deltf);
end
for ii=l+2:1:N/2
Gd(ii)=0;
x(ii)=0;
end
Gd(N/2+1)=0;
x(N/2+1)=0;
a=[];
xc=zeros(1,N);
for j=1:N/2
a(j)=2*pi*rand();
b(j)=cos(a(j))+i*sin(a(j));
xc(j)=b(j).*x(j);
end
xc(N/2+2:N)=fliplr(conj(xc(2:N/2))); %共轭对称
xifft=ifft(xc);
for j=1:N
t(j)=(j-1)*deltl;
end
xaver=mean(xifft)
xsd=std(xifft)
xmax=max(xifft)
subplot(2,1,1)
plot(t,xifft)
xlabel('沿道路走向长度(m)');
ylabel('路面不平度(mm)');
fx=[];
for i=1:N/2+1
fx(i)=(i-1)*deltf;
end
[power,fre]=psd(xifft,N,fs,boxcar(length(xifft)),0,0.97);
power(1:k)=0;
power(l+1:N/2)=0;
subplot(2,1,2);
loglog(fre,power,'r',fx,Gd(1:N/2+1),'b')
axis([10^(-2) 10 10^(-2) 10^6])
xlabel('空间频率n/m-1');
ylabel('功率谱密度Gd(n)*10^(-6)m2/(m-1)')
grid on
以上是程序,变换完后,我用power/Gd 二者之间总是错5倍。
计算结果:
2014-8-12 14:39 上传
matlab对频谱傅里叶逆变换,基于功率谱的傅里叶逆变换问题相关推荐
- 基于MATLAB的频谱、能量谱、三分之一倍频程分析
目录 1. 相关概念介绍 2. 代码实现 3. 场景仿真: [若觉文章质量良好且有用,请别忘了点赞收藏加关注,这将是我继续分享的动力,万分感谢!] 1. 相关概念介绍 以信号为例,信号在时域下的图形可 ...
- Matlab中频谱、功率谱、功率谱密度的关系以及如何求信号的功率
提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档 文章目录 Matlab中频谱.功率谱.功率谱密度的关系 一.频谱.功率谱.功率谱密度 二.MATLAB中求信号的功率 Matlab中频 ...
- 信号频域分析方法的理解(频谱、能量谱、功率谱、倒频谱、小波分析)
信号的频域分析方法多种多样,这里针对较为常见的(频谱.能量谱.功率谱.倒频谱.小波分析)集中进行说明.这些方法的MATLAB代码实现参见文章频域特征值提取的MATLAB代码实现(频谱.功率谱.倒频谱) ...
- 功率谱 魏凤英统计程序_频谱、能量谱、功率谱、功率谱估计
关于频谱.能量谱.功率谱 对于能量信号,常用能量谱来描述.所谓的能量谱,也称为能量谱密度. 是指用密度的概念表示信号能量在各频率点的分布情况.也即是说,对能量谱在频域上积分就可以得到信号的能量.能量谱 ...
- OPENCV计算机视觉图像处理频域傅里叶 DFT 变换低通滤波逆变换IDFT
OpenCV计算机视觉图像频域傅里叶 DFT 变换低通滤波逆变换IDFT 实验室做图像的,经常用到这部分,为了检测屏幕,看过好多博客,试用过许多代码,这个算是我找到的比较好用的,也容易改. 傅里叶变换 ...
- 功率谱有什么用_频谱、能量谱、功率谱、功率谱估计
关于频谱.能量谱.功率谱 对于能量信号,常用能量谱来描述.所谓的能量谱,也称为能量谱密度. 是指用密度的概念表示信号能量在各频率点的分布情况.也即是说,对能量谱在频域上积分就可以得到信号的能量.能量谱 ...
- 频域分析方法的理解(信号、频谱、能量谱、功率谱、倒频谱、小波分析)(终极教程)
信号频域分析方法的理解(频谱.能量谱.功率谱.倒频谱.小波分析) 参考文献,值得一读 https://blog.csdn.net/fengzhuqiaoqiu/article/details/1011 ...
- matlab 图片倒影_计算物理基于matlab方法研究水中倒影问题
计算物理基于matlab方法研究水中倒影问题 基于 Matlab 方法研究水中倒影问题[摘 要] 本文介绍了用 matlab 研究倒影问题的方法,利用 matlab 可视化的优点可 以直观的得出结果. ...
- matlab统计所有股票分析,MATLAB金融算法分析实战:基于机器学习的股票量化分析...
MATLAB金融算法分析实战:基于机器学习的股票量化分析 作者:吴婷;余胜威 编著 出版日期:2017年07月 文件大小:32.24M 支持设备: ¥50.00在线试读 适用客户端: 言商书局 iPa ...
- rls算法matlab实现,第5章基于RLS算法的数据预测与MATLAB实现MATLAB实现.PDF
第5章基于RLS算法的数据预测与MATLAB实现MATLAB实现 第 5章 基于 RLS算法的数据预测与 第5章 基于RLS算法的数据预测与MATLAB实现 MATLAB实现 RLS 1795 递归最 ...
最新文章
- 2016年度工作总结
- 甲小姐对话稚晖君:深度学习并非AI的终点
- 第十六周程序阅读(1)
- [Struts]Token 使用及原理
- 领域驱动设计-从贫血模型到充血模型
- Thinkphp js、css压缩类minify
- linux中如何复制文件并重命名_linux 下怎么复制一个文件到另外一个目录并且重命名...
- 车联网领域,传统TSP企业做错了什么 ?
- 商汤港理工提出基于聚类的联合建模时空关系的 GroupFormer 用于解决群体活动识别问题,性能SOTA...
- java实现pop客户机程序_java客户机如何通过JEE服务器(如weblogic、jboss等)进行身份验证的概念...
- 编写android驱动程序,Android 驱动编写LED-NDK程序
- c++函数可变参数的使用
- java CRC32
- 测试用例的几种常见设计方法
- 【MySQL数据库设计与应用(四)】视图
- yolov3损失函数分析
- 对抗攻击常见方法汇总
- 爬取豆瓣 TOP250 电影排行榜
- 他被称为“中国杀毒软件之父”,没想到身患残疾,38岁才接触电脑!
- 树莓派 zero 通过 max31865 连接 PT100 热电阻 测量温度
热门文章
- 从打造“兴趣电商”到“完全闭环”,抖音电商之路
- 仅在 localhost 中部署并使用 RSSHub
- itextpdf生成列表基本用法
- access 数据更新update
- 9张图看懂AI人工智能现状!从事AI职业女性稀缺!
- Oracle中TO_NUMBER()函数的用法
- 编码器/译码器(Verilog HDL)|计算机组成
- SQL 考勤打卡问题
- gyp ERR! stack Error: EACCES: permission denied,解决
- 建筑智能化集成系统及架构简介