各位大虾:
目前在matlab绘制等高线上遇到了一些障碍,我想绘制一个类似于下图的一个等高线图
主要是吸波材料中反射率的等高线图,x轴为频率,y轴为涂层厚度,现已有设定涂层厚度d下,反射率的matlab计算程序,现要求涂层厚度d在[0,10]内等间距的变化,间距为0.1,求最终作的反射率R的等高线图。
clear;
figure;
M=1,N=1,A=1,B=1;
DATA=load('F:\博士学位论文资料\datam\2012.3.12\datam\10.txt'); %数据导入
f=10^9*DATA(:,1);
epsilon_real=M*DATA(:,2);
epsilon_imag=N*DATA(:,3);
mu_real=A*DATA(:,4);
mu_imag=B*DATA(:,5);
mu=mu_real+i*mu_imag;
epsilon=epsilon_real+i*epsilon_imag;
c=3*10^8;
d=2*10^-3;
K=2*pi*f.*sqrt(conj(mu).*conj(epsilon))/c;
eta=sqrt(mu_real.*(1-i*tan(angle(mu)))./(epsilon_real.*(1-i*tan(angle(epsilon)))));
Z=eta.*tanh(i*K*d);
gamma=(Z-1)./(Z+1);
R=20*log10(abs(gamma));
a=f/10^9;
plot(a,R,'-');
xlabel('f/GHz');ylabel('R/dB');%2~18GHz
hold on;
grid on;
![]()
1111.jpg
![]()
11111111111.png |