第2章 MATLAB语言程序设计实例解析.pptVIP

【例2-1】隐函数的求导。 对于隐函数求偏导问题时，diff()函数并不能直接用于求解。考虑到隐函数偏导数求解的一般方法：假设隐函数的数学表达式为 则可以通过隐函数求出自变量之间的偏导数。具体可以用下面的公式 求出 ： 上式的MATLAB实现即为： F=-diff(f,xj)/diff(f,xi) 【例2-2】参数方程求导。 假设参数方程 ，则 可以由递推公式求出： 根据上述递推公式编写求解参数方程高阶导数的通用程序diff_para.m。 【例2-3】计算二重积分 ，其中 Step 1：首先绘制出其积分区域 fplot([x^2,2-x^2],[-2 2 -0.5 2.5]) % 绘制曲线图 text([-1.4,-1.4],[2,0.05],{\leftarrowy=x^2,\leftarrowy=2-x^2}) % 添加标注 A=double(solve(x^2-(2-x^2))) % 将符号型结果转换成浮点型 hold on xn=linspace(A(1),A(2),100); fill([xn,fliplr(xn)],[xn.^2,fliplr(2-xn.^2)],‘b’) % 填充积分区域 Step 2：计算积分值 syms x y; I=int(int(x^2,y,x^2,2-x^2),x,min(A),max(A)) % 嵌套调用求解重积分 运行结果： I =8/15 【例2-4】设 ，试求解 。 解：分析题意可知，首先应将u和v看成x,y的函数，故应先求解方程组得到u和v的表达式，再进行题目的相关问题的求解。 syms x y u v [u,v]=solve(x*u+y*v=0,y*u+x*v=1,u,v) % 求解多项式方程组 dudxdy=diff(diff(u,x),y) % 求二阶偏导数 输出结果： u =-y/(x^2-y^2) v =1/(x^2-y^2)*x dudxdy =2/(x^2-y^2)^2*x+8*y^2/(x^2-y^2)^3*x 【例2-5】if语句实例。 从键盘中输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出。 解：该问题是一个典型的条件判断语句形式，编写如下语句即可。 c=input(请输入一个字符:,s); % 提示输入字符串 if c=A c=Z disp(char(c+a-A)); % 输出c对应的小写字符 elseif c=a c=z disp(char(abs(c-a+A))); % 输出c对应的大写字符 elseif c=0 c=9 disp(abs(c-0)); % 输出c对应的数字 else disp(c); % 输出c end 运行结果： 请输入一个字符:A a 请输入一个字符:a A 请输入一个字符:9 9 请输入一个字符:? ? 【例2-6】switch语句实例（原题见书本【例2-40】）。 解：首先我们要提取出开关表达式，这里若直接选用price作为开关表达式，则表达式i则不好表示，因为price的取值理论上是无穷的，考虑到将price化为整数，这样图2-9中的表达式i也就不难表示了。下面给出具体的程序： price=input(请输入商品价格：); switch fix(price) case num2cell(0:199) % 价格小于200 rate=0; case num2cell(200:499) % 价格大于等于200但小于500 rate=3/100; case num2cell(500:999) % 价格大于等于500但小于1000 rate=5/100; case num2cell(1000:2499) % 价格大于等于1000但小于2500 rate=8/100; case num2cell(2500:4999) % 价格大于等于2500但小于5000 rate=10/100; otherwise % 价格大于等于5000 rate=14/100; end sale_price=price*(1-rate) % 输出商品实际销售价格 运行结果： 请输入商品价格：1225 sale_price =1127 【例2-7】循环语句实例。 求