声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

123
返回列表 发新帖
楼主: qijunshuai

[综合] 求助全息谱,瀑布图

[复制链接]
 楼主| 发表于 2007-7-16 11:31 | 显示全部楼层

回复 #28 miaomiao 的帖子

这段程序我也看过,不过没能理解!!而且这里面也没有用带内插计算的fft方法.而且你给的这段程序好像是三维全息谱吧!!
回复 支持 反对
分享到:

使用道具 举报

发表于 2007-7-18 10:06 | 显示全部楼层

回复 #3 qijunshuai 的帖子

本帖最后由 wdhd 于 2016-7-6 10:02 编辑

  能把图传上来分享一下吗?你说的瀑布图是三维谱阵吗?
  我前一段用waterfall编了一个三维谱阵图,发上来大家分享一下
  % waterfall program
  clc;
  clf;
  clear;
  close all;
  [filename path]=uigetfile('*.txt','please input you want to get file','file name');
  hectorload=[path filename];
  x=load(hectorload);
  [m,n]=size(x);
  a=x(:,1);
  b=x(:,2);
  m;
  k=fix(m/1024);
  fs=10000;
  x2=b;
  t=1/fs:1/fs:(m)/fs;
  t=t';
  fid=fopen(hectorload,'r');
  m1=fscanf(fid,'%d',1);
  s2=x2;
  a1=polyfit(t,s2,m1);
  s2=s2-polyval(a1,t);
  b=s2;
  %%%%%%%%%%%%%%
  n1=k;
  k1=50;
  ii1=1;
  %%%%%%%%%%%%%%%%%%%
  nn=80;
  n2=fix(n1/k1);
  nn2=n1-ii1;
  for i=ii1:1:n1;
  l1=1024*(i-1);
  l2=1024*i;
  j=1;
  if l1==0;
  l1=1;
  l2=1025;
  end
  t=l1:1:l2;
  Y=fft(b(t),1024)/1024;
  Pyy=sqrt(2*Y.*conj(Y));
  f=fs*(0:511)/1024;
  z=Pyy;
  Z(:,i)=z(1:nn);
  end
  Z=Z';
  T=Z;
  T1=T;
  T1(:,1)=0;
  for i=1:n2;
  T2(i,:)=T1(k1*i,:);
  end
  n=2*n2-1;
  f=1:0.0001:nn*0.00001;
  t=1:2:n;
  figure(1)
  waterfall(f,t,T2);
untitled.jpg

[ 本帖最后由 hector1982 于 2007-7-18 10:34 编辑 ]

waterfall.doc

45 KB, 下载次数: 89

评分

1

查看全部评分

 楼主| 发表于 2007-8-18 13:40 | 显示全部楼层
原帖由 hector1982 于 2007-7-18 10:06 发表
能把图传上来分享一下吗?你说的瀑布图是三维谱阵吗?
我前一段用waterfall编了一个三维谱阵图,发上来大家分享一下
% waterfall program
clc;
clf;
clear;
close all;

[filename path]=uigetfile('*. ...

waterfall.bmp
 楼主| 发表于 2007-8-18 13:41 | 显示全部楼层

回复 #33 qijunshuai 的帖子

这就是我画的瀑布图
发表于 2007-10-30 16:36 | 显示全部楼层
原帖由 qijunshuai 于 2007-5-13 21:51 发表

z=x+y*i;    fz=fft(z);
ZReal=real(fz)/512;ZImag=imag(fz)/512;    ZR1=ZReal(1:256);ZI1=ZImag(1:256);
for i=1:256
ZR2(i)=ZReal(512-i+1);   ZI2(i)=ZImag(512-i+1);
end;
XP=sqrt(ZR1.^2+ZI1.^2);  ...

请问各位高手,我用这段程序画出的图怎么是一行小点点,faiP中的P和faiR中的R是甚么意思啊
我对matlab不是太懂, qijunshuai 说把fai可以用angle替换,如果那样P,R怎么搞啊?
发表于 2015-4-17 22:21 | 显示全部楼层
qijunshuai 说把fai可以用angle替换,如果那样P,R怎么搞啊?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-5 02:34 , Processed in 0.070332 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表