控制系统matlab仿真实验报告2

09自动化张家明080312009081
控制系统仿真实验报告
实验2:MATLAB的数学运算、M文件及外部数据操作
班级：09自动化报告人：张家明完成时间：2012-10-22
1. M文件和函数的创建
1）创建一个计算阶乘的函数
2）创建一个M文件，并用它调用1）所创建的函数进行阶乘计算
3）创建一个能读取外部数据的函数，并建立一个有数据的外部文件，使用该函数读取文件中的数据并求数据之和、平均值。

1).M函数：
function result=jc(n)
result=1;
for i=1:n
result=result*i;
end
result
2).创建一个M文件，在文件中输入：
result=jc(6)
运行结果：
result =
720
3).M函数：
function [data,average,total]=cal(x1)
Fid=fopen(x1,'rt');
if Fid==-1
disp('file does not exist !')
else
data=fscanf(Fid, '%f',inf);
fclose(Fid);
data
average=mean(data)
total=sum(data)
end
运行结果：
>> cal('data.dat')
09自动化张家明080312009081
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
average =
10.5000
total =
210
2. MATLAB的程序设计
function result=find_zero(m,n)
a=rand(m,n)<0.7
result=zeros(m,1);
for i=1:m
for j=1:n
if ~a(i,j)
result(i)=j;
break;
end
end
if result(i)==0
result(i)=Inf;
end
end
result
运行结果：
>> find_zero(3,4)
a =
0 1 1 1
1 0 1 1
1 0 0 0
result =
1
2
2
3. 编写函数：判断一个点与三角形的位置关系，能够给出点在三角形内部，在三角形边上，还是在三角形外部的信息。

(可选误差为1.0e-3)
function cal(x,a,b,c)
if nargin<4
disp('the number of input variable is less than 4')
else
length_1=sqrt((a(1)-b(1))^2+(a(2)-b(2))^2);%计算三角形三边长
length_2=sqrt((b(1)-c(1))^2+(b(2)-c(2))^2);
length_3=sqrt((a(1)-c(1))^2+(a(2)-c(2))^2);
p=(length_1+length_2+length_3)/2; %半周长
s=sqrt(p*(p-length_1)*(p-length_2)*(p-length_3));%计算三角形面积
%计算x,a,b三点组成的三角形面积
len1=sqrt((a(1)-b(1))^2+(a(2)-b(2))^2);
len2=sqrt((b(1)-x(1))^2+(b(2)-x(2))^2);
len3=sqrt((a(1)-x(1))^2+(a(2)-x(2))^2);
p1=(len1+len2+len3)/2;
s1=sqrt(p1*(p1-len1)*(p1-len2)*(p1-len3));
%计算x,a,c三点组成的三角形面积
len4=sqrt((a(1)-c(1))^2+(a(2)-c(2))^2);
len5=sqrt((a(1)-x(1))^2+(a(2)-x(2))^2);
len6=sqrt((c(1)-x(1))^2+(c(2)-x(2))^2);
p2=(len4+len5+len6)/2;
s2=sqrt(p2*(p2-len4)*(p2-len5)*(p2-len6));
%计算x,c,b三点组成的三角形面积
len7=sqrt((c(1)-b(1))^2+(c(2)-b(2))^2);
len8=sqrt((b(1)-x(1))^2+(b(2)-x(2))^2);
len9=sqrt((c(1)-x(1))^2+(c(2)-x(2))^2);
p3=(len7+len8+len9)/2;
s3=sqrt(p3*(p3-len7)*(p3-len8)*(p3-len9));
%判断
if (s1==0||s2==0||s3==0)&&(abs(s1+s2+s3-s)<=0.001) disp('the point is on the boarder of triangle') elseif abs(s1+s2+s3-s)<=0.001
disp('the Point is inside the triangle ') elseif abs(s1+s2+s3-s)>0.001
disp('the point is outside the triangle') end
end
运行结果：
>> a=[1 1];b=[2 1];c=[2 2];x=[1.8 1.1];
>> cal(x,a,b,c)
the Point is inside the triangle .。