急寻happy，小波系数画图问题！！！

zhyuer

<v:shapetype> <v:stroke joinstyle="miter"></v:stroke><v:formulas><v:f eqn="if lineDrawn pixelLineWidth 0"></v:f><v:f eqn="sum @0 1 0"></v:f><v:f eqn="sum 0 0 @1"></v:f><v:f eqn="prod @2 1 2"></v:f><v:f eqn="prod @3 21600 pixelWidth"></v:f><v:f eqn="prod @3 21600 pixelHeight"></v:f><v:f eqn="sum @0 0 1"></v:f><v:f eqn="prod @6 1 2"></v:f><v:f eqn="prod @7 21600 pixelWidth"></v:f><v:f eqn="sum @8 21600 0"></v:f><v:f eqn="prod @7 21600 pixelHeight"></v:f><v:f eqn="sum @10 21600 0"></v:f></v:formulas><v:path connecttype="rect" gradientshapeok="t" extrusionok="f"></v:path><lock aspectratio="t" v:ext="edit"></lock></v:shapetype><v:shape><v:imagedata></v:imagedata></v:shape><BR>怎么得到这个时频等值线，信号是，在时间0到10秒频率为1Hz,在10到20秒频率为0.5Hz的正弦信号。

zhyuer

<v:shapetype> <v:shapetype> 图怎么上不去呀<BR><v:stroke joinstyle="miter"></v:stroke><v:formulas><v:f eqn="if lineDrawn pixelLineWidth 0"></v:f><v:f eqn="sum @0 1 0"></v:f><v:f eqn="sum 0 0 @1"></v:f><v:f eqn="prod @2 1 2"></v:f><v:f eqn="prod @3 21600 pixelWidth"></v:f><v:f eqn="prod @3 21600 pixelHeight"></v:f><v:f eqn="sum @0 0 1"></v:f><v:f eqn="prod @6 1 2"></v:f><v:f eqn="prod @7 21600 pixelWidth"></v:f><v:f eqn="sum @8 21600 0"></v:f><v:f eqn="prod @7 21600 pixelHeight"></v:f><v:f eqn="sum @10 21600 0"></v:f></v:formulas><v:path connecttype="rect" gradientshapeok="t" extrusionok="f"></v:path><lock aspectratio="t" v:ext="edit">我的程序是close all<BR>fs=100;<BR>t1=1/fs:1/fs:10;<BR>s1=sin(2*pi*t1);<BR>t2=10+1/fs:1/fs:20;<BR>s2=sin(2*pi*0.5*t2);<BR>n=2*length(t1);<BR>s=zeros(1,n);<BR>s(1:n/2)=s1;<BR>s(n/2+1:n)=s2;<BR>t=zeros(1,n);<BR>t(1:n/2)=t1;<BR>t(n/2+1:n)=t2;<BR>figure(1),plot(t,s)<BR>m=abs(fft(s));<BR>f=(0:length(s)-1)*fs/length(s);<BR>figure(2),plot(f,m)<BR>figure,specgram(s,[0:0.1:3],fs,512,200)<BR>最后的函数可能不是用specgram</lock></v:shapetype></v:shapetype>

zhyuer

《波在瞬态空化噪声分析中的应用》图4
帮我看一下，谢谢

aspen

happy最近在写文章，所以可能来得少点

zhyuer

你能帮我解决么？谢谢

happy

对时频等值线不怎么了解，只知道可以通过s变换，根据s矩阵来绘制<BR>等值线可以用contour命令来画<BR>你还是先研究研究s变换吧

zhyuer

谢谢！！！