MATLAB实验报告
学院：光电学院
班级：073-1
姓名：***
学号：************
实验五 函数文件
1.定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数文件。

程序设计：
function [e ln s c]=num(x) e=exp(x) ln=log(x) s=sin(x) c=cos(x) end
运行结果： >> num(5i) e =
0.2837 - 0.9589i ln =
1.6094 + 1.5708i s =
0 +74.2032i c =
74.2099 ans =
0.2837 - 0.9589i
2.一物理系统可用下列方程组来表示: ⎥⎥⎥
⎥
⎥⎦⎤
⎢⎢⎢⎢⎢⎣⎡=
⎥⎥⎥⎥⎥⎦
⎤⎢⎢
⎢⎢⎢
⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢
⎣⎡----g g m m N N a a m m m m 2121212111001cos 00
0sin 00cos 0sin 0sin cos θ
θθθθθ
从键盘输入
m
1
、
m
2
和θ的值，求
N a a 121、、和
N
2
的值。

其中g 取9.8，输入θ时以角度为单位。

程序设计： 函数文件in.m:
function [a1,a2,N1,N2]=in(m1,m2,t) g=9.8;
A=[m1*cos(t) -m1 -sin(t) 0;m1*sin(t) 0 cos(t) 0;0 m2 -sin(t) 0;0 0 -cos(t) 1]; C=[0;m1*g;0;m2*g]; B=inv(A)*C;
a1=B(1);
a2=B(2); N1=B(3); N2=B(4); end
调用in.m 的命令文件： >> m1=1;m2=2;t=30*pi/180; >> [a1,a2,N1,N2]=in(m1,m2,t) 运行结果： a1 =
6.5333 a2 =
1.8860 N1 =
7.5440 N2 =
26.1333 4.设 f(x)=
01
.01
1
.01
)
3()
2(4
2
+++--x x ，
编写一个MATLAB 函数文件fx.m ，使得调用f(x)时，x 可用矩阵代入，得出的f(x)为同阶矩阵。

程序设计： 函数文件fx.m:
function A=fx(x)
A=1./((x-2).^2+0.1)+1./(((x-3).^4)+0.01) end
调用fx.m 的命令文件：
>> A=fx([1 2;2 3;4 3]) 运行结果: A =
0.9716 10.9901 10.9901 100.9091 1.2340 100.9091 5.已知y=
)
20()30()
40(f f f +
（1）当f(n)=n+10ln(n 2+5)时，求y 的值。

程序设计：
函数文件fn.m:
function x=fn(n)
x=n+10*log(n^2+5)
end
调用fn.m的命令文件：
y=fn(40)/(fn(30)+fn(20)) 运行结果:
y =
0.6390
（2）当f(n)=1⨯2+2⨯3+3⨯4+…+n⨯(n+1)时，求y的值。

程序设计：
函数文件fn.m:
function x=fn(n)
if n<=0
x=0;
else
x=n*(n+1)+fn(n-1);
end
end
调用fn.m的命令文件：
y=fn(40)/(fn(30)+fn(20)) 运行结果:
y =
1.7662。