《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页