各位高手来看看。我不知道这个怎么弄

2019-07-18 12:25发布


我想把以下程序编为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 分
友情提示: 此问题已得到解决,问题已经关闭,关闭后问题禁止继续编辑,回答。