首页->【超声成像实验室】

8 1

% This Script can describe the probe properity;
E0=0;E1=0;
E2=0;E3=0;

f0=6.62e6;  %probe Center ferq
c=1540;    %Ultrasound Speed lamada=
k=2*pi*f0/c;    
d0=0.5/1000;  % 0.25 lamda=0.232mm
d1=0.05/1000;  % 0.25 lamda
angle = -pi:0.005:pi;

N0=16;
N1=32;

for n=1:N0
    E0=E0+exp(-1i*(n-1)*k*d0.*cos(angle));
    E1=E1+exp(-1i*(n-1)*k*d1.*cos(angle));
    %E=E+exp(1i*(n-1)*k*d.*cos(angle));
end

for n=1:N1
    E2=E2+exp(-1i*(n-1)*k*d0.*cos(angle));
    E3=E3+exp(-1i*(n-1)*k*d1.*cos(angle));
    %E=E+exp(1i*(n-1)*k*d.*cos(angle));
end

%E=sin(k*N*0.5*d*sin(angle))./sin(k*0.5*d*sin(angle));
ee0=abs(E0)/max(abs(E0));
ee1=abs(E1)/max(abs(E1));
ee2=abs(E2)/max(abs(E2));
ee3=abs(E3)/max(abs(E3));
fontsize=12;

subplot(1,2,1),p0=polarplot(angle,ee0,'b',angle,ee2,'r','LineWidth',1);...
    title('d=0.5mm,N=16/32线性阵元探头指向性','FontName','TIMES NEW ROMAN|宋体'...
    ,'FontSize',fontsize,'FontUnits','points');
ax=gca;
ax.FontName='TIMES NEW ROMAN |宋体';
ax.FontSize=fontsize;
ax.FontUnits='points';
ax.ThetaGrid='on';
ax.GridLineStyle='--';
ax.GridColor=[0 0 0];
ax.LineWidth=1;
legend('16阵元','32阵元','Location','northeast');
subplot(1,2,2),p1=polarplot(angle,ee1,'b',angle,ee3,'r','LineWidth',1);
    title('d=0.05mm,N=16/32线性阵元探头指向性','FontName','TIMES NEW ROMAN|宋体'...
    ,'FontSize',fontsize,'FontUnits','points');
ax=gca;
ax.FontName='TIMES NEW ROMAN|宋体';
ax.FontSize=fontsize;
ax.FontUnits='points';
ax.GridColor=[0 0 0];
ax.LineWidth=1;
ax.GridColorMode='auto';
ax.ThetaGrid='on';
ax.GridLineStyle='--';
legend('16阵元','32阵元','Location','northeast');


仿真结果

TIM截图20180402204355.png


楼主可见

  1. shwnyoo 1#

    医用超声探头指向性MATLAB例程,