curvefit10.m

数值分析与计算的PPT幻灯片与MAT源码

%curvefit09.m

q=[1.05 1.18 1.29 1.30 1.30 1.42 1.50 1.52 1.46 1.60 1.69 ...
   1.81 1.93 1.95 2.01 2.00 2.09 1.96 2.20 2.12 2.16 2.08 ...
   2.24 2.56 2.34 2.45 2.58];
   
 
k=[1.04 1.06 1.16 1.22 1.27 1.37 1.44 1.53 1.57 2.05 2.51 ... 
   2.63 2.74 2.82 3.34 3.24 3.61 4.10 4.36 4.77 4.75 4.54 ...
   4.54 4.58 4.58 4.58 4.54];
   
l=[1.05 1.08 1.18 1.22 1.17 1.30 1.39 1.47 1.31 1.43 1.58 ...
   1.59 1.66 1.68 1.65 1.62 1.86 1.93 1.96 1.95 1.95 1.58 ...
   1.67 1.82 1.60 1.61 1.64];

n=length(l);
for i=1:n
   z0=log(q);
   x0=log(k);
   y0=log(l);
end
a=[ones(n,1),x0',y0'];
b=z0';

k0=1:0.1:5;
l0=1:0.1:2;
n0=length(k0);
m0=length(l0);
s=0;
for i=1:n0
   for j=1:m0
      z=1.1790*(k0(i))^0.4193*(l0(j))^0.0467;
      s=s+1;
   end
end
s1=s
axis([0 5 0 2 0 5])
plot3(k,l,q)
grid on

hold on
mesh(k0,l0,z)