MATLAB全部实验及答案(5)

利润（万元） 70 122 144 152 174 196 202 要求用多项式拟合法预测出该企业2011年和2012年的利润，并绘出拟合图。 解：x=2004:2010;

y=[70 122 144 152 174 196 202]; k=polyfit(x,y,1); %得到系数 x1=2004:2013; y1=polyval(k,x1);

plot(x,y,'-ro',x1,y1,'-b*');

title('某乡镇企业2004-2013年的生产利润') legend('04-10年的利润','04-13年的利润趋势线') xlabel('年份'); ylabel('利润');

text(2011,polyval(k,2011),'<---2011年利润') text(2012,polyval(k,2012),'<---2012年利润') grid on

4. 先使用命令load加载census.mat得到一组数据（cdate,pop），该数据是美国自1790到1990年（以10年为一单位）的总人口。利用回归法拟合，预测出美国在2010年的人口总数，写出拟合模型y?a0?a1x?a2x2(其中y为pop,x为cdate)，并绘出拟合图。 解：clear all,clc load('census.mat'); nian=[1990:2010]; p=polyfit(cdate,pop,2);

z=polyval(p,cdate); zz=polyval(p,nian); my=mean(pop); mz=mean(z);

R=sum((pop-my).*(z-mz))./sqrt(sum((pop-my).^2).*sum((z-mz).^2)) plot(cdate,pop,'*',cdate,z,'r',nian,zz,'b-')

xlabel('年 份(十年)','FontName','隶书','FontSize',16) ylabel('人口（百万）','FontName','隶书','FontSize',16)

title('第四题 美国自1790-1990年的总人口','FontName','隶书','FontSize',20)

text(1990, 248.73,'<----之后为拟合数据') legend('原数据','拟合曲线')

实验九 MATLAB绘图

实验内容

1、 二维图形的绘制

1)设计程序，在一个图形窗体中画出函数y=ex和y=logx的特性曲线。其中y=ex中x取值为（-2:0.1:2），y=logx中x取值为（0.1:0.1:5）。 要求：要有图例、标题（'二维图'）、坐标轴标签（'x轴数据'和'y轴数据'）；每条曲线的clm自己设置。

解：clear,clc

x1=-2:0.1:2; y1=exp(x1); x2=0.1:0.1:5; y2=log(x2);

plot(x1,y1,'-r',x2,y2,'-b'); title('二维图');

legend('y=exp(x)','y=log(x)'); xlabel('x轴数据'); ylabel('y轴数据'); grid on;

2)绘制tan图形，其中x的取值为（0:pi/30:pi/2）。

要求：要有图例、标题、坐标轴标签（'x轴数据'和'y轴数据'）；曲线的自己设置。绘制出的图形应如下右图所示。

解

：x=0:pi/100:pi/2; y=tan(x);

plot(x,y,'-r+');

axis([0 pi/2 0 5]);

clm title('正切函数'); legend('y=tan(x)'); xlabel('x轴数据');

ylabel('y轴数据'); 3）绘制曲线 x=tcos(3t)

y=tsin2t 其中-π≤t≤π，步长取π/100。 要求：要有图例、标题、坐标轴标签；每条曲线的clm

自己设置。

解：clear,clc t=-pi:pi/100:pi; x=t.*cos(3*t); y=t.*sin(t).^2; plot(x,y); title('曲线图');

legend(strvcat('x=tcos(3t)','y=tsin^2t')); xlabel('x轴数据'); ylabel('y轴数据'); 2、 三维图形的绘制

分别用线框图和曲面图表现函数z=cos(x)sin(y)/y，其中[-1.5pi,1.5pi]，y=x，要求：要有标题、坐标轴标签

x的取值为

解：clear,clc

[x,y]=meshgrid(-1.5*pi:1.5*pi); z=cos(x).*sin(y)./y; subplot(1,2,1); surf(x,y,z)

title('三维曲线图01'); xlabel('x轴数据'); ylabel('y轴数据'); zlabel('z轴数据'); subplot(1,2,2); mesh(x,y,z);

title('三维曲线图02'); xlabel('x轴数据'); ylabel('y轴数据'); zlabel('z轴数据');