波束形成

波束形成1.问题描述:数字波束形成器是全数字化超声成像的基础,也是高性能彩超的保证。数字波束形成包括发射和接收两个部分。数字是接收波束形成的关键技术,它通过使用顺序储存器FIFO或随机存取存储器双端口RAM替代模拟式波束形成器中的LC延时线来实现波束聚焦,即以数字延时补偿替代模拟延时的补偿。数字延时不仅能实现精确延时补偿,实现所谓的逐点跟踪式动态聚焦,还能方便实现动态孔径、动态变迹控制,克服模拟式延时补偿存在的诸多固有缺点,通道数增加不受限制,是图像品质得以全面提高。2.部分程序:close..

大家好,又见面了,我是你们的朋友全栈君。

1.问题描述:

 数字波束形成器是全数字化超声成像的基础,也是高性能彩超的保证。

数字波束形成包括发射和接收两个部分。数字是接收波束形成的关键技术,它通过使用顺序储存器FIFO或随机存取存储器双端口RAM替代模拟式波束形成器中的LC延时线来实现波束聚焦,即以数字延时补偿替代模拟延时的补偿。数字延时不仅能实现精确延时补偿,实现所谓的逐点跟踪式动态聚焦,还能方便实现动态孔径、动态变迹控制,克服模拟式延时补偿存在的诸多固有缺点,通道数增加不受限制,是图像品质得以全面提高。

2.部分程序:

 

close all
clear all

%%波束形成
f=[4000,5000,6000];%信号频率,单位Hz
w=[1,1,1,1,1;%每一行表示一种加权
    1,2,4,2,1;
    4,2,1,2,4];
fs=4*f;%采样频率取信号频率的4倍
T=0.1;%信号长度,单位s
N=T*fs;%采样点数
c=1500;%声速,m/s
d=0.15;%阵元间距,单位m
jay=sqrt(-1);%复单位
theta=-90:1:90;%观测角度范围,单位度
tao=d*sin(pi*theta/180)/c;%相邻阵元时延
xx=zeros(3,2400);%准备产生三种频率下的参考信号

for i=1:3
    xx(i,1:N(i))=exp(jay*(2*pi*f(i)/fs(i)*(0:N(i)-1)));%三种频率下的参考信号
end

for j=1:3%三种频率循环
    for i=1:length(tao)%不同角度循环
        for k=(j-1)*5+1:(j-1)*5+5%第一种频率下的五个阵元的信号
            x(k,:)=xx(j,:)*exp(-jay*2*pi*f(j)*(k-1)*tao(i));
        end
        r=x((j-1)*5+1:(j-1)*5+5,:)*x((j-1)*5+1:(j-1)*5+5,:)’/N(j);%求某一角度下的波束图,用公式B=w*(∑x*x’)*w’/N.
        B((j-1)*3+1,i)=w(1,:)*r*w(1,:)’;%第j个频率下的第一种加权的波束图输出
        B((j-1)*3+2,i)=w(2,:)*r*w(2,:)’;%第j个频率下的第二种加权的波束图输出
        B((j-1)*3+3,i)=w(3,:)*r*w(3,:)’;%第j个频率下的第三种加权的波束图输出
    end
end
%c产生的举证B表示B=[B11;B12;B13;B21;B22;B23;B31;B32;B33];第一个下角标表示第几个频率,第二个下角标表示第几种加权而得到的波束
for k=1:9%分别对9个波束图进行归一化和取对数
    B(k,:)=real(B(k,:))/max(real(B(k,:)));%归一化
    B(k,:)=20*log10(B(k,:));%取对数
