| 查看: 3213 | 回复: 1 | |||
凝雪小小新虫 (小有名气)
|
[求助]
窗口傅里叶变换 怎么画时频图。。。
|
|
加窗傅里叶变换，加hamming窗，序列是离散的，存在数值A中可能有1000左右的数值在A中，怎么做它的加窗傅里叶变换。做傅里叶变换时不能利用fft快速傅里叶变换，因为不是二的整次幂。我写了个小程序，可是有问题，时频图也画不出来 clc;clear; fid=fopen('C:\Documents and Settings\Administrator\桌面\Matlab数据程序\必需基因\M65145.txt','r');%里面是1072个碱基序列。 b=fscanf(fid,'%s'); %l=length(a); %b=a(1:1638); N=length(b); fclose(fid); A=zeros(1,N); for i=1:1:N if b(i)=='a'; A(i)=0.1260; else if b(i)=='t'; A(i)=0.1335; else if b(i)=='g'; A(i)=0.0806; else if b(i)=='c'; A(i)=0.1340; end end end end end M=109; %Lap=M/2; Lap=90; Lshift=M-Lap; w=hamming(M); Nshift=floor((length(A)-Lap)/Lshift);% 计算窗函数移动次数 tf=zeros(M/2,Nshift);% 生成时频矩阵 for i=0:Nshift-1; xi=A((i*Lshift+1)  i*Lshift+M));% 截取数据段xw=xi.*w'; %数据相乘 n=0 M-1);k=0M-1);temp=xw*exp(-j*2*pi/M).^(n'*k);%离散傅里叶变换 tf(:,i+1)=temp(1:M/2);% 将频率为0-fs/2范围内的频率点数傅里叶值存入时频矩阵 end tf_t=(0:Nshift-1)*Lshift;% 计算时间 tf_f=(0:M/2-1)/M;% 计算频率 [T,F]=meshgrid(tf_t,tf_f); contour(T,F,abs(tf));  未命名.jpg |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有78人回复
-
filecode=后面第一个是大写字母
已经有8人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有13人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 傅里叶变换 频谱分析
已经有12人回复
凝雪小小
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 421.8
- 散金: 10
- 帖子: 105
- 在线: 61.7小时
- 虫号: 1394131
- 注册: 2011-09-07
- 性别: GG
- 专业: 控制理论与方法
2楼2012-09-27 20:36:34
