matlab符号运算实验报告.doc

实验报告实验名称： 数学建模课程设计1院 系 部: 数理 专业班级： 计科1102学生姓名： 张野 学 号: 1111080223指导教师： 石玉英 成 绩： 华北电力大学华北电力大学实习报告书学生姓名 张野性 别 男成 绩实习名称数学建模课程设计1 实习班级计科1102 实习周数1带队教师石玉英实习单位数理系地 点数理系机房 1.指出下面的 M1，M2，M3 分别是什么，并上机验证。 a=1; b=2; c=3;d=4; M1=a,b;c,d; M2=a,b;c,d; M3=sym(a,b;c,d);M1 = 1 2 3 4M2 =a,b;c,dM3 = a, b c, d2.下面语句计算出来的 c1，c2 相等吗，为什么？上机验证。 a1=1e10; b1=1e-10; c1=(a1+b1-a1)/b1; c1=0 a2=sym(a1); b2=sym(b1); c2=(a2+b2-a2)/b2; c2=1c1=0c2=1变量类型不一样3.设 A 是一个符号矩阵 (定义如下)，试指出 findsym(A,1) 的输出结果，并由此能得出什么结论？ syms a b t u v x y; A=a+b*x,sin(t)+u; x*exp(-t),log(y)+vA = a+b*x, sin(t)+u x*exp(-t), log(y)+v试指出下面两条命令的结果是否相同，如果不同，哪个是正确的？为什么？ factor(sym(12345678901234567890) factor(12345678901234567890) 错误ans =(2)*(3)2*(5)*(101)*(3803)*(3607)*(27961)*(3541) syms x;f=-x+2*x2+3*x3+4*x4;g=5*x5-6*x2+9; A=simple(f*g)B=simple(f/g)C=simple(fg)int(f)int(g,1,2)diff(f,x,2)D=f/glimit(D,x,0)用ezplot命令绘出 在-4,4之间的图像。ezplot(x2*sin(x),-4,4)试分别生成 5 阶的单位阵、8 阶均匀分布的随机矩阵。eye(5)rand(8)生成列向量 x=1, 3, 5, 7, 9, , 29。a=1:2:29x=a生成以题目2中向量 x 的元素为对角线的矩阵 A，并输出 A 的行数。A=diag(x)size(A,1)生成一个与 A 同阶的正态分布的随机矩阵 B。B=randn(15)生成由矩阵 A 与 B向量乘法得到的矩阵 C。C=A.*B生成由矩阵B 与 C矩阵乘积得到的矩阵 D。D=B*C生成一个由矩阵 D 的第 8、4、10、13 行和第 7、1、6、9、2 列组成的子矩阵 E。E=D(8 4 10 13,7 1 6 9 2)求出矩阵 E 的最大元素。max(max(E) 用ezplot3命令绘出单位球的图像 syms a b;ezsurf(sin(a)*cos(b),sin(a)*sin(b),cos(a),0,2*pi,0,2*pi)（1）建立M脚本文件求1到100的和。（2）建立M函数文件求 1+2+3+.+n并求n=50,100,150的值。（1）clear;s=0;for i=1:100 s=s+i;end（2）function sum1(n)s=0;while n=1 s=s+n; n=n-1;enddisp(s); sum1(50) 1275 sum1(100) 5050 sum1(150) 11325建立M函数文件作以下计算:1) 自然数n的阶乘;2) n中取m的组合。 （1）function f=mymutiply(n)if(n1) f=(2*n-1)*legrende(n-1)-(n-1)*legrende(n-2)/n; end endendfunction hanshu(x)if(x=1&x10) disp(2*x-1); else disp(3*x-11); endendclear;A=1 4 6 8;2 3 7 0;-2 -1 5 7;1 1 -1 0;B=-3 -7 6 -5;-2 1 3 -1;0 0 1 3;2 6 -1 0;A(A=-1)=-2;B(B0)=1;B(B=0)=A(B=