qijunshuai 发表于 2007-7-16 11:31

回复 #28 miaomiao 的帖子

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

hector1982 发表于 2007-7-18 10:06

回复 #3 qijunshuai 的帖子

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

  能把图传上来分享一下吗?你说的瀑布图是三维谱阵吗?
  我前一段用waterfall编了一个三维谱阵图,发上来大家分享一下
  % waterfall program
  clc;
  clf;
  clear;
  close all;
  =uigetfile('*.txt','please input you want to get file','file name');
  hectorload=;
  x=load(hectorload);
  =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);


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

qijunshuai 发表于 2007-8-18 13:40

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

=uigetfile('*. ...

qijunshuai 发表于 2007-8-18 13:41

回复 #33 qijunshuai 的帖子

这就是我画的瀑布图

bluelonger 发表于 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怎么搞啊?

osbertbovey 发表于 2015-4-17 22:21

qijunshuai 说把fai可以用angle替换,如果那样P,R怎么搞啊?
页: 1 2 [3]
查看完整版本: 求助全息谱,瀑布图