数学建模的资料(自己整理的).docxVIP

一、阶乘问题： 用 FOR 循环： clear n=1; s=0; for i=1:20 n=n*i; s=s+n; end  % 函数文件 factor.m function f=factor(n) if (n=1) f=1; else f=n*factor(n-1); end % main.m 用 while 循环：clear n=1; 用 while 循环：clear n=1; s=0; i=1; while i=20 n=n*i; s=s+n; i=i+1; %while 循环需要 end fprintf(the answer is %g \n,s) % clear; s=0; for i=1:5 s=s+factor(i); end fprintf( s=%g \n,s) fprintf(jiecheng of 20 de he is %g\n,s) 注：输出函数用 注：输出函数用 fprintf 用 disp 会出现错误。也可以直接调用阶乘函数 factorial(n):n 的阶乘 二：符号运算：在 MATLAB 中要进行符号的运算，必须定义符号对象。（syms alpha beta x y） MATLAB MATLAB 的符号运算，可以求方程的解、求导、定积分、还可以解决线性微分 方 程 （ eg:dsolve(‘D2y+D y+5=0’, 初始条件， 自变量)）。 符号对象的建立 数学实验 符号对象的建立 数学实验 符号对象的建立：sym 和 syms sym 函数用来建立单个符号变量，一般调用格式为： 符号变量 = sym(A) 参数 A 可以是一个数或数值矩阵，也可以是字符串 ?符号对象的建sy立m和：syms ?sym命s 令用来多建个立符号变量，一般调用格式为： sym符s号变1量符号变2量..符. 号变n量 例： a=sym(a) a 是符号变量 例：symsa b c b=sym(1/3) b 是符号常量 a=sy(ma); b=sy(mb;) c=sy(mc;) c=sym([1 ab; c d]) c 是符号矩阵 Matlab符号运算举例 数学实验求一元二次方a程x2+bx+ Matlab符号运算举例 数学实验 求一元二次方a程x2+bx+c = 0的根 solv(ea*x^2+b*x+c) 求的根f (x) = c(osx)2的一次导数 x=sy(mx); diff(cos(x)^2) 计算f (x) =x2 在区间[a, b] 上的定积分 symsa b x; int(x^2,a,b) 创建一个符号矩阵： s=sym([sin(x) 5*y]) s = [ sin(x), cos(y)] cos(y); tan(x) [ tan(x), 5*y] 符号表达式的替换 数学实验 subs 举例 数学实验 用给定的数据替换符号表达式中的指定的符号变量 subs(f,x,a) 用 a 替换字符函数 f 中的字符变量 x a 是可以是 数/数值变量/表达式或 字符变量/表达式 若 x 是一个由多个字符变量组成的数组或矩阵， 则 a 应该具有与 x 相同的形状的数组或矩阵。 例：指出下面各条语句的输出结果 f=sym(2*u); subs(f,u,2) f2=subs(f,u,u+2) a=3; subs(f2,u,a+2) subs(f2,u,a+2) syms x y f3=subs(f,u,x+y) subs(f3,[x,y],[1,2]) 下面的命令运行结果会是什么？  f=2*u ans=4 f2=2*(u+2) ans=14 ans=2*((a+2)+2) f3=2*x+2*y ans=6 cos(pi/6.0) ans = 0.8660 cos(pi\6.0) ans = -0.3326 注意：\和/的区别 注意：在绘制三维图形的时候，mesh(X,Y,Z,) 绘制由矩 阵 X,Y,Z 所确定 的曲面网格图。注意 x,y 这样写表示一个等宽的行矩阵。 注意矩阵的乘方最好用 x.^2; 矩阵的点除 也最好用 ./ ； x.^2 表示矩阵的每一个元素平方。而 x^2 表示 x*x;两者不同. 在画图中一般都是 用点乘  subs(f3,[x,y],[x+y,x+y]) 空间三维作图 数学实验 绘制由函数 z=z(x,y) 确定的曲面时，首先需产生一个网格矩阵，然后计算函数在各网格点上的值。 网格生成函数：meshgrid [X,Y]= meshgrid(x,y) x, y 为给定的向量， X, Y是网格划分后得到的网格矩阵 若 x = y, 则可简写为 [X,Y]= meshgrid(x) 例： x=[-8:0.5:8]; y=[-8:0.5:8];