我想把以下程序编为
labview的程序,我已经做好了三步了,但是后面的循环不会了,请指教
Y=fft(Sheet)
Pyy = abs(Y)
X = imag(Y)
Z = real(Y)
f = 100*(0:2048)/4096;
figure(1)
plot(f,Pyy(1:2049))
hold on
plot(f,X(1:2049))
hold on
plot(f,Z(1:2049))
title('Frequency content of Sheet')
xlabel('frequency (MHz)')
Y1=fft(Sheet1)
Pyy1 = abs(Y1)
X1 = imag(Y1)
Z1 = real(Y1)
f1 = 100*(0:2048)/4096;
figure(2)
plot(f1,Pyy1(1:2049))
hold on
plot(f1,X1(1:2049))
hold on
plot(f1,Z1(1:2049))
title('Frequency content of Sheet1')
xlabel('frequency (MHz)')
G=X.*X1+Z.*Z1
K=X.*Z1-Z.*X1
M=atan(-G./K)
figure(3)
plot(f1,M(1:2049))
b=3.1415926
c=150
for i=2:1:c
a=abs(M(i+1)-M(i))
if a>1
for(j=i+1:1:c+1)
M(j)=M(j)-b
end
end
end
figure(4)
plot(f1,M(1:2049))
t=150
for n=1:1:t
E=M(n+2)-M(n+1)
B=(f1(n+1)-f1(n))*2*3.1415926
Q(n)=abs(E./B)
end
figure(5)
plot(f(1:1:150),Q(1:1:150))
-
-
使用FFT分析信号实部和虚部(修改版).vi
下载积分: 积分 -1 分
24.4 KB, 下载次数: 2, 下载积分: 积分 -1 分
一周热门 更多>