Hi Bradley,

Thanks for your reply.After reading photoacoustic examples,I added source.p0 = smooth(source.p0,true) in my MATLAB program,and I tried to adjust the value of dt Nt to change the sampling frequency,the imaging results still aliasing.I don't know why this happens.I want to try to reproduce part speckle illumination imaging of this article(https://aip.scitation.org/doi/10.1063/1.5113476).Do you have some more detailed suggestions to solve imaging aliasing？

And I have another question.In linear photoacoustic,incident laser(Gaussian) and radiation area keep invariant，photoacoustic conversion efficiency is a constant.And now,I just change the incident laser to Gaussian speckle light，the results in simulation,I found that calculating the photoacoustic conversion efficiency of each different speckle is different,I am quite confused.

Sincerely look forward to your answer.

Lkk