MATLAB程序设计-.pptVIP

例7-5 计算二重定积分（1）建立一个函数文件fxy.m：function f=fxy(x,y)global ki;ki=ki+1; %ki用于统计被积函数的调用次数f=exp(-x.^2/2).*sin(x.^2+y);（2）调用dblquad函数求解。global ki;ki=0;I=dblquad(@fxy,-2,2,-1,1) kiI = 1.5745ki = 1038 如果使用inline函数，则命令如下：f=inline(exp(-x.^2/2).*sin(x.^2+y),x,y);I=dblquad(f,-2,2,-1,1)I = 1.5745 7.1 求定积分 例7-6 计算三重定积分命令如下：fxyz=inline(4*x.*z.*exp(-z.*z.*y-x.*x),x,y,z);triplequad(fxyz,0,pi,0,pi,0,1,1e-7)ans= 1.7328 补充示例 第六章: 实验内容：第1题，第4题 思考练习：第4题 第七章: 实验内容：第1题 思考练习：第3题，第4题 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 MATLAB基础 《通信电路原理》--无九 《通信电路原理》--无九 《通信电路原理》--无九 《通信电路原理》--无九 《通信电路原理》--无九 《通信电路原理》--无九 《通信电路原理》--无九 《通信电路原理》--无九 《通信电路原理》--无九 《通信电路原理》--无九 《通信电路原理》--无九 《通信电路原理》--无九 MATLAB程序设计 实验7 MATLAB方程求解与数值微积分 Dr Wang.Buyun - Lecture Notes * * dsolve y=dsolve(eq1,eq2, ... ,cond1,cond2, ... ,v) 其中 y 为输出的解， eq1、eq2、. . . 为微分方程，cond1、cond2、...为初值条件， v 为自变量 例 1：求微分方程 的通解，并验证。 y=dsolve(Dy+2*x*y=x*exp(-x^2),x) syms x; diff(y)+2*x*y - x*exp(-x^2) 微分方程求解 微分方程求解 * 例 2：求微分方程 满足初值条件  的特解，并画出解函数的图形。 y=dsolve(x*Dy+y-exp(x)=0, y(1)=2*exp(1), x) ezplot(y); 微分方程求解 * 例3：求微分方程组 在初值条件 下的特解，并画出解函数的图形。 [x,y]=dsolve(Dx+5*x+y=exp(t),Dy-x-3*y=0, ... x(0)=1, y(0)=0, t) ezplot(x,y,[0,1.3]); 注：解微分方程组时，如果所给的输出个数与方程个数相同，则方程组的解 按词典顺序 输出。 其它运算 * 反函数 finverse(f,v)：求 f 关于指定变量 v 的反函数 finverse(f)：求 f 关于默认变量的反函数 syms x t; f=x^2+2*t; g1=finverse(f,x) g2=finverse(f,t) 例：计算函数 的反函数 6.4 最优化问题求解 例6-13 求f(x)=x3-2x-5在[0,5]内的最小值点。 (1) 建立函数文件mymin.m。 function fx=mymin(x) fx=x.^3-2*x-5; (2) 调用fmin函数求最小值点。 x = fminbnd(mymin,0,5); x= 0.8165 例：求函数 在