数学实验之1.MATLAB软件简介.ppt

MATLAB 软件 例 作出由MATLAB的函数peaks产生的二元函数的曲面及其等值线图。 键入：[X,Y,Z]=peaks(30); surf(X,Y,Z); figure(2); contour(X,Y,Z,16); figure(3); contour3(X,Y,Z,16); 3）等高线 打开另一个图形窗口 返 回 5. M文件的编写 M文件 命令文件 函数文件 四种循环和选择控制结构 for循环 while循环 if-else-end结构 switch-case-end结构 返 回 M文件 在命令窗口进入程序编辑器输入、保存M文件，然后在命令窗口键入你建立的 M文件的文件名即可逐次执行其中的命令。 在此空白框键入M文件 保存按钮 命令M文件：一些命令行的组合。 函数M文件：定义函数的程序，具有特定格式。 例： x=[2 3]; f=cost(x) function f=cost(x) f=100*(x(2)-x(1)^2)^2… -(1-x(1))^2; M文件 函数M文件的形式为 function [输出变量列表]=函数名（输入变量列表） 函数体语句； M文件 返 回 逻辑运算符： （与运算）， |（或运算），~（非运算） 关系运算符： , , =, =, = =（等于），~=（不等于） 关系运算比较两个数值，当给出的关系式成立时，结果为1（表示真），否则为0（表示假）。 循环选择控制结构 返 回 如： s=0;n=1 while s=100 s=s+n^2; n=n+1; end n,s While循环的一般形式为： while 关系表达式 {statments} end 关系表达式 statments End的后续语句 Yes No while循环 i0=1-exp(-1)；n=0;in=i0; hold on plot(n,in,.) while (in=exp(-1)/(n+1))… (in=t) x=in; n=n+1；in=1-n*x plot(n,in,.) t= 1/(n+1); end hold off while循环 返 回 如, for k=1:4 x(k)=1/k; end format rat %设置输出格式为有理数 x For循环的一般形式 for x=x1:step:xn {statments} end x=x1 x=xn statments x=x+step End的后续语句 Yes No for循环 function f=jiecheng(n) f=1; for i=1:n f=f*i; end for循环 建立一个函数M文件求自然数n的阶乘 ? 注：循环可以嵌套，即循环体内又包含另一个完整的循环结构，内嵌的循环中还可以嵌套循环，形成多层循环。 for循环 返 回 最简单的形式： if 关系表达式 {statments} end 关系表达式 statments End的后续语句 Yes No if-else-end结构 当有三个或更多的选择时： if 关系表达式(1) {statments}(1) elseif 关系表达式(2) {statments}(2) ... elseif 关系表达式(n) {statments}(n) else {statments}(n+1) end 例如，可用以下程序得到符号函数。 function y=SIGN(x) if x0 y=-1; elseif x==0 y=0; else y=1; end 已知任意两个多项式（不一定同阶）的系数，求两个多项式的和。 if-else-end结构 ? if-else-end结构 function f=sumpoly(p1,p2) n=abs(length(p1)-length(p2)); if length(p1)length(p2) p12=[p1,zeros(1,n)]； p22=p2; elseif length(p1)length(p2) p22=[p2,zeros(1,n)]; p12=p1; end f=p12+p22; if-else-end结构 返 回 switch lower(NAME) case {zhanghua, lijiang}, disp( He comes from China.) case pet