上机实验4

实验五编程

一、实验目的：熟悉M文件的建立和编辑；三种程序控制结构，命令文件和函数文件

二、实验内容：

1，定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数文件。程序设计：

function [E,L,S,C] = fushu( x )

E=exp(x);

L=log(x);

S=sin(x);

C=cos(x);

end

2，某公司员工的工资计算方法如下：

（1）工作时数超过120小时者，超过部分加发15%。

（2）工作时数低于60小时者，扣发700元。

（3）其余按每小时84元计发。

试编程按输入的工号和该号员工的工时数，计算应发工资。

function [gz] = suangongzi( gh ,gs )

if gs>120

gz=gs*84+(gs-120)*84*0.15;

else

if gs>60

gz=gs*84;

else

gz=gs*84-700;

end

end

end

3,已知⎪⎪⎩⎪⎪⎨⎧>+-=======---3

,23,12,01,1321321

n n n n f f f f f f f n n n n 求f f 1001~中： （1）最大值、最小值、各数之和。

（2）正数、零、负数的个数。

f(1)=1;

f(2)=0;

f(3)=1;

for n=4:100

f(n)=f(n-1)-2*f(n-2)+f(n-3);

end

max=max(f)

min=min(f)

sum=sum(f)

l=0;

s=0;

z=0;

for n=1:100

if f(n)==0

l=l+1;

else

if f(n)>0

z=z+1;

else

s=s+1;

end

end

end

s

l

z

4，对于任意给出激光器腔的结构：L, R 1, R 2, 判断腔的类型。

function []=qiangxing(L,R1,R2)

g1=1-L/R1;

g2=1-L/R2;

if 0

fprintf('稳腔\n');

else

if g1*g2==0

fprintf('临界腔\n');

else

if g1==1&&g2==1

fprintf('平行平面腔\n');

else

if g1==-1&&g2==-1

fprintf('对称共心腔\n');

else

fprintf(' 非稳腔\n');

end

end

end

end

end

5，利用幂级数展开式

+-+-=!

7!5!3sin 7

53x x x x x 计算sin x 的近似值（要求误差小于10-8）。

function s=jsin(x)

s=0;

format long

for i=1:20

s=x^(2*i-1)/factorial(2*i-1)*(-1)^(i-1)+s;

end

6，利用定积分的定义，计算定积分1

20

I x dx =⎰的近似值 I=0;

for x=0:0.00000001:1

I=I+x.^2.*0.00000001;

end

I7，设函数32()41f x x x =-+，求方程()0f x =在[0,1]上根的近似值

x=0:.0001:1;

y=abs(x.^3-4.*x.^2+1);

[m,n]=min(y);

x(n)

8，设银行利率为3.6，将10000元人民币存入银行，问多长时间会连本带利翻一番？

s=10000;

i=0;

while s<20000

s=s+s*0.036;

i=i+1;

end

i s=10000;

i=0;

while s<20000

s=s+s*0.036;

i=i+1;

end

i

9，在共焦腔中，任意傍轴光线在其中两次往返自行闭合。

syms R1 R2 L r1 s1 ;

A=1-2*L/R2;

B=2*L-2*L^2/R2;

C=-2/R1-2/R2+4*L/R1/R2;

D=-2*L/R1+(1-2*L/R1)*(1-2*L/R2);

T=[A B; C D];

R1=L;

R2=L;

M=eval(T)^2;

N=M*[r1;s1]

10，求1！+2！+3！+4！+5！的值

sum(factorial(1:5))