MATLAB实验五 函数文件

光电学院 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.一物理系统可用下列方程组来表示: ?m1cos???m1sin??0??0?m1?sin?0cos?0??a1??0??????0??a2??m1g??0??N1??0??????g1?????m2N???2?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)=

1(x?2)2??0.11，

?0.01(x?3)4m02?sin??cos?

从键盘输入

编写一个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(30)?f(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)=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

3

本文来源：https://www.bwwdw.com/article/xnr8.html