实验1 Matlab程序设计基础

实验1 Matlab 程序设计基础
一、实验名称：Matlab 程序设计基础.
二、实验目的：掌握Matlab 语言，会用Matlab 进行一些简单的程序设计.
三、实验要求：
1、对以下问题，编写M 文件：
（1）用起泡法对10个数由小到大排序.即将相邻两个数比较，将小的调到前头.
（2）有一个4×5矩阵，编程求其最大值及其所处的位置.
（3）编程求∑=20
1!n n .
（4）一球从100米高度自由落下，每次落地后反跳回原高度的一半，再落下.求它在第10次落地时，共经过多少米？第10次反弹有多高？
（5）有一函数f(x,y)=x^2+sin(x*y)+2*y ，写一程序，输入自变量的值，输出函数值.
2、在同一平面中的两个窗口分别画出心形线和马鞍面.
要求：
（1）在图形上加格栅、图例和标注.
（2）定制坐标.
（3）以不同角度观察马鞍面.
3、请编写绘制以下图形的Matlab 命令，并展示绘得的图形：
（1）4,12222=+=+y x y x 分别是椭圆14/2
2=+y x 的内切圆和外切圆； （2）指数函数和对数函数x y ln =的图像关于直线x y =对称.
四、实验内容：
1、1）源程序：
function p=paixu(a)
a=[1,2,5,6,8,7,9,4,3,0];
for i=9:-1:1
for j=1:i
if a(j)>a(j+1)
t=a(j);
a(j)=a(j+1);
a(j+1)=t;
end
end
end
a
调试结果：a =
0 1 2 3 4 5 6 7 8 9
2）源程序：
function j=juzhen(a)
a=rand(4,5);
m=a(1);ind=[1 1];
for i=1:size(a,1)
for j=1:size(a,2)
if m<a(i,j)
m=a(i,j);
ind=[i j];
end
end
end
m
ind
a(ind(1),ind(2))
调试结果：m =0.9502 ind =3 5
ans =0.9502 ans =5
3）源程序：function r = fun(n)
n=20
result = 1;
r = 0;
for ii = 1 : n
result = result * ii;
re= r + result;
end
调试结果：ans = 2.5613e+018
4）源程序：function q=xiaoqiu(s)
s=zeros(1,10);
s(1)=100;
for i=1:9
s(i+1)=s(i)/2;
end
L=s(1)+2*sum(s(2:10))
h=s(10)
调试结果：L = 299.6094 h =0.1953
5）源程序：function f=fun(x,y)
f=x^2+sin(x*y)+2*y;
调试结果：在command window界面中敲入fun(1,2)，
得到ans = 5.9093
2、源程序：
clear
ax1=subplot(1,2,1);
t=0:.01*pi:2*pi;
f1=polar(t,1+cos(t));
title('心形线'); legend(f1,'r=1+cost')
ax2=subplot(1,2,2);
x1=linspace(-20,20,100);
y1=linspace(-20,20,100);
[xx1,yy1]=meshgrid(x1,y1);
zz1=xx1.*yy1;
f2=surf(xx1,yy1,zz1);
title('马鞍面'); legend(f2,'z=x^2-y^2','Location','NorthWest') xlabel('x ')
ylabel('y ')
zlabel('z')
axis([-20 20 -20 20 -400 400])
setAllowAxesRotate(rotate3d,ax1,false);
grid on
调试结果：
3、1）源程序：
t=linspace（-pi, pi,120);
x1=2*sin(t);
y1=2*cos(t);
x2=sin(t);
y2=cos(t);
x3=2*sin(t);
y3=cos(t);
plot(x1,y1,x2,y2,x3,y3)
调试结果：
2）源程序：x1=linspace(0,5,10);x2=linspace(0,5,10);
y1=exp(x1);
y2=log(x2);
y3=x1;
plot(x1,y1,x2,y2,x1,y3)
axis([0,5,0,5])
调试结果：。