output problem for penumbral blur. plzzzz check the program
1 次查看(过去 30 天)
显示 更早的评论
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
0 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Matrix Indexing 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!