“…*4./50,'k-'); end grid; xlabel('Frequency (kHz)','fontsize',18) ylabel('Velocity (mm/s)','fontsize',18) set(gca,'fontsize',16) end hold off eps=zeros(num/2,10); for n=1:1:num/2 for m=1:1:10 eps(n,m) = b(n,m)./5./c(n). *1e6; end end eps(13,:) = eps(13,:)./5; p=zeros(num,2); t=0:2:24; for n=1:1:num/2 p(n,:)=polyfit(eps(n,:),g(n,:),1); yfit(n,:) =polyval(p(n,:),eps(n,:)); yresid(n,:)=g(n,:)-yfit(n,:); SSresid=sum(yresid(n,:).^2); SStotal=(length(g(n,:))-1) * var(g(n,:)); rsq(n)=1-SSresid/SStotal; end for n=1:1:(num/2) pt(n)=p(n,1); fit(n,:)=(p(n,1)*(0:.2:5) + p(n,2)); end figure plot(t,abs(pt)) xlabel('Time (days)','fontsize', 18 ;w=2*pi*P(1,1,n-1); %young(n/2)=m*w*c(n/2)/tan(w*L/c(n/2))/S; rho=8030; young(n/2)=rho*c(n/2).^2; maxf(n/2)=(P(1,1,n-1)-P(1,37,n))/P(1,1,n-1); end…”