res=[100;20];
thick=1000;
u=4*pi*1e-7;
load r.dat
rr=length(r);
for i=1:rr
w=2*pi*r(i);
u=4*pi*1e-7;
k=complex(1,-1)*sqrt((u*w)/2.0);
Z=k*sqrt(res(nlayer));
for j=2:nlayer
j1=nlayer-j+1;
Z0=k*sqrt(res(j1));
e1=exp(-2*k*thick(j1)/sqrt(res(j1)));
e=(1-e1)/(1+e1);
Z=Z0*(Z/Z0+e)/(Z*e/Z0+1);
end
rho(i)=abs(Z)^2/(u*w)
end
%作图
axes(handles.axes1);
loglog(r,rho,'r-o');
set(gca,'xdir','reverse');
grid on;
axis([10^-4 10^3 1 10^-4 10^3]);
title('大地电磁测深视电阻率曲线图');
legend('原始数据');
xlabel('f/Hz');ylabel('RHO /ohm');
function zhengyanedit1_Callback(hObject, eventdata, handles)
% hObject handle to zhengyanedit1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of zhengyanedit1 as text
% str2double(get(hObject,'String')) returns contents of zhengyanedit1 as a double
% --- Executes during object creation, after setting all properties.
function zhengyanedit1_CreateFcn(hObject, eventdata, handles)
% hObject handle to zhengyanedit1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called