把三段程序 放在一起 第一段程序结束时加上 hold on 成为:(%scatter(t,exp(y))取消了数据点,想要数据点的话去掉%)
t=[ 6.855 9.854 17.14 22.14 27.99 33.13 37.99 41.13 45.99 54.13 54.98 60.13 64.7 69.84 76.12 82.98 88.83 ];
y=[ 2.558e+005 8.221e+005 4.275e+005 4.735e+005 1.705e+005 1.918e+005 1.866e+005 1.624e+005 1.49e+005 1.016e+005 1.107e+005 4.717e+004 5.386e+004 5.715e+004 1.576e+004 3.725e+004 3.319e+004 ];
y=log(y);
x=lsqcurvefit(@(x,t)2*log(t)+x(1)+x(2)*t+x(3)*t.^2,[1 1 -10],t,y)
%scatter(t,exp(y))
t=0:1:100;
y=2*log(t)+x(1)+x(2)*t+x(3)*t.^2;
plot(t,exp(y),'r-')
hold on
t=[2.999 7.798 13 17.99 22.99 27.99 33.13 37.99 42.99 47.99 53.13 58.13 63.27 68.27 72.98 77.98 82.98 86.98 92.97 97.97];
y=[1.372e+005 5.62e+005 3.266e+005 1.993e+005 5.324e+005 1.705e+005 1.918e+005 1.866e+005 2.409e+005 6.98e+004 3.323e+004 9505 8.387e+004 1.173e+005 4.554e+004 6.036e+004 3.725e+004 1.629e+004 5.337e+004 3.371e+004];
y=log(y);
x=lsqcurvefit(@(x,t)2*log(t)+x(1)+x(2)*t+x(3)*t.^2,[1 1 -10],t,y)
%scatter(t,exp(y))
t=3:1:98;
y=2*log(t)+x(1)+x(2)*t+x(3)*t.^2;
plot(t,exp(y),'g-')
t=[0.997 6.141 11 16 21.28 26.14 30.99 35.85 41.13 45.99 50.84 55.84 60.55 65.98 71.12 76.12 81.12 86.26 91.4 95.97];
y=[4.973e+004 1.261e+005 9.898e+005 4.182e+005 3.616e+005 9.781e+004 6.21e+004 5.806e+005 1.624e+005 1.49e+005 1.224e+005 1.244e+005 4.205e+004 5.005e+004 4.157e+004 1.576e+004 1.414e+004 4.875e+004 8531 5.323e+004];
y=log(y);
x=lsqcurvefit(@(x,t)2*log(t)+x(1)+x(2)*t+x(3)*t.^2,[1 1 10],t,y)
%scatter(t,exp(y))
t=1:1:96;
y=2*log(t)+x(1)+x(2)*t+x(3)*t.^2;
plot(t,exp(y),'b-')
xlabel('频率 ')
ylabel('振幅')