MATLAB实验五 函数文件
发布时间：2010-04-15 来源：文档文库
小
中
大
字号：
光电学院 073-1 刘颖
200713503117
MATLAB实验报告
学院：光电学院 班级：073-1 姓名：刘颖
学号：200713503117
1
光电学院 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.一物理系统可用下列方程组来表示: m1cosm1sin00m1sin0cos0a100a2m1g0N10g1m2N2a1=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=1(x220.11，0.01(x34m02sincos
从键盘输入编写一个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=m1、m2和θ的值，求a1、a2、N1和N的值。其中g取9.8，输2入θ时以角度为单位。 程序设计： 函数文件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;
2 f(40
f(30f(20（1）当f(n=n+10ln(n2+5时，求y的值。
光电学院 073-1 刘颖
200713503117 程序设计： 函数文件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 