死灵不归路
死灵不归路
关注数: 21
粉丝数: 18
发帖数: 75
关注贴吧数: 17
户外旅行,来涨知识,来攒勇气 混个脸熟~~~
本人小白,哪位大神可以教教我下面这段代码,怎么带入数据呢 function[tfr1,tfr2,tfr3,t,f]=ACS(x,Q,level,M,verbose,c,d,emits); %本程序需要调用DiscreteTFDs工具箱中的函数 %x 信号 %Q 自适应线性调频小波分解迭代次数 %level 线性调频小波参数估计方式:1/2/3,默认值为2 %M 分辨率参数,默认值为64 %verbose 自适应线性调频小波分解结果显示:0/1,默认值为0 %c 初始化线性调频率,默认值为0 %d 初始化持续时间长度,默认值为50 %emits EM算法的最大迭代次数,默认值为5 %tfr1 自适应线性调频小波谱图 %tfr2 归一化最优线性调频小波的Wigner-Ville分布 %tfr3 改进的自适应线性调频小波谱图 %t 时频平面中的时间坐标【1,N】,N为信号长度 %f 时频平面中的频率坐标【0,0.5】 N=length(x);%信号长度 x1=zeros(N,1);%初始化最优线性调频小波 x2=zeros(N,1);%初始化归一化的线性调频小波 tfr1=zeros(N,N);%初始化自适应线性调频小波谱图 tfr2=zeros(N,N);%初始化归一化的最优线性调频小波的Wigner-Ville分布 tfr3=zeros(N,N);%初始化改进的自适应线性调频小波谱图 [P,res]=find_chirplets(x,Q,level,M,verbose,c,d,emits); %自适应线性调频小波分解,p为最优线性调频小波参数,res为残余信号 [row,col]=size(P); PN=[ones(row,1) P(:,2:col)];%幅值归一化的最优线性调频小波参数 for loop=1:Q x1=chirplets(N,P(loop,:));%重构最优线性调频小波 [tfrx1,t,f]=tfrwv(x1);%最优线性调频小波的Wigner-Ville分布 tfr1=tfr1+tfrx1;%自适应线性调频小波谱图 x2=chirplets(N,PN(loop,:));%重构归一化的最优线性调频小波 [tfrx2,t,f]=tfrwv(x2);%归一化最优线性调频小波的Wigner-Ville分布 tfr2=tfr2+tfrx2;%归一化最优线性调频小波的Wigner-Ville分布 end threshold=mean(mean(tfr2));%时频滤波器阈值 tffilter=zeros(N,N);%初始化时频滤波器 for i=1:N for j=1:N if tfr2(i,j)>threshold tffilter(i,j)=1.0;%时频滤波器赋值 else tffilter(i,j)=0.0;%时频滤波器赋值 end end end [tfr,t,f]=tfrwv(x);%信号的Wigner-Ville分布 tfr3=tfr.*tffilter;%改进的自适应线性调频小波谱图
据说一个人的七夕是这样子滴
喝伊利金典牛奶 酸的还有点臭 目测那牛是吃肉的
青龙出什么装很无敌? 刚玩,求指导
1
下一页