matlab实验4MATLAB基本编程.doc

《matlab实验4MATLAB基本编程.doc》由会员分享，可在线阅读，更多相关《matlab实验4MATLAB基本编程.doc（21页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-datematlab实验4MATLAB基本编程实验1 MATLAB基本功能实验4 MATLAB基本编程一、 实验目的1. 熟悉MATLAB的脚本编写2. 熟悉MATLAB的函数编写3. 熟悉MATLAB的变量使用4. 熟悉MATLAB的程序控制结构二、 实验内容1. 分别选用if或switch结构实现下述函数表示。 -1 x-a f(x)= x/a -axa 1 xa2. 根

2、据ex=1+x+x2/2!+xn/n!+近似计算指数，当与指数函数的误差小于0.01时停止，分别用for和while结构实现。3. 记录上题的调试过程。刚开始对于指数直接用e，后来改成exp4. 提示用户输入1或2，如输入1时，执行第一题脚本；如输入2时，提示用户输入随机数序列长度，然后执行第二题的函数。5. 从键盘输入数值，迭代计算xn+1=3/(xn+2)。三、 思考题1. 创建脚本实现随机系列的各元素由大到小排列，其中随机数服从U(-5,9)的均匀分布，并且序列长度为10。2. 创建函数实现指定长度的随机数序列的各元素由大到小的排列，其中随机数服从N(3,9)的高斯分布。3. 分别用fo

3、r和while结构实现下述函数计算。(1) 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;endfunction output=function2(x,a)switch x&a case x(-a)&xa output=x/a; otherwise output=1;endfunction y=function3(x)y=0;for i=0:inf y=

4、y+(xi/factorial(i); if(exp(x)-y0.01) y=y+(xi/factorial(i); i=i+1;endfunction 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);endfunction y=function6()n=input(input a number :);y=1.5;if n=0 y=1.5;else for x=1:1:n y=3/

5、(y+2); endendfunction y=function7n=input(请输入n的值:);a=normrnd(3,9,1,n);y=fliplr(sort(a)function y=function8n=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 function9x=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);funct

6、ion y=function10A=1,2,3;0,1,2;0,0,1; n=input(n=); t=input(t=); y=0;for i=1:1:n x=Ai*exp(Ai*t); y=y+x;end disp(y);function y=function11A=1,2,3;0,1,2;0,0,1; n=input(n=); t=input(t=); y=0;while(n) x=An*exp(An*t); y=y+x; n=n-1;endy=y+exp(A*t);disp(y);function p=f(x,y)x=input(请输入x的值:)y=input(请输入y的值:)if y=1 p=sin(x);elseif y=2 p=cos(x);else p=sin(x)*cos(x); endfunction f=function13(x,y) switch y case 1 f=sin(x); case 2 f=cos(x); otherwise f=sin(x)*cos(x); end脚本a=unifrnd(-5,9,1,10);y=sort (a);y=rot90(y)-