嗨,伙计们,我的代码没有运行我无法调试,能有人帮我吗?
23日视图(30天)
显示旧的评论
%常数和参数
α= 2;%剪切波速
β= 1;%压缩波速度
频率= 0:15;%的频率范围
ε= 1平台以及;%小ε值
ρ= 1;%定义ro的价值
%定义源和接收器的位置
源= (0,0,0);%源位置(x, y, z)
站= [0.01,0,0;0.02、0、0;0.03、0、0;0.04、0、0;0.05,0,0);%接收机位置(x, y, z)
%计算不同频率和im G11站
结果= 0(长度尺寸(站,1),(频率));
为i = 1:尺寸(站,1)
接收机=站(我:);%接收机位置(x, y, z)
距离=规范(接收机-源);%源和接收器之间的欧氏距离
γ=(接收机-源)/距离;%单位向量从源到接收机
γ= repmat(伽马、长度(频率),1);%对每个频率重复γ
为j = 1:长度(频率)
ω= 2 *π*频率(j);
q =ω/α;
f1 =(β^ 2 /α^ 2)* (1 - 1 j * 2 . / (q *距离+ε)- 2。/(问^ ^ 2 + 2 *距离ε))* exp (1 j * q *距离);
f2 =(β^ 2 /α^ 2)* (1 j。/ (q *距离+ε)+ 1 /(问^ ^ 2 + 2 *距离ε))* exp (1 j * q *距离);
结果(i, j) = calculate_G11_imag (f1, f2,γ(j,:),ω,罗依,α,β);
结束
结束
%绘制
图;
为i = 1:尺寸(站,1)
情节(频率、图像放大(结果(我:))“DisplayName的”sprintf (“站:% s”mat2str(站(我,:))));
持有在;
结束
持有从;
包含(“频率”);
ylabel (“im G11”);
标题(“im G11不同电台的变化”);
传奇(“位置”,“东北”);
网格在;
ylim ([10]);
%定义函数来计算im G11
函数结果= calculate_G11_imag (f,γ,ω,罗依,α,β)
大小(f)
大小(γ)
结果=ω/(12 *π* ro) *(1 /(α^ 3)+ 2 /(β^ 3))* (f。*γ)总和;
结束