untitled334.m

为配套王能超版本数值分析简明教程编写,方便学生教师演示!

x=[179   168   175   176   166   166   190   188   173   186   171   168   172   165   161   170   168   170   171   178   178   170   171   173   168   175   156 158   158   161   159   162   159   156   163   160   158   164   162   158   150   164   156   166   174   163   162   160   160   154];
y=[67.0000   50.0000   63.0000   62.0000   57.0000   47.0000   70.0000   80.0000   60.0000   80.0000   59.0000   60.0000   52.0000   59.0000   67.0000   57.0000 48.0000   65.0000   85.0000   90.0000   60.0000   55.0000   57.0000   59.0000   53.0000   60.0000   44.0000   51.0000   45.0000   50.0000   47.0000   48.0000 50.0000   43.0000   48.0000   43.0000   44.0000   55.0000   52.0000   51.0000   40.0000   48.5000   48.0000   55.0000   51.0000   49.0000   51.0000   47.0000 50.0000   45.0000];
s=log(y);
y1=0.019711*x+0.695157;
e=s-y1;
%plot(y1,e,'*');
hold on
for i=50:-1:2
    for j=1:i-1
        if e(j+1)<e(j)
            t=e(j+1);
            e(j+1)=e(j);
            e(j)=t;
        end
    end
end
e;
for i=1:50
    p(i)=(i-0.375)/(50+0.25);
    K(i)=i;
end
K;
p;
S=1-p;
z=[-2.2400   -1.8500   -1.6200   -1.4600   -1.3300   -1.2200   -1.1200   -1.0300   -0.9500   -0.8700   -0.8000   -0.7300   -0.6700   -0.6100   -0.5500   -0.4900  -0.4400   -0.3800   -0.3300   -0.2800   -0.2300   -0.1700   -0.1200   -0.0700   -0.0200 0.0200    0.0700    0.1200    0.1700    0.2300    0.2800    0.3300    0.3800    0.4400    0.4900    0.5500    0.6100    0.6700    0.7300    0.8000    0.8700 0.9500    1.0300    1.1200    1.2200    1.3300    1.4600    1.6200    1.8500    2.2400];
%plot(z,e,'*g');
fprintf('%13.8f',max1(y1))
%y1
%ismember(y1,max1(y1))
h=(4.44024700-3.65180700)/10;
s;
l=1;
k=zeros(1,10);
for i=1:10
    for j=1:50
        if s(j)>=3.65180700+(i-1)*h & s(j)<=3.65180700+i*h
            %a=s(j);
            k(i)=k(i)+1;
            A(l)=s(j);
            l=l+1;
        end
    end
end
%ismember(A,4.0073);
A
K1=0;
for i=1:10
   K1=K1+k(i);
end
K1;
k
k(1)=k(1)+1;
k(10)=k(10)+1;
for i=1:10
    f(i)=k(i)/50;
end
f
h1=(e(50)-e(1))/10;
for i=1:10
    e1(i)=e(1)+(i-1)*h1;
end
e1
%plot(f,e1,'*')
bar(e1,f);
%stem(f,e1)