clear all
clc
A=readmatrix('Point_30mm.txt');
B=readmatrix('Point_20mm.txt');
x=0.01;
l=0.03;
for n=1:198
f=20+5n
p1=A(n,1)
p2=B(n,1)
h12=p2./p1
k=2pif/340
hi=exp(-1ikx)
hr=exp(1ikx)
c=(h12-hi)/(hr/h12)
r=cexp(1i2k*l)
R=abs(r)
a=1-R
figure(1)
plot (f,a,'.g',f,a,'k--')
xlim([25 1000])
end
画图为啥放在循环里了……
每次画一个点?那你不是应该用散点图画吗?而且如果你想做全部画在一张图片上,plot后面要加hold on