Question

用matlab或lingo解决矩阵的问题

用matlab或lingo解决矩阵，哦！~~

Answer 1

x=[7.0 10.5 13 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106.5 111.5 118 123.5 136.5 142 146 150 157 158];
y1=[44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68];
y2=[44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68];
m=7.0:0.001:158;
t1=interp1(x,y1,m,'spline');
plot(m,t1);
hold on
t2=interp1(x,y2,m,'spline');
plot(m,t2);
由此程序得到一个图形如何求面积 ？
%【1】

x=[7.0 10.5 13 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106.5 111.5 118 123.5 136.5 142 146 150 157 158];
y1=[44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68];
y2=[44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68];
m=7.0:0.001:158;
t1=interp1(x,y1,m,'spline');
plot(m,t1);
hold on
t2=interp1(x,y2,m,'spline');
plot(m,t2);
%面积
A=0.001*sum(t2-t1)

A = 8.5998e+003

%【2】可以用半径为1的圆的面积进行检验

x=-1:0.001:1;
y1=-sqrt(1-x.^2);
y2=sqrt(1-x.^2);
plot(x,y1,x,y2)
axis equal
aera=0.001*sum(y2-y1)

结果：
aera = 3.1416