资源描述：
《matlab符号运算实验报告》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、实验报告实验名称：数学建模课程设计1院系部:数理专业班级：计科1102学生姓名：张野学号:1111080223指导教师：石玉英成绩：华北电力大学9华北电力大学实习报告书学生姓名张野性别男成绩实习名称数学建模课程设计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=1234M2=[a,b;c,d]M3=[a,b][c,d]2.下面语句计算出来的c1，c2相等吗，为什么？
2、上机验证。>>a1=1e10;b1=1e-10;>>c1=(a1+b1-a1)/b1;c1=0>>a2=sym(a1);b2=sym(b1);9>>c2=(a2+b2-a2)/b2;c2=1c1=0c2=1变量类型不一样3.设A是一个符号矩阵(定义如下)，试指出findsym(A,1)的输出结果，并由此能得出什么结论？>>symsabtuvxy;>>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('1
3、2345678901234567890'))>>factor(‘12345678901234567890’)错误ans=(2)*(3)^2*(5)*(101)*(3803)*(3607)*(27961)*(3541)symsx;f=-x+2*x^2+3*x^3+4*x^4;g=5*x^5-6*x^2+9;A=simple(f*g)9B=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('x^2*sin(x)',[-4,4])试分别生成5阶的单位阵
4、、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.*B9生成由矩阵B与C矩阵乘积得到的矩阵D。D=B*C生成一个由矩阵D的第8、4、10、13行和第7、1、6、9、2列组成的子矩阵E。E=D([841013],[71692])求出矩阵E的最大元素。max(max(E))用ezplot3命令绘出单位球的图像symsa
5、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;fori=1:100s=s+i;end（2）functionsum1(n)s=0;9whilen>=1s=s+n;n=n-1;enddisp(s);>>sum1(50)1275>>sum1(100)5050>>sum1(150)11325建立M函数文件作以下计算:1)自然数n的阶乘;2)n中取m的组合。（1）functionf=mymu
6、tiply(n)if(n<=1)f=1;elsef=n*mymutiply(n-1);end（2）functionassemblage1(n,m)s=mymutiply(n)/mymutiply(n-m)/mymutiply(m);9disp(s);functionf=legrende(n)symsx;if(n==0)f=1;elseif(n==1)f=x;elseif(n>1)f=((2*n-1)*legrende(n-1)-(n-1)*legrende(n-2))/n;endendendfunctionhanshu(x)9if(x<1)disp(x)elseif(x>=1&x<10)di
7、sp(2*x-1);elsedisp(3*x-11);endendclear;A=[1468;2370;-2-157;11-10];B=[-3-76-5;-213-1;0013;26-10];A(A==-1)=-2;B(B<0)=1;B(B==0)=A(B==0);ABEnd用ezplot3命令绘出单位球的图像9symsab;ezsurf(sin(a)*cos(b),sin(a)*sin(b),cos(a),
