15年13级matlab实验报告(9)

《MATLAB程序设计》实验报告

m1=1; for i=1:10 m1=m1*i; if(m1>10100)

disp(['n=',num2str(i)]);

disp([num2str(i),'!=',num2str(m1)]); return; end end

>> factorial n=8 8!=40320 n=8 8!=40320

11、下列程序用来判断一个人的体温是否处于危险状态。调试程序是否正确，如果程序错误指出错误在哪里？并写出正确答案。

temp=input(‘请输入体温: temp = ’); if temp < 36.5 disp(‘体温偏低’); elseif temp > 36.5

%错误，这个条件划分错误，会导致下面的elseif无法执行

disp(‘体温正常’);

elseif temp > 38.0

disp(‘体温偏高！’);

elseif temp > 39

disp(‘体温高！！’);

end

错误改正之后： 15/11/2

temp=input('请输入体温: temp = '); if temp < 36.5

共64 页，第 41页

《MATLAB程序设计》实验报告

disp('体温偏低');

elseif temp <37.5 %将体温正常的条件满足区域定义在（36.5,37.5） disp('体温正常'); elseif temp <39.0 disp('体温偏高！'); elseif temp > 39 disp('体温高！！'); end

12、使用嵌套for循环命令创建下列矩阵：

?5?1?A??0??0??0源程序及实验结果： forC.m 15-11-2

%嵌套for循环创建一个矩阵

clear A for i=1:5 for j=1:5 if i==j

A(i,j)=5;%对角线赋值为5 elseif j==i+1||j==i-1 A(i,j)=1; else A(i,j)=0; end end end

A

1000?5100??1510?

?0151?0015?? >> forC A =

共64 页，第 42页

《MATLAB程序设计》实验报告

5 1 0 0 0 1 5 1 0 0 0 1 5 1 0 0 0 1 5 1 0 0 0 1 5

实验四 MATLAB绘图功能

一、实验目的

1、掌握MATLAB二维图形绘制命令及其图形控制； 2、熟悉MATLAB三维图形绘制命令及其图形控制； 3、熟悉特殊二维图形、三维图形的绘制方法。 二、实验内容

1、在同一坐标系下绘制下面三个函数在t?[0,40]的图象，要求三种曲线采用不同颜色以及不同的线型，同时对每条曲线利用函数进行标注，并添加图例。

y1?ty2?ty3?4?e?0.1tsin(t)源程序： 15-11-23 %画曲线图 %作者

t=0:0.1:40;

y1=t;

y2=sqrt(t);

y3=4*pi*exp(-0.1*t).*sin(t);

plot(t,y1,'.r',t,y2,'-b',t,y3,'--k') text(10,10,'y=t');

text(25,6,'y=sqrt(t)');

gtext('y=4*pi.*exp(-0.1*t).*sin(t)');

legend('y=t','y=sqrt(t)','y=4*pi.*exp(-0.1*t).*sin(t)')

共64 页，第 43页

《MATLAB程序设计》实验报告

2、编写程序，选择合适的步距，绘制下面函数在区间[-6,6]中的图象，并对分段的曲线进行标注，同时添加x轴和y轴的说明。

?sinx,?y(x)??x,??x?6,?源程序： x=-6:0.1:6; y=[]; for x0=x if x0<=0

y=[y,sin(x0)]; elseif x0<=3 y=[y,x0]; else y=[y,-x0+6]; end end

plot(x,y)

title('分段曲线函数'); xlabel('variable X'); ylabel('variable Y'); gtext('y=sin(x)'); text(2,2.1,'y=x'); gtext('y=-x+6') 实验结果：

x?00?x?3 x?3共64 页，第 44页

《MATLAB程序设计》实验报告

3、利用subplot函数在同一绘图窗口中用不同颜色和线型绘制以下两个函数在t?[-2?,2?]范围内的图象。

y1?20.5t源程序：

y2?2e?0.2t

clear

t=-2*pi:0.1*pi:2*pi; y1=2.^0.5*abs(t); y2=2*exp(-0.2*t); subplot(1,2,1)

plot(t,y1);title('y1=2^0.5*abs(t)'); subplot(1,2,2);

plot(t,y2);title('y2=2*exp(-0.2*t)'); 实验结果：

共64 页，第 45页