output problem for penumbral blur. plzzzz check the program

1 次查看（过去 30 天）

TULIKA 2015-3-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/182422-output-problem-for-penumbral-blur-plzzzz-check-the-program

    X =linspace(-8,8,10);
    g=.000080;
    Lambda=.000000001;
    R1=(1.-(X.^2));
    H=(exp(i*2*pi.*R1)).*g/Lambda;
    plot(X,H);
    R=X.^2;
    Lambda=.000000001;
    dis=10;
    F = (exp(i.*pi.*R))./(Lambda.*dis);
    plot(X,F);
    plot(X,H);
    l=15e-4;
    w=4e-4;
    Ts=1e-7;
    t0=0:Ts:w;
    t1=w+Ts:Ts:w+(5000*Ts);
    t2=w+(5001*Ts):Ts:l-Ts;
    s1=ones(size(t1));
    s0=zeros(size(t0));
    s2=zeros(size(t2));
    s=[s0 s1 s2];
    N=1; %How many repetitions of the wave
    ss=repmat(s,1,N);
    t=0:Ts:N*l-Ts;
    plot(t,ss);
    F=imresize(F,[1 15000]);
    H=imresize(H,[1 15000]);
    az=ss.*F;
    X=imresize(X,[1 15000]);
    plot(X,az);
    plot(t,az);
    az1=fftshift(fft2(fftshift(az)));
    az2=fftshift(fft2(fftshift(H)));
    az3=az1.*az2;
    az4=fftshift(ifft2(fftshift(az3)));
    plot(X,az4);
    plot(t,az4);
    there is no error but the output is unacceptable.output should be like
    <<http://in.mathworks.com/help/matlab/examples/square-wave-from-sine-waves.html>>
     last image of this link