日日行,不怕千万里;时时学,不怕千万卷
分类:
2008-05-27 15:48:54
%使用fft检测相位 tic; clear all; clc; sample=4096; t=linspace(0,4*pi,sample); %基于fft求相位的算法 A=sin(t+pi/3); B=sin(t+3*pi/4); N=length(t); tp=5; Ws=N/tp; R=0.0001*rand(1,N); P=0.0001*rand(1,N); Z=A+R; O=B+P; a=fft(A); [F,I]=max(abs(a)); %对A和B信号做fft,并且求相位 b=angle(i*real(a(I))-imag(a(I)))*180/pi a=fft(B); [F,I]=max(abs(a)); c=angle(i*real(a(I))-imag(a(I)))*180/pi a=fft(Z); [F,I]=max(abs(a)); %对A和B信号分别加了躁声后的实际信号仿真做fft,并且求相位 d=angle(i*real(a(I))-imag(a(I)))*180/pi a=fft(O); [F,I]=max(abs(a)); e=angle(i*real(a(I))-imag(a(I)))*180/pi toc |
|