利用FFT计算非平稳随机信号WVD分布
clc
clear all
NN=100;
w1=0.001;
a1=0.0005;
w2=0.001;
a2=-0.005;
j=sqrt(-1);
nn=-3*NN-1:1:3*NN-1; %构造信号
%n=nn-3*NN;
x=exp(j*(2*pi*(w1*nn+1/2*a1*nn.^2))); %s1=ss1(nn-3*NN)
%ss1c=conj(ss1);
L=64; %加矩形窗
N=2*L;
w=zeros(N);
w(1)=0;
l=-(L-1):1:(L-1);
w(2:length(l))=1; %2L-1窗
%s1=ss1(3*NN-(L-1)+1:3*NN+(L-1)+1);
Wx=zeros(N,N);
G=zeros(N,N);
GG=zeros(N,N); %conj(G)
f=zeros(N,N);
for n=0:N-1 % n=nn-3*NN
% l=-(L-1):1:(L-1)
% G=w.*x(3*NN+n+l) %G(n,l)=w(l).*x(n+1)
%for
l=0:1:L-1;
% ll=L+l;
% nn=3*NN+n;
G(n+1,l+1)=w(L+1+l).*x(3*NN+n+L+1+l);
GG(n+1,l+1)=conj(w(L+1-l).*x(3*NN+n+L+1-l));
f(n+1,l+1)=G(n+1,l+1).*GG(n+1,l+1);
l=L:1:2*L-1;
G(n+1,l+1)=w(L+1+l-2*L).*x(3*NN+n+L+1+l-2*L);
GG(n+1,l+1)=conj(w(L+1-l+2*L).*x(3*NN+n+L+1-l+2*L));
f(n+1,l+1)=G(n+1,l+1).*GG(n+1,l+1);
end
F=f'; %FFT按列计算矩阵的傅立叶变换
Wx=fft(F);
Wx=Wx';
n=0:N-1;
k=0:N-1;
figure;
grid;
mesh(n,k,Wx);
figure;
contour(Wx);
D156
利用FFT计算非平稳随机信号WVD分布相关推荐
- 利用FFT计算非平稳随机信号的WVD分布
up目录 一.理论基础 二.核心程序 三.测试结果 一.理论基础 fft: 快速傅里叶变换 (fast Fourier transform), 即利用计算机计算离散傅里叶变换(DFT)的高效.快速计算 ...
- 平稳与非平稳随机信号的理解
最近在做stft的相关算法,发现stft的前提条件是在窗函数内信号必须被认为是平稳随机信号,随机信号的平稳和非平稳有什么区别么?总结如下: 1.首先理解随机信号概念,随机信号理论上讲是不能利用公式复现 ...
- matlab平稳随机信号,平稳随机信号处理.ppt
平稳随机信号处理 图c是用STFT求出的x(n)的联合时频分布后,再求幅平方得到的谱图.该图是三维图形的二维投影,一个轴是时间,一个轴是频率.由该图可以清楚地看到x(n)的时间与频率的关系. (c) ...
- 【20220629】【信号处理】(平稳随机信号)自相关函数性质的证明过程
目录 1. 偶函数 2. tau=0 处取得最大值 3. 周期函数的自相关函数也是周期函数,且周期和原函数相同 4. 对于非周期信号,当 tau 趋于无穷大时,自相关函数趋于信号平均值的平方 5. 两 ...
- 现代信号处理——平稳随机信号
一.平稳随机信号的定义 狭义(严)平稳随机序列:随机信号的统计特性(均值,方差,自相关,自协方差)不随时间平移而变化,这样的随机信号在实际生活中很难找到. 广义(宽)平稳随机序列:随机信号的均值和方差 ...
- 现代信号处理——平稳随机信号的功率谱密度
一.相关函数的z变换 平稳随机序列是非周期函数,且是能量无限信号,无法直接利用傅里叶变换进行分析. 由前面对自相关函数和自协方差函数的讨论可知: 随机信号的自相关和自协方差收敛了以后,便满足了能够去做 ...
- 现代信号处理——平稳随机信号通过线性系统
一.输出信号的数字特征 确定性信号通过线性时不变系统 如果平稳随机信号x(n)通过上述线性时不变系统的输出是y(n),问: y(n)是否为平稳随机信号? 如果x(n)的统计特性已知,如何确定y(n)的 ...
- 利用 FFT 运算恢复原来的信号,PSD和原来信号之间的关系
明:利用FFT运算实现信号的重构 一. 信号建模 % Use Fourier transforms to find the frequency components of a signal burie ...
- matlab对非平稳一维信号的小波变换
matlab实验代码如下 clc cleart=linspace(0,2*pi,1024);%信号采样时间点 signal=sin(10*t).*stepfun(t,pi);%实验信号%进行一维离散小 ...
最新文章
- three20 阅读到的代码
- golang 代码格式化工具 gofmt
- 任务与中断共享资源冲突示例
- asynchronous-logging-with-log4j-2--转
- 如何跟程序员谈一场没有Bug的恋爱
- 对Python参数类型详解以及学习中遇到的坑
- java中换行符号_Java中的回车换行符/n /r /t分析
- fastcgi java_FastCGI
- 开源 免费 java CMS - FreeCMS1.3-数据对象-mail
- 工程师如何“神还原”用户问题?闲鱼回放技术揭秘
- 积分上下限无穷_数学分析|第九章 定积分利用等价无穷小量和定积分定义解决数列极限问题总结...
- 我的HEVC码流分析工具MFC小笔记:树形控件使用及窗口缩放
- uva 1291 - Dance Dance Revolution ( dp )
- centos服务端ftp的搭建_centos下搭建ftp服务器
- Thinking in Java 4th(Java编程思想第四版)文档、源码、习题答案
- Spring源码学习笔记:经典设计模式之代理模式
- R语言maps包绘制世界地图并存为矢量图 超基础!
- poi 3.17合并单元格报错 java.lang.IllegalArgumentException: Merged region A1 must contain 2 or more cells
- API开发手册在线中文版
- Echarts图列legend动态设置选中或者不选中 selected