matlab实验4MATLAB基本编程

实验4 MATLAB基本编程

实验目的

熟悉MATLAB的脚本编写

熟悉MATLAB的函数编写

熟悉MATLAB的变量使用

熟悉MATLAB的程序控制结构

实验内容

分别选用if或switch结构实现下述函数表示。

-1 x≤-a

f(x)= x/a -a

1 x≥a

根据ex=1+x+x2/2!+…+xn/n!+…近似计算指数，当与指数函数的误差小于0.01时停止，分别用for和while结构实现。

记录上题的调试过程。

刚开始对于指数直接用e，后来改成exp

提示用户输入1或2，如输入1时，执行第一题脚本；如输入2时，提示用户输入随机数序列长度，然后执行第二题的函数。

从键盘输入数值，迭代计算xn+1=3/(xn+2)。

思考题

创建脚本实现随机系列的各元素由大到小排列，其中随机数服从U(-5,9)的均匀分布，并且序列长度为10。

创建函数实现指定长度的随机数序列的各元素由大到小的排列，其中随机数服从N(3,9)的高斯分布。

分别用for和while结构实现下述函数计算。sin(x)-cos(x)+sin(2x)+cos(2x)+……+sin(nx)+(-1)ncos(nx)

(2)

分别if 或switch结构实现下列函数表示

function y=function1(x,a)

if (x<=-a)

y=-1;

elseif (x-a)

y=x/a;

else

y=1;

end

function output=function2(x,a)

switch x&a

case x<=(-a)

output=-1;

case x>(-a)&x

output=x/a;

otherwise

output=1;

end

function y=function3(x)

y=0;

for i=0:inf

y=y+(x^i/factorial(i));

if(exp(x)-y<0.01)

break

end

end

function y=function4(x)

y=0;

i=0;

while(exp(x)-y>0.01)

y=y+(x^i/factorial(i));

i=i+1;

end

function n=function5()

r=input('please input 1 or 2:' );

switch r

case 1

disp(function1(1,2));

case 2

x=input('input a number:');

disp(function3(x));

end

function y=function6()

n=input('input a number :');

y=1.5;

if n==0

y=1.5;

else

for x=1:1:n

y=3/(y+2);

end

end

function y=function7

n=input('请输入n的值:');

a=normrnd(3,9,1,n);

y=fliplr(sort(a))

function y=function8

n=input('n=');

x=input('x=');

y=0;

for i=1:n

s=sin(i*x)+(-1)^i*cos(i*x);

y=y+s;

end

disp(y);

function function9

x=input('x=');

n=1;

y=0;

while n<=10

s=sin(n*x)+(-1)^n*cos(n*x);

y=y+s;

n=n+1;

end

disp(y);

function y=function10

A=[1,2,3;0,1,2;0,0,1];

n=input('n=');

t=input('t=');

y=0;

for i=1:1:n

x=A^i*exp(A^i*t);

y=y+x;

end

disp(y);

function y=function11

A=[1,2,3;0,1,2;0,0,1];

n=input('n=');

t=input('t=');

y=0;

while(n)

x=A^