实验五函数文件

1.1实验目的

1.2实验内容

1.3流程图

1.4程序清单

1.5运行结果及分析

1.6实验的收获与体会

1.1实验目的

理解函数文件的概念。
掌握matlab中定义和调用函数的方法。

1.2实验内容

1.3流程图

1.4程序清单

%第一题函数

function [ zhi,dui,sin1,cos1] = fu1( s )

%UNTITLED3 此处显示有关此函数的摘要

%   此处显示详细说明

zhi=exp(s);

dui=log(s);

sin1=sin(s);

cos1=cos(s);

end

%第二题函数

function [ a1,a2,N1,N2 ] = wlxt( m1,m2,sai)

%UNTITLED4 此处显示有关此函数的摘要

%   此处显示详细说明

g=9.8;

a=[m1*cos(sai),-m1,-sin(sai),0;m1*sin(sai),0,cos(sai),0;0,m2,-sin(sai),0;0,0,-cos(sai),1]^(-1)*[0;m1*g;0;m2*g];

a1=a(1);

a2=a(2);

N1=a(3);

N2=a(4);

end

%第三题函数

function a=sushu(b)

x=fix(b/10);

y=rem(b,10);

a=10*y+x;

sig=0;

for m=2:sqrt(b)

    if rem(b,m)==0

        sig=sig+1;

    end

end

for n=2:sqrt(a)

    if rem(a,n)==0

        sig=sig+1;

    end

end

if sig==0

    a=b;

else

    a=0;

end

%第四题函数

function f= fx( x )

%UNTITLED6 此处显示有关此函数的摘要

%   此处显示详细说明

f=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);

end

%第5题第1问函数

function fn=f51(n)

fn=n+10*log(n*n+5);

end

%第5题第2问函数

function fn=f52(n)

fn=0;

for k=1:n

    fn=fn+k*(k+1);

end

%% 实验五

clc

clear

%% 1

fprintf('\n-----------1------------\n');

s1=input('输入一个复数的实部:');

s2=input('输入一个复数的虚部:');

s=s1+s2*i;

[ zhi,dui,sin1,cos1] = fu1( s )

%fprintf('\n对于%f+%fi\n指数为:%f%fi，\n对数为:%f%fi，\n正弦为:%f%fi，\n余弦为:%f%fi\n',s1,s2,real(zhi),imag(zhi),real(dui),imag(dui),real(sin1),imag(sin1),real(cos1),imag(cos1));

%% 2

fprintf('\n-----------1.2------------\n');

m1=input('m1=');

m2=input('m2=');

sai=input('角度为：');

[ a1,a2,N1,N2 ]=wlxt(m1,m2,sai)

%% 3

fprintf('\n-----------3------------\n');

lwss=[];

for n=10:99

    [a]=sushu(n);

    if a>0

        lwss=[lwss,a];

    end

end

disp(['两位绝对素数有：',num2str(lwss)])

fprintf('\n\n');

%% 4

fprintf('\n------------4---------------\n');

x=rand(4,5)

f=fx(x)

%% 5

fprintf('\n-------------5--------------\n');

y1=f51(40)/(f51(30)+f51(20));

y2=f52(40)/(f52(30)+f52(20));

fprintf('\n-------------5.1--------------\n');

fprintf('\n      f(40)\ny = ————— = %f\n   f(30)+f(20)\n',y1);

fprintf('\n-------------5.2--------------\n');

fprintf('\n      f(40)\ny = ————— = %f\n   f(30)+f(20)\n',y2);