curve02.m

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

% curve02.m
x=[160.0 160.0 160.0 160.0 160.0 160.0 160.1 160.1 160.1 160.1...
      160.1 160.2 160.2 160.2 160.3 160.3 160.3 160.3 160.4 160.4...
      160.4 160.4 160.3 160.3 160.2 160.0 159.9 159.7 159.4 159.0...
      158.6 158.1 157.5 156.8 156.0 155.0 153.9 152.6 151.2 149.6...
      147.7 145.8 143.5 141.0 138.2 135.2 132.0 128.4 124.6 120.5...
      116.0 111.3 106.4 100.9 95.0 89.3 82.9 76.4 69.4 62.4 ...
      55.0 47.3 39.4 31.5 23.7 14.3 5.50 -2.5 -11.0 -20.5...
      -28.7 -38.0 -46.5 -55.4 -64.3 -72.7 -81.2 -89.0 -97.5 -104.0...
      -112.6 -119.7 -126.7 -133.0 -139.7 -145.5 -151.1 -156.1 -160.7 -164.9...
      -168.4 -172.5 -175.4 -178.3 -181.1 -183.2 -184.6 -184.8 -185.0 -185.0];
y=[0.0 0.0 0.0 0.1 0.2 0.3 0.4 0.6 0.8 1.2 1.6 1.9 2.4 3.1 3.7...
      4.5 5.4 6.4 7.5 8.8 10.1 11.7 13.3 15.1 16.9 19.0 21.2 23.6 26.1 28.8...
      31.6 34.5 37.7 41.0 44.4 48.0 51.7 55.7 59.6 63.7 68.1 72.2 76.8 81.3 86.0...
      90.7 95.3 100.1 104.8 109.4 114.1 118.7 123.1 127.6 132.0 136.0 140.0 143.7...
      147.3 150.5 153.5 156.2 158.6 160.6 162.2 163.7 164.6 165.1 165.3 165.0 ...
      164.3 163.1 161.6 159.6 157.1 154.3 151.0 147.5 143.2 139.0 134.1 129.1 ...
      123.6 118.1 111.7 105.0 99.1 92.7 86.1 79.5 73.5 65.7 59.3 52.3 44.6 37.7 ...
      32.5 31.7 31.0 31.0];
z=[0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 ...
   31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 ...
   59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 ...
   87 88 89 90 91 92 93 94 95 96 97 98 99];  

n=3;

axy=polyfit(x,y,n)
y0=polyval(axy,x);
plot(x,y0,'.r','markersize',20)
hold on
plot(x,y,'-b','markersize',40)
xlabel('x');
ylabel('y');
title('X-Y curve');

figure(2)
plot(y,z,'-b','markersize',40)
xlabel('y');
ylabel('z');
title('Y-Z curve');
axis([0 180 -20 120]);

figure(3)
plot(z,x,'-b','markersize',40)
xlabel('z');
ylabel('x');
title('z-x curve');
axis([-20 120 -200 200 ]);