function NU=carpenteri(beta,d)a=beta(1);b=beta(2);NU=(a+b./d).^0.5;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%d=[1,2,3,4];NU=[40,30,25,22];[A,res]=lsqcurvefit('carpenteri',ones(1,2),d,NU);a=A(1);b=A(2);x=[1:0.1:100];y=(a+b./x).^0.5;plot(x,y,'r-');hold on;scatter(d,NU);axis([0 100 0 100]);xlabel('试件特征尺寸L','FontSize',16)ylabel('名义断裂能','FontSize',20)title('carpenteri分形','FontName','隶书' ,'FontSize',20)legend('fenxing','shuju')grid on现在想接着求给出的数据，与拟合出来的公式之间的相关系数，应该怎么编程？？
|
