我的MATLAB程序点RUN运行后没任何反应。既不出图,workspace里面也没有任何矩阵量,也不报错。随便改个错误出来就会报错,workspace里面也有值了.运行其他程序也都正常,就是这个没任何反应。换到别人的电脑上运行也是一样。附程序:
clear;
clc;
A=[10.25,13.9,17.3,20.7,24.1,27.5,30.9,34.3];%A矩阵为8个阵元的平均半径;
c=5.6e6;%C为声速;
li=2;%ri为阵元宽度;
y0=75;%y0为预期聚焦深度;
Ti=zeros(1,8);
for n=1:8
Ti(1,n)=(sqrtm(A(1,n)^2+y0^2)-y0)/c;%延时Ti;
end
[x,y]=meshgrid(0:0.1:40,60:0.1:100);%声场范围,这里取X方向0~40mm,Y方向60~100mm;
z=zeros(401);
for m=1:400;
for n=1:400
for p=1:8
R0i=sqrtm(y0^2+A(1,p)^2);
Ri=sqrtm(y(m,n)^2+A(1,p)^2);
z0=quad('sin(x)',0,2*3.1415926);
%z0=SmartSimpson('A(1,p)*li*exp(-j*(w/c)*(Ri-R0i))/Ri',0,2*3.1415926);
z(m,n)=z(m,n)+z0;
% Pi=li*A(1,p)*exp((-j)*(x/c)*((sqrtm(y0^2+(0.1*(Bx-1bx))^2+A(1,p)^2-0.2*(BX-1)cosx)-R0i)/ Ri=sqrtm(y0^2+(0.1*(Bx-1))^2+A(1,p)^2-0.2*(BX-1)cosx)% Ri=sqrtm(y0^2+(0.1*(Bx-1))^2+A(1,p)^2-0.2*(BX-1)cosx)
end
end
end
mesh(x,y,z);