MATLAB实验五函数文件.docVIP

MATLAB实验报告 学院：光电学院 班级：073-1 姓名：刘颖 学号：200713503117 实验五 函数文件 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取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)=，编写一个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= （1）当f(n)=n+10ln(+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)=12+23+34+…+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 光电学院 073-1 刘颖 200713503117 3