matlab符号运算实验报告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 4 M2 = [a,b;c,d] M3 = [ a, b] [ c, d] 2.下面语句计算出来的 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=1 c1=0 c2=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)+v] A = [ 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*x^2+3*x^3+4*x^4;g=5*x^5-6*x^2+9; A=simple(f*g) B=simple(f/g) C=simple(f\g) int(f) int(g,1,2) diff(f,x,2) D=f/g limit(D,x,0) 用ezplot命令绘出 在[-4,4]之间的图像。 ezplot(x^2*sin(x),[-4,4]) 试分别生成 5 阶的单位阵、8 阶均匀分布的随机矩阵。 eye(5) rand(8) 生成列向量 x=[1, 3, 5, 7, 9, … , 29]。 a=[1:2:29] x=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; end disp(s); sum1(50) 1275 sum1(100) 5050 sum1(150) 11325 建立M函数文件作以下计算: 1) 自然数n的阶乘; 2) n中取m的组合。 （1）function f=mymutiply(n) if(n=1) f=1; else f=n*mymutiply(n-1); end （2）function assemblage1(n