wigner程序
请问wigner分布函数tfrwv的源代码是不是matlab工具箱自带的?怎么才能调用?还是需要自己找源代码的m文件?望有高手指点。非常感谢!:@)wigner程序
谁能给一份tfrwv.m的源程序代码啊?还有wigner-vill的三维图函数代码。迫切需要。谢谢。 tfrwv.m是时频分析工具箱里的,需要下载,在论坛里找找吧,应该有至于三维图代码,这里有个简单的例子,你可以参考
%本程序用来演示WVD
sig=fmlin(128,0.1,0.4);%线性调频信号
tfrwv(sig)%产生菜单
=tfrwv(sig);
figure(1)%绘制contour图
contour(t,f,abs(tfr))
xlabel('time/s')
ylabel('frequency/Hz')
title(' contour')
figure(2)%绘制mesh三维图
mesh(t,f,abs(tfr))
xlabel('time/s')
ylabel('frequency/Hz')
zlabel('Ampititude')
title(' mesh')
回复 沙发 若菱 的帖子
??? Undefined command/function 'fmlin'这个是新版本的函数吗?还是自己装的呢?我的7.0找不到。
回复 板凳 sogooda 的帖子
Matlab R2008a也找不到此函数Reference Page Not Found回复 沙发 若菱 的帖子
为什么我运行这个程序有错误无法运行呢?总是提示:??? Error using ==> tfrview
unsupported matlab version. please send an email.
Error in ==> tfrqview at 176
tfrview(tfr,sig,t,method,param);
Error in ==> tfrwv at 73
tfrqview(tfr,x,t,'tfrwv');
Error in ==> Untitled4 at 3
tfrwv(sig)%产生菜单 我下载的tfrwv。m好像也不能用。我的程序是:
N=1024;
fs=1024;
dt=1/fs;
df=fs/N;
t=0:dt:(N-1)*dt;
x=cos(2*pi*30*t+0.5*sin(2*pi*15*t))+sin(2*pi*120*t);
figure(1)
plot(t,x)
figure(2)
=tfrwv(x)
plot(tfr,f)
运行就出错,系统提示:
??? Error using ==> tfrwv
X must have one or two columns
Error in ==> WVDceshi at 15
=tfrwv(x)
这是什么问题?迫切需要高人解答。 本帖最后由 牛小贱 于 2015-1-4 09:32 编辑
试试这个
N=1024;
fs=1024;
dt=1/fs;
df=fs/N;
t=0:dt:(N-1)*dt;
x=cos(2*pi*30*t+0.5*sin(2*pi*15*t))+sin(2*pi*120*t);
x=x';
figure(1)
plot(t,x)
figure(2)
=tfrwv(x)
plot(tfr,f)
增加个转置 x=x'; 有什么目的呢 ?
回复 8楼 graduate 的帖子
根据错误提示,将行矩阵转化为列矩阵 谢谢你的回答。但请问这个程序最后是不是还要加上坐标的变换才能得到正确的频率和时间坐标轴呢?
谢谢了。:@)
回复 10楼 zhuxiaoxun 的帖子
不好意思,这个就不清楚啦! 谢谢唠嗑啊。我觉得时频工具箱里很多函数都经过归一化处理,所以有些函数出来的结果不对。在其他方法中也出现了这种现象。比如hilbert变换等,这些好像最后都用到了这样一段程序:
=toimage(A,f,t,length(t));
for k=1:size(E,1)
bjp(k)=sum(E(k,:))*1/fs;
end
plot(f,bjp);
不知道是不是这个意思。等待大家指教啊。 这个工具箱需要下载吗
页:
[1]