《实验5_函数文件_参考解答》由会员分享，可在线阅读，更多相关《实验5_函数文件_参考解答（4页珍藏版）》请在金锄头文库上搜索。
1、高教社刘卫国高教社刘卫国MATLAB 程序设计与应用 （第二版）实验参考解答程序设计与应用 （第二版）实验参考解答 实验五：函数文件实验五：函数文件 实验五实验五 函数文件函数文件 教材教材 P365 一、实验目的一、实验目的 1. 理解函数文件的概念。 2. 掌握定义和调用 MATLAB 函数的方法。 二、实验内容二、实验内容 1. 定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用 该函数文件。 答：程序设计：答：程序设计： 函数文件：函数文件：complexnum.m function expval,logval,sinval,cosval=complexnum(x
2、) %COMPLEXNUM compute complex number value expval=exp(x); logval=log(x); sinval=sin(x); cosval=cos(x); 命令文件：命令文件：exp0501.m clear all; close all; clc; x=input(请输入一个复数：); expval,logval,sinval,cosval=complexnum(x); disp(给定复数的指数是：, num2str(expval) disp(给定复数的对数是：, num2str(logval) disp(给定复数的正弦是：, num2str(
3、sinval) disp(给定复数的余弦是：, num2str(cosval) 运行结果：运行结果： 请输入一个复数：1+2i 给定复数的指数是：-1.1312+2.4717i 给定复数的对数是：0.80472+1.1071i 给定复数的正弦是：3.1658+1.9596i 给定复数的余弦是：2.0327-3.0519i 2. 一物理系统可用下列方程组来表示： 11121112220cossin0sin0cos0 00sin000cos1mm m gmNm Nm g = 从键盘输入、和1m2m的值，求1、2、和的值。其中 g 取 9.8，输入1N2N时以角度为单位。 要求：定义一个求解线性方程
4、组=AXB的函数文件，然后在命令文件中调用该函数文MATLAB 语言语言 WangChengyou 山东大学威海分校山东大学威海分校 1高教社刘卫国高教社刘卫国MATLAB 程序设计与应用 （第二版）实验参考解答程序设计与应用 （第二版）实验参考解答 实验五：函数文件实验五：函数文件 件。 答：程序设计：答：程序设计： 函数文件：函数文件：equationsolution.m function J=equationsolution(m1,m2,theta) H=m1*cos(theta*pi/180), -m1, -sin(theta*pi/180), 0; m1*sin(theta*pi/1
5、80), 0, cos(theta*pi/180), 0; 0, m2, -sin(theta*pi/180), 0; 0, 0, -cos(theta*pi/180), 1; K=0;m1*9.8;0;m2*9.8; J=inv(H)*K; 命令文件：命令文件：exp0502.m clear all; close all; clc; m1=input(m1=); m2=input(m2=); theta=input(theta=); J=equationsolution(m1, m2, theta); disp(方程组的解alpha1：, num2str(J(1) disp(方程组的解alp
6、ha2：, num2str(J(2) disp(方程组的解N1：, num2str(J(3) disp(方程组的解N2：, num2str(J(4) 运行结果：运行结果： m1=2 m2=3 theta=45 方程组的解alpha1：8.6621 方程组的解alpha2：2.45 方程组的解N1：10.3945 方程组的解N2：36.75 3. 一个自然数是素数，且它的数字位置经过任意对换后仍为素数，则称是绝对素数， 例如 13 是绝对素数。试求所有两位绝对素数。 要求：定义一个判断素数的函数文件。 答：程序设计：答：程序设计： 函数文件：函数文件：prime.m function a=pri
7、me(b) %absolute prime x=fix(b/10); y=rem(b,10); a=10*y+x; signab=0; for m=2:sqrt(b) if rem(b,m)=0 MATLAB 语言语言 WangChengyou 山东大学威海分校山东大学威海分校 2高教社刘卫国高教社刘卫国MATLAB 程序设计与应用 （第二版）实验参考解答程序设计与应用 （第二版）实验参考解答 实验五：函数文件实验五：函数文件 signab=signab+1; end end for n=2:sqrt(a) if rem(a,n)=0 signab=signab+1; end end if s
8、ignab=0 a=b; else a=0; %这里赋值0，表示不是素数 end 命令文件：命令文件：exp0503.m clear all; close all; clc; twoprime=; for n=10:99 a=prime(n); if a0 twoprime=twoprime,a; end end disp(所有的两位绝对素数是：, num2str(twoprime) 运行结果：运行结果： 所有的两位绝对素数是：11 13 17 31 37 71 73 79 97 4. 设2411( )(2)0.1(3)0.01f xxx=+，编写一个 MATLAB 函数文件 fx.m，使得调
9、用( )f x时，x可用矩阵代入，得出的( )f x为同阶矩阵。 答：程序设计：答：程序设计： 函数文件：函数文件：fx.m function f=fx(x) f=1./(x-2).2+0.1)+1./(x-3).4+0.01); 命令文件：命令文件：exp0504.m clear all; close all; clc; x=1,2;3+4i,5 f=fx(x) 运行结果：运行结果： x = 1.0000 2.0000 3.0000 + 4.0000i 5.0000 MATLAB 语言语言 WangChengyou 山东大学威海分校山东大学威海分校 3高教社刘卫国高教社刘卫国MATLAB 程
10、序设计与应用 （第二版）实验参考解答程序设计与应用 （第二版）实验参考解答 实验五：函数文件实验五：函数文件 f = 0.9716 10.9901 -0.0482 - 0.0280i 0.1724 5. 已知(40) (30)(20)fyff=+(1) 当时，求y的值。 2( )10ln(5)f nnn=+ (2) 当时，求y的值。 ( )1 22 33 4(1)f nnn= + + + +? 答：答： (1) 程序设计：程序设计： 函数文件：函数文件：f.m function fn=f(n) fn=n+10*log(n*n+5); 命令文件：命令文件：exp05051.m clear all; close all; clc; y=f(40)/(f(30)+f(20) 运行结果：运行结果： y = 0.6390 (2) 程序设计：程序设计： 函数文件：函数文件：f1.m function fn=f1(n) fn=0; for k=1:n fn=fn+k*(k+1); end 命令文件：命令文件：exp05052.m clear all; close all; clc; y=f1(40)/(f1(30)+f1(20) 运行结果：运行结果： y = 1.7662 MATLAB 语言语言 WangChengyou 山东大学威海分校山东大学威海分校 4
