实验五 函数文件答案.docVIP

实验五 函数文件 1. 定义一个函数文件，求给定复数的指数、对数、正弦、余弦，并在命令文件中调用该函数文件。 函数文件： function [a,b,c,d] fun1 x a exp x ; b log x ; c sin x ; d cos x ; 命令文件： x input x ? ; [a,b,c,d] fun1 x 2. 一物理系统可用下列方程组来表示： 从键盘输入和的值，求和的值。其中g取9.8，输入时以角度为单位。 要求：定义一个求解线性方程组的函数文件，然后在命令文件中调用该函数文件。 函数文件： function [a1,a2,N1,N2] fun2 m1,m2,theta g 9.8; A [m1*cos theta -m1 -sin theta 0;m1*sin theta 0 cos theta 0;0 m2 -sin theta 0;0 0 -cos theta 1]; B [0 m1*g 0 m2*g]; X inv A *B; a1 X 1 ; a2 X 2 ; N1 X 3 ; N2 X 4 ; 命令文件： m1 input m1 ? ; m2 input m2 ? ; theta input theta ? ; [a1,a2,N1,N2] fun2 m1,m2,theta 3. 一个自然数是素数，且它的数字位置经过任意对换后仍为素数，则称是绝对素数。例如13是绝对素数。试求所有两位绝对素数。 要求：定义一个判断素数的函数文件。 函数文件： function y fun3 x p isprime x ; if p 1 y x else y []; end 命令文件： for n 10:99 m fun3 n ; a fix m/10 ; b rem m,10 ; x 10*b+a; y fun3 x ; end 4. 设，编写一个MATLAB函数文件fx.m，使得调用时，x可用矩阵代入，得出的为同阶矩阵。 函数文件： function y fx x y 1./ x-2 .^2+0.1 +1./ x-3 .^2+0.01 ; 命令文件： x input x ? ; y fx x 5. 已知 （1）当时，求y的值。 （2）当时，求y的值。 1 函数文件： function y fun5 n y n+10*log n^2+5 ; 命令窗口中调用： y fun5 40 / fun5 30 +fun5 20 2 函数文件： function y fun6 n y 0; for m 1:n y y+m* m+1 ; end 命令窗口中调用： y fun6 40 / fun6 30 +fun6 20 数学软件课程实验报告 1