I'm trying to figure out why after i run this code I'm not getting any messages and no plots are popping up

Question

Cheryl 2013-12-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/109880-i-m-trying-to-figure-out-why-after-i-run-this-code-i-m-not-getting-any-messages-and-no-plots-are-pop

回答： Image Analyst 2013-12-16

if true
  %close all
clear all
num_fft=1024;
dt=1;

f= [ 0 0 0 0 1 0 0.5 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0]

h=[ 2.76077e-006 0.000746586 0.0407636 0.449702 1.02038 0.673993 0.540762,... 0.225414 0.0211304 0.0411355 0.44933 1 0.449329 0.040762 0.000746586,... 2.76077e-006 0 0 0 0 0 0 0 0]

N= numel(h);

for i = 1:N

    x(i)=(i-1);
    H=fft(h,num_fft); %Take the FFT of h
    for j=1:num_fft;
        dfh(j)=(((1/num_fft*dt))*(j-1))-0.5;%Calculate frequency bins for H
    end
    HR=real(H);%Calculate real part of H;
    HI=image(H);%Calculate imaginary part of H;
    for i=1:num_fft;
        HM(i)=sqrt((HR(i)^2)+((HI(i)^2)));%calculate magnitude of H components
    end
    figure(1)
    plot((dfh),fftshift(HM),'k');
    title('Magnitude of H');
    xlabel('freq axis')
    ylabel('Magnitude axis');
    Htheta = fftshift(atan2(HI,HR));% Calculate the phase for H
    figure (2)
    plot(dfh,unwrap((Htheta),.1),'k');
    title('Phase Angle of H');
    xlabel('freq axis')
    ylabel('Phase axis')
    F=fft(f,num_fft); %Take the FFT of f
    for j=1:num_fft
        dfg(j)=(((1/num_fft*dt))*(j-1))-0.5;%Calculate frequency bins for F
    end
    FR=real(F);
    FI=image(F);
    for i=1:num_fft;
        FM(i)=sqrt((FR(i)^2)+((FI(i)^2)));%Calculate magnitude of F components
    end
    figure(3)
    plot((dfh),fftshift(FM),'k');
    title('Magnitude of F');
    xlabel('Freq axis');
    ylabel('Phase axis');
    FH=H
    for i=1:num_fft;
        if FM== 0;
            FDF(i)=0;
        else
            FDF(i)=FH(i)./F(i);
        end
    end
    for j=1:num_fft;
        dfg(j)=(((1/num_fft*dt))*(j-1))-0.5; 
    end
    GDFR=real(GDF);%calculate the real part of F;
    GDFI=image(GDF);%calculate the imaginary part of F;
    for i=1:num_fft;
        FDFM(i)=sqrt((GDFR(i)^2)+((GDFI(i)^2)));%Calculate magnitude of F components
    end
    figure (4)
    plot((dfh),fftshift(GDFM),'k');
    title('Magnitude of GDF');
    xlabel('freq axis')
    ylabel('Magnitude axis')
    FDFtheta=fftshift(atan2(GDFI,GDFR));%Calculate the phase for F
    figure (5)
    plot (dfh,unwrap((GDFtheta),.1),'k');
    title('Phase Angle of GDF');
    xlabel('freq axis');
    ylabel('Phase axis');
    Fdf=ifft(GDF,num_fft);
    Plot(GDF,'k');
    title('unknown g')
    xlabel('freq axis');
    ylabel('Phase axis');
end