matlab实验五答案详解

实验内容：

1.函数文件：function [y]=mat1(x)

y=[exp(x),log(x),sin(x),cos(x)];

函数调用：[y]=mat1(1+i)

y =

1.4687 + 2.2874i 0.3466 + 0.7854i 1.2985 + 0.6350i 0.8337 - 0.98i

2.函数文件：function [a,b,N,M]=shiyanwu2(m,n,t)

A=[m*cos(t*pi/180),-m,-sin(t*pi/180),0;m*sin(t*pi/180),0,cos(t*pi/180),0;0,n,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1]; B=[0,9.8*m,0,9.8*n]; C=inv(A)*B'; a=C(1); b=C(2); N=C(3); M=C(4);

函数调用：m1=input('m1='); m2=input('m2='); theta=input('theta=');

[a1,a2,N1,N2]=shiyanwu2(m1,m2,theta) m1=54 m2=65 theta=56 a1 =

9.4680 a2 =

2.4025 N1 =

188.3684 N2 =

742.3342

3、函数文件：function flag=mat3(x)

flag=1;

for i=2:sqrt(x) if rem(x,i)==0 flag=0; break; end end

函数调用：for i=10:99 j=10*rem(i,10)+fix(i/10); if mat3(i)&mat3(j) disp(i) end end 11

13

17

31

37

71

73

79

97

4. 程序代码

fx.m

function f= fx(x)

%fx fx求算x矩阵下的f(x)的函数值 A=0.1+(x-2).^2; B=0.01+(x-3).^4; f=1./A+1./B; 命令

>> x=input('输入矩阵x='); f=fx(x) 结果 f =

0.0211 0.0745 0.0101 1.2340 0.1724 0.0012

0.0016 0.0005 0.0003 5. （1） 程序代码 f.m

function f=f(x)

f=x+10*log(x^2+5); 命令

n1=input('n1='); n2=input('n2='); n3=input('n3='); y1=f(n1); y2=f(n2); y3=f(n3); y=y1/(y2+y3) 结果 n1=40 n2=30 n3=20 y =

0.6390 （2） 程序代码 g.m

function s= g(n) for i=1:n g(i)=i*(i+1); end

s=sum(g); 命令

>> n1=input('n1='); n2=input('n2='); n3=input('n3='); y1=g(n1); y2=g(n2); y3=g(n3); y=y1/(y2+y3) 结果 n1=40 n2=30 n3=20 y =

1.7662