end
figure(1);%频率f=4000,加权不同
subplot(3,1,1);
plot(theta,B(1,:));axis([-90 90 -60 0]);title(‘频率f=4000,加权[1,1,1,1,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,2);
plot(theta,B(2,:));axis([-90 90 -60 0]);title(‘频率f=4000,加权[1,2,4,2,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,3);
plot(theta,B(3,:));axis([-90 90 -60 0]);title(‘频率f=4000,加权[4,2,1,2,4]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);

figure(2);%频率f=5000,加权不同
subplot(3,1,1);
plot(theta,B(4,:));axis([-90 90 -60 0]);title(‘频率f=5000,加权[1,1,1,1,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,2);
plot(theta,B(5,:));axis([-90 90 -60 0]);title(‘频率f=5000,加权[1,2,4,2,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,3);
plot(theta,B(6,:));axis([-90 90 -60 0]);title(‘频率f=5000,加权[4,2,1,2,4]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);

figure(3);%频率f=6000,加权不同
subplot(3,1,1);
plot(theta,B(7,:));axis([-90 90 -60 0]);title(‘频率f=6000,加权[1,1,1,1,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,2);
plot(theta,B(8,:));axis([-90 90 -60 0]);title(‘频率f=6000,加权[1,2,4,2,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,3);
plot(theta,B(9,:));axis([-90 90 -60 0]);title(‘频率f=6000,加权[4,2,1,2,4]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);

figure(4);%加权[1,1,1,1,1],频率不同
subplot(3,1,1);
plot(theta,B(1,:));axis([-90 90 -60 0]);title(‘频率f=4000,加权[1,1,1,1,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,2);
plot(theta,B(4,:));axis([-90 90 -60 0]);title(‘频率f=5000,加权[1,1,1,1,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,3);
plot(theta,B(7,:));axis([-90 90 -60 0]);title(‘频率f=6000,加权[1,1,1,1,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);

figure(5);%加权[1,2,4,2,1],频率不同
subplot(3,1,1);
plot(theta,B(2,:));axis([-90 90 -60 0]);title(‘频率f=4000,加权[1,2,4,2,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,2);
plot(theta,B(5,:));axis([-90 90 -60 0]);title(‘频率f=5000,加权[1,2,4,2,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,3);
plot(theta,B(8,:));axis([-90 90 -60 0]);title(‘频率f=6000,加权[1,2,4,2,1]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);

figure(6);%加权[4,2,1,2,4],频率不同
subplot(3,1,1);
plot(theta,B(3,:));axis([-90 90 -60 0]);title(‘频率f=4000,加权[4,2,1,2,4]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,2);
plot(theta,B(6,:));axis([-90 90 -60 0]);title(‘频率f=5000,加权[4,2,1,2,4]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);
subplot(3,1,3);
plot(theta,B(9,:));axis([-90 90 -60 0]);title(‘频率f=6000,加权[4,2,1,2,4]’);xlabel(‘角度/°’);ylabel(‘波束图/dB’);

%计算主瓣宽度D和旁瓣级P
for k=1:9
    [i,j]=find(abs(B(k,:)+3)<0.5);%误差限设在0.5dB。
    D(k)=(j(length(j))-j(1))*1;%每一个点之间是一度
end    %根据计算D的结果可以发现,同一种频率下三种加权分辨率最高的是第三种,其次是第一种。同一加权下,频率越高,分辨率越高
disp(‘9个波束图的参数依次为,其中分别表示:’);
disp(‘f=4000,w=[1,1,1,1,1];f=4000,w=[1,2,4,2,1];f=4000,w=[4,2,1,2,4]’);
disp(‘f=5000,w=[1,1,1,1,1];f=5000,w=[1,2,4,2,1];f=5000,w=[4,2,1,2,4]’);
disp(‘f=6000,w=[1,1,1,1,1];f=6000,w=[1,2,4,2,1];f=6000,w=[4,2,1,2,4]’);
disp(‘波束宽度为:’);
ss=sprintf(‘%d ‘,D);
disp(ss);
NN=length(theta);%每一个波束图离散的点数
for k=1:9
    for i=(NN-1)/2+1:NN-1%从主波束开始往后找到第一个次级大的位置
        if B(k,i)>=B(k,i+1)&&B(k,i)>=B(k,i-1)
            CC(k)=B(k,i);%记录第一个次级大的分贝数
        end
    end
end
CC(2)=B(2,1);CC(5)=B(5,1);%由于在[-90,90]度的范围内,第二种和第五种情况下没有找到次级大,故用它在-90°的地方分贝数近似。
disp(‘旁絆级为:’);
s=sprintf(‘%f ‘,CC);
disp(s);

3.仿真结论:

波束形成

D00003

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/147650.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • windows搭建php环境_如何安装php环境

    windows搭建php环境_如何安装php环境使用apache+php+mysql的方式一、安装mysql1、下载mysql,下载地址:https://dev.mysql.com/downloads/installer/2、下载完成后解压我解压的路径是:D:\ProgramFiles\mysql-8.0.253、配置文件首先在解压的路径下查看是否含有my.ini的文件,如果没有则新建一个,内容如下:[mysqld]#设置3306端口port=3306#设置mysql的安装目录basedir=D:\ProgramFiles\

  • oracle未明确定义列[通俗易懂]

    oracle未明确定义列[通俗易懂]分类:SQL使用技巧2012-04-0616:121332人阅读评论(0)收藏举报运行环境:Oracle10gsqlplus环境下。 在查询语句中,经常会出现一个错误: SQL基础:ORA-00918:未明确定义列的错误。 当前遇到有两种情况。原因为:当查询语句中,查询的表(数据集)中有相同的字

  • Shell脚本备忘录

    Shell脚本备忘录目录1.jq1.1安装1.2几个常用例子1.2.1取出数组index=0的内容1.2.2取出数组index=0的name的内容1.2.3以key-value的格式取出数组index=0的name和city1.2.4以key-value的格式取出所有数组的name和city1.2.5以key-value的格式取出数组index=0的name和arrayBrowser的index=1的url1.2.6以key-value的格式取出所有数组的name和city并放在一个数组里(前后加上[])1.

  • 研招网的准考证还能下载吗_考研准考证还能下载吗

    研招网的准考证还能下载吗_考研准考证还能下载吗下载打印方式往下看,有手把手教学。【考生请注意】据了解,《准考证》由考生使用A4幅面白纸在规定时间内(2021年12月10日开始至12月27日)上网自行下载打印。《准考证》正、反两面在使用期间均不得涂改或书写。考生凭下载打印的《准考证》及居民身份证参加考研初试和复试。考生凭本人《准考证》及有效居民身份证按规定时间进入考场,对号入座。入座后将上述证件放在桌面左上角,以便核验。请考生下载《准考证》前认真阅读报考点所在地省(区、市)考生防疫有关要求,并认真关注《准考证》上“省级招生考试…

  • unity安装到一半失败

    unity安装到一半失败从unityhub安装unity总是到一半就提示失败,之前试过用取消同时安装visualstudio的选项,但后期重新安装很麻烦。从unity官方论坛找到了这些解决办法,一一排查过了,发现去掉Documentation可选项安装成功,visualstudio也正常安装了。…

  • RSA加密算法c++实现

    RSA加密算法c++实现RSA加密算法c++实现 先码个代码:/*************************************************************************&gt;FileName:RSA.cpp&gt;Author:YB&gt;E-mail:134…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号