Matlab中的最佳拟合线

Question

我需要一条最适合此图中显示的每个数据点的线。 到目前为止，这是我的代码。

load xsarea
load width
load depth
load da

figure;

loglog(da,width,'s',da,depth,'o',da,xsarea,'^');
xlabel('drainage area');
ylabel('depth/width/area');
hold on;

P = polyfit(log(da),log(width),1);
yfit = exp(P(1)+P(2).*log10(da));
plot(yfit,da,'ko');
label = ['log(depth) = 'num2str(p(1)) 'log(width) + 'num2str(p(2))];

你可以看到它是对数的。完成后，图表看起来应该像this document中的图2一样。一切都很好，除了我需要显示实际的行。