用MATLAB2010編寫的程序
第壹題:
clear;clc;?
x0=rand(10,1);
y0=rand(10,1);
figure
plot(x0,y0,'o')
hold?on
p=polyfit(x0,y0,2);
x1=0:0.01:1;
y1=polyval(p,x1);
plot(x1,y1,'-r')
axis([0,1,0,1])
axis?equal
第二題:
clear;clc
x0=rand(10,1);
y0=rand(10,1);
z0=rand(10,1);
figure
plot3(x0,y0,z0,'o')
hold?on
n=length(z0);
X=[ones(n,1)?x0?y0];
b=X\z0;
x2=0:0.01:1;y2=x2;
[X,Y]=meshgrid(x2,y2);
z2=b(1)+b(2)*X+b(3)*Y;
mesh(X,Y,z2)
axis([0,1,0,1,0,1])
axis?equal