|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
code 如下:
clc
clear
x1=linspace(0.5,0.78,1e2);
m1=[15.5 14.9 14.4 13.9 13.4 13.0 12.6 12.2 11.9 11.5 11.2 10.9 10.6 10.3 10.1];
m2=[1.205 1.13 1.063 1.002 0.974 0.896 0.850 0.807 0.767 0.730 0.696 0.665 0.636 0.609 0.583 ];
x2=[0.5 0.52 0.54 0.56 0.58 0.60 0.62 0.64 0.66 0.68 0.70 0.72 0.74 0.76 0.78];
a1=m2.^2-m1.^2;
a2=2.*m1.*m2;
a3=a1+1i.*a2;
a3c=interp1(x2,a3,x1,'spline');
plot(x1,real(a3c))
hold on
plot(x1,imag(a3c))
插值完还是分离的数值点,如果我得到x1=0.643的值,这个要怎么取出来?
|
|