求助：如何用小波分析实现这个图形？

kingsword1

各位高手，我现在有一组采集数据，想用小波分析来处理。我看了别人分析出来的图形是这样的（如附图），我不知道用什么命令来实现？我用了一个c=cwt(s,[0 5 10：5:25],'morl','3Dlvl');的命令，结果不对，望请赐教！谢谢！

kingsword1

是我说的不够详细吗？没有人能回答我的问题？

hyl2323

你画出来的是什么样子的呢？

hyl2323

改成
c=cwt(s,[0 5 10：5:25],'morl','3Dplot');
或者
c=cwt(s,[0 5 10：5:25],'morl','3Dabsglb');
画出来的就是三维的图，至于颜色方案，cwt默认的是pink，上图中用的是jet，你加一句
colormap(jet);
就可以了。

hyl2323

c=cwt(x,2:50,'morl','3Dplot');
colormap(jet);
x是我自己定义的一个信号。

hyl2323

c=cwt(x,2:50,'morl','3Dplot');
colormap(pink);


好像jet更好看，和ansys里面的差不多

kingsword1

我试试看呀，不过在这先谢谢了。有什么疑问，希望再给我指出！谢谢

kingsword1

哦，对了，我这个是时间-频率分析图，我怎么样才能将我的信号转换为频率到这里面去呢？我现在现有的就只有传感器采集来的数据（加速度）。能不能告诉我？

kingsword1

我直接将我采集的数据作为x带入到cwt（x，2:50,'morl','3Dplot'），计算后出来的图形是这个（附图1），与正确的图形（附图2）有天壤之别，可能是我的频率还没有转换进去，但不知道如何转入呢？

hyl2323

采样频率除以尺度才是实际的以Hz为单位的频率轴，你的实际信号是低频，出现上面的不相同的原因估计是尺度的范围选择不对，将尺度的下限选大一些（>2），这样频率轴的最高频率就会减小，频率范围集中在你的信号频率范围段（25Hz以内），那样就对了。你再试试

kingsword1

我的程序如下，得出的图形不对呀！
x=Data(1500*500:1505*500);
c=cwt(x,1:25,'morl','3Dplot');
colormap(pink);
还有就是cwt里面的频率参数为什么不能从0开始呀？

kingsword1

我照着本论坛上面的一个程序自己编了一下，程序如下：
figure(1);
t=5;fs=500;
Ts=1/fs;N=t/Ts;
x=Data(1500*500:1505*500);
% x=awgn(x,-3,'measured');
subplot(221)
plot(0:Ts:t,real(x));xlabel('time-sec');ylabel('cable acceleration-m/s^2');axis([0,5,-1,1]);
n=0:N-1;
Xk=fft(x);
mXk=abs(Xk(1:N/2));
k=(0:N/2-1*fs/N);
subplot(212),plot(k,mXk);
xlabel('Frequency in Hz units');title('Magnitude of Data');axis([0,100,0,100]);
figure(2);
Nw=500;
L=Nw/2;
Tn=(N-Nw)/L+1;
nfft=25;
TF=zeros(Tn,nfft);
for i=1:Tn
    xw=x((i-1)*10+1:i*10+10);
    temp=fft(xw,nfft);
    temp=fftshift(temp);
    TF(i,:)=temp;
end
fnew=((1:nfft)-nfft/2)*fs/nfft;
tnew=(1:Tn)*L*Ts;
[F,T]=meshgrid(fnew,tnew);
mesh(F,T,abs(TF));
其中的Nw的含义是什么呀？画出的三维时频图的频率是对称的，我如果要单边的如何改动呢？
我画出的图1是原信号的时程和功率谱，其中原信号中的横坐标为时间，幅值只有0.5～1左右，而如果用采样次数直接画出，幅值有4左右，这为什么会不同，有什么区别吗？功率谱也是，我用自己编的程序算出这段时间的功率谱在4Hz左右才有一个高峰值，其他频率处都很小，这里确出现了几个峰值，为什么？
图2是时频三维图，图像与我贴出的样板相差也很远，不知道如何改动这个程序才能达到我想要的样本那样的图形？

hyl2323

由尺度到频率的转换关系应该是采样频率/尺度*小波中心频率。我上面的帖子写错了，特此更正。向学术精神道歉！！！！！！！！！！！！！！！！！！！11

hyl2323

cwt里的是尺度参数，不是频率参数，两者成倒数关系。
自己多研究研究。

xjxu

领教，多谢！