clc
clear
close all
aa2=0;
ka=[5,10,15];
ka=ka./10;
for i=1:3
ka(i)=10^ka(i);
end
aa2=1;
A=sqrt(ka.*2*aa2);
r=0:0.05:9;
for k=1:3
for i=1:length(r)
fr(k,i)=r(i)*2/aa2*exp(-r(i)^2/aa2-ka(k))*besseli(0,2*r(i)*sqrt(ka(k)
/aa2));
end
end
figure
plot(r,fr(1,:),'m-');
hold on;
plot(r,fr(2,:),'r-.');
hold on;
plot(r,fr(3,:),'--');
hold on
fm=150;
f2=0:1:149;
lenf2=length(f2);
for i=1:lenf2
Hw(i)=(1-(f2(i)/fm)^2)^(-0.5);
end
w=zeros(1,lenf2*4);
w(1:lenf2)=Hw;
for i=1:lenf2*2
w(lenf2*2+i)=w(lenf2*2-i+1);
end
hn=ifft(w,lenf2*4);
hn=real(hn);
hw=fft(hn);
figure
plot(Hw,'ro');
评论2