大学数学实验MATLAB简介2.ppt

（3）求导数 diff(S,v) 求表达式S对变量v的一阶导数. diff(S,v,n) 求表达式S对变量v的n阶导数. 例如：设A= ,求键入命令: syms a b x; A= [1/(1+a),(b+x)/cos(x);1,exp(x^2)]; diff(A,x) ans = [0,1/cos(x)+(b+x)/cos(x)^2*sin(x)] [0, 2*x*exp(x^2)] 例：求y=sinx+ex的三阶导数,键入命令: diff(sin(x)+x*exp(x),3) 得 ans = -cos(x)+3*exp(x)+x*exp(x) 例：求y=(lnx)x的导数.可键入命令： p=(log(x))^x; p1=diff(p,x) 得：p1 = log(x)^x*(log(log(x))+1/log(x)) 例：求y=xf(x2)的导数.可键入命令： p=x*f(x^2);p1=diff(p,x) 得：p1 = f(x^2)+2*x^2*D(f)(x^2) （4）求积分 int(P) 对表达式P进行不定积分. int(P,v) 以v为积分变量对P进行不定积分. int(P,v,a,b) 以v为积分变量,以a为下限,b为上限对P进行定积分. 例：求 ,可键入 int(-2*x/(1+x^2)^2) 得 ans = 1/(1+x^2) 例：求 ,可键入 int(x*log(1+x),0,1) 得 ans = 1/4 例：求 可键入： int(2*x,sin(t),log(t)) 得： ans = log(t)^2-sin(t)^2 对（符号）矩阵积分 （5）符号函数作图 ezplot(f) 在默认区间[-2π,2π]绘制y=f(x)的函数图 ezplot(f,[a,b]) 在区间[a,b]绘制y=f(x)的函数图 ezplot(x,y,[tmin, tmax]) 绘制由参数议程x=x(t),y=y(t), tmin=t=tmax表示的曲线 ezpolar(f,[a,b]) 绘制由极坐标议程r=f(θ),a=θ=b表示的曲线 ezplot3,ezmesh,ezsurf,ezcontore类似。 A.6 程序设计—M文件的编写 6.1 M文件简介 Matlab的程序文件，以.m为扩展名 脚本式M文件 函数式M文件 1.脚本式M文件 脚本式M文件中的命令访问MATLAB工作空间中的所有变量，其结果也保留在工作空间 运行： 在MATLAB命令窗口中输入文件名，注意加入有哪些信誉好的足球投注网站路径（File-Set Path） 直接RUN运行M文件。 例1： x=4;y=6;z=2; items=x+y+z; const=x*25+y*22+z*99 average_cost=cost/items. 2.函数式M文件及其调用 function [输出变量列表]=函数名（输入变量列表） 函数体语句； 函数式M文件的文件名最好与其函数名相同。 函数式M文件中能访问输入变量和自身中变量，不能访问MATLAB工作空间中的变量。 例2： function f=fun(x) f=100*(x(2)-x(1)^2)-(1-x(1))^2; 调用： x=[2 3];f=fun(x) 例3.多个输出 function[F,G]=fun2(x) F=2*x(1)^2+2*x(2)^2-2*x(1)*x(2)-4*x(1)-6*x(2); G=[x(1)+5*x(2)-5 2*x(1)^2-x(2); -x(1) -x(2)]; 调用： x1=[4 5];[F1,G1]=fun2(x1) 6.2 运算符 算术运算符: + - * / ^ 关系运算符: = = == ~= 逻辑运算符: | ~ 6.3 循环结构 1.for for x=表达式 语句; end 例子 for k=1:4 x(k)=1/k; end format rat x 2.while循环 while 关系表达式 语句; end 例子： n=0;EPS=1; while(1+EPS)1 EPS=EPS/2;n=n+1; end n,EPS=EPS*2 6.4 分支结构 1.if-else-end结构（elseif） if 关系表达式 语句 else 语句 end 2.switch-case-end结构 switch 表达式（标量或字符串） case 值1 语句体1; case {值2.1，值2.2，…} 语句体2; … otherwise 语句体n end 例A.7 操练 操练三 2）3