求大神看看 这种时频三维波形要用什么函数画

flydonghao

求大神回复

猫头鹰先生

这是最简单的WVD！楼主这个图很简单的，有专门工具箱可以搞定。

粤语残片

猫头鹰先生 发表于 2014-4-1 00:07
这是最简单的WVD！楼主这个图很简单的，有专门工具箱可以搞定。

请问可以用小波做出这样的图吗

猫头鹰先生

粤语残片 发表于 2014-4-1 08:57
请问可以用小波做出这样的图吗

当然可以的，小波变换三维视频图都可以得到。

猫头鹰先生

……………………

粤语残片

猫头鹰先生 发表于 2014-4-1 10:09
……………………

咦！好漂亮的图，请问可以请教一下绘图方法吗，就是具体的步骤，谢谢~PWV不是伪魏格纳分布吗，不是小波啊，我是菜鸟，见笑了啊

kkkttt

subplot(4,4,[2,3,4])
绘制原始信号图
subplot(4,4,[5,9,13])
绘制功率谱
subplot(4,4,[6,7,8,10,11,12,14,15,16])
绘制视频图

大概应该是这样一个程序流程，三维时频图可以用surf或者mesh等命令实现

kkkttt

一个示例子，功率谱部分没有绘制，自己完成吧

1. % cwt_test
   
2. clear all
   
3. close all
   
4. clc
   
5. 
   
6. %%
   
7. Ts = 0.001;
   
8. Fs = 1/Ts;
   
9. f1 = 20;
   
10. f2 = 50;
    
11. f3 = 100;
    
12. 
    
13. dt = 0.2;
    
14. t1 = (0:Ts:dt-Ts) + 0;
    
15. t2 = (0:Ts:dt-Ts) + dt;
    
16. t3 = (0:Ts:dt-Ts) + 2*dt;
    
17. 
    
18. y1 = sin(2*pi*f1*t1);
    
19. y2 = sin(2*pi*f2*t2);
    
20. y3 = sin(2*pi*f3*t3);
    
21. 
    
22. t = [t1 t2 t3];
    
23. y = [y1 y2 y3];
    
24. 
    
25. subplot(4,4,[2,3,4])
    
26. plot(t,y)
    
27. xlim([t(1) t(end)])
    
28. ylim([min(y) max(y)])
    
29. xlabel('时间t')
    
30. ylabel('信号y(t)')
    
31. title('原始信号')
    
32. 
    
33. %%
    
34. scale = 1:50;
    
35. cw2 = cwt(y,scale,'morl');
    
36. 
    
37. %%
    
38. subplot(4,4,[5,9,13])
    
39. 
    
40. 
    
41. 
    
42. subplot(4,4,[6,7,8,10,11,12,14,15,16])
    
43. [X,Y] = meshgrid(t,5/(2*pi)./scale*Fs);
    
44. mesh(X,Y,abs(cw2))
    
45. title('时频图')
    
46. xlabel('时间')
    
47. ylabel('频率')
    
48. xlim([t(1) t(end)])
    
49. set(gca,'ylim',[0,max(max(Y))])
    
50. set(gca,'YScale','log')

复制代码

粤语残片

kkkttt 发表于 2014-4-2 10:28
一个示例子，功率谱部分没有绘制，自己完成吧

出现错误Maximum variable size allowed by the program is exceeded.，网上搜说要改算法，我外行，不是，什么都不行，怎么破，大侠，我是一维数据，怎么改成符合你的程序的矩阵啊，谢谢！

flydonghao

谢谢大神回复  确实是用工具箱画出来的  STFT、PWD、WD都有~