[发明专利]基于介质有限元叠加的磁声耦合成像声信号求解方法

权利要求书

1.一种基于介质有限元叠加的磁声耦合成像声信号求解方法，包括有基于有限元平台的空间求解和基于时间叠加的时间求解两部分，其特征在于，具体是：首先进行磁声耦合成像介质的空间剖分求解；然后根据得到的空间剖分求解结果，进行磁声耦合时间域求解。

2.根据权利要求1所述的基于介质有限元叠加的磁声耦合成像声信号求解方法，其特征在于，在进行磁声耦合成像介质的空间剖分求解中，对磁声成像介质进行多物理场有限元剖分求解，将各剖分单元编号为1,2，....n，记录剖分单元的空间位置r_i，其中i=1,2，...，n，记录求解的各剖分单元电流密度幅值a_i，所述的n取100以上的整数。

3.根据权利要求1所述的基于介质有限元叠加的磁声耦合成像声信号求解方法，其特征在于，所述的进行磁声耦合时间域求解包括如下步骤：

（1）根据各剖分单元的空间位置r_i计算各剖分单元到超声换能器的距离li和各剖分单元到超声换能器的距离l_i与换能器轴线的角度θ_i；

（2）设介质声学均匀，由公式delay_i=l_i/v计算由各剖分单元到超声换能器的距离l_i形成的时间延迟delay_i，其中v为超声在介质中的传播速度；

（3）设置磁声耦合激励信号的时间波形f(t)；

（4）根据激励信号的时间波形f(t)计算各剖分单元的电流密度时域信号波形J_i(t)=a_if(t)；

（5）将步骤（4）求得的各个剖分单元的电流密度时域信号波形J_i(t)，以及步骤（1）求得的距离l_i和角度θ_i，步骤（2）求得的传播时间延迟delay_i代入公式：

p(r,t)=-14πΣicosθi▿·[Ji(t+delayi)·B0]]]>

求解，即得到磁声信号时域波形解，其中B₀为静磁场。