下载文档原格式
第2章MATLAB概论1、与其他计算机语言相比较，MA TLAB语言突出的特点是什么？答：起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富，可扩展性强.2、MA TLAB系统由那些部分组成？答：开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装.第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可.4、MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？答：在MATLAB操作桌面上有五个窗口，在每个窗口的右下角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口称为独立的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面称为独立窗口，在独立窗口的view菜单中选择Dock，菜单项就可以将独立的窗口重新防止的桌面上.5、如何启动M文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动.在命令窗口中键入edit命令时也可以启动M文件编辑/调试器.6、存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可.7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中.8、如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成.在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上.9、在MA TLAB中有几种获得帮助的途径？答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器.（2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息.（3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数.（4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数.注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息.第3章 MATLAB 数值运算3.1在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？ 答：在Command Window 窗口输入操作:>> a=[5 7 3;4 9 1]3.2有几种建立矩阵的方法？各有什么优点？ 答：（1）直接输入法，如a=[1 2 3 4]，优点是输入方法方便简捷；（2）通过M 文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改； （3）由函数建立，如y=sin(x)，可以由MATLAB 的内部函数建立一些特殊矩阵； （4）通过数据文件建立，该方法可以调用由其他软件产生数据.3.3在进行算术运算时，数组运算和矩阵运算各有什么要求？答：进行数组运算的两个数组必须有相同的尺寸.进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数.渗釤呛俨匀谔鱉调硯錦。
MATLAB基础实验指导书漳州师范学院物电系2010年10月目录实验一MATLAB环境的熟悉与基本运算 (2)实验二MATLAB数值运算 (8)实验三MATLAB语言的程序设计 (12)实验四MATLAB的图形绘制 (16)实验五采用SIMULINK的系统仿真 (20)实验六MATLAB在电路中的应用 (25)实验七MATLAB在信号与系统中的应用 (27)实验八MATLAB在控制理论中的应用 (29)实验一 MATLAB环境的熟悉与基本运算一、实验目的：1．熟悉MATLAB开发环境2．掌握矩阵、变量、表达式的各种基本运算二、实验基本知识：1.熟悉MATLAB环境:MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器文件和搜索路径浏览器。
2.掌握MATLAB常用命令3.MATLAB变量与运算符变量命名规则如下：（1）变量名可以由英语字母、数字和下划线组成（2）变量名应以英文字母开头（3）长度不大于31个（4）区分大小写MATLAB中设置了一些特殊的变量与常量，列于下表。
MATLAB运算符，通过下面几个表来说明MATLAB的各种常用运算符表2 MATLAB算术运算符表3 MATLAB关系运算符表4 MATLAB逻辑运算符表5 MATLAB特殊运算4.MATLAB的一维、二维数组的寻访表6 子数组访问与赋值常用的相关指令格式5.MATLAB的基本运算表7 两种运算指令形式和实质内涵的异同表6.MATLAB的常用函数表8 标准数组生成函数表9 数组操作函数三、实验内容1、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推）2、学习使用clc、clear，观察command window、command history和workspace等窗口的变化结果。
3、初步程序的编写练习，新建M-file，保存（自己设定文件名，例如exerc1、exerc2、exerc3……），学习使用MATLAB的基本运算符、数组寻访指令、标准数组生成函数和数组操作函数。
1.判断三角形clca=input('');b=input('');c=input('');if a+b>c&&b+c>a&&a+c>aswitch 1case a==b==cdisp(1)case a==b|b==c|a==cdisp(2)otherwise disp(3)endelsedisp(0)end2、费波纳切数列：1、1、2、3、5、8、13……，其递推公式为：a(1)=a(2)=1，a(n)=a(n-1)+a(n-2)(n>=3)。
编程实现求n项的费波纳切数列。
要求：用M语言脚本文件实现；对输入的错误数据格式（如字符串、小数等）能够给出错误提示。
（函数）function y=feibonaqie(n)%求费波纳切数列a(1)=1;a(2)=1;%n=input('请输入项数n：');if(n==1)a(n)=1;elseif(n==2)a(n)=1;elseif(n>=3)i=3;while(i<=n)a(i)=a(i-1)+a(i-2);i=i+1;endendy=a(n);（脚本）clc;a(1)=1;n=input('请输入项数n ：');if(n==1)a(n)=1;elseif(n==2)a(n)=1;elseif(n>=3)i=3;while(i<=n)a(i)=a(i-1)+a(i-2);i=i+1;endendsprintf('a(%u)=%u',n,a(n))9X9乘法表for i=1:9;A=[];for j=1:i;A=[A,sprintf('%d ×%d=%2d ',j,i,i*j)];enddisp(A)end3、古代的数学家祖冲之利用正多边形逼近的割圆法计算了常数pi 的值，请编程实现，要求精确到小数点后的第10位。
第4章 MATLAB程序流程控制习题4一、选择题1．下列关于脚本文件和函数文件的描述中不正确的是（)。
A A．函数文件可以在命令行窗口直接运行B．去掉函数文件第一行的定义行可转变成脚本文件C．脚本文件可以调用函数文件D．函数文件中的第一行必须以function开始2．下列程序的输出结果是( ）。
Dy=10；if y==10y=20；elseif y>0y=30enddisp（y)A．1 B．30 C．10 D．203．有以下语句：a=eye（5)；for n=a(2:end,：)for循环的循环次数是（）。
CA．3 B．4 C．5 D．104．设有程序段k=10;whilekk=k—1end则下面描述中正确的是（）。
AA．while循环执行10次B．循环是无限循环C．循环体语句一次也不执行D．循环体语句执行一次5．有以下程序段：x=reshape（1:12,3,4）;m=0;n=0；for k=1：4if x（：，k)<=6m=m+1;elsen=n+1；endend则m和n的值分别是（）。
CA．6 6 B．2 1 C．2 2 D．1 26．调用函数时，如果函数文件名与函数名不一致,则使用（)。
A A．函数文件名B．函数名C．函数文件名或函数名均可D．＠函数名7．如果有函数声明行为“function [x,y,z］=f1（a，b,c）”,则下述函数调用格式中错误的是（）。
BA．x=f1（a,b,c) B．[x,y，z，w]=f1(a,b,c)C．[x，b，z］=f1（a，y,c) D．[a，b]=f1(x,y，z）8．执行语句“fn=@(x）10＊x;"，则fn是（）。
AA．匿名函数B．函数句柄C．字符串D．普通函数9．执行下列语句后，变量A的值是( ）。
D>> f=@（x，y) log（exp（x+y）);〉〉 A=f(22,3);A．22，3B．22 C．3 D．2510．程序调试时用于设置断点的函数是（）。
第4章 MATLAB程序流程控制习题4一、选择题1．下列关于脚本文件和函数文件的描述中不正确的是（）。
AA．函数文件可以在命令行窗口直接运行B．去掉函数文件第一行的定义行可转变成脚本文件C．脚本文件可以调用函数文件D．函数文件中的第一行必须以function开始2．下列程序的输出结果是（）。
Dy=10;if y==10y=20;elseif y>0y=30enddisp(y)A．1 B．30 C．10 D．203．有以下语句：a=eye(5);for n=a(2:end,:)for循环的循环次数是（）。
CA．3 B．4 C．5 D．104．设有程序段k=10;while kk=k-1end则下面描述中正确的是（）。
AA．while循环执行10次B．循环是无限循环C．循环体语句一次也不执行D．循环体语句执行一次5．有以下程序段：x=reshape(1:12,3,4);m=0;n=0;for k=1:4if x(:,k)<=6m=m+1;elsen=n+1;endend则m和n的值分别是（）。
CA．6 6 B．2 1 C．2 2 D．1 26．调用函数时，如果函数文件名与函数名不一致，则使用（）。
A A．函数文件名B．函数名C．函数文件名或函数名均可D．@函数名7．如果有函数声明行为“function [x,y,z]=f1(a,b,c)”，则下述函数调用格式中错误的是（）。
BA．x=f1(a,b,c) B．[x,y,z,w]=f1(a,b,c)C．[x,b,z]=f1(a,y,c) D．[a,b]=f1(x,y,z)8．执行语句“fn=@(x) 10*x;”，则fn是（）。
AA．匿名函数B．函数句柄C．字符串D．普通函数9．执行下列语句后，变量A的值是（）。
D>> f=@(x,y) log(exp(x+y));>> A=f(22,3);A．22,3B．22 C．3 D．2510．程序调试时用于设置断点的函数是（）。
%Exerc ise 1(1)r oots([1 11])%Exer cise1(2)roots([3 0 -4 0 2 -1])%Exerc ise 1(3)p=zero s(1,24);p([1 17 1822])=[5 -6 8 -5];ro ots(p)%E xerci se 1(4)p1=[2 3];p2=conv(p1,p1);p3=co nv(p1, p2);p3(end)=p3(en d)-4; %原p3最后一个分量-4r oots(p3)%Exer cise2fun=inli ne('x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x');fzer o(fun,2)%Exer cise3fun=inli ne('x^4-2^x');fplot(fun,[-2 2]);gr id on;fze ro(fu n,-1),fzer o(fun,1),f minbn d(fun,0.5,1.5)%Exe rcise 4fu n=inl ine('x*sin(1/x)','x');fp lot(f un, [-0.10.1]);x=z eros(1,10);fori=1:10, x(i)=fz ero(f un,(i-0.5)*0.01);end;x=[x,-x]%Ex ercis e 5f un=in line('[9*x(1)^2+36*x(2)^2+4*x(3)^2-36;x(1)^2-2*x(2)^2-20*x(3);16*x(1)-x(1)^3-2*x(2)^2-16*x(3)^2]','x');[a,b,c]=fso lve(f un,[0 0 0])%E xerci se 6fun=@(x)[x(1)-0.7*si n(x(1))-0.2*cos(x(2)),x(2)-0.7*cos(x(1))+0.2*sin(x(2))]; [a,b,c]=fsolv e(fun,[0.5 0.5])%E xerci se 7clear; clo se; t=0:p i/100:2*pi;x1=2+sqr t(5)*cos(t); y1=3-2*x1+sq rt(5)*sin(t);x2=3+s qrt(2)*cos(t);y2=6*sin(t);pl ot(x1,y1,x2,y2); gri d on; %作图发现4个解的大致位置，然后分别求解y1=fsolv e('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[1.5,2])y2=fsolv e('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[1.8,-2])y3=fsol ve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[3.5,-5])y4=fso lve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[4,-4])%Exerc ise 8(1)c lear;fun=inlin e('x.^2.*(x.^2-x-2)');fp lot(f un,[-2 2]);grid on;%作图观察x(1)=-2;x(2)=fminb nd(fu n,-1,-0.5);x(4)=fmi nbnd(fun,1,2);fun2=inlin e('-x.^2.*(x.^2-x-2)');x(3)=f minbn d(fun2,-0.5,0.5);x(5)=2feval(fun,x)%答案: 以上x(2)(4)是局部极小，x(1)(3)(5)是局部极大，从最后一句知道x(1)全局最大， x(4)最小。
《MATLAB 原理及应用》实验报告三.课后练习题答案1.为 ⎪⎩⎪⎨⎧<-=>+=01001x x x xx x y 编写赋值程序。
程序如下：①建立如下的M 文件：x=input('x=');%让用户通过键盘输入数值、字符串或表达式if x>0y=x+1;elseif x==0y=x;else x<0y=x-1;e nd程序执行结果如下>> kh1 %在当前工作目录下，文件名为“kh.1.m ” x=1>> yy =22.使用for ... end循环的array向量编程求出1+3+5...+100 的值程序如下：sum=0;>> for k=1:2:100sum=sum+k;end>> sumsum =25003.计算1+3+5...+100 的值，当和大于1000时终止计算。
程序如下：sum=0;for m=1:2:100; %建立1 3 5….100的向量if sum<=1000 %如果sum小于1000则可以继续加sum=sum+m; %累加elsebreak; %若sum的结果不符合条件就跳出整个循环endend结果为：sum =1024k =653.1计算从1开始多少个自然数之和超过100。
程序如下：>> sum=0;n=0;>> while sum<=100n=n+1;sum=sum+n;end结果为：n =14sum =1054.求1！+2！+3！+……+8！的值程序如下：n=1;sum=1;for m=2:8; %循环7次使得得到各次阶乘n=n*m;sum=sum+n; %累加end结果为:sum =462335.写程序，判断一年是否为闰年，符合下面两条件之一：（1990~2014）A、能被4整除，不能被100整除B、能被400整除程序如下count=0;for y=1990:2014;if((rem(y,4)==0&rem(y,100)~=0)|(rem(y,4)==0&rem(y,400)~=0));count=count+1;endend结果为:count =5。
欢迎共阅习题二1.如何理解“矩阵是MATLAB 最基本的数据对象”？ 答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
(2) 删除矩阵A 的第7号元素。
答：A(7)=[](3) 将矩阵A 的每个元素值加30。
答：A=A+30;(4) 求矩阵A 的大小和维数。
答：size(A);ndims(A);(5) 将向量 t 的0元素用机器零来代替。
答：t(find(t==0))=eps; (6) 将含有12个元素的向量 x 转换成34⨯矩阵。
答：reshape(x,3,4); (7) 求一个字符串的ASCII码。
答：abs(‘123’); 或E 。
答：B=A(1:3,:); C=A(:,1:2); D=A(2:4,3:4);E=B*C;(2) 分别求E<D 、E&D 、E|D 、~E|~D 和find(A>=10&A<25)。
答：E<D=010001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，E&D=110111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，11⎡⎤⎢⎥答：student(1).id='0001';student(1).name='Tom'; student(1).major='computer';student(1).grade=[89,78,67,90,86,85]; 8.建立单元矩阵B 并回答有关问题。
B{1,1}=1;B{1,2}='Brenden';B{2,1}=reshape(1:9,3,3); B{2,2}={12,34,2;54,21,3;4 ,23,67};(1)size(B)和ndims(B)的值分别是多少？答：size(B) 的值为2, 2。
ndims(B) 的值为2。
(2)B(2)和B(4)的值分别是(2)建立5×6随机矩阵A，其元素为[100,200]范围内的随机整数。
第4章M A T L A B程序流程控制-习题答案第4章MATLAB程序流程控制习题4一、选择题1．下列关于脚本文件和函数文件的描述中不正确的是（）。
AA．函数文件可以在命令行窗口直接运行B．去掉函数文件第一行的定义行可转变成脚本文件C．脚本文件可以调用函数文件D．函数文件中的第一行必须以function开始2．下列程序的输出结果是（）。
Dy=10;ify==10y=20;elseify>0y=30enddisp(y)A．1B．30C．10D．203．有以下语句：a=eye(5);forn=a(2:end,:)for循环的循环次数是（）。
CA．3B．4C．5D．104．设有程序段k=10;whilekk=k-1end则下面描述中正确的是（）。
AA．while循环执行10次B．循环是无限循环C．循环体语句一次也不执行D．循环体语句执行一次5．有以下程序段：x=reshape(1:12,3,4);m=0;n=0; fork=1:4ifx(:,k)<=6m=m+1; elsen=n+1; end end则m 和n 的值分别是（）。
CA ．66B ．21C ．22D ．126．调用函数时，如果函数文件名与函数名不一致用（）。
A A ．函数文件名B ．函数名 C ．函数文件名或函数名均可D ．@函数名 7．如果有函数声明行为“f unction[x,y,z]=f1(a,b,c)”，则下述函数调用格式中错误的是 （）。
B A ．x=f1(a,b,c)B ．[x,y,z,w]=f1(a,b,c) C ．[x,b,z]=f1(a,y,c)D ．[a,b]=f1(x,y,z) 8．执行语句“f n=@(x)10*x;”，则fn 是（）。
A A ．匿名函数B ．函数句柄C ．字符串D ．普通函数 9．执行下列语句后，变量A 的值是（）。
D >>f=@(x,y)log(exp(x+y)); >>A=f(22,3); A ．22,3B ．22C ．3D ．25 10．程序调试时用于设置断点的函数是（）。
A A ．dbstopB ．dbclearC ．dbcontD ．dbstack 二、填空题 1．将有关M A T L A B 命令编成程序存储在展名为.m 的文件中，该文件称 为。
M 文件2．有语句“f ork=[12;34]”引导的循环结构，其循环体执行的次数为。
1 3．M A T L A B 中用于控制不确定重复次数的循环中需要终止该循用的语句为。
while ⋯end ，break 4．函数文件由语句引导。
在函数定义时，函数数。
在调用函数时，输入输出为参数。
function ，形式，形参，实际，实参 5．在MATLAB 中，函数文件中的变量是变量。
定义变量是函数间 传递信息的一种手段用命令定义。
局部，全局，global 6．应用程序的错误有两类，一类是错误，另一类是运行时的错误，即23错误。
MATLAB程序调试方法有两种，一是利用进行程序调试，二是利用进行程序调试。
语法，逻辑，调试函数，调试工具三、应用题1．写出下列程序的输出结果。
s=0;a=[12,13,14;15,16,17;18,19,20;21,22,23];fork=aforj=1:4ifrem(k(j),2)~=0s=s+k(j);endendends2．分别用if语句和switch语句实现以下计算，其中a、b、c的值从键盘输入。
2axbxc,3.5 5.5语句:a=input('a=');b=input('b=');c=input('c=');x=input('x=');<=x&x<y=a*x^2+b*x+c;<=x&x<y=a*(sin(b))^c+x;<=x&x<y=log(abs(b+c/x));endyswitch语句：a=input('a=');b=input('b=');3c=input('c=');x=input('x=');switchfix(x*10)casenum2cell(5:14)y=a*x^2+b*x+c;casenum2cell(15:34)y=a*(sin(b))^c+x;casenum2cell(35:54)y=log(abs(b+c/x));otherwisedisp('nsrdxcw')endy3．产生20个两位随机整数，输出其中小于平均值的偶数。
x=fix(10+rand(20)*89);a=sum(x)/20;fori=1:20ifx(i)<aifrem(x(i),2)==0disp(x(i))endendend4．输入20个数，求其中最大数和最小数。
要求分别用循环结构和调用MATLAB的max函数、min函数来实现。
循环结构:fora=1:20A(a)=input('请输入20个数：');endAmi=A(1);ma=A(1);fori=1:20ifmi>A(i)mi=A(i);endifma<A(i)ma=A(i);45end end disp(ma) disp(mi) max 、min 函数: fora=1:20A(a)=input('请输入20个数：'); end Ami=min(A),ma=max(A) 5．已知s=1+2+22+23+⋯+263分别用循环结用MATLAB 的sum 函数求s 的值。
循环结构: s=0;forx=0:63 s=2^x+s; enddisp(s) sum 函数: i=0:63; f=2.^i; sum(f)6．当n 分别取100、1000、10000时，求下列各式的值：1111n1(1)1(1)(ln2) 234n111π(2)1()357411111 (3)()n4166443224466(2n)(2n)π (4)()133557(2n1)(2n1)2 要求分别用循环结构和(1): 循环结构: n =i n p u t5s=0;fora=1:ns=s+(-1)^(a+1)/a;ends向量运算:n=input('请输入n的值：');a=1:n;f=(-1).^(a+1)./a;s=sum(f)(2):循环结构:n=input('请输入n的值：');s=0;fora=1:n;s=s+(-1)^(a+1)/(2*a-1);ends向量运算:n=input('请输入n的值：');a=1:n;f=(-1).^(a+1)./(2*a-1);s=sum(f)(3):循环结构:n=input('请输入n的值：');s=0;fora=1:ns=s+(1/4)^a;ends向量运算:n=input('请输入n的值：');a=1:n;f=(1/4).^a;s=sum(f)(4):67循环结构:n=input('请输入n的值：');s=1;fora=1:nf=(2*a)^2/((2*a-1)*(2*a+1));s=s*f;ends向量运算:n=input('请输入n的值：');a=1:n;f=(2*a).^2./((2*a-1).*(2*a+1));s=prod(f)7．编写一个函数文件，求小于任意自然数n的斐波那契（Fibnacci）数列各项。
斐波那契数列定义如下：1f1f12fff(n2)nn1n2functionf=Fibnacci(n)fori=1:nifi==1f=1;elseifi==2f=1;elsef=Fibnacci(i-1)+Fibnacci(i-2);endendn=input('n=');f=Fibnacci(n)8．编写一个函数文件，用于求两个矩阵的乘积和点乘，然后在脚本文件中调用该函数。
function[C,D]=CJ(A,B)C=A*B;D=A.*B;7A=[123;456;789]; B=[111;111;111]; [C,D]=CJ(A,B)9．先用函数的递归调用定义一个函数文件求nmi，然后调用该函数文件求i1100501012kk。
kk1k1kfunctionf=xt4_10(m,n)ifn<1f=0;elsef=n^m+xt4_10(m,n-1);endxt4_10(1,100)+xt4_10(2,50)+xt4_10(-1,10)10．写出下列程序的输出结果。
脚本文件：globalxx=1:2:5;y=2:2:6;sub(y);xy函数文件：functionfun=sub(z)globalxz=3*x;x=x+z;8。
