MATLAB程序设计

1、实验名称： matlab 程序设计一、实验目的与要求：matlab作为一种高级语言，不仅可以以命令方式完成操作，也具有数据结构、控制流、输入输出和面向对象编程的能力，本次实验通过熟悉和掌握m脚本文件和m函数文件的建立与使用、程序控制流程语句的使用，掌握基本的matlab 编程和调试能力。二、实验原理、内容和步骤1. 编写 m脚本文件， t 的范围为 0 ，2 ，步长取 0.05 ，计算函数y1=5e-2tsin(4t),y2=5e-2tcos(4t) 的值，并将变量t ，y1 和 y2 放在同一矩阵的a 的三行中。2fibonacci数列的各元素为：1,1,2,3,5,8,满足关系：f1=1；

2、f2=1；fn=fn-1+fn-2 用 m函数文件实现，数列的元素个数为输入变量。要求有完整的文件开头注释。3编写一个m函数文件，实现三类字符的统计。如下的两段文章，任选择一个，要求根据调用函数时传入的实参值的不同，对应的统计出英文大写字母的个数、英文小写字母个数或空格的个数。4编写一个m函数文件，以存储多个欲排序的数据向量为输入参数，用冒泡法排序后输出。5编写 m函数文件，实现分段绘制曲面，绘制每段曲面为一个子函数。要求在一个窗口的四个子窗口中分别用 plot3,mesh,surf,contour,四个绘图命令，绘制该曲面并观察区别。15457.0117575.015457.0),(5.17

3、5.375.065.175.375.0222222yxeyxeyxeyxzxxyxyxxy6根据输入参数个数实现当没有输入参数时，显示信息；当有一个参数时，则以该参数为边长绘制红色填充的正方形；当有两个参数，以两个参数为长和宽绘制绿色填充的矩形。三、结论分析与实验体会1代码： t=0:0.05*pi:2*pi; y1=5.*exp(-2).*t).*sin(4.*t); y2=5.*exp(-2).*t).*cos(4.*t); y=t;y1;y2 2m函数文件代码：function f=fibonacci(n) %fibonacci 数列%n 元素个数%f1=1 ；f2=1；fn=fn-1+

4、fn-2 函数元素满足的关系n=input(输入个数 ); f(1)=1;f(2)=1; for i=3:1:n 第二题结果示例： f(i)=f(i-1)+f(i-2); end 图 1 3m函数文件代码：function f(str,x); %big 大写字母数第三题结果示例：%small 小写字母数%space 空格数%str 字符串m=length(str); big=0;small=0;space=0; for i=1:1:m if str(i)= space=space+1; 图 2 elseif (str(i)=a) big=big+1; elseif (str(i)=a) sma

5、ll=small+1; end end if x=1 disp(big); elseif x=2 disp(space); elseif x=3 disp(small); end 4m函数文件代码： function t(str) %str 要排列的数字str=input(请输入数 ); 第四题结果示例：m=length(str); for i=1:m for j=2:m-i+1 if (str(j-1)=str(j) 图 3 temp=str(j-1); str(j-1)=str(j); str(j)=temp; end end end disp(str); 5m函数文件代码：functio

6、n g() x1=-2:0.1:2; x2=-2:0.1:2; xx1,xx2=meshgrid(x1,x2); n1,m1=size(xx1); z=zeros(n1,m1); for n=1:n1 for m=1:m1 if xx1(n,m)+xx2(n,m)1 z(n,m)=calp1(xx1(n,m),xx2(n,m); elseif (xx1(n,m)+xx2(n,m)-1)&(xx1(n,m)+xx2(n,m)=1) z(n,m)=calp2(xx1(n,m),xx2(n,m); else z(n,m)=calp3(xx1(n,m),xx2(n,m); 第五题结果示例： e

7、nd end end subplot (2,2,1) plot3(xx1,xx2,z);subplot (2,2,2) mesh (xx1,xx2,z); subplot (2,2,3) surf (xx1,xx2,z); subplot (2,2,4) contour (xx1,xx2,z); function z=calp1(x,y) z=0.5457*exp(-0.75*y2-3.75*x2-1.5*x); 图 5 function z=calp2(x,y) z=0.7575*exp(-y2-6*x2); function z=calp3(x,y) z=0.5457*exp(-0.75*y2-3.75*x2+1.5*x); 6m函数文件代码：function h(varargin) %nargin 输入变量个数%varargin 各变量元素内容if nargin=0 disp(无输入参数 ) elseif nargin=1 r=varargin1; rectangle(position,0 0 r r